def slowtest3(self):
print "Testing FallbackStellarEvolution: evolve 3 stars at the same time"
# results of original code (not really check the numbers, tests have been relaxed because
# different timestepping of evolve_model and evolve_one_step)
results = dict()
results["10.0 MSun"] = dict(sse_age=3644655.52487 | units.yr,
sse_mass=9.99289747724 | units.MSun,
sse_rad=4.35949272485 | units.RSun,
sse_L=6180.22077675 | units.LSun,
evtwin_age=2999999.99732 | units.yr,
evtwin_mass=9.99832250931 | units.MSun,
evtwin_rad=4.1727904434 | units.RSun,
evtwin_L=6175.5214329 | units.LSun)
results["20.0 MSun"] = dict(sse_age=2609167.84395 | units.yr,
sse_mass=19.8895071502 | units.MSun,
sse_rad=7.30432858615 | units.RSun,
sse_L=51232.436257 | units.LSun,
evtwin_age=2999999.99792 | units.yr,
evtwin_mass=19.877338008 | units.MSun,
evtwin_rad=6.77369109826 | units.RSun,
evtwin_L=54123.7589452 | units.LSun)
results["40.0 MSun"] = dict(
sse_age=2869735.64001
| units.yr, # for some reason this is not 2617093.59096
sse_mass=37.9172423183 | units.MSun,
sse_rad=15.1388511451 | units.RSun,
sse_L=322674.980598 | units.LSun,
evtwin_age=2999999.99923 | units.yr,
evtwin_mass=37.0862227121 | units.MSun,
evtwin_rad=13.5982316053 | units.RSun,
evtwin_L=320099.608846 | units.LSun)
instance = FallbackStellarEvolution()
instance._main_se.parameters.max_age_stop_condition = 3. | units.Myr
stars = Particles(3, mass=[10., 20., 40.] | units.MSun)
stars = instance.particles.add_particles(stars)
stars.initial_mass = stars.mass
instance.commit_particles()
while instance.model_time < .05 | units.Myr:
instance.evolve_model()
for star in stars:
self.assertEqual(instance.ActiveModel[star].__class__.__name__,
"EVtwin")
while instance.model_time <= 3. | units.Myr:
instance.evolve_model()
for star in stars:
self.assertEqual(instance.ActiveModel[star].__class__.__name__,
"SSE")
for star in stars:
m = star.initial_mass
self.assertAlmostRelativeEqual(
results[str(m)]["sse_age"],
instance._FBTimeseries[star].particles[0].SSEAgeAtSwitch, 3)
evstar = star.as_particle_in_set(instance._main_se.particles)
self.assertAlmostRelativeEqual(results[str(m)]["evtwin_age"],
evstar.age, 7)
self.assertAlmostRelativeEqual(results[str(m)]["evtwin_mass"],
evstar.mass, 2)
self.assertAlmostRelativeEqual(results[str(m)]["evtwin_rad"],
evstar.radius, 2)
self.assertAlmostRelativeEqual(results[str(m)]["evtwin_L"],
evstar.luminosity, 2)
instance.stop()
def slowtest3(self):
print "Testing FallbackStellarEvolution: evolve 3 stars at the same time"
# results of original code (not really check the numbers, tests have been relaxed because
# different timestepping of evolve_model and evolve_one_step)
results=dict()
results["10.0 MSun"]=dict(sse_age=3644655.52487 | units.yr,
sse_mass=9.99289747724 | units.MSun,
sse_rad=4.35949272485 | units.RSun,
sse_L=6180.22077675| units.LSun,
evtwin_age=2999999.99732 | units.yr,
evtwin_mass=9.99832250931 | units.MSun,
evtwin_rad=4.1727904434 | units.RSun,
evtwin_L=6175.5214329 | units.LSun)
results["20.0 MSun"]=dict(sse_age= 2609167.84395| units.yr,
sse_mass= 19.8895071502| units.MSun,
sse_rad=7.30432858615| units.RSun,
sse_L=51232.436257| units.LSun,
evtwin_age= 2999999.99792| units.yr,
evtwin_mass=19.877338008| units.MSun,
evtwin_rad=6.77369109826| units.RSun,
evtwin_L= 54123.7589452| units.LSun)
results["40.0 MSun"]=dict(sse_age= 2869735.64001| units.yr, # for some reason this is not 2617093.59096
sse_mass=37.9172423183| units.MSun,
sse_rad=15.1388511451| units.RSun,
sse_L=322674.980598| units.LSun,
evtwin_age= 2999999.99923| units.yr,
evtwin_mass=37.0862227121 | units.MSun,
evtwin_rad= 13.5982316053| units.RSun,
evtwin_L=320099.608846| units.LSun)
instance = FallbackStellarEvolution()
instance._main_se.parameters.max_age_stop_condition=3.| units.Myr
stars = Particles(3,mass=[10.,20.,40.] | units.MSun)
stars = instance.particles.add_particles(stars)
stars.initial_mass=stars.mass
instance.commit_particles()
while instance.model_time<.05| units.Myr:
instance.evolve_model()
for star in stars:
self.assertEqual(instance.ActiveModel[star].__class__.__name__,"EVtwin")
while instance.model_time<=3.| units.Myr:
instance.evolve_model()
for star in stars:
self.assertEqual(instance.ActiveModel[star].__class__.__name__,"SSE")
for star in stars:
m=star.initial_mass
self.assertAlmostRelativeEqual(results[str(m)]["sse_age"],instance._FBTimeseries[star].particles[0].SSEAgeAtSwitch,3)
evstar=star.as_particle_in_set(instance._main_se.particles)
self.assertAlmostRelativeEqual(results[str(m)]["evtwin_age"],evstar.age,7)
self.assertAlmostRelativeEqual(results[str(m)]["evtwin_mass"],evstar.mass,2)
self.assertAlmostRelativeEqual(results[str(m)]["evtwin_rad"],evstar.radius,2)
self.assertAlmostRelativeEqual(results[str(m)]["evtwin_L"],evstar.luminosity,2)
instance.stop()
