Fix embarrassing bug that broke BIMBAM for larger datasets - found it when porting Plink to the new Analyze function

1 files changed, 1 insertions, 1 deletions

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);