Geodesolver for the Geodesic Equation

NASA has a need to determine the distribution of dark matter around planets like Earth. Doing so necessitates the ability to repeatedly solve the geodesic equation in simulations of weakly interacting particles streaming through bodies described by a Schwarzschild metric. Geodesolver numerically solves the geodesic equation orders of magnitudes faster than Runge-Kutta for arbitrary density profiles.