diff --git a/test/harmonicenergytest.cc b/test/harmonicenergytest.cc
index 402b86e6a55a71ed09a28e5c0edcab7d6c10505b..7a15f469c98d3898fe79dc9a9ee8c3a3c3095f2d 100644
--- a/test/harmonicenergytest.cc
+++ b/test/harmonicenergytest.cc
@@ -85,4 +85,15 @@ int main(int argc, char** argv)
     coefficients[2] = uv;
     
     testEnergy<GridType>(grid, coefficients);
+
+    // //////////////////////////////////////////////////////////
+    //   Test the approximation to the Hesse matrix
+    // //////////////////////////////////////////////////////////
+    
+    HarmonicEnergyLocalStiffness<GridType::LeafGridView,TargetSpace> assembler;
+
+    assembler.assemble(*grid->leafbegin<0>(), coefficients);
+
+    std::cout << "Hessian: \n" << assembler.A[0][0] << std::endl;
+
 }