hamers / mse Goto Github PK

configuration = "[2,2]"
N_bodies = 4
masses = [1.57769899738381, 8.154962292095899, 1.0343146279646283, 0.18580975749217013]
metallicities = [0.0212712771643236, 0.0212712771643236, 0.0212712771643236, 0.0212712771643236]
semimajor_axes = [1.5440469715166245, 0.20541527337354934, 3533.170180677806]
eccentricities = [0.7035748259529817, 0.5248321204665087, 0.14734933489357616]
inclinations = [0.8487787595371594, 0.7101780895642751, 1.3124128651467024]
arguments_of_pericentre = [2.3129289604683665, 3.4062372205116516, 0.8662035988977094]
longitudes_of_ascending_node = [3.28413868542158, 0.7185962342019538, 0.021155103046979377]
tend = 1.000000E+08
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [31.997069636649776, 28.53940494801786, 6.9278001267971145, 1.2520098815628715]
metallicities = [0.016362325519875612, 0.016362325519875612, 0.016362325519875612, 0.016362325519875612]
semimajor_axes = [8.895146947970506, 0.3470533934010974, 8177.992315216426]
eccentricities = [0.3951159692596551, 0.3573953911246622, 0.9029334691318619]
inclinations = [1.496563735256394, 1.611801200734346, 2.4233252268815257]
arguments_of_pericentre = [4.197566048960249, 4.78663949033743, 4.569441274290801]
longitudes_of_ascending_node = [3.293170967703992, 0.05172382002093702, 1.368786108311765]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

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
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.

configuration = "2"
N_bodies = 2
masses = [4.02704075397777, 4.379012755929251]
metallicities = [0.027766048578118464, 0.027766048578118464]
semimajor_axes = [0.10025035833229334]
eccentricities = [0.5714652662243841]
inclinations = [2.835480367248477]
arguments_of_pericentre = [5.451797077856235]
longitudes_of_ascending_node = [5.838150796777557]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1]
object_types = [1, 1]

configuration = "2"
N_bodies = 2
masses = [4.5079693703370385, 5.674085796223652]
metallicities = [0.018763436864135542, 0.018763436864135542]
semimajor_axes = [0.12266441097847466]
eccentricities = [0.46544978345249416]
inclinations = [2.317283864797076]
arguments_of_pericentre = [2.0248133366672803]
longitudes_of_ascending_node = [6.03295667773142]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1]
object_types = [1, 1]

configuration = "2"
N_bodies = 2
masses = [3.9550374262920838, 4.32011032535225]
metallicities = [0.02812278814862606, 0.02812278814862606]
semimajor_axes = [0.15784038080509885]
eccentricities = [0.5787429789505243]
inclinations = [1.470918421189333]
arguments_of_pericentre = [1.0856897921808524]
longitudes_of_ascending_node = [3.8745510443739994]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1]
object_types = [1, 1]

configuration = "2"
N_bodies = 2
masses = [4.256086321614061, 3.8849988196712117]
metallicities = [0.01264658763353502, 0.01264658763353502]
semimajor_axes = [0.10359400348173937]
eccentricities = [0.47165459989818936]
inclinations = [1.5427614144716344]
arguments_of_pericentre = [5.825841253572492]
longitudes_of_ascending_node = [3.8046920904743344]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1]
object_types = [1, 1]

configuration = "2"
N_bodies = 2
masses = [3.859754117852526, 4.011959825109891]
metallicities = [0.01923419736817146, 0.01923419736817146]
semimajor_axes = [0.17721185333622858]
eccentricities = [0.4422824723597685]
inclinations = [2.620776275744103]
arguments_of_pericentre = [2.1475615587259593]
longitudes_of_ascending_node = [0.36744133224443276]
tend = 1.400000E+10
N_steps = 1000  

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
ODE_mass_changes.cpp -- ERROR in determine_E_0; E_0 = nan x  = nan e = nan (1.0/e)*(1.0 - 1.0/x) nan
ODE_mass_changes.cpp -- ERROR in determine_E_0; E_0 = nan x  = nan e = nan (1.0/e)*(1.0 - 1.0/x) nan

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

configuration = "[2,2]"
N_bodies = 4
masses = [7.540948688753799, 5.542806479907622, 1.6449369158111955, 8.523708895122976]
metallicities = [0.021331238724547827, 0.021331238724547827, 0.021331238724547827, 0.021331238724547827]
semimajor_axes = [193.12588279845278, 18.119725906827203, 4256.999687317027]
eccentricities = [0.7524873566774626, 0.15263088650184092, 0.3869766276151553]
inclinations = [0.8612693368583286, 2.5342304447719104, 1.048379593931638]
arguments_of_pericentre = [0.5027362011630949, 3.6066003077451625, 1.0431612857126726]
longitudes_of_ascending_node = [6.264325417514426, 4.867578242273096, 3.7973749864625987]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [3.762254469849947, 7.349744147945759, 8.579556659374338, 47.04141214391311]
metallicities = [0.02728435020513672, 0.02728435020513672, 0.02728435020513672, 0.02728435020513672]
semimajor_axes = [6.20045452202843, 4.1103742181516845, 277.63581001067365]
eccentricities = [0.8265914585781401, 0.7442378738496561, 0.2934071857082655]
inclinations = [1.8436231207845388, 2.9116722325910405, 2.5066085538633645]
arguments_of_pericentre = [4.576486222775116, 1.7128248807379545, 1.0029244850244559]
longitudes_of_ascending_node = [6.160384217581603, 5.154441680206168, 4.286053523206743]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

[CVODE ERROR]  CVode
  At t = 7.34677e+06 and h = 0.000473526, the error test failed repeatedly or with |h| = hmin.

[CVODE ERROR]  CVodeGetDky
 Illegal value for t.t = 7.96553e+06 is not between tcur - hu = 7.96553e+06 and tcur = 7.96553e+06.

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]

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.

