diff options
| author | Pjotr Prins | 2017-10-13 13:07:40 +0000 |
|---|---|---|
| committer | Pjotr Prins | 2017-10-13 15:27:24 +0000 |
| commit | a5cee35c058b48f725a9eceabc205eb0cd1cbb07 (patch) | |
| tree | 47c3d1f0b332490c9919033c7219f44010b2b886 /src/lmm.cpp | |
| parent | a610dd723a233aed1abe31aa32e3137b23b5f983 (diff) | |
| download | pangemma-a5cee35c058b48f725a9eceabc205eb0cd1cbb07.tar.gz | |
Replacing first dgemm - tests fail
Diffstat (limited to 'src/lmm.cpp')
| -rw-r--r-- | src/lmm.cpp | 3 |
1 files changed, 2 insertions, 1 deletions
diff --git a/src/lmm.cpp b/src/lmm.cpp index 71aa184..7b32330 100644 --- a/src/lmm.cpp +++ b/src/lmm.cpp @@ -42,6 +42,7 @@ #include "gzstream.h" #include "io.h" +#include "fastblas.h" #include "lapack.h" #include "lmm.h" @@ -1315,7 +1316,7 @@ void LMM::Analyze(std::function< SnpNameValues(size_t) >& fetch_snp, gsl_matrix_submatrix(UtXlarge, 0, 0, inds, 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); |