mse.py -- Internal error with code 10 occurred -- stopping simulation.

[pavanvyn]

Version 93c1a62cf452ccacb9eb6f980c5a376a755c532c

configuration = "[2,2]"
N_bodies = 4
masses = [1.356199705391964, 3.699593611492698, 7.898460782002915, 10.061139971442456]
metallicities = [0.023261271697795446, 0.023261271697795446, 0.023261271697795446, 0.023261271697795446]
semimajor_axes = [0.31067780415984153, 0.2239681079857236, 6.587690513916891]
eccentricities = [0.45020820700243996, 0.3787241714302784, 0.48009146224834015]
inclinations = [1.5537506705381328, 1.3403534424487569, 1.6354386296853225]
arguments_of_pericentre = [5.942264725130637, 5.080407386264119, 5.290640482035836]
longitudes_of_ascending_node = [1.924628206941482, 0.5288133897436507, 2.014825166740314]
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 -- outer_binary 5 donor 3 inner_binary 4
mass_changes.cpp -- ODE_handle_RLOF_triple_mass_transfer -- Omega_donor_dot nan

[CVODE ERROR]  CVode
At t = 1.81046e+08 and h = 1.09347e-63, the corrector convergence test failed repeatedly or with |h| = hmin.

mse.py -- Internal error with code  10 occurred -- stopping simulation.

[hamers]

There was a problem related to the triple CE routine; should now be fixed in a422cc7.