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. #4 #32

Closed pavanvyn closed 3 years ago

pavanvyn commented 3 years ago

Version 7ec7a30659cbf87c09cd8ccc3fffedf1e529d1f4

configuration = "[2,2]"
N_bodies = 4
masses = [8.395638165566435, 9.683102968524427, 2.4079043649334837, 10.196160604508613]
metallicities = [0.017834708328831494, 0.017834708328831494, 0.017834708328831494, 0.017834708328831494]
semimajor_axes = [0.2937518022746324, 0.3631286282157685, 26.922506682980465]
eccentricities = [0.04849264493622715, 0.4220476144564787, 0.8049918930616927]
inclinations = [2.7172138113497986, 1.775144912962339, 0.24740766324612773]
arguments_of_pericentre = [5.289502192487941, 0.981029983263301, 2.9788044266407474]
longitudes_of_ascending_node = [0.11867389430683109, 0.6137180302196914, 1.8481696174701707]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

Error:

tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity spin_vec[0] is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity spin_vec[1] is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity spin_vec[2] is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity h_vec_unit[0] is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity h_vec_unit[1] is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity h_vec_unit[2] is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity e_vec_unit[0] is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity e_vec_unit[1] is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity e_vec_unit[2] is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity q_vec_unit[0] is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity q_vec_unit[1] is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity q_vec_unit[2] is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity spin_vec_dot_e_vec_unit is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity spin_vec_dot_h_vec_unit is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity spin_vec_dot_q_vec_unit is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity C is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity j_p4_inv is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity X_rot is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity Y_rot is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity Z_rot is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity Z_TB is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity h is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity n is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity binary->dh_vec_dt[i] is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity binary->de_vec_dt[i] is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity star->dspin_vec_dt[i] is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity binary->dh_vec_dt[i] is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity binary->de_vec_dt[i] is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity star->dspin_vec_dt[i] is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity binary->dh_vec_dt[i] is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity binary->de_vec_dt[i] is NaN
tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity star->dspin_vec_dt[i] is NaN
mass_changes.cpp -- ODE_handle_RLOF_triple_mass_transfer -- Omega_donor_dot nan

[CVODE ERROR]  CVode
  At t = 2.23957e+07 and h = 5.8391e-64, the corrector convergence test failed repeatedly or with |h| = hmin.

structure.cpp -- set_positions_and_velocities -- ERROR: quantity true_anomaly is NaN
structure.cpp -- set_positions_and_velocities -- ERROR: quantity r[0] is NaN
structure.cpp -- set_positions_and_velocities -- ERROR: quantity r[1] is NaN
structure.cpp -- set_positions_and_velocities -- ERROR: quantity r[2] is NaN
structure.cpp -- set_positions_and_velocities -- ERROR: quantity true_anomaly is NaN
structure.cpp -- set_positions_and_velocities -- ERROR: quantity r[0] is NaN
structure.cpp -- set_positions_and_velocities -- ERROR: quantity r[1] is NaN
structure.cpp -- set_positions_and_velocities -- ERROR: quantity r[2] is NaN
mse.py -- Internal error with code  1 occurred -- stopping simulation.
hamers commented 3 years ago

The triple MT routines could lead to NaNs when the MT rate was too low; should be fixed in 67fa481.