about summary refs log tree commit diff
path: root/src/param.cpp
diff options
context:
space:
mode:
authorPjotr Prins2025-11-30 08:32:33 +0100
committerPjotr Prins2025-11-30 08:32:33 +0100
commit254113f4e02d5b423af302e9facb0ed580749ef5 (patch)
treee989992d84d210842d02936470d4f7f559dfd96e /src/param.cpp
parent882a426c70370e50a521435e1b2175598e0fe411 (diff)
downloadpangemma-254113f4e02d5b423af302e9facb0ed580749ef5.tar.gz
K: reads geno mdb
Diffstat (limited to 'src/param.cpp')
-rw-r--r--src/param.cpp3
1 files changed, 2 insertions, 1 deletions
diff --git a/src/param.cpp b/src/param.cpp
index 05e4a7b..3b74e53 100644
--- a/src/param.cpp
+++ b/src/param.cpp
@@ -1294,7 +1294,7 @@ void PARAM::ReadBIMBAMGenotypes(gsl_matrix *UtX, gsl_matrix *K, const bool calc_
 }
 
 void PARAM::MdbCalcKin() {
-  error = !mdb_calc_kin(file_geno, is_loco(), indicator_snp, a_mode - 20);
+  error = !mdb_calc_kin(file_geno, is_loco(), a_mode - 20, indicator_snp);
   return;
 }
 
@@ -1310,6 +1310,7 @@ void PARAM::CalcKin(gsl_matrix *matrix_kin) {
       error = true;
     }
   } else {
+    // indicator_snp is an array of bool
     error = !BimbamKin(file_geno, setKSnps, indicator_snp, a_mode - 20, d_pace,
                        matrix_kin, ni_max == 0);
   }