configuration = "[2,2]"
N_bodies = 4
masses = [15.323692042564135, 1.6540947669787935, 5.241738203733809, 3.5086521703113407]
metallicities = [0.023706677179601062, 0.023706677179601062, 0.023706677179601062, 0.023706677179601062]
semimajor_axes = [0.2723205994887504, 6.670910534497935, 1494.5689757397226]
eccentricities = [0.9486004207915687, 0.7598476935402387, 0.3575651064444511]
inclinations = [2.097232229847697, 1.0032437330254709, 1.1710219687310373]
arguments_of_pericentre = [3.013109209250083, 5.304241335167159, 2.5010940601444904]
longitudes_of_ascending_node = [0.7535527705006922, 0.2052269088280549, 2.088198932501593]
tend = 1.000000E+08
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [2.1520091801507344, 5.554758696239336, 9.212293366739805, 1.6459606034011065]
metallicities = [0.027386795625515235, 0.027386795625515235, 0.027386795625515235, 0.027386795625515235]
semimajor_axes = [2.382875060027812, 0.7544912181403461, 811.252248356948]
eccentricities = [0.8093054948682695, 0.31201634649125515, 0.6729912139877858]
inclinations = [0.6649404868857367, 1.7232537730978637, 2.159940140983064]
arguments_of_pericentre = [4.36295901814875, 1.432200347389523, 5.704844019422479]
longitudes_of_ascending_node = [2.6315280576610904, 3.064845943730461, 4.36370731366349]
tend = 1.000000E+08
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [47.643992189057975, 50.445366858771266, 45.25627339750874, 17.62837673193206]
metallicities = [0.012862801475893733, 0.012862801475893733, 0.012862801475893733, 0.012862801475893733]
semimajor_axes = [0.6689433925633086, 0.28122681885568385, 3867.268837446437]
eccentricities = [0.46369985674959907, 0.1768538859917899, 0.6226246106696653]
inclinations = [0.9837918622467343, 2.0849834819627473, 1.5386344118860225]
arguments_of_pericentre = [6.050730734464304, 1.5180479341436717, 0.5739194600286861]
longitudes_of_ascending_node = [3.1026518087318626, 2.3138754908474626, 2.1752334688954122]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

stellar_evolution.cpp -- evolve_stars -- ERROR: quantity sse_main_sequence_timescale is infinite
stellar_evolution.cpp -- evolve_stars -- ERROR: quantity spin_vec[0] is NaN
stellar_evolution.cpp -- evolve_stars -- ERROR: quantity spin_vec[1] is NaN
stellar_evolution.cpp -- evolve_stars -- ERROR: quantity spin_vec[2] is NaN
mse.py -- Internal error with code  1 occurred -- stopping simulation.

configuration = "[2,2]"
N_bodies = 4
masses = [17.70951475073602, 1.4958771186273825, 12.296863337016326, 9.762121983758005]
metallicities = [0.01271026805068375, 0.01271026805068375, 0.01271026805068375, 0.01271026805068375]
semimajor_axes = [1.2526386291001999, 0.7567784142837688, 10.648335006122476]
eccentricities = [0.0899656679182423, 0.4019667882222402, 0.30648264211183285]
inclinations = [1.9375972000564443, 2.8475517094939904, 1.7992830237120063]
arguments_of_pericentre = [0.41898377406774956, 0.21106201244258357, 6.240703876311492]
longitudes_of_ascending_node = [2.8247485391533376, 0.44179318199817735, 1.700820500240725]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

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]

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.

configuration = "[2,2]"
N_bodies = 4
masses = [2.755442754007535, 2.0332518165757953, 2.8774179125584496, 8.108252480444714]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [32.97933822901889, 42.980491285413216, 572.5055242696573]
eccentricities = [0.9188877553282521, 0.7596664642379991, 0.5577001633296457]
inclinations = [1.2265908662498253, 0.46121911223032525, 1.6047317212982855]
arguments_of_pericentre = [4.410327569281705, 3.674837786568595, 3.1890961286489623]
longitudes_of_ascending_node = [1.9388515780323803, 3.8140139924068066, 0.24438883452436594]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [7.7652942105962754, 10.112395173547398, 0.2843108569399502, 1.7796145333159648]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [2.3308716518022834, 0.6052503586975648, 445.2335224138562]
eccentricities = [0.21170415809675638, 0.4634892414029424, 0.7570790839985683]
inclinations = [1.4584412919632603, 1.200859224863009, 1.6843651335054188]
arguments_of_pericentre = [0.05565943570831119, 5.038530769485884, 3.3255979179661614]
longitudes_of_ascending_node = [3.9651442726276085, 5.37931691352577, 4.068974312830167]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [4.844524218671435, 8.11891899747903, 7.833078005839005, 0.46753945289745924]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [12.62371561001661, 6.004199751613622, 6736.579264887057]
eccentricities = [0.764259851004382, 0.383047422054649, 0.8409918894045305]
inclinations = [1.011816183280177, 2.0220865953864893, 2.022218711344736]
arguments_of_pericentre = [2.0056590053492602, 2.7923321965459107, 2.240399302614785]
longitudes_of_ascending_node = [3.362167998015483, 0.7757501917277166, 3.640222539487186]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

N_bodies = 4
masses = [35.48451711269334, 3.754291259792952, 1.5698448495850958, 2.458995707154814]
metallicities = [0.02856425149565036, 0.02856425149565036, 0.02856425149565036, 0.02856425149565036]
semimajor_axes = [13.3262155657255, 3.4979312630635655, 6609.517083248268]
eccentricities = [0.23737870375549824, 0.23197438121205205, 0.49232074675514736]
inclinations = [1.6635484797295157, 0.6920134796825302, 1.4541731092472017]
arguments_of_pericentre = [1.6311984421749086, 5.668718789942289, 3.0727957388875455]
longitudes_of_ascending_node = [0.8760625104853672, 3.237978932106665, 0.8929969114254929]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [22.004427064590327, 10.67637249995399, 21.350801343398434, 17.83059913226916]
metallicities = [0.027777104472852476, 0.016972739261366773, 0.006195420811250992, 0.025374648115150562]
semimajor_axes = [16.24594083770805, 194.98606734659958, 1652.989574959285]
eccentricities = [0.22516602530571117, 0.05854204586705469, 0.16390500070201974, 0.2066390631287358]
inclinations = [1.5488332269228329, 1.8906322113481837, 0.8630347336065893]
arguments_of_pericentre = [0.8330545554295615, 2.481413211077939, 2.5436993516042894]
longitudes_of_ascending_node = [3.9221880631496515, 3.980372509445738, 2.7136189206018697]
tend = 1.000000E+08
N_steps = 5000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]
plot_filename = "results/test"

