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)