Completely conservative and covariant numerical methodology for N-body problems with distance-dependent potentials