Centre Blaise Pascal » HPL sur GPU

#ifdef HPL_CALL_CUBLAS

   double                    alpha = ALPHA, beta = BETA;

   int                       IONE = 1;

#define CUBLASM                 M

#define CUBLASN                 N

#define CUBLASK                 K

#define CUBLASlda               LDA

#define CUBLASldb               LDB

#define CUBLASldc               LDC

   char                      ctransa, ctransb;

   if(      TRANSA == HplNoTrans ) ctransa = 'N';

   else if( TRANSA == HplTrans   ) ctransa = 'T';

   else                            ctransa = 'C';

   if(      TRANSB == HplNoTrans ) ctransb = 'N';

   else if( TRANSB == HplTrans   ) ctransb = 'T';

   else                            ctransb = 'C';

   if( ORDER == HplColumnMajor )

   {

      CUBLAS_DGEMM( &ctransa, &ctransb, &CUBLASM, &CUBLASN, &CUBLASK, 

		    &alpha, A, &CUBLASlda, B, &CUBLASldb, &beta, C, &CUBLASldc,

		    &IONE, &IONE );

   }

   else

   {

      CUBLAS_DGEMM( &ctransb, &ctransa, &CUBLASN, &CUBLASM, &CUBLASK, 

		    &alpha, B, &CUBLASldb, A, &CUBLASlda, &beta, C, &CUBLASldc,

		    &IONE, &IONE );

   }

#endif