Absolute positions on the axis linking particles
The distance on the axis linking two cells is currently computed with respect to the coordinate origin of the whole volume not of the cell. That could lead to round-up errors as we only use floating point precision.
Should doubles be used instead ? Or can we store the position along the axis using the cell's corner as origin ? If so can the merger-sort still be used ?