diff options
| author | Pjotr Prins | 2017-10-13 15:23:00 +0000 |
|---|---|---|
| committer | Pjotr Prins | 2017-10-13 15:27:24 +0000 |
| commit | fdb48997ee3ed2b326a92f8e0cc7f72a4b38d8c8 (patch) | |
| tree | 5f62e06dbacdec2d4ee60da9112508615d42fc1f /src/mvlmm.cpp | |
| parent | 7eca3c49b7790007a4190b73209cab9ffb2bb117 (diff) | |
| download | pangemma-fdb48997ee3ed2b326a92f8e0cc7f72a4b38d8c8.tar.gz | |
Refactored debug settings
Replaced eigenlib_dgemm with fast_dgemm - 10-30% speed gain for GEMMA
Diffstat (limited to 'src/mvlmm.cpp')
| -rw-r--r-- | src/mvlmm.cpp | 7 |
1 files changed, 4 insertions, 3 deletions
diff --git a/src/mvlmm.cpp b/src/mvlmm.cpp index c0dc143..95ef14b 100644 --- a/src/mvlmm.cpp +++ b/src/mvlmm.cpp @@ -39,6 +39,7 @@ #include "gsl/gsl_vector.h" #include "eigenlib.h" +#include "fastblas.h" #include "gzstream.h" #include "io.h" #include "lapack.h" @@ -3250,8 +3251,8 @@ void MVLMM::AnalyzeBimbam(const gsl_matrix *U, const gsl_vector *eval, gsl_matrix_submatrix(UtXlarge, 0, 0, UtXlarge->size1, l); time_start = clock(); - eigenlib_dgemm("T", "N", 1.0, U, &Xlarge_sub.matrix, 0.0, - &UtXlarge_sub.matrix); + fast_dgemm("T", "N", 1.0, U, &Xlarge_sub.matrix, 0.0, + &UtXlarge_sub.matrix); time_UtX += (clock() - time_start) / (double(CLOCKS_PER_SEC) * 60.0); gsl_matrix_set_zero(Xlarge); @@ -3717,7 +3718,7 @@ void MVLMM::AnalyzePlink(const gsl_matrix *U, const gsl_vector *eval, gsl_matrix_submatrix(UtXlarge, 0, 0, UtXlarge->size1, l); time_start = clock(); - eigenlib_dgemm("T", "N", 1.0, U, &Xlarge_sub.matrix, 0.0, + fast_dgemm("T", "N", 1.0, U, &Xlarge_sub.matrix, 0.0, &UtXlarge_sub.matrix); time_UtX += (clock() - time_start) / (double(CLOCKS_PER_SEC) * 60.0); |