diff --git a/dune/gfe/localgeodesicfefunction.hh b/dune/gfe/localgeodesicfefunction.hh
index 83187d0dded7cfac401315d2952b33102683f9a5..56a565551cec87155291042c6e41d7faffd27c17 100644
--- a/dune/gfe/localgeodesicfefunction.hh
+++ b/dune/gfe/localgeodesicfefunction.hh
@@ -301,11 +301,11 @@ evaluateDerivativeOfGradientWRTCoefficient(const Dune::FieldVector<ctype, dim>&
     assembler.assembleHessian(q,dFdq);
     
    
-	Dune::FieldMatrix<ctype,embeddedDim,embeddedDim> mixedDerivative = TargetSpace::secondDerivativeOfDistanceSquaredWRTFirstAndSecondArgument(coefficients_[coefficient], q);
+    Dune::FieldMatrix<ctype,embeddedDim,embeddedDim> mixedDerivative = TargetSpace::secondDerivativeOfDistanceSquaredWRTFirstAndSecondArgument(coefficients_[coefficient], q);
     Tensor3<double,embeddedDim,embeddedDim,dim+1> dpDwF(0);
-	for (int i=0; i<embeddedDim; i++)
-		for (int j=0; j<embeddedDim; j++)
-			dpDwF[i][j][coefficient] = mixedDerivative[i][j];
+    for (int i=0; i<embeddedDim; i++)
+        for (int j=0; j<embeddedDim; j++)
+            dpDwF[i][j][coefficient] = mixedDerivative[i][j];
     
     
     // dFDq is not invertible, if the target space is embedded into a higher-dimensional