This allows to merge the points of unstructured grids
using a double precision tolerance. The previous
implementation only allowed to merge points using
single precision tolerance which might not be accurate
enough for some applications.

Moreover, added a test to ensure that this new functionality
works as expected.