from amuse.community import * import sys sys.path.append('../../twobodymp') from interface import twobodympInterface from interface import twobodymp import numpy if __name__ == '__main__': instance = twobodympInterface() instance.initialization(20) result = instance.new_particle([0.99, 0.01], [0.1, 0.1], [0.0, .2], [0.0, 0.0], [0.0, 0.0], [0.0, 0.0], [0.0, 1.0], [0.0, 0.0]) instance.commit_particles() for i in numpy.arange(0.01, 5.0, 0.1): instance.evolve_system(i) pos = instance.get_position(1) vel = instance.get_velocity(1) print "{0} {1} {2} {3} {4} {5}".format(pos.x_, pos.y_, pos.z_, vel.vx_, vel.vy_, vel.vz_) instance.stop()
vz0 += dt_kick / 2 * az0 vx1 += dt_kick / 2 * ax1 vy1 += dt_kick / 2 * ay1 vz1 += dt_kick / 2 * az1 instance.set_velocity_mp(0, nstr(vx0, 100), nstr(vy0, 100), nstr(vz0, 100)) instance.set_velocity_mp(1, nstr(vx1, 100), nstr(vy1, 100), nstr(vz1, 100)) print x0, y0, z0, x1, y1, z1 if __name__ == '__main__': perturber = tidal = tidal_field_from_file('tideseq-1013801.dat') instance = twobodympInterface(redirection='none') instance.initialization(PRECISION) mp.dps = PRECISION mu = 1. rinit = .2 vinit = 2.95 energy = 0.5 * vinit**2 - mu / rinit a = -mu / (2 * energy) ecc = 1 - rinit / a rmax = a * (1 + ecc) print a, ecc, rmax tperiod = 2 * numpy.pi / numpy.sqrt(mu) * a**(3. / 2) tidal = tidal_field_from_file('tideseq-1013801.dat')