Remove unused code

1 files changed, 0 insertions, 39 deletions

diff --git a/src/ldr.cpp b/src/ldr.cpp
index f70eb85..a4ef40a 100644
--- a/src/ldr.cpp
+++ b/src/ldr.cpp
@@ -29,7 +29,6 @@
 #include <stdio.h>
 #include <stdlib.h>
 
-// #include "Eigen/Dense"
 #include "gsl/gsl_blas.h"
 #include "gsl/gsl_cdf.h"
 #include "gsl/gsl_eigen.h"
@@ -46,7 +45,6 @@
 #include "param.h"
 
 using namespace std;
-// using namespace Eigen;
 
 void LDR::CopyFromParam(PARAM &cPar) {
   a_mode = cPar.a_mode;
@@ -73,40 +71,3 @@ void LDR::CopyFromParam(PARAM &cPar) {
 
 void LDR::CopyToParam(PARAM &cPar) { return; }
 
-/*
-// X is a p by n matrix.
-void LDR::VB(const vector<vector<unsigned char>> &Xt, const gsl_matrix *W_gsl,
-             const gsl_vector *y_gsl) {
-
-  // Save gsl_vector and gsl_matrix into Eigen library formats.
-  MatrixXd W(W_gsl->size1, W_gsl->size2);
-  VectorXd y(y_gsl->size);
-  VectorXd x_col(y_gsl->size);
-
-  double d;
-  for (size_t i = 0; i < W_gsl->size1; i++) {
-    d = gsl_vector_get(y_gsl, i);
-    y(i) = d;
-    for (size_t j = 0; j < W_gsl->size2; j++) {
-      W(i, j) = gsl_matrix_get(W_gsl, i, j);
-    }
-  }
-
-  // Initial VB values by lm.
-  cout << indicator_snp[0] << " " << indicator_snp[1] << " " << indicator_snp[2]
-       << endl;
-  uchar_matrix_get_row(Xt, 0, x_col);
-
-  for (size_t j = 0; j < 10; j++) {
-    cout << x_col(j) << endl;
-  }
-
-  // Run VB iterations.
-  // TO DO.
-
-  // Save results.
-  // TO DO.
-
-  return;
-}
-*/