diff --git a/dune/gfe/localgeodesicfefunction.hh b/dune/gfe/localgeodesicfefunction.hh
index 7b554c8822311d7d4887f2f2cfa230c4dc7554ac..f79f2c3b43094667119d79c0a83b3b9a274451d0 100644
--- a/dune/gfe/localgeodesicfefunction.hh
+++ b/dune/gfe/localgeodesicfefunction.hh
@@ -317,11 +317,11 @@ evaluateDerivativeOfGradientWRTCoefficient(const Dune::FieldVector<ctype, dim>&
     for (size_t i=0; i<result.size(); i++) {
         
         //
-        Tensor3<double,embeddedDim,embeddedDim,embeddedDim> dcDqF
+        Tensor3<double,embeddedDim,embeddedDim,embeddedDim> dvDqF
            =  TargetSpace::thirdDerivativeOfDistanceSquaredWRTFirst1AndSecond2Argument(coefficients_[i], q);
     
     
-        result[i] = dFdqPseudoInv * ( dcDqF[i] * dFdqPseudoInv * dFdw - dpDwF[i]) * B;   
+        result[i] = dFdqPseudoInv * ( dvDqF[i] * dFdqPseudoInv * dFdw - dpDwF[i]) * B;   
      
     }