From cc840ea408a6dea80f41c308cdcbc0e9f146ccc1 Mon Sep 17 00:00:00 2001
From: Pjotr Prins
Date: Wed, 11 Oct 2017 11:00:17 +0000
Subject: Fix embarrassing bug that broke BIMBAM for larger datasets - found it
 when porting Plink to the new Analyze function

---
 src/lmm.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

(limited to 'src/lmm.cpp')

diff --git a/src/lmm.cpp b/src/lmm.cpp
index 1193700..134fbf9 100644
--- a/src/lmm.cpp
+++ b/src/lmm.cpp
@@ -1414,7 +1414,7 @@ void LMM::AnalyzeBimbam(const gsl_matrix *U, const gsl_vector *eval,
     gsl_vector_memcpy(&Xlarge_col.vector, x);
     c++; // count SNPs going in
 
-    if (c == msize)
+    if (c % msize == 0)
       batch_compute(msize);
   }
   batch_compute(c % msize);
-- 
cgit v1.2.3