GRHD: also impose EoS limits when computing final primitives
This went unnoticed so far, because we have been running only with polytropes and there epsilon is basically unused.
This went unnoticed so far, because we have been running only with polytropes and there epsilon is basically unused.