configuration = "[2,2]"
N_bodies = 4
masses = [2.8133821822469973, 12.037556492139071, 5.149535129780804, 3.2424237734711907]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [0.18443531422416343, 9.427516288686817, 931.5276795166421]
eccentricities = [0.39437497094845303, 0.9315504971650672, 0.6032846058387555]
inclinations = [2.7348886113270523, 1.463945003952422, 2.8466968152999352]
arguments_of_pericentre = [0.25002794691340974, 6.046526651354287, 5.641895170882967]
longitudes_of_ascending_node = [1.817969185732822, 3.3471582774311806, 3.610853397758905]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [6.358273438605692, 14.27935838907645, 13.698266857550399, 15.055584377109724]
metallicities = [0.02312100278834555, 0.02312100278834555, 0.02312100278834555, 0.02312100278834555]
semimajor_axes = [0.17330036459059014, 5.0296529663235665, 7009.032006154533]
eccentricities = [0.2921602295459753, 0.4052202950794419, 0.44846009559641115]
inclinations = [2.0332354658512686, 1.5546687625255196, 2.4161044234182047]
arguments_of_pericentre = [5.615224200633939, 6.098070470598655, 5.365432975087523]
longitudes_of_ascending_node = [4.5944239685851755, 3.329072562051022, 0.8348800010005164]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

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.

configuration = "2"
N_bodies = 2
masses = [3.0570774984253286, 3.4335505060992015]
metallicities = [0.02, 0.02]
semimajor_axes = [0.10308355683896057]
eccentricities = [0.4772310671310243]
inclinations = [1.2968326028476609]
arguments_of_pericentre = [6.203836066580958]
longitudes_of_ascending_node = [4.0540468980697595]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1]
object_types = [1, 1]

configuration = "2"
N_bodies = 2
masses = [1.465998881349993, 1.2309546317228084]
metallicities = [0.02, 0.02]
semimajor_axes = [0.20439504179243606]
eccentricities = [0.887249092698013]
inclinations = [1.8833950553687608]
arguments_of_pericentre = [3.260033210427437]
longitudes_of_ascending_node = [5.828964710110773]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1]
object_types = [1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [1.172601355296052, 0.9705824261600782, 3.9456487519118917, 8.292214646979753]
metallicities = [0.017433728781051934, 0.017433728781051934, 0.017433728781051934, 0.017433728781051934]
semimajor_axes = [1.3326804114590938, 2.4482621313522195, 3539.4224364683173]
eccentricities = [0.5992544499080654, 0.43727958575359643, 0.9127659988692236]
inclinations = [1.8921877658757038, 1.845711838692201, 1.1913708833707037]
arguments_of_pericentre = [4.5914563893933185, 3.3145582976268524, 2.3573766578639206]
longitudes_of_ascending_node = [3.598051398180878, 2.7318439241437784, 3.556589478138784]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [62.159704663171716, 59.774910076680904, 31.417961828319992, 2.03289118796106]
metallicities = [0.02493711162848645, 0.02493711162848645, 0.02493711162848645, 0.02493711162848645]
semimajor_axes = [5.0602172728575505, 8.333273949324111, 1225.6916909975089]
eccentricities = [0.4623847217492819, 0.732755835950576, 0.7914329440080197]
inclinations = [2.324101373467737, 0.9582927328595943, 1.3051385381800236]
arguments_of_pericentre = [3.895245648188573, 4.459100942021893, 0.1100232142371]
longitudes_of_ascending_node = [1.985013688142783, 1.000873058513418, 4.103131684156803]
tend = 1.000000E+08
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [7.755298110394679, 13.02455302992831, 1.7306612927312375, 1.7045592918749755]
metallicities = [0.027715890415101217, 0.027715890415101217, 0.027715890415101217, 0.027715890415101217]
semimajor_axes = [0.9629082312726244, 3.162434677535477, 1800.1953209157884]
eccentricities = [0.33701500691806957, 0.5524906988041355, 0.7528957425541132]
inclinations = [0.7291109296646838, 1.7463778481022854, 0.9665018746484156]
arguments_of_pericentre = [2.49910590507181, 2.9581175250276113, 2.5587319122352388]
longitudes_of_ascending_node = [4.960759633082254, 5.883604722095698, 4.674119055391521]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

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]

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.

configuration = "[2,2]"
N_bodies = 4
masses = [1.645565025632771, 10.348235262980317, 9.55095296584496, 0.3480582406322891]
metallicities = [0.014968509519472145, 0.014968509519472145, 0.014968509519472145, 0.014968509519472145]
semimajor_axes = [2.42088693886329, 0.4146994168496801, 20.664357544851985]
eccentricities = [0.9092548265625824, 0.776923012142269, 0.25593059458018647]
inclinations = [2.9932024686188776, 2.0347693879441078, 1.8489847085043605]
arguments_of_pericentre = [4.880462107911308, 6.075874899920778, 2.435156369599932]
longitudes_of_ascending_node = [4.314600705932014, 6.251153046365085, 4.685161323858285]
tend = 1.000000E+08
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]


configuration = "[2,2]"
N_bodies = 4
masses = [4.383480384532576, 7.160840701715892, 20.851071525285477, 3.6247548403600547]
metallicities = [0.028496476268530253, 0.028496476268530253, 0.028496476268530253, 0.028496476268530253]
semimajor_axes = [0.1788656870848023, 2.291532219553283, 46.44540686369012]
eccentricities = [0.8173456420460182, 0.4537599738970811, 0.5207020213932612]
inclinations = [2.2280224331462257, 1.6849197839593588, 1.5621041180066848]
arguments_of_pericentre = [5.650009690126222, 4.842792598611465, 4.494645511518449]
longitudes_of_ascending_node = [3.5981902103200354, 2.0645303010088742, 0.4403064017173386]
tend = 1.000000E+08
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [11.545724854578776, 1.0881783040872224, 3.3115343779347834, 6.010473778620222]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [218.78017428049276, 153.5405019638936, 703.9382128960965]
eccentricities = [0.8217645445206639, 0.6461807039125479, 0.07027389562515694]
inclinations = [1.897725151813237, 1.482110956942365, 1.672984517289217]
arguments_of_pericentre = [5.215497922485309, 5.6057896113067995, 1.5135481989561066]
longitudes_of_ascending_node = [2.855244987772052, 3.9616874484784317, 0.1728142629395301]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

index_system = 7909
configuration = "[2,2]"
N_bodies = 4
masses = [48.723227950163476, 50.55586184999055, 0.6663968367064421, 14.291396349689649]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [121.29162525631203, 0.9079747821422273, 1818.947985533156]
eccentricities = [0.6091125996566156, 0.5409414732050581, 0.6442229873070678]
inclinations = [2.4878784849388733, 2.3424456314611697, 1.4797802972432195]
arguments_of_pericentre = [4.547909648545598, 4.293142200399169, 1.7682030312785086]
longitudes_of_ascending_node = [4.4724965520039754, 0.6355538540261445, 3.3518132859648975]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]


