diff --git a/src/mpi_interface/MPI_Ops.cpp b/src/mpi_interface/MPI_Ops.cpp
index 91dbd36fe5ec81f14aa8bec868ebcf834c9af18c..93dc1f9d4ecec371443981518158fa68e6f586b0 100644
--- a/src/mpi_interface/MPI_Ops.cpp
+++ b/src/mpi_interface/MPI_Ops.cpp
@@ -59,18 +59,6 @@ void mpi_reduce_op::set_op(std::string project_type, double cross_cor_max, int n
 
 std::vector<node_sc_pair> mpi_reduce_op::select_top_feats(std::vector<node_sc_pair> in_vec_1, std::vector<node_sc_pair> in_vec_2)
 {
-    for(int ff = 0; ff < in_vec_1.size(); ++ff)
-    {
-        std::cout << ff << '\t' << in_vec_1[ff]._feat->feat(0)->arr_ind() << '\t' << in_vec_1[ff]._feat->expr() << std::endl;
-    }
-    std::cout << std::endl;
-
-    for(int ff = 0; ff < in_vec_2.size(); ++ff)
-    {
-        std::cout << ff << '\t' << in_vec_2[ff]._feat->feat(0)->arr_ind() << '\t' << in_vec_2[ff]._feat->expr() << std::endl;
-    }
-    std::cout << std::endl;
-
     // Set up an output vector
     std::vector<node_sc_pair> out_vec;
     out_vec.reserve(N_SIS_SELECT);
@@ -92,11 +80,5 @@ std::vector<node_sc_pair> mpi_reduce_op::select_top_feats(std::vector<node_sc_pa
         ++ff;
     }
 
-    for(int ff = 0; ff < out_vec.size(); ++ff)
-    {
-        std::cout << ff << '\t' << out_vec[ff]._feat->feat(0)->arr_ind() << '\t' << out_vec[ff]._feat->expr() << std::endl;
-    }
-    std::cout << std::endl;
-
     return out_vec;
 }