  RelativisticDynamics.jl is an open-source Julia package for relativistic
spin-orbital dynamics in the gravitational strong-field of a Kerr spacetime.
Existing codes for modelling the dynamics of spinning objects like pulsars in
the strong-field regime are generally lacking, since such systems occupy an
intermediate regime that is generally overlooked. At the "low" end of this
regime there are post-Newtonian descriptions which neglect the influence of the
pulsar spin on the underlying spacetime metric ("spin-curvature" coupling). At
the "high" end there are the full numerical relativity solutions which are
primarily applicable toe two black holes with a mass ratio $\mathcal{O}(1)$,
and are computationally intractable for pulsar systems observed over a large
number of orbital cycles. RelativisticDynamics.jl aims to bridge this gap by
providing a modern, fast code for accurate numerical evolution of spinning
relativistic systems via the Mathisson-Papetrou-Dixon formalism. Julia is a
modern language that solves the "two language problem", enabling fast dynamic
typing and JIT compilation on conjunction with petaflop performance, comparable
with numerical languages that are better known in the astrophysics community
such as C or Fortran. RelativisticDynamics.jl is written to be fully type
flexible, being able to support arbitrary number formats, and fully
differentiable via automatic differentiation.
