mse.py -- Internal error with code 1 occurred -- stopping simulation. #12 (ERROR: quantity binary->dh_vec_dt[i] is infinite)

[pavanvyn]

Version bf1154c2efa9e47a68c9a613954a19396f4ad5e6

error_code = 1
configuration = "[2,2]"
N_bodies = 4
masses = [8.206047120001081, 0.444161072756691, 2.6660612887736503, 0.8227416666787222]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [0.21671790930287946, 0.23461307773500406, 1.1549400998669097]
eccentricities = [0.16740717501968838, 0.7217742010530627, 0.06183471994993741]
inclinations = [2.08829834904585, 1.2947270599667524, 1.7897735592924715]
arguments_of_pericentre = [2.753986993051959, 4.751985971525096, 3.356920200647366]
longitudes_of_ascending_node = [1.335398333382004, 3.2593494989656526, 2.5957654583003795]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]
kick_distribution_sigma_km_s_NS = 265.0
kick_distribution_sigma_km_s_BH = 50.0
common_envelope_timescale = 1e2

Error

tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity binary->dh_vec_dt[i] is infinite

[hamers]

Issues #52-#54 might be related to the moment of inertia being zero in stellar_evolution.cpp -- compute_moment_of_inertia(), but this is hard to verify because the issues are nearly impossible to reproduce. Version e27838d now checks for zero moment of inertia and will give a separate error code (38), which should help narrowing down the problem if it occurs again. I will close these issues for now.