@@ -931,22 +931,15 @@ testdict = Dict([:name => "test"])
931931@named sys = ODESystem (eqs, t, metadata = testdict)
932932@test get_metadata (sys) == testdict
933933
934- @variables P (t)= 0 Q (t)= 2
935- ∂t = D
936-
937- eqs = [∂t (Q) ~ 1 / sin (P)
938- ∂t (P) ~ log (- cos (Q))]
939- @named sys = ODESystem (eqs, t, [P, Q], [])
940- sys = complete (debug_system (sys));
941- prob = ODEProblem (sys, [], (0 , 1.0 ));
942- du = zero (prob. u0);
943- if VERSION < v " 1.8"
944- @test_throws DomainError prob. f (du, [1 , 0 ], prob. p, 0.0 )
945- @test_throws DomainError prob. f (du, [0 , 2 ], prob. p, 0.0 )
946- else
947- @test_throws " -cos(Q(t))" prob. f (du, [1 , 0 ], prob. p, 0.0 )
948- @test_throws " sin(P(t))" prob. f (du, [0 , 2 ], prob. p, 0.0 )
949- end
934+ @variables P (t) = NaN Q (t) = NaN
935+ @named sys = ODESystem ([
936+ D (Q) ~ 1 / sin (P)
937+ D (P) ~ log (- cos (Q))
938+ ], t, [P, Q], [])
939+ sys = complete (debug_system (sys))
940+ prob = ODEProblem (sys, [], (0.0 , 1.0 ))
941+ @test_throws " log(-cos(Q(t))) errors" prob. f ([1 , 0 ], prob. p, 0.0 )
942+ @test prob. f ([0 , 2 ], prob. p, 0.0 )[1 ] == 1 / 0
950943
951944let
952945 @variables x (t) = 1
0 commit comments