diff --git a/dune/gfe/localgeodesicfefunction.hh b/dune/gfe/localgeodesicfefunction.hh
index f79f2c3b43094667119d79c0a83b3b9a274451d0..bb2098ef903f6078e7ac2cef83da3d6d1d7bd26f 100644
--- a/dune/gfe/localgeodesicfefunction.hh
+++ b/dune/gfe/localgeodesicfefunction.hh
@@ -313,17 +313,13 @@ evaluateDerivativeOfGradientWRTCoefficient(const Dune::FieldVector<ctype, dim>&
     // best way, though.
     Dune::FieldMatrix<ctype,embeddedDim,embeddedDim> dFdqPseudoInv = pseudoInverse(dFdq);
     
-    // Put it all together
-    for (size_t i=0; i<result.size(); i++) {
-        
-        //
-        Tensor3<double,embeddedDim,embeddedDim,embeddedDim> dvDqF
-           =  TargetSpace::thirdDerivativeOfDistanceSquaredWRTFirst1AndSecond2Argument(coefficients_[i], q);
-    
+    //
+    Tensor3<double,embeddedDim,embeddedDim,embeddedDim> dvDqF
+       =  TargetSpace::thirdDerivativeOfDistanceSquaredWRTFirst1AndSecond2Argument(coefficients_[coefficient], q);
     
+    // Put it all together
+    for (size_t i=0; i<result.size(); i++)
         result[i] = dFdqPseudoInv * ( dvDqF[i] * dFdqPseudoInv * dFdw - dpDwF[i]) * B;   
-     
-    }
     
 }
 #endif