diff options
Diffstat (limited to 'src')
| -rw-r--r-- | src/fastblas.cpp | 26 |
1 files changed, 11 insertions, 15 deletions
diff --git a/src/fastblas.cpp b/src/fastblas.cpp index c15e34f..e971b22 100644 --- a/src/fastblas.cpp +++ b/src/fastblas.cpp @@ -208,25 +208,21 @@ static void fast_cblas_dgemm(const char *TransA, const char *TransB, const doubl /* - Use the fasted/supported way to call BLAS dgemm + Use the fast/supported way to call BLAS dgemm */ void fast_dgemm(const char *TransA, const char *TransB, const double alpha, const gsl_matrix *A, const gsl_matrix *B, const double beta, gsl_matrix *C) { - if (is_legacy_mode()) { - eigenlib_dgemm(TransA,TransB,alpha,A,B,beta,C); - } else { - fast_cblas_dgemm(TransA,TransB,alpha,A,B,beta,C); - - #ifdef DISABLE - if (is_check_mode()) { - // ---- validate with original implementation - gsl_matrix *C1 = gsl_matrix_alloc(C->size1,C->size2); - eigenlib_dgemm(TransA,TransB,alpha,A,B,beta,C1); - enforce_msg(gsl_matrix_equal(C,C1),"dgemm outcomes are not equal for fast & eigenlib"); - gsl_matrix_free(C1); - } - #endif + fast_cblas_dgemm(TransA,TransB,alpha,A,B,beta,C); + +#ifdef DISABLE + if (is_check_mode()) { + // ---- validate with original implementation + gsl_matrix *C1 = gsl_matrix_alloc(C->size1,C->size2); + eigenlib_dgemm(TransA,TransB,alpha,A,B,beta,C1); + enforce_msg(gsl_matrix_equal(C,C1),"dgemm outcomes are not equal for fast & eigenlib"); + gsl_matrix_free(C1); } +#endif } |