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

mst.c -- ERROR: U not finite. #40

Closed pavanvyn closed 3 years ago

pavanvyn commented 3 years ago

Version 3ae9006280819f1bc190ce06e5191bb669a90876

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

There was a bug in the triple CE routine where the age of the donor star was not properly passed onto hrdiag(), leading to NaNs. Fixed in 8c5e2fb.