about summary refs log tree commit diff
path: root/src
diff options
context:
space:
mode:
authorPjotr Prins2021-08-16 10:56:02 +0200
committerPjotr Prins2021-08-16 10:56:02 +0200
commit0dbe0670fe437ddadd713f8abc404a4be215ac7f (patch)
tree75d425d21560fa8580b2c4ec0331feca9e3a5351 /src
parent42de3667c90964206c3764be1867527d8eeb8cfc (diff)
downloadpangemma-0dbe0670fe437ddadd713f8abc404a4be215ac7f.tar.gz
README.md
README
Diffstat (limited to 'src')
-rw-r--r--src/mvlmm.cpp6
1 files changed, 6 insertions, 0 deletions
diff --git a/src/mvlmm.cpp b/src/mvlmm.cpp
index 4b8db05..51ed901 100644
--- a/src/mvlmm.cpp
+++ b/src/mvlmm.cpp
@@ -3666,6 +3666,8 @@ void MVLMM::AnalyzePlink(const gsl_matrix *U, const gsl_vector *eval,
     t_last++;
   }
   for (vector<SNPINFO>::size_type t = 0; t < snpInfo.size(); ++t) {
+    cout << snpInfo[t].rs_number << endl;
+
     if (t % d_pace == 0 || t == snpInfo.size() - 1) {
       ProgressBar("Reading SNPs", t, snpInfo.size() - 1);
     }
@@ -3752,6 +3754,10 @@ void MVLMM::AnalyzePlink(const gsl_matrix *U, const gsl_vector *eval,
       gsl_matrix_set_zero(Xlarge);
 
       for (size_t i = 0; i < l; i++) {
+        double t = (double(clock()-time_start)/(double(CLOCKS_PER_SEC)));
+        printf("%.2f ",t);
+        if (i>10) exit(1);
+        cout << i << "," << l  << endl;
         gsl_vector_view UtXlarge_col = gsl_matrix_column(UtXlarge, i);
         gsl_vector_memcpy(&X_row.vector, &UtXlarge_col.vector);