@@ -22,16 +22,18 @@ ps = [p => zeros(3, 3),
22
22
q => 1.0 ]
23
23
tspan = (0.0 , 10.0 )
24
24
@mtkcompile sys = System (eqs, t)
25
- prob = ODEProblem (sys, u0, tspan, ps )
25
+ prob = ODEProblem (sys, [u0; ps], tspan )
26
26
sol = solve (prob, Tsit5 ())
27
27
28
28
mtkparams = parameter_values (prob)
29
29
new_p = rand (14 )
30
- gs = gradient (new_p) do new_p
31
- new_params = SciMLStructures. replace (SciMLStructures. Tunable (), mtkparams, new_p)
32
- new_prob = remake (prob, p = new_params)
33
- new_sol = solve (new_prob, Tsit5 ())
34
- sum (new_sol)
30
+ @test_broken begin
31
+ gs = gradient (new_p) do new_p
32
+ new_params = SciMLStructures. replace (SciMLStructures. Tunable (), mtkparams, new_p)
33
+ new_prob = remake (prob, p = new_params)
34
+ new_sol = solve (new_prob, Tsit5 ())
35
+ sum (new_sol)
36
+ end
35
37
end
36
38
37
39
@testset " Issue#2997" begin
50
52
sys = mtkcompile (sys)
51
53
52
54
function x_at_0 (θ)
53
- prob = ODEProblem (sys, [sys. x => 1.0 ], ( 0.0 , 1.0 ), [ sys. ργ0 => θ[1 ], sys. h => θ[2 ]])
55
+ prob = ODEProblem (sys, [sys. x => 1.0 , sys. ργ0 => θ[1 ], sys. h => θ[2 ]], ( 0.0 , 1.0 ) )
54
56
return prob. u0[1 ]
55
57
end
56
58
61
63
@named sys = System (
62
64
Equation[], t, [], [a, b, c, d, e, f, g, h],
63
65
continuous_events = [ModelingToolkit. SymbolicContinuousCallback (
64
- [a ~ 0 ] => [c ~ 0 ], discrete_parameters = c)])
66
+ [a ~ 0 ] => [c ~ 0 ], discrete_parameters = c, iv = t )])
65
67
sys = complete (sys)
66
68
67
69
ivs = Dict (c => 3 a, b => ones (3 ), a => 1.0 , d => 4 , e => [5.0 , 6.0 , 7.0 ],
@@ -116,7 +118,7 @@ fwd, back = ChainRulesCore.rrule(remake_buffer, sys, ps, idxs, vals)
116
118
sys = mtkcompile (sys)
117
119
118
120
# Find initial throw velocity that reaches exactly 10 m after 1 s
119
- dprob0 = ODEProblem (sys, [D (y) => NaN ], (0.0 , 1.0 ), [] ; guesses = [y => 0.0 ])
121
+ dprob0 = ODEProblem (sys, [D (y) => NaN ], (0.0 , 1.0 ); guesses = [y => 0.0 ])
120
122
function f (ics, _)
121
123
dprob = remake (dprob0, u0 = Dict (D (y) => ics[1 ]))
122
124
dsol = solve (dprob, Tsit5 ())
0 commit comments