Approximation-based methods, such as the cubetree algorithm, have proven to be significantly faster than traditionalmethods for complex force evaluations near small irregular bodies. Suchmethods also hold the promise of simplifying the inclusion of experimental data to update the force model. However, the cubetree algorithm does not preserve intrinsic properties of the gravitational force such ...