diff --git a/test/localgeodesicfestiffnesstest.cc b/test/localgeodesicfestiffnesstest.cc
index dbaf6ee68fbde875e0351cdc9ae295d945c333eb..e1106a9c4f9646dc9872c6538a94aa7b0689cbf3 100644
--- a/test/localgeodesicfestiffnesstest.cc
+++ b/test/localgeodesicfestiffnesstest.cc
@@ -130,7 +130,7 @@ void testUnitVector3d()
         for (size_t j=0; j<nDofs; j++)
             for (size_t k=0; k<nDofs; k++)
                 embeddedHessian[j][k] = TargetSpace::secondDerivativeOfDistanceSquaredWRTSecondArgument(coefficients[j],
-                                                                                                        coefficients[j]);
+                                                                                                        coefficients[k]);
         Matrix<FieldMatrix<double,2,2> > hessian(nDofs,nDofs);
         
         // transform to local tangent space bases