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. #2 #25

Closed pavanvyn closed 3 years ago

pavanvyn commented 3 years ago

Version 93c1a62cf452ccacb9eb6f980c5a376a755c532c

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]

Error

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.
hamers commented 3 years ago

SSE for some reason turns a He star into a CHeB star with zero core mass, causing problems for CE evolution; workaround has been added in 8f1a333.