stellar_evolution.cpp -- evolve_stars -- ERROR: quantity sse_main_sequence_timescale is infinite

[pavanvyn]

Version 01880c7ea445664e127b775eb6adbf25a8c2615e

System:

configuration = "[2,2]"
N_bodies = 4
masses = [47.643992189057975, 50.445366858771266, 45.25627339750874, 17.62837673193206]
metallicities = [0.012862801475893733, 0.012862801475893733, 0.012862801475893733, 0.012862801475893733]
semimajor_axes = [0.6689433925633086, 0.28122681885568385, 3867.268837446437]
eccentricities = [0.46369985674959907, 0.1768538859917899, 0.6226246106696653]
inclinations = [0.9837918622467343, 2.0849834819627473, 1.5386344118860225]
arguments_of_pericentre = [6.050730734464304, 1.5180479341436717, 0.5739194600286861]
longitudes_of_ascending_node = [3.1026518087318626, 2.3138754908474626, 2.1752334688954122]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

Error:

stellar_evolution.cpp -- evolve_stars -- ERROR: quantity sse_main_sequence_timescale is infinite
stellar_evolution.cpp -- evolve_stars -- ERROR: quantity spin_vec[0] is NaN
stellar_evolution.cpp -- evolve_stars -- ERROR: quantity spin_vec[1] is NaN
stellar_evolution.cpp -- evolve_stars -- ERROR: quantity spin_vec[2] is NaN
mse.py -- Internal error with code  1 occurred -- stopping simulation.

I have not been encountering many more errors (other than 'Killed' terminations) recently. Hopefully, this is among the final ones.

[hamers]

The problem arises from a problem in CE evolution, which could be called after contact evolution when the donor was not a giant. This has been fixed in 00eb2ca (now collision is invoked in the case of contact and when neither donor and accretor are giants).