The current implementation of an accurate distance calculation is not accurate for distances less than 1m. Since the haversine function is more robust against rounding error of floating point arithmetic, the great circle distance should be replaced by an accurate implementation of the haversine function.

Related to #4498