index_system = 8752
configuration = "[2,2]"
N_bodies = 4
masses = [45.14307817214505, 11.374085884881296, 13.649374604859617, 2.8441445240656473]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [50.272222946543, 1.7421167030987932, 723.9109160983443]
eccentricities = [0.2653987139971449, 0.9207728286617188, 0.6233839885058325]
inclinations = [0.9559966649182329, 2.2097000017077995, 1.2783637914899817]
arguments_of_pericentre = [1.536014917333147, 5.2193136310201265, 5.723101792252989]
longitudes_of_ascending_node = [0.07171332655419127, 3.167398716289923, 0.28751130541481895]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]


index_system = 86146
configuration = "[2,2]"
N_bodies = 4
masses = [8.018034511510015, 17.292270297474218, 75.62450534711071, 56.75109107398421]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [0.3472542357103747, 6.422814495917538, 30.683676666048434]
eccentricities = [0.37292476017542925, 0.4119348325795392, 0.26490743523065124]
inclinations = [2.655225572716907, 2.3983294395351242, 2.45618057056905]
arguments_of_pericentre = [5.7358755736208415, 5.377547214351855, 4.193824944719226]
longitudes_of_ascending_node = [3.012021595002201, 1.636517908311905, 2.8414774897562167]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [8.161200003208195, 8.140763068381933, 10.007763242074555, 2.1034835491014237]
metallicities = [0.015708744421724558, 0.015708744421724558, 0.015708744421724558, 0.015708744421724558]
semimajor_axes = [1.0222922886794827, 0.4107728806816176, 11.240011114483286]
eccentricities = [0.5677485083081685, 0.34906630495318375, 0.4589338159127716]
inclinations = [1.0280068750594389, 1.4389260994929471, 1.8113949640485276]
arguments_of_pericentre = [5.935484660745404, 0.6222960663547583, 4.873459527988544]
longitudes_of_ascending_node = [3.4790712090292044, 1.0115987157490396, 0.408639005084607]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

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.79223e+07 and h = 1.30343e-64, the corrector convergence test failed repeatedly or with |h| = hmin.

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

configuration = "[2,2]"
N_bodies = 4
masses = [5.454617915708907, 5.261190930216584, 8.740613846022839, 6.112137453098097]
metallicities = [0.019345190451989645, 0.019345190451989645, 0.019345190451989645, 0.019345190451989645]
semimajor_axes = [2.2761090983089503, 0.8429224364730526, 28.184848414403664]
eccentricities = [0.43540600066038854, 0.8102339449906077, 0.5878448254861636]
inclinations = [1.4165173412384477, 0.5303433685865508, 1.9108843710441161]
arguments_of_pericentre = [4.654042896241758, 1.2820386595450664, 5.0249043846692105]
longitudes_of_ascending_node = [2.0064752362645937, 3.994763020982937, 4.786445201421792]
tend = 1.000000E+08
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [3.077097805465111, 15.839702359488674, 5.517510578935515, 2.3828930195118176]
metallicities = [0.01630303006890956, 0.01630303006890956, 0.01630303006890956, 0.01630303006890956]
semimajor_axes = [0.22431991892662603, 31.90022634843763, 6576.536305682837]
eccentricities = [0.7430836382128856, 0.5204313220275024, 0.9289400466733657]
inclinations = [1.7941748834200397, 1.1646898681822524, 1.2992456528220557]
arguments_of_pericentre = [5.987315439449751, 5.74661821195682, 0.9299369122614957]
longitudes_of_ascending_node = [0.5383802291327379, 5.087132778096285, 5.68592590161848]
tend = 1.000000E+08
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

ODE_tides.cpp -- FATAL ERROR: t_V is NaN
binary_index 4 star_index 0 companion_index 1
t_V -nan 
st 1
pr 1
M 3.1982
star->convective_envelope_mass 1e-10
m 15.6046
a -0.0556011469472139
e 1.76966042663074
rp 0.042794
R 0.00936525
star->convective_envelope_radius 4.64913e-13
star->luminosity 0.0292151
k_AM 0.00590012
=============================
Printing system; N=7; integration_flag=0 size = 7
index 0 -- body -- parent 4 m 3.198203953142868 r 0.00936525 st 1 mc 0 minit 15.6046 menv 1e-10 epoch 5.9257e+06 age 274305 rc 0 renv 4.64913e-13 lum 0.0292151 Spin_freq 12429.8 Omega_crit 12398.1
index 1 -- body -- parent 4 m 15.604590298397904 r 0.0362814 st 1 mc 0 minit 15.6046 menv 1e-10 epoch -1.93714e+06 age 8.13714e+06 rc 0 renv 4.64913e-13 lum 9.84024 Spin_freq 1834.33 Omega_crit 3591.54
index 2 -- body -- parent 5 m 5.517510578935515 r 0.0131999 st 1 mc 0 minit 5.51751 menv 1e-10 epoch -1.81365e+06 age 8.01366e+06 rc 0 renv 4.64913e-13 lum 0.217252 Spin_freq 3948.53 Omega_crit 9731.83
index 3 -- body -- parent 5 m 2.382893019511818 r 0.00795339 st 1 mc 0 minit 2.38289 menv 1e-10 epoch -1.81365e+06 age 8.01366e+06 rc 0 renv 4.64913e-13 lum 0.0093551 Spin_freq 4389.84 Omega_crit 13674.2
index 4 -- binary -- parent 6 child1 0 child2 1 m 18.8028 a -0.0556011 e 1.76966 rp 0.042794 
index 5 -- binary -- parent 6 child1 2 child2 3 m 7.9004 a 31.9015 e 0.512455 rp 15.5534 
index 6 -- binary -- parent -1 child1 4 child2 5 m 26.7032 a 6557.74 e 0.929059 rp 465.212

configuration = "[2,2]"
N_bodies = 4
masses = [14.224950339314399, 0.8897139647426828, 13.94306058151649, 25.529234346782907]
metallicities = [0.016552486981039515, 0.016552486981039515, 0.016552486981039515, 0.016552486981039515]
semimajor_axes = [321.64856782898676, 13.781481292380843, 6610.827792364814]
eccentricities = [0.8341440727777816, 0.5671605191088954, 0.126947856450332]
inclinations = [1.2502977473560803, 2.219627718811098, 2.428346067703022]
arguments_of_pericentre = [5.854524108992268, 6.282864677583723, 1.6957614719421723]
longitudes_of_ascending_node = [1.8960677331378781, 1.0366560788684072, 2.8757182594407933]
tend = 1.000000E+08
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [10.605765840906358, 12.808546742652963, 1.5973515434842434, 8.893001415954778]
metallicities = [0.014420010133190024, 0.014420010133190024, 0.014420010133190024, 0.014420010133190024]
semimajor_axes = [1.5036199933179712, 0.5811001991526038, 291.46590663564245]
eccentricities = [0.4438059792712014, 0.46105663564734384, 0.4221648585910908]
inclinations = [1.7748346856690747, 0.40983594087210423, 1.8404934917875604]
arguments_of_pericentre = [1.8700278727270907, 0.43102394164651825, 2.2149958071719147]
longitudes_of_ascending_node = [1.4589307864827417, 4.793613840703381, 0.7013812935222983]
tend = 1.000000E+08
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

