def testNumericReplicator0(self):
"""Numeric QuickReplicator should yield the same result as QuickReplicator"""
for i in xrange(50):
p = Dynamics.RandomDistribution(20)
f = array([random.random() for k in xrange(20)])
r = Dynamics._numeric_QuickReplicator(p, f)
r2 = Dynamics._numeric_QuickReplicator(p, f)
self.failUnless(AlmostEqual(r, r2))
def testNumericReplicator4(self):
"""Little noise should only matter a little"""
for i in xrange(50):
p = Dynamics.RandomDistribution(100)
f = array([random.random() for k in xrange(100)])
r = Dynamics._numeric_QuickReplicator(p, f)
r2 = Dynamics._numeric_Replicator(p, f, 0.01)
self.failUnless(AlmostEqual(r, r2, 0.0101))
def testNumericReplicator3(self):
"""Noise should matter"""
for i in xrange(50):
p = Dynamics.RandomDistribution(100)
f = array([random.random() for k in xrange(100)])
r = Dynamics._numeric_QuickReplicator(p, f)
r2 = Dynamics._numeric_Replicator(p, f, 0.5)
self.failIf(AlmostEqual(r, r2))
