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

Closed pavanvyn closed 3 years ago

pavanvyn commented 3 years ago

Version bbfb219c4c121dc1a2498728fa6b1f32a04763de

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]

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
mse.py -- Internal error with code  1 occurred -- stopping simulation.
hamers commented 3 years ago

There was a problem with mass accreting naked He GB stars, causing SSE to produce NaNs that propagated further. Should be resolved in 93c1a62.