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

U not finite #17

Closed pavanvyn closed 3 years ago

pavanvyn commented 3 years ago

Version 4b4ba38ba77004b458c9d76f61eecb7710450eae

configuration = "[2,2]"
N_bodies = 4
masses = [7.540948688753799, 5.542806479907622, 1.6449369158111955, 8.523708895122976]
metallicities = [0.021331238724547827, 0.021331238724547827, 0.021331238724547827, 0.021331238724547827]
semimajor_axes = [193.12588279845278, 18.119725906827203, 4256.999687317027]
eccentricities = [0.7524873566774626, 0.15263088650184092, 0.3869766276151553]
inclinations = [0.8612693368583286, 2.5342304447719104, 1.048379593931638]
arguments_of_pericentre = [0.5027362011630949, 3.6066003077451625, 1.0431612857126726]
longitudes_of_ascending_node = [6.264325417514426, 4.867578242273096, 3.7973749864625987]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

Errror : U not finite

Edit : Another system with the same error

configuration = "[2,2]"
N_bodies = 4
masses = [3.762254469849947, 7.349744147945759, 8.579556659374338, 47.04141214391311]
metallicities = [0.02728435020513672, 0.02728435020513672, 0.02728435020513672, 0.02728435020513672]
semimajor_axes = [6.20045452202843, 4.1103742181516845, 277.63581001067365]
eccentricities = [0.8265914585781401, 0.7442378738496561, 0.2934071857082655]
inclinations = [1.8436231207845388, 2.9116722325910405, 2.5066085538633645]
arguments_of_pericentre = [4.576486222775116, 1.7128248807379545, 1.0029244850244559]
longitudes_of_ascending_node = [6.160384217581603, 5.154441680206168, 4.286053523206743]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]
hamers commented 3 years ago

This was due to a bug in MSTAR; it has been resolved in a4289e0d257ec808437c5e1566ca08c707caa70b.