def test_DrawR(self):
D = 1e-12
v = -3e-8
kf = 1e-8
L = 2e-7
r0 = L / 2
gf = mod.GreensFunction1DRadAbs(D, kf, v, r0, L)
t = 1e-3
r = gf.drawR(0.5, t)
self.failIf(r < 0 or r > L)
r1 = gf.drawR(0.0, t)
r2 = gf.drawR(0.999999999999, t)
self.failIf(r1 < 0 or r1 > L)
self.failIf(r2 < 0 or r2 > L)
self.assertAlmostEqual(r1, 0)
self.assertAlmostEqual(r2, L)
print " "
print "GreensFunction1DRadAbs_test.py : test_DrawR : r1 =", r1
print "GreensFunction1DRadAbs_test.py : test_DrawR : r2 =", r2
def test_DrawR_squeezed(self):
D = 1e-12
v = -3e-8
kf = 1e-8
L = 0.02e-8
r0 = L / 2
gf = mod.GreensFunction1DRadAbs(D, kf, v, r0, L)
t = 1e-6
r0 = 0
gf.setr0(r0)
r = gf.drawR(0.5, t)
self.failIf(r < 0 or r > L)
# near s
r0 = 0.0001e-8
gf.setr0(r0)
r = gf.drawR(0.5, t)
self.failIf(r < 0 or r > L)
print " "
print "GreensFunction1DRadAbs_test.py : test_DrawR_squeezed : 1st drawn r =", r
# near a
r0 = L - 0.0001e-8
gf.setr0(r0)
r = gf.drawR(0.5, t)
self.failIf(r < 0 or r > L)
self.failIf(r < 0 or r > L)
print "GreensFunction1DRadAbs_test.py : test_DrawR_squeezed : 2nd drawn r =", r
def test_Instantiation(self):
D = 1e-12
v = -3e-8
kf = 1e-8
L = 2e-7
r0 = L / 2
gf = mod.GreensFunction1DRadAbs(D, kf, v, r0, L)
self.failIf(gf == None)
def test_DrawTime_a_equal_sigma(self):
D = 1e-12
v = -3e-8
kf = 1e-8
L = 0
r0 = L
gf = mod.GreensFunction1DRadAbs(D, kf, v, r0, L)
t = gf.drawTime(0.5)
self.assertEqual(0.0, t)
print " "
print "GreensFunction1DRadAbs_test.py : test_DrawTime_a_equal_sigma : t =", t
def test_DrawTime_r0_equal_sigma_kf_zero(self):
D = 1e-12
v = -3e-8
kf = 0.0 # note this
L = 1e-7
r0 = 0
gf = mod.GreensFunction1DRadAbs(D, kf, v, r0, L)
t = gf.drawTime(0.5) # This does not converge for negative v ~< -1e-16
self.failIf(t < 0.0 or t >= numpy.inf)
print " "
print "GreensFunction1DRadAbs_test.py : test_DrawTime_r0_equal_sigma_kf_zero : t =", t
def test_DrawR_r0_equal_sigma(self):
D = 1e-12
v = -3e-8
kf = 1e-8
L = 2e-7
r0 = 0
t = 1e-3
gf = mod.GreensFunction1DRadAbs(D, kf, v, r0, L)
r = gf.drawR(0.5, t)
self.failIf(r < 0 or r > L)
def no_test_DrawTime_r0_equal_sigma_kf_large(self):
D = 1e-12
v = -3e-8
kf = 1e-8
L = 20e-7
r0 = 1e-12
gf = mod.GreensFunction1DRadAbs(D, kf, v, r0, L)
t = gf.drawTime(0.5)
self.failIf(t < 0.0 or t >= numpy.inf)
print " "
print "GreensFunction1DRadAbs_test.py : test_DrawTime_r0_equal_sigma_kf_large : t =", t
def test_DrawTime_a_near_sigma(self):
D = 1e-12
v = -3e-8
kf = 1e-8
L = 2e-14
r0 = L / 2
gf = mod.GreensFunction1DRadAbs(D, kf, v, r0, L)
print "drawTime...a near sigma"
t = gf.drawTime(0.5)
print "done"
self.failIf(t <= 0.0 or t >= numpy.inf)
print " "
print "GreensFunction1DRadAbs_test.py : test_DrawTime_a_near_sigma : t =", t
def test_DrawR_zerot(self):
D = 1e-12
v = -3e-8
kf = 1e-8
L = 1e-7
r0 = L / 2
gf = mod.GreensFunction1DRadAbs(D, kf, v, r0, L)
t = 0.0
r = gf.drawR(0.5, t)
self.assertEqual(r0, r)
print " "
print "GreensFunction1DRadAbs_test.py : test_DrawR_zerot : r =", r
def test_DrawEventType(self):
D = 1e-12
v = -3e-8
kf = 1e-8
L = 2e-7
r0 = L / 2
gf = mod.GreensFunction1DRadAbs(D, kf, v, r0, L)
t = gf.drawTime(0.5)
eventType = gf.drawEventType(0.5, t)
self.failIf(eventType != 12 and eventType != 13 and eventType != 14)
print " "
print "GreensFunction1DRadAbs_test.py : test_DrawEventType : eventType =", eventType
eventType = gf.drawEventType(0.999999, t)
self.assertEqual(eventType, 13) # ESCAPE
eventType = gf.drawEventType(0.0, t)
self.assertEqual(eventType, 14) # REACTION
def test_DrawTime(self):
D = 1e-12
v = -3e-8
kf = 1e-8
L = 2e-7
r0 = 5e-8
gf = mod.GreensFunction1DRadAbs(D, kf, v, r0, L)
t = gf.drawTime(0.5)
self.failIf(t <= 0.0 or t >= numpy.inf)
print " "
print "GreensFunction1DRadAbs_test.py : test_DrawTime : t =", t
t = gf.drawTime(0.0)
self.failIf(t < 0.0 or t >= numpy.inf)
print "GreensFunction1DRadAbs_test.py : test_DrawTime : t =", t
t = gf.drawTime(1 - 1e-16)
self.failIf(t <= 0.0 or t >= numpy.inf)
print "GreensFunction1DRadAbs_test.py : test_DrawTime : t =", t
