|
1 | 1 | using OrdinaryDiffEq, DiffEqDevTools
|
2 |
| -probArr = Vector{ODETestProblem}(2) |
3 |
| -bigprobArr = Vector{ODETestProblem}(2) |
| 2 | +probArr = Vector{ODEProblem}(2) |
| 3 | +bigprobArr = Vector{ODEProblem}(2) |
4 | 4 |
|
5 | 5 | const linear_bigα = parse(BigFloat,"1.01")
|
6 | 6 | f = (t,u) -> (linear_bigα*u)
|
7 |
| -analytic = (t,u0) -> u0*exp(linear_bigα*t) |
| 7 | +(p::typeof(f))(::Type{Val{:analytic}},t,u0) = u0*exp(linear_bigα*t) |
8 | 8 | """Linear ODE on Float64"""
|
9 |
| -prob_ode_bigfloatlinear = ODETestProblem(f,parse(BigFloat,"0.5"),analytic) |
| 9 | +prob_ode_bigfloatlinear = ODEProblem(f,parse(BigFloat,"0.5"),(0.0,1.0)) |
10 | 10 |
|
11 | 11 | f = (t,u,du) -> begin
|
12 | 12 | for i in 1:length(u)
|
13 | 13 | du[i] = linear_bigα*u[i]
|
14 | 14 | end
|
15 | 15 | end
|
| 16 | +(p::typeof(f))(::Type{Val{:analytic}},t,u0) = u0*exp.(linear_bigα*t) |
16 | 17 | """2D Linear ODE, bigfloats"""
|
17 |
| -prob_ode_bigfloat2Dlinear = ODETestProblem(f,map(BigFloat,rand(4,2)).*ones(4,2)/2,analytic) |
| 18 | +prob_ode_bigfloat2Dlinear = ODEProblem(f,map(BigFloat,rand(4,2)).*ones(4,2)/2, (0.0,1.0)) |
18 | 19 |
|
19 | 20 | probArr[1] = prob_ode_linear
|
20 | 21 | probArr[2] = prob_ode_2Dlinear
|
|
0 commit comments