Révision 9 src/blas/HPL_dtrsm.c

HPL_dtrsm.c (revision 9)
871 871
                  TRANS, DIAG, N, M, ALPHA, A, LDA, B, LDB );
872 872
   }
873 873
#endif
874

  
874 875
#ifdef HPL_CALL_FBLAS
875 876
   double                    alpha = ALPHA;
876 877
#ifdef StringSunStyle
......
969 970
#endif
970 971
   }
971 972
#endif
973

  
974
#ifdef HPL_CALL_CUBLAS
975
   double                    alpha = ALPHA;
976

  
977
   int                       IONE = 1;
978

  
979
#define  CUBLASM                M
980
#define  CUBLASN                N
981
#define  CUBLASlda              LDA
982
#define  CUBLASldb              LDB
983

  
984
   char                      cside, cuplo, ctran, cdiag;
985

  
986
   if(      TRANS == HplNoTrans ) ctran = 'N';
987
   else if( TRANS == HplTrans   ) ctran = 'T';
988
   else                           ctran = 'C';
989
   cdiag = ( DIAG == HplUnit  ? 'U' : 'N' );
990

  
991
   if( ORDER == HplColumnMajor )
992
   {
993
      cside = ( SIDE == HplRight ? 'R' : 'L' );
994
      cuplo = ( UPLO == HplLower ? 'L' : 'U' );
995

  
996
      CUBLAS_DTRSM( &cside, &cuplo, &ctran, &cdiag, &CUBLASM, &CUBLASN, &alpha,
997
		    A, &CUBLASlda, B, &CUBLASldb, &IONE, &IONE, &IONE, &IONE );
998
   }
999
   else
1000
   {
1001
      cside = ( SIDE == HplRight ? 'L' : 'R' );
1002
      cuplo = ( UPLO == HplLower ? 'U' : 'L' );
1003

  
1004
      CUBLAS_DTRSM( &cside, &cuplo, &ctran, &cdiag, &CUBLASN, &CUBLASM, &alpha,
1005
		    A, &CUBLASlda, B, &CUBLASldb, &IONE, &IONE, &IONE, &IONE );
1006
   }
1007
#endif
972 1008
/*
973 1009
 * End of HPL_dtrsm
974 1010
 */

Formats disponibles : Unified diff