binary_evolution.cpp -- handle_mass_transfer_cases -- ERROR: t_MT = nan; a 1.51987 e 0.46048 R_Lc 0.599562 x 1.8535 E_0 nan 

configuration = "[2,2]"
N_bodies = 4
masses = [24.028002450633135, 5.424563409803646, 2.761589046240006, 12.75282280725989]
metallicities = [0.013879815952003946, 0.013879815952003946, 0.013879815952003946, 0.013879815952003946]
semimajor_axes = [25.422147872812605, 3.0990529280790575, 2054.9262925007934]
eccentricities = [0.727431546566552, 0.5610129004052001, 0.8752639158022109]
inclinations = [0.915042492257424, 0.08114248791422286, 0.8029112225963015]
arguments_of_pericentre = [2.6041072099556115, 0.8010991006190188, 5.281902603019346]
longitudes_of_ascending_node = [0.37547007211639055, 2.2008189035172423, 5.778884167237088]
tend = 1.000000E+08
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [4.24275623639273, 5.0571180390992, 25.798725559974713, 62.622689057887115]
metallicities = [0.016382018773432913, 0.016382018773432913, 0.016382018773432913, 0.016382018773432913]
semimajor_axes = [24.085555880866927, 5.830229657929401, 3394.4411880980797]
eccentricities = [0.8367931310467791, 0.9033769390241188, 0.7861739694189925]
inclinations = [0.404309096333188, 1.2617307012094159, 1.4292098376420659]
arguments_of_pericentre = [1.4121119633443193, 2.818529409637241, 3.885826308846624]
longitudes_of_ascending_node = [2.0623091725410294, 3.352629275054458, 0.3810936148264323]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [25.472672442547363, 29.840431054608715, 1.578442917811624, 0.1701283229267224]
metallicities = [0.027177148943916173, 0.027177148943916173, 0.027177148943916173, 0.027177148943916173]
semimajor_axes = [73.96419622249603, 0.5468550519414457, 331.69787410952017]
eccentricities = [0.9234048754427818, 0.8599075001749081, 0.26298598653151656]
inclinations = [2.84029791686314, 1.6496060566125157, 0.5987913752691788]
arguments_of_pericentre = [4.538739296301567, 6.27932481197671, 3.9414333484223065]
longitudes_of_ascending_node = [5.543214461262542, 5.102103794424898, 3.384179592192939]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [0.9858421344726858, 2.2197813915714204, 32.624086808030626, 12.941314253710717]
metallicities = [0.02045373811168327, 0.02045373811168327, 0.02045373811168327, 0.02045373811168327]
semimajor_axes = [0.37311328592805276, 0.10395244202071213, 73.67556982848069]
eccentricities = [0.6440307309975877, 0.4457454703638011, 0.9024697950467355]
inclinations = [1.6153302307185153, 1.1938660511146522, 2.7900756114919387]
arguments_of_pericentre = [1.8797746705934975, 0.504944076085212, 1.881296488666065]
longitudes_of_ascending_node = [3.626528823982504, 3.1835675943716324, 5.849273516037526]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

stellar_evolution.cpp -- evolve_stars -- ERROR: quantity spin_vec[0] is NaN
stellar_evolution.cpp -- evolve_stars -- ERROR: quantity spin_vec[1] is NaN
stellar_evolution.cpp -- evolve_stars -- ERROR: quantity spin_vec[2] is NaN
stellar_evolution.cpp -- evolve_stars -- ERROR: quantity spin_vec[0] is NaN
stellar_evolution.cpp -- evolve_stars -- ERROR: quantity spin_vec[1] is NaN
stellar_evolution.cpp -- evolve_stars -- ERROR: quantity spin_vec[2] is NaN

configuration = "2"
N_bodies = 2
masses = [1.4033902660164514, 1.4497387409123248]
metallicities = [0.01083138031566544, 0.01083138031566544]
semimajor_axes = [1.3374445598964206]
eccentricities = [0.3174495897133802]
inclinations = [1.1381511990978366]
arguments_of_pericentre = [0.2904225417557361]
longitudes_of_ascending_node = [2.1074651848995956]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1]
object_types = [1, 1]

