Refactor GRHD project

Major rework of GRHD project to include new broadcasting interface. The goal is to get a TOV evolution with the AV method.

PR contains various refactors:

• initialdata/: We promote this to a separate env in order to add CairoMakie for plotting.
• refactored and improve initialdata/SphericalAccretion.jl solver so that it agrees with various references: solver from bugner's thesis (extract from bamps), plots from two papers from Papadoulos.
• refactored initialdata/TOV.jl to not depend on dg1d.jl anymore, instead use OrindaryDiffEq and CairoMakie for plots
• copy of cons2prim from SRHD. these two versions should be merged eventually, but will do in a separate PR
• implement spherical1d formulation: this is a rewrite of the evolution GRHD equations to work with non-zero radial shift vector. this is needed to evolve the spherical accretion test in cowling approximation where there should be no dynamics.
• implement rescaled_spherical1d formulation: this is a rewrite of the previous formulation based on the generalized Valencia formulation (cf. Montero, Baumgarte paper from 2017); the auxiliary metric is chosen such that coordinate singularities at r=0 in the source terms are cancelled; doesn't work atm
• allow the spherical1d to include r=0 in the domain by using symmetry conditions to evaluate source terms at r=0 (workaround for rescaled_spherical1d not working right now)
• new ref tests: bondi accretion and TOV star
• implement a differentiate! method on Mesh in order to compute derivatives by hand, if needed

added 1d GRHD equation initialdata refactor labels

added 2 commits

• 2fc2ef07 - refactor TOV solver to use OrdinaryDiffEq
• 6e00d156 - wrap SphericalAccretion solver into module; recompute bondi_accretion.h5 with less radial points

Compare with previous version

changed the description

added 12 commits

• 6e00d156...1b4ed3ed - 2 earlier commits
• d14ab9d5 - implement initialdata setup for tov case
• 6c245a13 - add specialized equations for the rescaled GRHD version
• 96328fae - specialize timespeed computation based on equations
• b3e57ccd - implement rhs! for rescaled_spherical1d formulation
• 3fad66fe - add atmosphering to cons2prim
• f848d425 - rename cons2prim routines
• 05ed0660 - derivative computation needs jacobi
• 691dd6ed - fixup initialdata setup
• 4c721c2e - fix typos
• 2027484a - fix missing lapse factor in source terms

Compare with previous version

added 1 commit

• 46aebd69 - update notebooks

Compare with previous version

added 1 commit

• a961f68b - add grhd example par files

Compare with previous version

added 9 commits

• c83dc738 - mesh: implement differentiate for FVElement
• c9469a1d - implement equations for FVElement and rescaled_spherical1d
• b16652f0 - implement rhs for FVElement and rescaled_spherical1d
• 47d8fe53 - add formulation parameter
• e2dba042 - implement tov initial data for spherical1d formulation
• c9e7fd29 - fix: set atm density in bondi_accretion initial data
• 3a777c26 - cons2prim: use non-zero values for rho, eps thresholds
• f019c514 - wip: change maxspeed implementation for tov star initial data
• 304a8e96 - add tov example that runs with spherical1d formulation

Compare with previous version

added 2 commits

• fbc98d10 - make impose_symmetry! dispatch on project
• efe0920c - make spherical1d formulation work when r=0 is included by removing potential...

Compare with previous version

added 4 commits

• d21e08de - cleanup
• 5d166d3f - spherical1d: fix missing shift from maxspeed computation
• df95ad12 - IntegrationTests: remove reltol from grhd_bondi_accretion test
• 9ad3049d - IntegartionTests: add grhd_tov_spherical1d ref test

Compare with previous version

changed the description

added 3 commits

• 29d34fee - cleanup: remove codegen for now
• bff4db28 - remove comments
• 757f358b - update examples once more

Compare with previous version

changed the description

enabled an automatic merge when the pipeline for 757f358b succeeds

merged

mentioned in commit 50ac5d57