def get_component_binary_elements(comp1, comp2):
kep = Kepler(redirection = "none")
kep.initialize_code()
mass = comp1.mass + comp2.mass
pos = comp2.position - comp1.position
vel = comp2.velocity - comp1.velocity
kep.initialize_from_dyn(mass, pos[0], pos[1], pos[2],
vel[0], vel[1], vel[2])
a,e = kep.get_elements()
r = kep.get_separation()
E,J = kep.get_integrals() # per unit reduced mass, note
kep.stop()
return mass,a,e,r,E
def get_component_binary_elements(comp1, comp2):
kep = Kepler(redirection="none")
kep.initialize_code()
mass = comp1.mass + comp2.mass
pos = comp2.position - comp1.position
vel = comp2.velocity - comp1.velocity
kep.initialize_from_dyn(mass, pos[0], pos[1], pos[2], vel[0], vel[1],
vel[2])
a, e = kep.get_elements()
r = kep.get_separation()
E, J = kep.get_integrals() # per unit reduced mass, note
kep.stop()
return mass, a, e, r, E
