diff options
author | Pjotr Prins | 2020-05-22 11:21:45 -0500 |
---|---|---|
committer | Pjotr Prins | 2020-05-22 11:21:45 -0500 |
commit | f1cd914e6f20c9a162e16d7283477c1b98d005d1 (patch) | |
tree | 0a23068f9b06525ded025450c209b3a1dcf94b38 /src/lmm.cpp | |
parent | 862ace03212ed17bdc1ab349bfab55543720a980 (diff) | |
parent | b309569fe9497befa008ac2d2cbc04f2e861ce76 (diff) | |
download | pangemma-f1cd914e6f20c9a162e16d7283477c1b98d005d1.tar.gz |
Merge branch 'master' of github.com:genenetwork/GEMMA
Diffstat (limited to 'src/lmm.cpp')
-rw-r--r-- | src/lmm.cpp | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/lmm.cpp b/src/lmm.cpp index 16fc3c8..b811594 100644 --- a/src/lmm.cpp +++ b/src/lmm.cpp @@ -1582,7 +1582,7 @@ void LMM::Analyze(std::function< SnpNameValues(size_t) >& fetch_snp, continue; double geno = gs[i]; - if (std::isnan(geno)) { + if (isnan(geno)) { gsl_vector_set(x_miss, pos, 1.0); n_miss++; } else { @@ -1829,7 +1829,7 @@ void LMM::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); |