def test_potential_integer_division2():
"""Choose variables that can trigger integer division."""
theta = 1
a = 1
I = 1
dt = 2
Nt = 4
u, t = decay.solver(I=I, a=a, T=Nt * dt, dt=dt, theta=theta)
u_de = np.array([decay.u_discrete_exact(n, I, a, theta, dt) for n in range(Nt + 1)])
diff = np.abs(u_de - u).max()
tol = 1e-14
nt.assert_almost_equal(diff, 0, delta=1e-14)
def test_potential_integer_division2():
"""Choose variables that can trigger integer division."""
theta = 1
a = 1
I = 1
dt = 2
Nt = 4
u, t = decay.solver(I=I, a=a, T=Nt * dt, dt=dt, theta=theta)
u_de = np.array(
[decay.u_discrete_exact(n, I, a, theta, dt) for n in range(Nt + 1)])
diff = np.abs(u_de - u).max()
tol = 1E-14
nt.assert_almost_equal(diff, 0, delta=1E-14)
def test_u_discrete_exact2():
"""
Compare result from solver against
formula for the discrete solution.
"""
theta = 0.8
a = 2
I = 0.1
dt = 0.8
Nt = int(8 / dt) # no of steps
u, t = decay.solver(I=I, a=a, T=Nt * dt, dt=dt, theta=theta)
u_de = np.array([decay.u_discrete_exact(n, I, a, theta, dt) for n in range(Nt + 1)])
diff = np.abs(u_de - u).max()
nt.assert_almost_equal(diff, 0, delta=1e-14)
tol = 1e-14
success = diff < tol
assert success
def test_u_discrete_exact2():
"""
Compare result from solver against
formula for the discrete solution.
"""
theta = 0.8
a = 2
I = 0.1
dt = 0.8
Nt = int(8 / dt) # no of steps
u, t = decay.solver(I=I, a=a, T=Nt * dt, dt=dt, theta=theta)
u_de = np.array(
[decay.u_discrete_exact(n, I, a, theta, dt) for n in range(Nt + 1)])
diff = np.abs(u_de - u).max()
nt.assert_almost_equal(diff, 0, delta=1E-14)
tol = 1E-14
success = diff < tol
assert success
