fixed Mesh::eval, surface plots look better now

* Mesh::eval is more tolerant (up to EQ_TOL to be precise) to the
barycentric coordinates so that if we evaluate corners we get proper
results displayed
* helmholtz uses the proper surface plots now for visualization