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

ODE_tides.cpp -- FATAL ERROR #3 #15

Closed pavanvyn closed 3 years ago

pavanvyn commented 3 years ago

Version f4f9bf2b2bd52ece4af30ec149ade1a40f78585e

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]

The full error -

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

This should now have been resolved in 4b4ba38.