configuration = "2"
N_bodies = 2
masses = [1.332040740163885, 1.3068911975703883]
metallicities = [0.023009827126615208, 0.023009827126615208]
semimajor_axes = [0.10297972574554977]
eccentricities = [0.5102998299876382]
inclinations = [0.8405107912056219]
arguments_of_pericentre = [5.970126102009413]
longitudes_of_ascending_node = [1.6207052542414224]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1]
object_types = [1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [2.638583237816104, 11.959695596521696, 6.358649581584067, 1.7713076776714505]
metallicities = [0.01936543324410305, 0.01936543324410305, 0.01936543324410305, 0.01936543324410305]
semimajor_axes = [0.4814231000939026, 0.17852091882440735, 9.751992806474941]
eccentricities = [0.18575870217165172, 0.3469796694644274, 0.6486515558764904]
inclinations = [1.1094613994856506, 1.9612175944915886, 1.2983684023065827]
arguments_of_pericentre = [0.9243483647991751, 6.0020126733204195, 2.0123820399751877]
longitudes_of_ascending_node = [4.511247208668151, 4.796830595614424, 2.781652923782494]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [9.330796046420812, 12.31874606447617, 16.34011180640447, 13.733130727748412]
metallicities = [0.011882243485479766, 0.011882243485479766, 0.011882243485479766, 0.011882243485479766]
semimajor_axes = [1.4629582700386954, 1.4271541401366976, 25.780532498621106]
eccentricities = [0.566902289493805, 0.5943991368249637, 0.08250802769876078]
inclinations = [0.9443444795543554, 1.3009692248388145, 1.0869803781110217]
arguments_of_pericentre = [0.06487086198048052, 5.471328815619844, 5.03384927758664]
longitudes_of_ascending_node = [2.589633119365464, 4.3262187011442235, 1.0469814067237144]
tend = 1.000000E+08
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

ODE_tides.cpp -- FATAL ERROR: t_V is NaN
binary_index 4 star_index 2 companion_index 3
t_V -nan 
st 4
pr 1
M 15.9889
star->convective_envelope_mass 1e-10
m 13.618
a -11.454118012798
spin -41.7167 -3.24768 22.0532 norm 47.2988
e 1.10257229779413
rp 1.17488
R 0.655889
rg 0.08
star->convective_envelope_radius 4.64913e-13
star->luminosity 23.2044
k_AM 0.0366236

configuration = "[2,2]"
N_bodies = 4
masses = [7.692620887411417, 42.76798396758461, 40.64553793350696, 12.057569043812478]
metallicities = [0.02649954172569116, 0.012814905599314315, 0.002473080642553742, 0.007022191275437818]
semimajor_axes = [5.050106477176275, 144.261944580885, 1215.5796316956648]
eccentricities =  [0.23987555486249884, 0.49304674164813916, 0.11993551412336234, 0.26816414144104966]
inclinations = [2.643398245158917, 1.219303456138767, 1.6285303742502966]
longitudes_of_ascending_node = [0.31748334346942136, 4.233797061352631, 2.474914893060231]
arguments_of_pericentre = [2.643398245158917, 1.219303456138767, 1.6285303742502966]
tend = 100000000.0
N_steps = 5000
stellar_types = []
object_types=[]

configuration = "[2,2]"
N_bodies = 4
masses = [5.093058447760394, 12.648237722921836, 6.443580107515282, 12.342633840803302]
metallicities = [0.028430542650177776, 0.028430542650177776, 0.028430542650177776, 0.028430542650177776]
semimajor_axes = [0.12401049479541143, 0.3497077963137055, 3.247264834545161]
eccentricities = [0.47201203636119604, 0.6515837152119581, 0.47068967710714166]
inclinations = [3.0839014161281986, 2.088592665666578, 0.3044507497210534]
arguments_of_pericentre = [1.8014437875856235, 3.704441107488157, 3.336967919674607]
longitudes_of_ascending_node = [4.133943340139425, 0.790055599595232, 0.5130315279982524]
tend = 1.000000E+08
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

ODE_tides.cpp -- FATAL ERROR: t_V is NaN
binary_index 5 star_index 2 companion_index 3
t_V -nan 
st 1
pr 1
M 6.44358
star->convective_envelope_mass 1e-10
m 12.3426
a -1.15259944261353
spin 2289.98 -2031.26 -1575.68 norm 3442.79
e 1.01358019277815
rp 0.0156525
R 0.0149387
rg 0.08
star->convective_envelope_radius 4.64913e-13
star->luminosity 0.349338
k_AM 0.00928437
=============================
Printing system; N=7; integration_flag=0 size = 7
index 0 -- body -- parent 4 m 5.093058448107017 r 0.0130446 st 1 mc 0 minit 5.09306 menv 1e-10 epoch 0 age 99140.3 rc 0 renv 4.64913e-13 lum 0.151223 Spin_freq 3964.5 Omega_crit 9517.53
index 1 -- body -- parent 4 m 12.648237709082121 r 0.0222775 st 1 mc 0 minit 12.6476 menv 1e-10 epoch -8.82727e-05 age 99140.3 rc 0 renv 4.64913e-13 lum 3.2361 Spin_freq 2130.09 Omega_crit 6720.42
index 2 -- body -- parent 5 m 6.443580107669352 r 0.0149387 st 1 mc 0 minit 6.44358 menv 1e-10 epoch 0 age 99140.3 rc 0 renv 4.64913e-13 lum 0.349338 Spin_freq 3442.79 Omega_crit 8735.2
index 3 -- body -- parent 5 m 12.342633828396878 r 0.0219502 st 1 mc 0 minit 12.342 menv 1e-10 epoch -8.23321e-05 age 99140.3 rc 0 renv 4.64913e-13 lum 3.00086 Spin_freq 2169.85 Omega_crit 6787.76
index 4 -- binary -- parent 6 child1 0 child2 1 m 17.7413 a 0.124757 e 0.479357 rp 0.0649538 
index 5 -- binary -- parent 6 child1 2 child2 3 m 18.7862 a -1.1526 e 1.01358 rp 0.0156525 
index 6 -- binary -- parent -1 child1 4 child2 5 m 36.5275 a 3.25139 e 0.47049 rp 1.72164

configuration = "[2,2]"
N_bodies = 4
masses = [2.637289392039294, 30.422508761203126, 18.524433804979367, 21.904410732825582]
metallicities = [0.0287060880866142, 0.0287060880866142, 0.0287060880866142, 0.0287060880866142]
semimajor_axes = [0.11679976142288281, 17.002803212706528, 1925.0810470138933]
eccentricities = [0.9001032171589015, 0.6784999547217095, 0.8617681964238854]
inclinations = [1.7255463919491807, 2.4556618697569945, 0.693026467836522]
arguments_of_pericentre = [1.7802895541513786, 0.8728659577824223, 5.644177993305278]
longitudes_of_ascending_node = [1.7266035085075135, 5.927816691601421, 2.683845230261317]
tend = 1.000000E+08
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

ODE_tides.cpp -- FATAL ERROR: t_V is NaN
binary_index 4 star_index 2 companion_index 3
t_V -nan 
st 1
pr 1
M 24.764
star->convective_envelope_mass 1e-10
m 12.293
a -1730.86043128734
spin -985.4 -845.446 -1016.02 norm 1648.66
e 1.00453245734712
rp 7.84505
R 0.0539933
rg 0.08
star->convective_envelope_radius 4.64913e-13
star->luminosity 33.3296
k_AM 0.00645829
=============================
Printing system; N=4; integration_flag=0 size = 4
index 0 -- body -- parent -1 m 8.521199150612446 r 1.67973e-07 st 14 mc 8.5212 minit 9.74572 menv 1e-10 epoch 5.94826e+06 age 1.77814e+06 rc 1.67973e-07 renv 4.64913e-13 lum 2.7104e-14 Spin_freq 2e+08 Omega_crit 2.66424e+11
index 2 -- body -- parent 4 m 24.764019992532162 r 0.0539933 st 1 mc 0 minit 24.8727 menv 1e-10 epoch 3.59554e+06 age 4.13086e+06 rc 0 renv 4.64913e-13 lum 33.3296 Spin_freq 1648.66 Omega_crit 2492.19
index 3 -- body -- parent 4 m 12.292965233098277 r 6.00564 st 4 mc 6.15582 minit 20.6449 menv 3.24701 epoch -331599 age 8.058e+06 rc 0.00335852 renv 5.0921 lum 39.2977 Spin_freq 0.000745338 Omega_crit 1.49682
index 4 -- binary -- parent -1 child1 2 child2 3 m 37.057 a -1730.86 e 1.00453 rp 7.84505

configuration = "[2,2]"
N_bodies = 4
masses = [5.7460828788155265, 8.357801422681348, 1.4648323385777384, 0.975317498667244]
metallicities = [0.02997867474838753, 0.02997867474838753, 0.02997867474838753, 0.02997867474838753]
semimajor_axes = [0.2818201949248253, 0.12278006722826543, 10.450979239812405]
eccentricities = [0.23547098492944143, 0.5780702156414276, 0.5829107575734138]
inclinations = [1.0951374930164222, 1.9667019164797899, 2.269031514165921]
arguments_of_pericentre = [2.959631886126624, 2.7891867643123223, 2.5693711806875448]
longitudes_of_ascending_node = [2.073510581755221, 4.465879312867095, 3.865192290024434]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

stellar_evolution.cpp -- compute_spin_angular_momentum_from_spin_frequency -- ERROR: quantity S is NaN
nbody_evolution.cpp -- create_mstar_instance_of_system -- R->gbs_tolerance 1e-10 R->stopping_condition_tolerance 1e-10 R->output_time_tolerance 1e-06
nbody_evolution.cpp -- analyze_mstar_system -- starting integration
U not finite.

configuration = "[2,2]"
N_bodies = 4
masses = [1.0053773093746776, 41.44700146086816, 0.23477380962735328, 33.8908268398115]
metallicities = [0.008100239195764945, 0.022055820663677846, 0.028865656353523145, 0.00746259430559874]
semimajor_axes = [115.23146688356738, 118.4083862543678, 1144.5038115817467]
eccentricities = [0.11154081632030916, 0.4763745057584925, 0.22356268930881368, 0.4232043362355639]
inclinations = [2.197478952675784, 0.9344257397126362, 2.5566212518846934]
arguments_of_pericentre = [2.4913190451021308, 5.536134662197809, 3.652245172807698]
longitudes_of_ascending_node = [5.540106670427087, 4.351304311534335, 4.556907034931878]
tend = 1.000000E+08
N_steps = 5000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]
plot_filename = "results/test"

