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. #5 #36

Closed pavanvyn closed 3 years ago

pavanvyn commented 3 years ago

Version 9aafcda90c2016b50bc646a3157045f5fb620901

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

The function compute_eccentric_anomaly_from_mean_anomaly() would not converge in this case, finally giving a NaN; this is resolved in 7a0b943 by checking for non-convergence and in the latter case simply setting the sampled eccentric anomaly equal to mean anomaly