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

Possible case of stalling #6

Closed pavanvyn closed 3 years ago

pavanvyn commented 3 years ago

I don't know if this is a problem, but some configurations run for a lot of time, with no end in sight. One such set of parameters is this - configuration = "[2,2]" N_bodies = 4 masses = [5.678484806683704, 17.734276507507087, 0.3077184448926684, 16.96659279614795] metallicities = [0.01568104521664872, 0.025068510046561455, 0.006705614803401992, 0.01546545088235576] semimajor_axes = [193.90635623544225, 193.16896749992142, 997.2798950542823] eccentricities = [0.3696338795586719, 0.10695168146863654, 0.24230746455432545, 0.06681659983986149] inclinations = [0.3442294115191764, 1.5444893364050116, 1.7832183025348438] longitudes_of_ascending_node = [1.3927608569078445, 6.0273416705134615, 1.194482337762545] arguments_of_pericentre = [0.3442294115191764, 1.5444893364050116, 1.7832183025348438] tend = 1.0e8 N_steps = 5000 stellar_types = [1,1,1,1] object_types=[]

hamers commented 3 years ago

In this particular system the semisecular regime was entered, and the problem was that the subsequent time steps imposed on MSTAR were very short. This should now have been resolved (ae2a5567aba2224531671d6bb0214e2ec5008617).