configuration = "[2,2]"
N_bodies = 4
masses = [8.932593431368586, 9.61377558081899, 10.389622413272619, 11.381783486114147]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [2.9955454269732003, 2.565151620608462, 536.6209259454847]
eccentricities = [0.3366875609652914, 0.1713143428375245, 0.7263038653127867]
inclinations = [2.4586523682892047, 1.5943050547546478, 0.6791199653107841]
arguments_of_pericentre = [3.923983927349352, 4.340222121054781, 3.832341057029352]
longitudes_of_ascending_node = [0.284556890984554, 3.4728930794221284, 2.3879167043330285]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [11.936308470789099, 7.547039803563019, 10.485550503630007, 14.9226023204833]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [0.7373304083674703, 19.468351016464595, 301.10897805970336]
eccentricities = [0.07776443321165889, 0.6847039014680483, 0.6227576939518841]
inclinations = [0.444389012062038, 1.725357029474652, 2.620807358435491]
arguments_of_pericentre = [3.523291155163777, 4.4182808111681044, 0.4903957377732229]
longitudes_of_ascending_node = [1.3104762524702338, 3.235687730517438, 4.081576882875382]
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,kick_distribution_sigma_km_s_BH = 0.0,0.0


configuration = "[2,2]"
N_bodies = 4
masses = [1.6337969506029681, 4.2750927137387515, 11.877354130165823, 10.288911109538986]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [24.278312077726156, 15.8970375523649, 1793.9673042722384]
eccentricities = [0.5521380989079157, 0.041501194627405105, 0.03098032451290733]
inclinations = [2.744543967647522, 0.8542249583240987, 2.114534648076657]
arguments_of_pericentre = [1.5503620961209863, 1.3694357099696468, 3.042271831766447]
longitudes_of_ascending_node = [4.519472526396828, 2.0803372971183856, 2.2168615570984365]
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,kick_distribution_sigma_km_s_BH = 0.0,0.0


configuration = "[2,2]"
N_bodies = 4
masses = [10.155717437998788, 10.694291523444676, 9.377971441266535, 9.393160984603881]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [326.68621136397473, 55.380787636848346, 9718.700044479629]
eccentricities = [0.7676894566629486, 0.7164696698248222, 0.07778716556753808]
inclinations = [0.8839497330900175, 2.097147399211663, 2.6356261816685924]
arguments_of_pericentre = [3.561688172958032, 4.457275773557397, 4.630435476178859]
longitudes_of_ascending_node = [1.3847490327754812, 2.2650396692046706, 1.9912735600647105]
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,kick_distribution_sigma_km_s_BH = 0.0,0.0

WARNING -- mse.py -- Internal error with code  4 occurred -- -- stopping the simulation but saving data/making plots if specified in command line arguments.

index_system = 34800
error_code = 39
configuration = "2"
N_bodies = 2
masses = [5.772395761720277, 0.5785211527442907]
metallicities = [0.02, 0.02]
semimajor_axes = [0.1010678585666266]
eccentricities = [0.43458777309961805]
inclinations = [1.0153103975311522]
arguments_of_pericentre = [2.942442673342705]
longitudes_of_ascending_node = [4.122304273623329]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1]
object_types = [1, 1]

common_envelope_evolution.cpp -- binary_common_envelope_evolution -- no core present (MC1 = 0)! M1 = 5.75362; skipping



index_system = 27424
error_code = 39
configuration = "2"
N_bodies = 2
masses = [4.907418150967316, 0.6481168394048793]
metallicities = [0.02, 0.02]
semimajor_axes = [0.35195824001524556]
eccentricities = [0.8528390249574894]
inclinations = [1.6787408441193545]
arguments_of_pericentre = [4.199371241058291]
longitudes_of_ascending_node = [3.4766074853019773]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1]
object_types = [1, 1]

common_envelope_evolution.cpp -- binary_common_envelope_evolution -- no core present (MC1 = 0)! M1 = 4.85786; skipping



index_system = 92268
error_code = 39
configuration = "2"
N_bodies = 2
masses = [4.307161203158041, 0.7377170266034233]
metallicities = [0.02, 0.02]
semimajor_axes = [0.18065798231294078]
eccentricities = [0.7158621563074218]
inclinations = [2.855349660347169]
arguments_of_pericentre = [0.5221935033068231]
longitudes_of_ascending_node = [5.508429160847378]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1]
object_types = [1, 1]

