def run_kepler(mass, semi, ecc, time):
kep = Kepler(redirection='none')
kep.initialize_code()
kep.set_longitudinal_unit_vector(1.0, 1.0,
0.0)
kep.initialize_from_elements(mass, semi, ecc)
a,e = kep.get_elements()
p = kep.get_periastron()
print "elements:", a, e, p
kep.transform_to_time(time)
x,y,z = kep.get_separation_vector()
print "separation:", x,y,z
x,y,z = kep.get_longitudinal_unit_vector()
print "longitudinal:", x,y,z
pos = [1, 0, 0] | nbody_system.length
vel = [0, 0.5, 0] | nbody_system.speed
kep.initialize_from_dyn(mass, pos[0], pos[1], pos[2],
vel[0], vel[1], vel[2])
a,e = kep.get_elements()
p = kep.get_periastron()
print "elements:", a, e, p
kep.transform_to_time(time)
x,y,z = kep.get_separation_vector()
print "separation:", x,y,z
x,y,z = kep.get_velocity_vector()
print "velocity:", x,y,z
x,y,z = kep.get_longitudinal_unit_vector()
print "longitudinal:", x,y,z
kep.set_random(42)
kep.make_binary_scattering(0.5 | nbody_system.mass,
0.5,
0.5 | nbody_system.mass,
0.0 | nbody_system.speed,
0.0 | nbody_system.length,
1.e-6,
0)
kep.stop()
