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

ZeroDivisionError: 0.0 cannot be raised to a negative power (no error code or hard exit) #22

Closed pavanvyn closed 3 years ago

pavanvyn commented 3 years ago

Version 0cf1cab2b0e0ca163f94a5f0ec63bfc5aa04bbd5

configuration = "[2,2]"
N_bodies = 4
masses = [7.755298110394679, 13.02455302992831, 1.7306612927312375, 1.7045592918749755]
metallicities = [0.027715890415101217, 0.027715890415101217, 0.027715890415101217, 0.027715890415101217]
semimajor_axes = [0.9629082312726244, 3.162434677535477, 1800.1953209157884]
eccentricities = [0.33701500691806957, 0.5524906988041355, 0.7528957425541132]
inclinations = [0.7291109296646838, 1.7463778481022854, 0.9665018746484156]
arguments_of_pericentre = [2.49910590507181, 2.9581175250276113, 2.5587319122352388]
longitudes_of_ascending_node = [4.960759633082254, 5.883604722095698, 4.674119055391521]
tend = 1.400000E+10
N_steps = 1000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

This was run on ada.

hamers commented 3 years ago

During a CE, dgcore.f determined that the two cores would explode in a SN and leave a massless remnant; the latter was not properly handled by binary_common_envelope_evolution(); now fixed (53d977e).