diff --git a/src/localgeodesicfefunction.hh b/src/localgeodesicfefunction.hh
index 4842e7d1127239cf2c5e8cc8dfc05d9cbc3dc48d..aaf8af6e10ff051834c1cf6df842ef23f353d654 100644
--- a/src/localgeodesicfefunction.hh
+++ b/src/localgeodesicfefunction.hh
@@ -117,7 +117,7 @@ evaluateDerivativeFD(const Dune::FieldVector<ctype, dim>& local)
         forward[i]  += eps;
         backward[i] -= eps;
         
-        EmbeddedTangentVector fdDer = evaluate(forward) - evaluate(backward);
+        EmbeddedTangentVector fdDer = evaluate(forward).globalCoordinates() - evaluate(backward).globalCoordinates();
         fdDer /= 2*eps;
         
         for (int j=0; j<EmbeddedTangentVector::size; j++)