Skip to content

GRHD: Fixes to c2p

Florian Atteneder requested to merge fa/c2p-fix into main

This addresses two issues:

  • Utilize cold EoS when we limit velocity. This should avoid artificial spikes in epsilon (internal energy density) due to the limiting, because the derived relations used in the Kastaun paper no longer hold.
  • Make atm_threshold_factor = 1.0 the default. In some tests I noticed that having this greater than 1.0 an artificial plateau in the evolved density appears, which I did not like. This seems to fix that.

Merge request reports