aboutsummaryrefslogtreecommitdiff
path: root/src/lmm.cpp
diff options
context:
space:
mode:
authorPjotr Prins2020-05-22 11:21:45 -0500
committerPjotr Prins2020-05-22 11:21:45 -0500
commitf1cd914e6f20c9a162e16d7283477c1b98d005d1 (patch)
tree0a23068f9b06525ded025450c209b3a1dcf94b38 /src/lmm.cpp
parent862ace03212ed17bdc1ab349bfab55543720a980 (diff)
parentb309569fe9497befa008ac2d2cbc04f2e861ce76 (diff)
downloadpangemma-f1cd914e6f20c9a162e16d7283477c1b98d005d1.tar.gz
Merge branch 'master' of github.com:genenetwork/GEMMA
Diffstat (limited to 'src/lmm.cpp')
-rw-r--r--src/lmm.cpp4
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);