def test0(self): centre, orbiter = self.sun_and_earth() mercury = MercuryWayWard() #,debugger='xterm') mercury.initialize_code() mercury.central_particle.add_particles(centre) mercury.orbiters.add_particles(orbiter) mercury.commit_particles() self.assertAlmostEqual(mercury.central_particle.j4, .0 | units.AU**4) self.assertAlmostEqual(mercury.central_particle.mass, 1.98892e+30 | units.kg, 3) self.assertAlmostEqual(mercury.central_particle.mass, 1.0 | units.MSun, 3) self.assertEquals(mercury.get_number_of_orbiters(), 1) self.assertEquals(mercury.orbiters.position, [[1, 0, 0]] | units.AU) self.assertEquals(mercury.orbiters.density, 1.0 | units.g / units.cm**3) self.assertEquals(mercury.orbiters.angular_momentum, [[1.0, 0.0, 0.0]] | units.MSun * units.AU**2 / units.day) mercury.evolve_model(365.24 | units.day) self.assertAlmostEqual(mercury.orbiters.position, [[1.0, 0.0, 0.0]] | units.AU, 1) self.assertAlmostEqual( mercury.kinetic_energy + mercury.potential_energy, mercury.total_energy, 3) mercury.stop()
def test0(self): centre, orbiter = self.sun_and_earth() mercury = MercuryWayWard()#,debugger='xterm') mercury.initialize_code() mercury.central_particle.add_particles(centre) mercury.orbiters.add_particles(orbiter) mercury.commit_particles() self.assertAlmostEqual(mercury.central_particle.j4, .0|units.AU**4) self.assertAlmostEqual(mercury.central_particle.mass, 1.98892e+30 |units.kg, 3) self.assertAlmostEqual(mercury.central_particle.mass, 1.0 |units.MSun, 3) self.assertEquals(mercury.get_number_of_orbiters(),1) self.assertEquals(mercury.orbiters.position, [[1,0,0]] | units.AU) self.assertEquals(mercury.orbiters.density, 1.0|units.g/units.cm**3 ) self.assertEquals(mercury.orbiters.angular_momentum, [[1.0, 0.0, 0.0]] | units.MSun*units.AU**2/units.day) mercury.evolve_model(365.24 | units.day) self.assertAlmostEqual(mercury.orbiters.position, [[1.0, 0.0, 0.0]] | units.AU, 1) self.assertAlmostEqual(mercury.kinetic_energy+mercury.potential_energy,mercury.total_energy,3) mercury.stop()