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)
