SRHD: Fix FV slope limiter method to use up-to-date primitives in evolution
Due to the cons2prim
methods being tied to D, S, tau
registered as dynamic_variables
, we need some juggling to get the primitives of the current step.
This inconvenience is due to
- not (yet) separating the time stepping from the flux computation for the FV methods,
- and the 'hard coding' of the primitives to the
cons2prim
methods.
Side note: We face the same issue with the GRHD
version.
This PR also adds some more reftests to cover the blastwave1
exercise.
To this end we also updated the FV solver to work with non-periodic bdry conditions with bdry values taken from initial data.