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()