simpler evaluate of the function value.  Now only works with in the triangle, because there are some singularities lurking outside

[[Imported from SVN: r4076]]