def test2(self): instance = self.new_instance_of_an_optional_code(SeBa) p = Particle() p.mass = 5 | units.MSun p.metallicity = 0.02 p = instance.particles.add_particle(p) instance.evolve_model(130 | units.Myr) print(p) self.assertAlmostRelativeEqual(p.mass, 0.9906 | units.MSun, 4)
def test2(self): instance = self.new_instance_of_an_optional_code(SeBa) p = Particle() p.mass = 5 | units.MSun p.metallicity = 0.02 p = instance.particles.add_particle(p) instance.evolve_model(130 | units.Myr) print p self.assertAlmostRelativeEqual(p.mass, 0.9906 | units.MSun, 4)
def test10(self): """ Test supernova stopping condition """ instance = self.new_instance_of_an_optional_code(SeBa) instance.stopping_conditions.supernova_detection.enable() p = Particle() p.mass = 10 | units.MSun p.metallicity = 0.02 p = instance.particles.add_particle(p) instance.set_supernova_kick_velocity(0.0|units.kms) instance.evolve_model(30 | units.Myr) self.assertEquals(instance.stopping_conditions.supernova_detection.is_set(), True) self.assertEquals(instance.stopping_conditions.supernova_detection.particles(0)[0].key, p.key) self.assertAlmostRelativeEqual(p.age, 27.35866 | units.Myr, 4) self.assertAlmostRelativeEqual(p.mass, 1.22632 | units.MSun, 4) self.assertAlmostRelativeEqual(p.natal_kick_velocity, [0,0,0] | units.kms, 4)
def test10(self): """ Test supernova stopping condition """ instance = self.new_instance_of_an_optional_code(SeBa) instance.stopping_conditions.supernova_detection.enable() p = Particle() p.mass = 10 | units.MSun p.metallicity = 0.02 p = instance.particles.add_particle(p) instance.set_supernova_kick_velocity(0.0 | units.kms) instance.evolve_model(30 | units.Myr) self.assertEqual( instance.stopping_conditions.supernova_detection.is_set(), True) self.assertEqual( instance.stopping_conditions.supernova_detection.particles(0) [0].key, p.key) self.assertAlmostRelativeEqual(p.age, 27.35866 | units.Myr, 4) self.assertAlmostRelativeEqual(p.mass, 1.22632 | units.MSun, 4) self.assertAlmostRelativeEqual(p.natal_kick_velocity, [0, 0, 0] | units.kms, 4)