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

stellar_evolution.cpp -- compute_spin_angular_momentum_from_spin_frequency -- ERROR: quantity S is NaN #11

Closed pavanvyn closed 3 years ago

pavanvyn commented 3 years ago

Version 816bfbf151389ae330ba9315f2e796d1c0c88e97

configuration = "[2,2]"
N_bodies = 4
masses = [62.159704663171716, 59.774910076680904, 31.417961828319992, 2.03289118796106]
metallicities = [0.02493711162848645, 0.02493711162848645, 0.02493711162848645, 0.02493711162848645]
semimajor_axes = [5.0602172728575505, 8.333273949324111, 1225.6916909975089]
eccentricities = [0.4623847217492819, 0.732755835950576, 0.7914329440080197]
inclinations = [2.324101373467737, 0.9582927328595943, 1.3051385381800236]
arguments_of_pericentre = [3.895245648188573, 4.459100942021893, 0.1100232142371]
longitudes_of_ascending_node = [1.985013688142783, 1.000873058513418, 4.103131684156803]
tend = 1.000000E+08
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

This system gives an error - stellar_evolution.cpp -- compute_spin_angular_momentum_from_spin_frequency -- ERROR: quantity S is NaN

hamers commented 3 years ago

The changes that were made to resolve Issue #10 should have also resolved this issue (new version baa6d17c57b9cef377a911e58f6d93312e864423).