@@ -3,32 +3,6 @@ using ModelingToolkit: t_nounits as t, D_nounits as D
3
3
4
4
@testset " CheckInit" begin
5
5
abstol = 1e-10
6
- @testset " Sundials + ODEProblem" begin
7
- function rhs (u, p, t)
8
- return [u[1 ] * t, u[1 ]^ 2 - u[2 ]^ 2 ]
9
- end
10
- function rhs! (du, u, p, t)
11
- du[1 ] = u[1 ] * t
12
- du[2 ] = u[1 ]^ 2 - u[2 ]^ 2
13
- end
14
-
15
- oopfn = ODEFunction {false} (rhs, mass_matrix = [1 0 ; 0 0 ])
16
- iipfn = ODEFunction {true} (rhs!, mass_matrix = [1 0 ; 0 0 ])
17
-
18
- @testset " Inplace = $(SciMLBase. isinplace (f)) " for f in [oopfn, iipfn]
19
- prob = ODEProblem (f, [1.0 , 1.0 ], (0.0 , 1.0 ))
20
- integ = init (prob, Sundials. ARKODE ())
21
- u0, _, success = SciMLBase. get_initial_values (
22
- prob, integ, f, SciMLBase. CheckInit (), Val (SciMLBase. isinplace (f)); abstol)
23
- @test success
24
- @test u0 == prob. u0
25
-
26
- integ. u[2 ] = 2.0
27
- @test_throws SciMLBase. CheckInitFailureError SciMLBase. get_initial_values (
28
- prob, integ, f, SciMLBase. CheckInit (), Val (SciMLBase. isinplace (f)); abstol)
29
- end
30
- end
31
-
32
6
@testset " Sundials + DAEProblem" begin
33
7
function daerhs (du, u, p, t)
34
8
return [du[1 ] - u[1 ] * t - p, u[1 ]^ 2 - u[2 ]^ 2 ]
0 commit comments