.

src/data-structures/program_state.hh

@@ -315,7 +315,7 @@ template <class VectorTEMPLATE, class ScalarVectorTEMPLATE> class ProgramState {
     solveLinearProblem(dirichletNodes, contactNetwork.matrices().elasticity, ell0, u,
                        parset.sub("u0.solver"));
 
-    //print(u, "initial u:");
+    print(u, "initial u:");
 
     // Initial acceleration: Computed in agreement with Ma = ell0 - Au
     // (without Dirichlet constraints), again assuming dPhi(v = 0) = 0
@@ -348,7 +348,7 @@ template <class VectorTEMPLATE, class ScalarVectorTEMPLATE> class ProgramState {
     solveLinearProblem(noNodes, contactNetwork.matrices().mass, accelerationRHS, a,
                        parset.sub("a0.solver"));
 
-    //print(a, "initial a:");
+    print(a, "initial a:");
   }
 
 private:

src/utils/debugutils.hh

@@ -118,7 +118,8 @@
    void print(const std::vector<T>& x, std::string message){
       std::cout << message << std::endl;
       for (size_t i=0; i<x.size(); i++) {
-          std::cout << x[i] << std::endl;
+          //std::cout << x[i] << std::endl;
+          print(x[i], "");
       }
       std::cout << std::endl << std::endl;
    }