diff --git a/test/localgeodesicfestiffnesstest.cc b/test/localgeodesicfestiffnesstest.cc
index 73dabd801a985b3a5cfa55a63fe199d754f517a5..08f0b7fcac029d9656bd6a0711114e4266d24a4f 100644
--- a/test/localgeodesicfestiffnesstest.cc
+++ b/test/localgeodesicfestiffnesstest.cc
@@ -179,6 +179,15 @@ void testHessian()
         std::cout << "hessian:" << std::endl;
         printmatrix(std::cout, hessian, "hessian", "--");
         
+        ///////////////////////////////////////////////////////////////////////////////////
+        //  Abort if there is a difference
+        ///////////////////////////////////////////////////////////////////////////////////
+        Matrix<FieldMatrix<double,spaceDim,spaceDim> > difference = hessian;
+        difference -= fdHessian;
+        
+        if (difference.infinity_norm() > 1e-4)
+            assert(false);
+        
     }
 
 }