tolerance cannot be choosen that large; in case of non-connected meshes we could end up interpolating too far outside the (wrong) element, leading to bad values

useful for remapping a mesh onto a (slightly modified) CAD, for sensitivity analyses

boundary layers in 3D are advancing (filtering intersecting layers using rtrees)
randomness of 2D meshes is now close to be ok (we still have problems with OCC ...)

improve reproductibility of 2D mesh algorithm (compare vertices and triangles by number, not by pointer)

when reparamOnFace fails, try to find the closest point before resorting to simply using the linear interpolation

fullMatrix::gemm

Allowed the periodical voronoi diagram using the -microstructure option to deal with rectangle of arbitrary size. Xmax, Ymax and Zmax are to be specified as entry 3 to 5 of the input file.

TODO : avoid topological unconsistencies for (e.g.) cones when
delaunay meshing is involved.