search

hamers / mse

Multiple Stellar Evolution -- A Population Synthesis Code for Multiple-Star Systems
GNU General Public License v3.0
9 stars 4 forks source link

mse.py -- Internal error with code 1 occurred -- stopping simulation. #13 (ERROR: quantity star->dspin_vec_dt[i] is infinite) #54

Closed pavanvyn closed 3 years ago

pavanvyn commented 3 years ago

Version bf1154c2efa9e47a68c9a613954a19396f4ad5e6

error_code = 1
configuration = "[2,2]"
N_bodies = 4
masses = [0.40861479984518045, 2.659738237884771, 10.448274790463161, 7.193906466262495]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [0.02216333783479387, 0.8355036354698007, 14.050528605772518]
eccentricities = [0.06605294043639066, 0.3586354128939225, 0.6963407484891593]
inclinations = [1.6051269334624083, 1.7958261441386474, 1.0391056042837474]
arguments_of_pericentre = [2.8339204378082146, 5.506120483974405, 1.483643493338533]
longitudes_of_ascending_node = [5.671699434031563, 4.881965555627381, 1.2134387231070838]
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 = 1e4

Error

tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity star->dspin_vec_dt[i] is infinite
hamers commented 3 years ago

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.