https://hgpu.org/?p=1338
Integrating Post-Newtonian Equations on Graphics Processing Units