|
59 | 59 | tspan = (0.0, 6.0)
|
60 | 60 |
|
61 | 61 | op = ODEProblem(pend, u0map, tspan, parammap)
|
62 |
| - osol = solve(op) |
| 62 | + osol = solve(op, Vern9()) |
63 | 63 |
|
64 | 64 | bvp = SciMLBase.BVProblem{true, SciMLBase.AutoSpecialize}(pend, u0map, tspan, parammap)
|
65 | 65 | for solver in solvers
|
|
111 | 111 | end
|
112 | 112 |
|
113 | 113 | u0 = [1., 2.]; p = [1.5, 1., 1., 3.]
|
114 |
| - genbc_iip = ModelingToolkit.generate_function_bc(lksys, u0, [1, 2], tspan, true) |
115 |
| - genbc_oop = ModelingToolkit.generate_function_bc(lksys, u0, [1, 2], tspan, false) |
| 114 | + fns = ModelingToolkit.generate_function_bc(lksys, u0, [1, 2], tspan) |
| 115 | + genbc_oop, genbc_iip = ModelingToolkit.eval_or_rgf.(fns) |
116 | 116 |
|
117 | 117 | bvpi1 = SciMLBase.BVProblem(lotkavolterra!, bc!, [1.,2.], tspan, p)
|
118 | 118 | bvpi2 = SciMLBase.BVProblem(lotkavolterra!, genbc_iip, [1.,2.], tspan, p)
|
|
141 | 141 | end
|
142 | 142 |
|
143 | 143 | u0 = [1, 1.]
|
144 |
| - genbc_iip = ModelingToolkit.generate_function_bc(lksys, u0, [1], tspan, true) |
145 |
| - genbc_oop = ModelingToolkit.generate_function_bc(lksys, u0, [1], tspan, false) |
| 144 | + fns = ModelingToolkit.generate_function_bc(lksys, u0, [1], tspan) |
| 145 | + genbc_oop, genbc_iip = ModelingToolkit.eval_or_rgf.(fns) |
146 | 146 |
|
147 | 147 | bvpi1 = SciMLBase.BVProblem(lotkavolterra!, bc!, u0, tspan, p)
|
148 | 148 | bvpi2 = SciMLBase.BVProblem(lotkavolterra!, genbc_iip, u0, tspan, p)
|
|
0 commit comments