pub fn compute_pressure_from_eos( density: f64, target_density: f64, eos_scale: f64, eos_exponent: f64, negative_pressure_scale: f64, ) -> f64