From ec8a139cb07b46bd0b9e5de2ea8db1f7f335a56c 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 c0ecb2a..4bca3fb 100644
--- a/src/lmm.cpp
+++ b/src/lmm.cpp
@@ -1424,7 +1424,7 @@ void LMM::Analyze(std::function< SnpNameValues(size_t) >& fetch_snp,
     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