@@ -386,9 +386,9 @@ errno_t GPU_SVD_computeControlMatrix(
386
386
387
387
save_fits (ID_VTmatrix_name , "!matVT.fits" );
388
388
389
- cudaStat = cublasSgemm (cublasH , CUBLAS_OP_T , CUBLAS_OP_T , n , m , n , & alpha , d_VT ,
389
+ cublasStatus_t cublasStat = cublasSgemm (cublasH , CUBLAS_OP_T , CUBLAS_OP_T , n , m , n , & alpha , d_VT ,
390
390
n , d_U , m , & beta , d_M , n );
391
- if (cudaStat != cudaSuccess )
391
+ if (cublasStat != CUBLAS_STATUS_SUCCESS )
392
392
{
393
393
printf ("cublasSgemm returned error code %d, line(%d)\n" , cudaStat , __LINE__ );
394
394
exit (EXIT_FAILURE );
@@ -411,11 +411,16 @@ errno_t GPU_SVD_computeControlMatrix(
411
411
}
412
412
413
413
414
- ID_Cmatrix = create_image_ID (ID_Cmatrix_name ,
415
- data .image [ID_Rmatrix ].md [0 ].naxis , arraysizetmp , _DATATYPE_FLOAT , 0 , 0 );
414
+ ID_Cmatrix = create_image_ID (
415
+ ID_Cmatrix_name ,
416
+ data .image [ID_Rmatrix ].md [0 ].naxis ,
417
+ arraysizetmp ,
418
+ _DATATYPE_FLOAT ,
419
+ 0 ,
420
+ 0 );
416
421
417
422
418
- // cudaStat = cudaMemcpy(data.image[ID_Cmatrix].array.F, d_M, sizeof(float)*m*n, cudaMemcpyDeviceToHost);
423
+ // cudaStat = cudaMemcpy(data.image[ID_Cmatrix].array.F, d_M, sizeof(float)*m*n, cudaMemcpyDeviceToHost);
419
424
420
425
h_M = (float * ) malloc (sizeof (float ) * m * n );
421
426
cudaStat = cudaMemcpy (h_M , d_M , sizeof (float ) * m * n , cudaMemcpyDeviceToHost );
0 commit comments