diff --git a/dune/gfe/localgeodesicfefunction.hh b/dune/gfe/localgeodesicfefunction.hh
index 8cff20deaedb0ae2d28f949a84ea09c6ccde28e6..10f794906cb7fdd2d9a32959a8d3a090c9cdd392 100644
--- a/dune/gfe/localgeodesicfefunction.hh
+++ b/dune/gfe/localgeodesicfefunction.hh
@@ -490,7 +490,7 @@ evaluateDerivativeOfGradientWRTCoefficient(const Dune::FieldVector<ctype, dim>&
     for (size_t k=0; k<coefficients_.size(); k++) {
         Dune::SymmetricMatrix<RT,embeddedDim> hesse = TargetSpace::secondDerivativeOfDistanceSquaredWRTSecondArgument(coefficients_[k], q);
         for (int i=0; i<embeddedDim; i++)
-            for (int j=0; j<embeddedDim; j++)
+            for (int j=0; j<i; j++)
                 dqdwF(i, j, k) = dqdwF(j, i, k) = hesse(i,j);
 
     }