diff --git a/src/localgeodesicfefunction.hh b/src/localgeodesicfefunction.hh
index bb7bbd82891cdc96f18217e2a2757a0641d6e88e..06d3f944b11963e5e3a5402cbc52eb77a44f7226 100644
--- a/src/localgeodesicfefunction.hh
+++ b/src/localgeodesicfefunction.hh
@@ -191,14 +191,18 @@ template <int dim, class ctype, class TargetSpace>
 TargetSpace LocalGeodesicFEFunction<dim,ctype,TargetSpace>::
 evaluate(const Dune::FieldVector<ctype, dim>& local)
 {
-    assert(dim==2);
-    ctype extraCoord = 1-local[0]-local[1];
-    ctype factor = extraCoord+local[0];
+    ctype extraCoord = 1;
+    for (int i=0; i<dim; i++)
+        extraCoord -= local[i];
 
-    TargetSpace result = TargetSpace::interpolate(coefficients_[0], coefficients_[1], local[0]/factor);
+    ctype normalizingFactor = extraCoord+local[0];
 
-    for (int i=1; i<dim; i++)
-        result = TargetSpace::interpolate(result, coefficients_[i+1], local[i]);
+    TargetSpace result = TargetSpace::interpolate(coefficients_[0], coefficients_[1], local[0]/normalizingFactor);
+
+    for (int i=1; i<dim; i++) {
+        normalizingFactor += local[i];
+        result = TargetSpace::interpolate(result, coefficients_[i+1], local[i] / normalizingFactor);
+    }
 
     return result;
 }