void OpenMPBLAS3::symm( const enum CBLAS_ORDER order, const enum CBLAS_SIDE Side, const enum CBLAS_UPLO Uplo, const IndexType m, const IndexType n, const float alpha, const float* A, const IndexType lda, const float* B, const IndexType ldb, const float beta, float* C, const IndexType ldc, SyncToken* syncToken ) { if ( syncToken ) { LAMA_LOG_WARN( logger, "no asynchronous execution for openmp possible at this level." ) } char SD, UL; #ifdef F77_CHAR F77_CHAR F77_SD, F77_UL; #else #define F77_SD &SD #define F77_UL &UL #endif #ifdef F77_INT F77_INT F77_M = m, F77_N = n, F77_lda = lda, F77_ldb = ldb; F77_INT F77_ldc = ldc; #else #define F77_M m #define F77_N n #define F77_lda lda #define F77_ldb ldb #define F77_ldc ldc #endif IndexType RowMajorStrg; RowMajorStrg = 0; if ( order == CblasColMajor ) { if ( Side == CblasRight ) { SD = 'R'; } else if ( Side == CblasLeft ) { SD = 'L'; } else { BLASHelper::XERBLA_cpu( RowMajorStrg, 2, "cblas_ssymm", "Illegal Side setting, %d\n", Side ); RowMajorStrg = 0; return; } if ( Uplo == CblasUpper ) { UL = 'U'; } else if ( Uplo == CblasLower ) { UL = 'L'; } else { BLASHelper::XERBLA_cpu( RowMajorStrg, 3, "cblas_ssymm", "Illegal Uplo setting, %d\n", Uplo ); RowMajorStrg = 0; return; } #ifdef F77_CHAR F77_UL = C2F_CHAR( &UL ); F77_SD = C2F_CHAR( &SD ); #endif F77_ssymm( F77_SD, F77_UL, &F77_M, &F77_N, &alpha, A, &F77_lda, B, &F77_ldb, &beta, C, &F77_ldc ); } else if ( order == CblasRowMajor ) { RowMajorStrg = 1; if ( Side == CblasRight ) { SD = 'L'; } else if ( Side == CblasLeft ) { SD = 'R'; } else { BLASHelper::XERBLA_cpu( RowMajorStrg, 2, "cblas_ssymm", "Illegal Side setting, %d\n", Side ); RowMajorStrg = 0; return; } if ( Uplo == CblasUpper ) { UL = 'L'; } else if ( Uplo == CblasLower ) { UL = 'U'; } else { BLASHelper::XERBLA_cpu( RowMajorStrg, 3, "cblas_ssymm", "Illegal Uplo setting, %d\n", Uplo ); RowMajorStrg = 0; return; } #ifdef F77_CHAR F77_UL = C2F_CHAR( &UL ); F77_SD = C2F_CHAR( &SD ); #endif F77_ssymm( F77_SD, F77_UL, &F77_N, &F77_M, &alpha, A, &F77_lda, B, &F77_ldb, &beta, C, &F77_ldc ); } else { BLASHelper::XERBLA_cpu( RowMajorStrg, 1, "cblas_ssymm", "Illegal order setting, %d\n", order ); } RowMajorStrg = 0; return; }
void cblas_ssymm(const enum CBLAS_ORDER Order, const enum CBLAS_SIDE Side, const enum CBLAS_UPLO Uplo, const int M, const int N, const float alpha, const float *A, const int lda, const float *B, const int ldb, const float beta, float *C, const int ldc) { char SD, UL; #ifdef F77_CHAR F77_CHAR F77_SD, F77_UL; #else #define F77_SD &SD #define F77_UL &UL #endif #ifdef F77_INT F77_INT F77_M=M, F77_N=N, F77_lda=lda, F77_ldb=ldb; F77_INT F77_ldc=ldc; #else #define F77_M M #define F77_N N #define F77_lda lda #define F77_ldb ldb #define F77_ldc ldc #endif extern int CBLAS_CallFromC; extern int RowMajorStrg; RowMajorStrg = 0; CBLAS_CallFromC = 1; if( Order == CblasColMajor ) { if( Side == CblasRight) SD='R'; else if ( Side == CblasLeft ) SD='L'; else { cblas_xerbla(2, "cblas_ssymm", "Illegal Side setting, %d\n", Side); CBLAS_CallFromC = 0; RowMajorStrg = 0; return; } if( Uplo == CblasUpper) UL='U'; else if ( Uplo == CblasLower ) UL='L'; else { cblas_xerbla(3, "cblas_ssymm", "Illegal Uplo setting, %d\n", Uplo); CBLAS_CallFromC = 0; RowMajorStrg = 0; return; } #ifdef F77_CHAR F77_UL = C2F_CHAR(&UL); F77_SD = C2F_CHAR(&SD); #endif F77_ssymm(F77_SD, F77_UL, &F77_M, &F77_N, &alpha, A, &F77_lda, B, &F77_ldb, &beta, C, &F77_ldc); } else if (Order == CblasRowMajor) { RowMajorStrg = 1; if( Side == CblasRight) SD='L'; else if ( Side == CblasLeft ) SD='R'; else { cblas_xerbla(2, "cblas_ssymm", "Illegal Side setting, %d\n", Side); CBLAS_CallFromC = 0; RowMajorStrg = 0; return; } if( Uplo == CblasUpper) UL='L'; else if ( Uplo == CblasLower ) UL='U'; else { cblas_xerbla(3, "cblas_ssymm", "Illegal Uplo setting, %d\n", Uplo); CBLAS_CallFromC = 0; RowMajorStrg = 0; return; } #ifdef F77_CHAR F77_UL = C2F_CHAR(&UL); F77_SD = C2F_CHAR(&SD); #endif F77_ssymm(F77_SD, F77_UL, &F77_N, &F77_M, &alpha, A, &F77_lda, B, &F77_ldb, &beta, C, &F77_ldc); } else cblas_xerbla(1, "cblas_ssymm", "Illegal Order setting, %d\n", Order); CBLAS_CallFromC = 0; RowMajorStrg = 0; return; }