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.