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

IndexError: list index out of range (no error code or hard exit) #23

Closed pavanvyn closed 3 years ago

pavanvyn commented 3 years ago

Version 0cf1cab2b0e0ca163f94a5f0ec63bfc5aa04bbd5

configuration = "[2,2]"
N_bodies = 4
masses = [1.172601355296052, 0.9705824261600782, 3.9456487519118917, 8.292214646979753]
metallicities = [0.017433728781051934, 0.017433728781051934, 0.017433728781051934, 0.017433728781051934]
semimajor_axes = [1.3326804114590938, 2.4482621313522195, 3539.4224364683173]
eccentricities = [0.5992544499080654, 0.43727958575359643, 0.9127659988692236]
inclinations = [1.8921877658757038, 1.845711838692201, 1.1913708833707037]
arguments_of_pericentre = [4.5914563893933185, 3.3145582976268524, 2.3573766578639206]
longitudes_of_ascending_node = [3.598051398180878, 2.7318439241437784, 3.556589478138784]
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

Two He WDs collide in this system (interesting!), and they are assumed to leave no remnant; the latter was improperly handled by collision_product(). Now fixed in 53d977e.