common_envelope_evolution.cpp -- binary_common_envelope_evolution -- no core present (MC1 = 0)! M1 = 4.2377; skipping

configuration = "[2,2]"
N_bodies = 4
masses = [3.8185810756007745, 8.961776206419307, 8.875922471380392, 0.0818910846179762]
metallicities = [0.019090329816641598, 0.019090329816641598, 0.019090329816641598, 0.019090329816641598]
semimajor_axes = [111.82872269402648, 37.83261352955516, 3634.403287332339]
eccentricities = [0.4883658209526968, 0.1948852647292638, 0.43906171089340085]
inclinations = [0.676684821448654, 1.5698512422448085, 1.9830045124192537]
arguments_of_pericentre = [4.8665607051995075, 1.1664849935341755, 0.7689294976218852]
longitudes_of_ascending_node = [4.383333555209011, 5.71251976775718, 1.7944128191722075]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

/opt/python-3.7.2/lib/python3.7/site-packages/matplotlib/collections.py:874: RuntimeWarning: invalid value encountered in sqrt
  scale = np.sqrt(self._sizes) * dpi / 72.0 * self._factor 

configuration = "[2,2]"
N_bodies = 4
masses = [3.548986035679437, 17.021172984446476, 8.12233857680884, 4.428527949432436]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [0.2015380813653181, 5.604271111382169, 53.586871789366896]
eccentricities = [0.3527075005652079, 0.9102901776980962, 0.47208484125535405]
inclinations = [1.3308787041010495, 1.9977131616727142, 1.6271690533278615]
arguments_of_pericentre = [6.227825330298463, 4.316756296360188, 4.147507183112485]
longitudes_of_ascending_node = [0.9923061086882445, 3.28162665237633, 0.7969508437343856]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [9.135347337893288, 4.1538867387257685, 0.8912380511169684, 6.179700070165928]
metallicities = [0.02854221674403863, 0.02854221674403863, 0.02854221674403863, 0.02854221674403863]
semimajor_axes = [0.46692544582780987, 25.21133488196536, 641.2916952582702]
eccentricities = [0.900214408003911, 0.9150568422407039, 0.31807705119021806]
inclinations = [1.4896174214872016, 2.1857845559212397, 2.272687154058045]
arguments_of_pericentre = [5.438711084842693, 2.318957850426181, 0.6664173295467685]
longitudes_of_ascending_node = [1.3713311698453714, 2.393325291426528, 4.08187847634804]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "[2,2]"
N_bodies = 4
masses = [15.382118687282885, 26.572723273311784, 11.027003732710213, 34.323766359074185]
metallicities = [0.011397589597065867, 0.011397589597065867, 0.011397589597065867, 0.011397589597065867]
semimajor_axes = [0.3302039618418993, 0.21560824114268048, 56.40763781359777]
eccentricities = [0.5263312638217957, 0.3874900116035225, 0.36246057783040386]
inclinations = [0.8246393085598901, 2.8401141135767025, 1.5802817679541676]
arguments_of_pericentre = [0.4570475683556849, 0.917048455620907, 2.9235281959997717]
longitudes_of_ascending_node = [1.2999779467170967, 4.00047047922905, 2.381220676836103]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

index_system = 5068
configuration = "[2,2]"
N_bodies = 4
masses = [11.480709507969538, 3.6457095442634584, 7.160267781685773, 34.51804481639648]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [1.9077232939054682, 0.29143687635844623, 39.58910363532721]
eccentricities = [0.9244502029730212, 0.3037430079738928, 0.10758301005433966]
inclinations = [0.5408156367311555, 2.0549600465725244, 1.2114510361461213]
arguments_of_pericentre = [5.0265337306641475, 4.832210273222782, 0.664336346748185]
longitudes_of_ascending_node = [4.193160664920405, 5.666932986590698, 1.2585027759919258]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]


index_system = 42187
configuration = "[2,2]"
N_bodies = 4
masses = [17.955099856415227, 25.43861143816699, 70.44729020382495, 70.54221254128075]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [2.2720248981610838, 3.4103402126941993, 30.21357327009001]
eccentricities = [0.7710879979782606, 0.7846935884744111, 0.35153389358305764]
inclinations = [1.5892470192609987, 1.5349852150676497, 1.9012925213102971]
arguments_of_pericentre = [2.6606954235213, 4.887883871569553, 3.689148814395094]
longitudes_of_ascending_node = [5.710763958698144, 2.0541324114926, 5.296973911249262]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]


index_system = 63284
configuration = "[2,2]"
N_bodies = 4
masses = [33.10094308475692, 5.529951813777178, 65.87393448664203, 39.860439925592296]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [5.019474555774855, 0.4384528664467132, 72.94295687449474]
eccentricities = [0.8062429673991279, 0.5007754835232592, 0.4880394392603675]
inclinations = [1.356931270359145, 2.0678048806045957, 1.3406356779928115]
arguments_of_pericentre = [3.247602606101761, 0.6490004544830226, 0.7808334231283096]
longitudes_of_ascending_node = [6.244496779767114, 4.638393053125136, 1.4325320669797317]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

configuration = "2"
N_bodies = 2
masses = [4.039663119247391, 9.528885714625995]
metallicities = [0.02, 0.02]
semimajor_axes = [14.558792717733997]
eccentricities = [0.7585923334122942]
inclinations = [1.2772339580887269]
arguments_of_pericentre = [4.410128770679282]
longitudes_of_ascending_node = [4.546122970738234]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1]
object_types = [1, 1]

tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity star->dspin_vec_dt[i] is infinite
==================================================
WARNING -- mse.py -- Internal error with code  1 occurred -- -- stopping the simulation but saving data/making plots if specified in command line arguments.

configuration = "[2,2]"
N_bodies = 4
masses = [5.061273813723217, 1.7892620395569252, 4.492164764857658, 9.312305062530513]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [7.361062303816427, 254.97604032440273, 1217.126882689642]
eccentricities = [0.6851465744813263, 0.4994877639138737, 0.2868688940630904]
inclinations = [1.8893561079033927, 0.6103322171552832, 1.197393715083195]
arguments_of_pericentre = [5.667626286819018, 4.429486502961903, 2.264566425049202]
longitudes_of_ascending_node = [4.776006925094928, 1.123648903947437, 4.748757830311203]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

tides.cpp -- compute_EOM_equilibrium_tide_BO_full -- ERROR: quantity star->dspin_vec_dt[i] is infinite
==================================================
WARNING -- mse.py -- Internal error with code  1 occurred -- -- stopping the simulation but saving data/making plots if specified in command line arguments.