def test_replace_Species_in_ComplexSpecies(self):
c1 = Complex([self.s1, self.s_old])
c2 = Complex([self.s1, c1])
self.assertTrue(c1.replace_species(self.s2, self.s_new) == c1)
self.assertTrue(
c1.replace_species(self.s_old, self.s_new) == Complex(
[self.s1, self.s_new]))
self.assertTrue(c2 == Complex([self.s1, c1]))
self.assertTrue(
c2.replace_species(self.s_new, self.s_old) == Complex(
[self.s1, Complex([self.s1, self.s_old])]))
def test_replace_species_with_a_non_massaction_reaction(self):
c1 = Complex([self.s1, self.s_old])
prop_hill_old = ProportionalHillPositive(k=1.,
s1=self.s1,
K=10,
d=self.s_old,
n=2)
r1 = Reaction([self.s1, self.s_old], [c1],
propensity_type=prop_hill_old)
prop_hill_new = ProportionalHillPositive(k=1.,
s1=self.s1,
K=10,
d=self.s_new,
n=2)
r1_new = Reaction([self.s1, self.s_new],
[c1.replace_species(self.s_old, self.s_new)],
propensity_type=prop_hill_new)
self.assertTrue(r1.replace_species(self.s_old, self.s_new) == r1_new)
def test_replace_Species_in_OrderedComplexSpecies(self):
oc1 = Complex([self.s1, self.s_old], ordered=True)
self.assertTrue(
oc1.replace_species(self.s_old, self.s_new) == Complex(
[self.s1, self.s_new], ordered=True))