def setup_method(self):
""" this is run before each test """
self.p = rates.Nucleus("p")
self.h1 = rates.Nucleus("H1")
self.d = rates.Nucleus("d")
self.he4 = rates.Nucleus("He4")
self.c12 = rates.Nucleus("C12")
self.o16 = rates.Nucleus("O16")
self.ni56 = rates.Nucleus("Ni56")
self.u238 = rates.Nucleus("U238")
self.he4_also = rates.Nucleus("he4")
def test_solar(self):
self.comp.set_solar_like()
sum = 0.0
for k in self.comp.X:
sum += self.comp.X[k]
assert sum == approx(1.0)
assert self.comp.X[rates.Nucleus("h1")] == approx(0.7)
def setup_method(self):
""" this is run before each test """
files = [
"c12-pg-n13-ls09", "c13-pg-n14-nacr", "n13--c13-wc12",
"n13-pg-o14-lg06", "n14-pg-o15-im05", "n15-pa-c12-nacr",
"o14--n14-wc12", "o15--n15-wc12"
]
self.rc = networks.RateCollection(files)
self.p = rates.Nucleus("p")
self.he4 = rates.Nucleus("he4")
self.c12 = rates.Nucleus("c12")
self.c13 = rates.Nucleus("c13")
self.n13 = rates.Nucleus("n13")
self.n14 = rates.Nucleus("n14")
self.n15 = rates.Nucleus("n15")
self.o14 = rates.Nucleus("o14")
self.o15 = rates.Nucleus("o15")
def setup_method(self):
""" this is run before each test """
self.nuclei = [
rates.Nucleus("h1"),
rates.Nucleus("he4"),
rates.Nucleus("c12"),
rates.Nucleus("o16"),
rates.Nucleus("n14"),
rates.Nucleus("ca40")
]
self.comp = networks.Composition(self.nuclei)
def setup_method(self):
""" this is run before each test """
files = [
"b17-nnn-c14-wc12", "he3-he3pp-he4-nacr", "he4-aag-c12-fy05",
"he4-npahe3-li7-mafo", "he4-pphe3-he3-nacr", "he6-gnn-he4-cf88",
"li7-tnna-he4-mafo", "n--p-wc12", "p-ng-d-an06", "t-gn-d-nk06",
"t-pn-he3-de04"
]
self.fn = networks.StarKillerNetwork(files)
self.n = rates.Nucleus("n")
self.p = rates.Nucleus("p")
self.d = rates.Nucleus("d")
self.t = rates.Nucleus("t")
self.he3 = rates.Nucleus("he3")
self.he4 = rates.Nucleus("he4")
self.he6 = rates.Nucleus("he6")
self.li7 = rates.Nucleus("li7")
self.b17 = rates.Nucleus("b17")
self.c12 = rates.Nucleus("c12")
self.c14 = rates.Nucleus("c14")
def setup_method(self):
""" this is run before each test """
# chapter-1
self.rate1 = rates.Rate("o15--n15-wc12")
# chapter-2
self.rate2 = rates.Rate("t-gn-d-nk06")
# chapter-3
self.rate3 = rates.Rate("he6-gnn-he4-cf88")
# chapter-4
self.rate4 = rates.Rate("c12-ag-o16-nac2")
# chapter-5
self.rate5 = rates.Rate("n15-pa-c12-nacr")
# chapter-6
self.rate6 = rates.Rate("he3-he3pp-he4-nacr")
# chapter-7
self.rate7 = rates.Rate("li7-tnna-he4-mafo")
# chapter-8
self.rate8 = rates.Rate("he4-aag-c12-fy05")
# chapter-9
self.rate9 = rates.Rate("he4-pphe3-he3-nacr")
# chapter-10
self.rate10 = rates.Rate("he4-npahe3-li7-mafo")
# chapter-11
self.rate11 = rates.Rate("b17-nnn-c14-wc12")
self.n = rates.Nucleus("n")
self.p = rates.Nucleus("p")
self.h1 = rates.Nucleus("H1")
self.d = rates.Nucleus("d")
self.h3 = rates.Nucleus("H3")
self.he3 = rates.Nucleus("He3")
self.he4 = rates.Nucleus("He4")
self.he6 = rates.Nucleus("He6")
self.li7 = rates.Nucleus("Li7")
self.b17 = rates.Nucleus("B17")
self.c12 = rates.Nucleus("C12")
self.c14 = rates.Nucleus("C14")
self.n15 = rates.Nucleus("N15")
self.o15 = rates.Nucleus("O15")
self.o16 = rates.Nucleus("O16")
self.ni56 = rates.Nucleus("Ni56")
self.u238 = rates.Nucleus("U238")
self.he4_also = rates.Nucleus("he4")
def test_get_molar(self):
self.comp.set_solar_like(Z=0.02)
molar = self.comp.get_molar()
assert molar[rates.Nucleus("he4")] == approx((0.3 - 0.02) / 4.0)
