diff --git a/dune/gfe/localgeodesicfefunction.hh b/dune/gfe/localgeodesicfefunction.hh
index eba6558bb9346ab09c18416d8ce92d9eae464816..e0d5cd6b19d35e74abc6a17fdb43e1895d483546 100644
--- a/dune/gfe/localgeodesicfefunction.hh
+++ b/dune/gfe/localgeodesicfefunction.hh
@@ -202,7 +202,7 @@ evaluateDerivative(const Dune::FieldVector<ctype, dim>& local) const
     TargetSpace q = evaluate(local);
 
     // Actually compute the derivative
-    evaluateDerivative(local,q);    
+    return evaluateDerivative(local,q);    
 }
 
 template <int dim, class ctype, class LocalFiniteElement, class TargetSpace>