def test3(self):
instance = AarsethZare(
nbody_system.nbody_to_si(1.0 | units.yr, 1.0 | units.MSun))
particles = self.new_system_of_sun_and_earth_and_moon()
instance.particles.add_particles(particles)
earth = instance.particles[1]
position_at_start = earth.position[0]
instance.evolve_model(365.0 | units.day)
position_after_full_rotation = earth.position[0]
self.assertAlmostRelativeEquals(position_at_start,
position_after_full_rotation, 4)
instance.evolve_model(1.5 | units.yr)
position_after_half_a_rotation = earth.position[0]
self.assertAlmostRelativeEquals(-position_at_start,
position_after_half_a_rotation, 2)
instance.evolve_model(1.75 | units.yr)
position_after_quarter_a_rotation = earth.position[1]
self.assertAlmostRelativeEquals(position_at_start,
position_after_quarter_a_rotation, 2)
instance.stop()
def test3(self):
instance = AarsethZare(nbody_system.nbody_to_si(1.0 | units.yr, 1.0 | units.MSun))
particles = self.new_system_of_sun_and_earth_and_moon()
instance.particles.add_particles(particles)
earth = instance.particles[1]
position_at_start = earth.position[0]
instance.evolve_model(365.0 | units.day)
position_after_full_rotation = earth.position[0]
self.assertAlmostRelativeEquals(position_at_start, position_after_full_rotation, 4)
instance.evolve_model(1.5 | units.yr)
position_after_half_a_rotation = earth.position[0]
self.assertAlmostRelativeEquals(-position_at_start, position_after_half_a_rotation, 2)
instance.evolve_model(1.75 | units.yr)
position_after_quarter_a_rotation = earth.position[1]
self.assertAlmostRelativeEquals(position_at_start, position_after_quarter_a_rotation, 2)
instance.stop()
def test0(self):
instance = AarsethZare()
instance.stop()
def test0(self):
instance = AarsethZare()
instance.stop()