Skip to content

Crash in planetary SESAME equation of state

The function SESAME_pressure_from_internal_energy() can return infinite values which then propagate in the code and lead to a crash (eventually by having particles with NaN positions as a result of infinite accelerations). In some cases the denominator in the equation for intp_u_2 turns out to be zero, meaning there are two consecutive values in table_log_u_rho_T which are equal. For instance the entries 59241 and 59242.

If the tables are correct then that function should be made robust to this. Or something needs to be done to prevent crashes down the line.

To upload designs, you'll need to enable LFS and have an admin enable hashed storage. More information