BUGFIX: coordinates on subsimplices need to be normalized (only 2d, so far)

[[Imported from SVN: r4067]]

src/localgeodesicfefunction.hh

@@ -191,7 +191,11 @@ template <int dim, class ctype, class TargetSpace>
 TargetSpace LocalGeodesicFEFunction<dim,ctype,TargetSpace>::
 evaluate(const Dune::FieldVector<ctype, dim>& local)
 {
-    TargetSpace result = TargetSpace::interpolate(coefficients_[0], coefficients_[1], local[0]);
+    assert(dim==2);
+    ctype extraCoord = 1-local[0]-local[1];
+    ctype factor = extraCoord+local[0];
+
+    TargetSpace result = TargetSpace::interpolate(coefficients_[0], coefficients_[1], local[0]/factor);
 
     for (int i=1; i<dim; i++)
         result = TargetSpace::interpolate(result, coefficients_[i+1], local[i]);