def test2(self):
print "Test backtrack_particles"
colliders = self.new_colliders()
colliders.radius = [1, 2] | units.RSun
collision = StellarEncounterInHydrodynamics(None, None, verbose=True)
relative_position = colliders[1].position - colliders[0].position
self.assertAlmostRelativeEqual(relative_position.length(),
1.0 | units.RSun, 7)
total_energy_before = colliders.kinetic_energy(
) + colliders.potential_energy()
self.assertTrue(total_energy_before > zero)
collision.backtrack_particles(colliders)
collision.kepler.stop()
relative_position = colliders[1].position - colliders[0].position
self.assertAlmostRelativeEqual(relative_position.length(),
15 | units.RSun, 3)
total_energy_after = colliders.kinetic_energy(
) + colliders.potential_energy()
self.assertAlmostRelativeEqual(total_energy_after, total_energy_before,
7)
def test2(self):
print "Test backtrack_particles"
colliders = self.new_colliders()
colliders.radius = [1, 2] | units.RSun
collision = StellarEncounterInHydrodynamics(None, None, verbose=True)
relative_position = colliders[1].position - colliders[0].position
self.assertAlmostRelativeEqual(relative_position.length(), 1.0 | units.RSun, 7)
total_energy_before = colliders.kinetic_energy() + colliders.potential_energy()
self.assertTrue(total_energy_before > zero)
collision.backtrack_particles(colliders)
collision.kepler.stop()
relative_position = colliders[1].position - colliders[0].position
self.assertAlmostRelativeEqual(relative_position.length(), 15 | units.RSun, 3)
total_energy_after = colliders.kinetic_energy() + colliders.potential_energy()
self.assertAlmostRelativeEqual(total_energy_after, total_energy_before, 7)
