TODO self-gravity

TODO using Tobias Moldenhauer’s correction to Baruteau’s implementation