jesterTOV.utils.locate_lowest_non_causal_point

jesterTOV.utils.locate_lowest_non_causal_point#

locate_lowest_non_causal_point(cs2)[source]#

Find the first point where the equation of state becomes non-causal.

The speed of sound squared \(c_s^2 = dp/d\varepsilon\) must satisfy \(c_s^2 \leq 1\) (in units where \(c = 1\)) for causality. This function locates the first density where this condition is violated.

Parameters:

cs2 (Float[Array, "n"]) – Speed of sound squared values

Return type:

Int[Array, '']

Returns:

Int[Array, “”] – Scalar array: Index of first non-causal point, or -1 if EOS is everywhere causal

Notes

This function is used to determine the maximum valid central pressure for TOV equation integration.