gemma,lmm,mvlmm: Add data debug printing everywhere.

1 files changed, 2 insertions, 0 deletions

diff --git a/src/mvlmm.cpp b/src/mvlmm.cpp
index 1d05bf1..f9cb0a2 100644
--- a/src/mvlmm.cpp
+++ b/src/mvlmm.cpp
@@ -2768,6 +2768,7 @@ void MphInitial(const size_t em_iter, const double em_prec,
                 gsl_matrix *B) {
 
   debug_msg("MphInitial");
+  write(Y, "Y in MphInitial");
   gsl_matrix_set_zero(V_g);
   gsl_matrix_set_zero(V_e);
   gsl_matrix_set_zero(B);
@@ -2971,6 +2972,7 @@ double PCRT(const size_t mode, const size_t d_size, const double p_value,
 void MVLMM::AnalyzeBimbam(const gsl_matrix *U, const gsl_vector *eval,
                           const gsl_matrix *UtW, const gsl_matrix *UtY) {
   debug_msg("entering");
+  write(UtY, "UtY in AnalyzeBimbam");
   igzstream infile(file_geno.c_str(), igzstream::in);
   if (!infile) {
     cout << "error reading genotype file:" << file_geno << endl;