def test6(self): print("Test 6: testing user interface (nbody units)") numpy.random.seed(345672) n = 10000 masses = new_flat_mass_distribution_nbody(n) self.assertEqual(len(masses), n) self.assertAlmostEqual(masses.sum(), 1.0 | nbody_system.mass) self.assertAlmostEqual(masses.mean(), 1.0 / n | nbody_system.mass) self.assertAlmostEqual(masses.amin(), 0.100145673289 / 17.5145247111 / n | nbody_system.mass) self.assertAlmostEqual(masses.amax(), 124.94980234 / 17.5145247111 / n | nbody_system.mass)
def test6(self): print "Test 6: testing user interface (nbody units)" numpy.random.seed(345672) n = 10000 masses = new_flat_mass_distribution_nbody(n) self.assertEqual(len(masses), n) self.assertAlmostEqual(masses.sum(), 1.0 | nbody_system.mass) self.assertAlmostEqual(masses.mean(), 1.0 / n | nbody_system.mass) self.assertAlmostEqual(masses.amin(), 0.100145673289 / 17.5145247111 / n | nbody_system.mass) self.assertAlmostEqual(masses.amax(), 124.94980234 / 17.5145247111 / n | nbody_system.mass)
def test7(self): print("Test 7: testing user interface (nbody units), optional args") numpy.random.seed(345672) n = 10000 masses = new_flat_mass_distribution_nbody(n, mass_min=10.0|nbody_system.mass, mass_max=100.0|nbody_system.mass) self.assertEqual(len(masses), n) self.assertAlmostEqual(masses.sum(), 1.0 | nbody_system.mass) self.assertAlmostEqual(masses.mean(), 1.0 / n | nbody_system.mass) self.assertAlmostEqual(masses.amin(), 10.0047015091 / 39.1111546565 / n | nbody_system.mass) self.assertAlmostEqual(masses.amax(), 99.9870310764 / 39.1111546565 / n | nbody_system.mass)
def test7(self): print "Test 7: testing user interface (nbody units), optional args" numpy.random.seed(345672) n = 10000 masses = new_flat_mass_distribution_nbody(n, mass_min=10.0|nbody_system.mass, mass_max=100.0|nbody_system.mass) self.assertEqual(len(masses), n) self.assertAlmostEqual(masses.sum(), 1.0 | nbody_system.mass) self.assertAlmostEqual(masses.mean(), 1.0 / n | nbody_system.mass) self.assertAlmostEqual(masses.amin(), 10.0047015091 / 39.1111546565 / n | nbody_system.mass) self.assertAlmostEqual(masses.amax(), 99.9870310764 / 39.1111546565 / n | nbody_system.mass)