apparently did the trick

src/OpenVolumeMesh/DifferentialOperators/EdgeWeightEvaluators.hh

@@ -74,8 +74,6 @@ public:
         double first_alpha = alpha;
         VecT first_nijk = nijk;
 
-        std::cout<<" first alpha : "<<first_alpha<<std::endl;
-
         hehf_iter++;
 
         int tet_count(0);
@@ -99,6 +97,7 @@ public:
 
             weight += tet_weight(alpha, beta, theta);
 
+#if 0
             std::cout<<" ===== TET "<<tet_count<<std::endl;
             std::cout<<"  xl : "<<xl<<std::endl;
             std::cout<<"  xk : "<<xk<<std::endl;
@@ -114,6 +113,7 @@ public:
             double tet_w = tet_weight(alpha, beta, theta);
             std::cout<<"  tet weight : "<<tet_w<<std::endl;
             std::cout<<"  --> cummulated weight : "<<weight<<std::endl;
+#endif
 
             //setting-up next iteration
             xk = xl;
@@ -131,7 +131,7 @@ public:
 
         weight += tet_weight(alpha, beta, theta);
 
-
+#if 0
         std::cout<<" ===== TET "<<tet_count<<std::endl;
         std::cout<<"  alpha : "<<alpha<<std::endl;
         std::cout<<"   beta : "<<beta<<std::endl;
@@ -142,6 +142,7 @@ public:
         double tet_w = tet_weight(alpha, beta, theta);
         std::cout<<"  tet weight : "<<tet_w<<std::endl;
         std::cout<<"  --> cummulated weight : "<<weight<<std::endl;
+#endif
 
         //multiply by Vol(i,j) (i.e. edge length)
         weight *= (mesh.vertex(xi) - mesh.vertex(xj)).norm() / 8.;
@@ -171,6 +172,7 @@ private:
             //weight computation following the paper's formula
             return cot(theta) * (2. * (cotan_alpha * cotan_beta)/cos(theta) - cotan_alpha*cotan_alpha - cotan_beta*cotan_beta);;
         }else{
+#if 0
             std::cout<<" can't compute tet weight : "<<std::endl;
             std::cout<<"  epsilon   : "<<std::numeric_limits<Scalar>::epsilon() <<std::endl;
             std::cout<<"  sin alpha : "<<std::setprecision(20)<<sin(alpha)<<std::endl;
@@ -186,7 +188,7 @@ private:
                        (fabs(sin(beta))  > std::numeric_limits<Scalar>::epsilon()) << " "<<
                        (fabs(sin(theta)) > std::numeric_limits<Scalar>::epsilon()) << " "<<
                        (fabs(cos(theta)) > std::numeric_limits<Scalar>::epsilon()) << std::endl;
-
+#endif
             return 0.f;
         }
     }

src/Unittests/unittests_laplacians.cc

@@ -250,7 +250,6 @@ TEST_F(DualLaplacianTest, GetPerHalfedgeWeight){
 }
 
 
-#if 0
 
 TEST_F(DualLaplacianTest, GetPerVertexLaplacian){
 
@@ -277,5 +276,4 @@ TEST_F(DualLaplacianTest, PrecomputedDualLaplacianGivesSameResults){
         auto on_the_fly = on_demand_laplacian[vertex];
     }
 }
-#endif