@@ -22,14 +22,14 @@ prs = exp.(2 .* rand(2))
2222ps = [ prs[1 ], prs[1 ] / exp (4 * rand ()),
2323 prs[2 ], prs[2 ] / exp (4 * rand ()) ]
2424
25- Nmax = 130
25+ Nmax = 100
2626u0 = zeros (Nmax+ 1 , Nmax+ 1 )
2727u0[1 ] = 1.0
2828
2929tt = [ 0.25 , 1.0 , 10.0 ]
3030
3131prob = convert (ODEProblem, sys, u0, 10.0 , ps)
32- sol = solve (prob, Vern7 (), abstol= 1e-9 , reltol = 1e- 6 , saveat= tt)
32+ sol = solve (prob, Vern7 (), abstol= 1e-6 , saveat= tt)
3333
3434@test marg (sol. u[1 ], dims= 2 ) ≈ pdf .(Poisson (ps[1 ] / ps[2 ] * (1 - exp (- ps[2 ] * tt[1 ]))), 0 : Nmax) atol= 1e-4
3535@test marg (sol. u[1 ], dims= 1 ) ≈ pdf .(Poisson (ps[3 ] / ps[4 ] * (1 - exp (- ps[4 ] * tt[1 ]))), 0 : Nmax) atol= 1e-4
@@ -44,7 +44,7 @@ A = convert(SparseMatrixCSC, sys, (Nmax+1, Nmax+1), ps, 0)
4444f = (du,u,t) -> mul! (du, A, u)
4545
4646probA = ODEProblem (f, u0, 10.0 )
47- solA = solve (prob, Vern7 (), abstol= 1e-9 , reltol = 1e- 6 , saveat= tt)
47+ solA = solve (prob, Vern7 (), abstol= 1e-6 , saveat= tt)
4848
4949@test sol. u[1 ] ≈ solA. u[1 ] atol= 1e-4
5050@test sol. u[2 ] ≈ solA. u[2 ] atol= 1e-4
0 commit comments