Introduce first test

1 files changed, 4 insertions, 4 deletions

diff --git a/src/param.cpp b/src/param.cpp
index 9e55693..82121b9 100644
--- a/src/param.cpp
+++ b/src/param.cpp
@@ -318,7 +318,7 @@ void PARAM::ReadFiles(void) {
 
     trim_individuals(indicator_idv, ni_max);
     trim_individuals(indicator_cvt, ni_max);
-    // The following reads the geno file to get the SNPs
+    // The following reads the geno file to get the SNPs - only for BIMBAM
     if (is_bimbam && ReadFile_bimbam_geno(file_geno, setSnps, W2, indicator_idv, indicator_snp,
                                           maf_level, miss_level, hwe_level, r2_level, mapRS2chr,
                                           mapRS2bp, mapRS2cM, snpInfo, ns_test) == false) {
@@ -392,7 +392,7 @@ void PARAM::ReadFiles(void) {
     infile.clear();
   }
 
-  // Read genotype and phenotype file for multiple BIMBAM files.
+  // Read genotype and phenotype file for multiple BIMBAM files. Note this goes untested.
   if (!file_mgeno.empty()) {
 
     // Annotation file before genotype file.
@@ -1273,8 +1273,8 @@ void PARAM::ReadBIMBAMGenotypes(gsl_matrix *UtX, gsl_matrix *K, const bool calc_
       error = true;
     }
   } else {
-      // Read BIMBAM. Not supposed to be mdb
-      if (is_mdb || ReadFile_geno(file_geno, indicator_idv, indicator_snp, UtX, K, calc_K) == false) {
+      // Read BIMBAM
+      if (is_bimbam && ReadFile_geno(file_geno, indicator_idv, indicator_snp, UtX, K, calc_K) == false) {
       error = true;
     }
   }