This fully implements self-gravity bar the periodic task over MPI.
Could you please check that is does not break anything on the hydro side ?
If you are interested in running with self-gravity I'd recommend using the updated EAGLE_12 case and run with
My next steps is to work on the formal requirement of using
-D alongside the
-G. I also need to finish working out the automated scripts to generate code for the higher-order multipoles but that is orthogonal to the integration of the code into the master.