def test_stratSRS2_R74(self, exact_solution_R74):
(dW, I, J, f, f_strat, G, G_separate, y0, tspan,
y) = exact_solution_R74
ySRS2 = sdeint.stratSRS2(f_strat, G_separate, y0, tspan, dW=dW, J=J)
_assert_close(ySRS2, y, 1e-2, 1e-2)
return ySRS2
def test_stratSRS2_R74(self, exact_solution_R74):
(dW, I, J, f, f_strat, G, G_separate, y0, tspan,y) = exact_solution_R74
ySRS2 = sdeint.stratSRS2(f_strat, G_separate, y0, tspan, dW=dW, J=J)
_assert_close(ySRS2, y, 1e-2, 1e-2)
return ySRS2
def test_stratSRS2_KPS445(self, exact_solution_KPS445):
(dW, I, J, f, f_strat, G, y0, tspan, y) = exact_solution_KPS445
ySRS2 = sdeint.stratSRS2(f_strat, G, y0, tspan, dW=dW, J=J)[:, 0]
_assert_close(ySRS2, y, 1e-2, 1e-2)
return ySRS2
def test_stratSRS2_KPS445(self, exact_solution_KPS445):
(dW, I, J, f, f_strat, G, y0, tspan, y) = exact_solution_KPS445
ySRS2 = sdeint.stratSRS2(f_strat, G, y0, tspan, dW=dW, J=J)[:,0]
_assert_close(ySRS2, y, 1e-2, 1e-2)
return ySRS2