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

U not finite #2 #29

Closed pavanvyn closed 3 years ago

pavanvyn commented 3 years ago

Version b1aa18219634329d0c3092c5a9b7d9487ed81f89

configuration = "[2,2]"
N_bodies = 4
masses = [5.7460828788155265, 8.357801422681348, 1.4648323385777384, 0.975317498667244]
metallicities = [0.02997867474838753, 0.02997867474838753, 0.02997867474838753, 0.02997867474838753]
semimajor_axes = [0.2818201949248253, 0.12278006722826543, 10.450979239812405]
eccentricities = [0.23547098492944143, 0.5780702156414276, 0.5829107575734138]
inclinations = [1.0951374930164222, 1.9667019164797899, 2.269031514165921]
arguments_of_pericentre = [2.959631886126624, 2.7891867643123223, 2.5693711806875448]
longitudes_of_ascending_node = [2.073510581755221, 4.465879312867095, 3.865192290024434]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

Error:

stellar_evolution.cpp -- compute_spin_angular_momentum_from_spin_frequency -- ERROR: quantity S is NaN
nbody_evolution.cpp -- create_mstar_instance_of_system -- R->gbs_tolerance 1e-10 R->stopping_condition_tolerance 1e-10 R->output_time_tolerance 1e-06
nbody_evolution.cpp -- analyze_mstar_system -- starting integration
U not finite.

The program does not return control to the terminal, but stalls after this point.

hamers commented 3 years ago

The problem was related to invalid spins being assigned in collision.cpp as a result of formally unbound initial collision/CE orbits that could occur in N-body mode. Should be fixed in 768c69b.