diff --git a/test/localgeodesicfestiffnesstest.cc b/test/localgeodesicfestiffnesstest.cc
index 08f0b7fcac029d9656bd6a0711114e4266d24a4f..426a4eb480606cf31c3ffed2b9c1481c910f74a9 100644
--- a/test/localgeodesicfestiffnesstest.cc
+++ b/test/localgeodesicfestiffnesstest.cc
@@ -146,11 +146,11 @@ void testHessian()
         embeddedHessian[1][1] = TargetSpace::secondDerivativeOfDistanceSquaredWRTSecondArgument(coefficients[0],
                                                                                                 coefficients[1]);
 
-        embeddedHessian[0][0] = TargetSpace::secondDerivativeOfDistanceSquaredWRTFirstAndSecondArgument(coefficients[0],
-                                                                                                coefficients[1]);
+        embeddedHessian[0][1] = TargetSpace::secondDerivativeOfDistanceSquaredWRTFirstAndSecondArgument(coefficients[0],
+                                                                                                        coefficients[1]);
 
-        embeddedHessian[0][0] = TargetSpace::secondDerivativeOfDistanceSquaredWRTFirstAndSecondArgument(coefficients[0],
-                                                                                                coefficients[1]);
+        embeddedHessian[1][0] = TargetSpace::secondDerivativeOfDistanceSquaredWRTFirstAndSecondArgument(coefficients[1],
+                                                                                                        coefficients[0]);
 
         Matrix<FieldMatrix<double,spaceDim,spaceDim> > hessian(nDofs,nDofs);