Reverted on using eigenlib_dgemm because of hanging K compute with dataset

1 files changed, 6 insertions, 0 deletions

diff --git a/src/fastblas.cpp b/src/fastblas.cpp
index 2d9b6e5..81705a6 100644
--- a/src/fastblas.cpp
+++ b/src/fastblas.cpp
@@ -228,3 +228,9 @@ void fast_dgemm(const char *TransA, const char *TransB, const double alpha,
   }
 #endif
 }
+
+void fast_eigen_dgemm(const char *TransA, const char *TransB, const double alpha,
+                      const gsl_matrix *A, const gsl_matrix *B, const double beta,
+                      gsl_matrix *C) {
+  eigenlib_dgemm(TransA,TransB,alpha,A,B,beta,C);
+}