305305 @variables A (t) B (t) C (t)
306306 @parameters k
307307 vrj = VariableRateJump (k * (sin (t) + 1 ), [A ~ A + 1 , C ~ C + 2 ])
308- js = complete (JumpSystem ([vrj], t, [A,C], [k]; name = :js , observed = [B ~ C* A]))
308+ js = complete (JumpSystem ([vrj], t, [A, C], [k]; name = :js , observed = [B ~ C * A]))
309309 oprob = ODEProblem (js, [A => 0 , C => 0 ], (0.0 , 10.0 ), [k => 1.0 ])
310310 jprob = JumpProblem (js, oprob, Direct (); rng)
311311 sol = solve (jprob, Tsit5 ())
329329 nothing
330330 end
331331 vrj2 = VariableRateJump (vrjrate, vrjaffect!)
332- oprob2 = ODEProblem ((du,u,p, t) -> (du .= 0 ; nothing ), [0 , 0 ], (0.0 , 10.0 ), (1.0 ,))
332+ oprob2 = ODEProblem ((du, u, p, t) -> (du .= 0 ; nothing ), [0 , 0 ], (0.0 , 10.0 ), (1.0 ,))
333333 jprob2 = JumpProblem (oprob2, Direct (), vrj2; rng)
334334 cmean2 = zeros (11 )
335335 for n in 1 : N
338338 end
339339 cmean2 ./= N
340340
341- @test all ( abs .(cmean .- cmean2) .<= .05 .* cmean)
342- end
341+ @test all (abs .(cmean .- cmean2) .<= 0 .05 .* cmean)
342+ end
0 commit comments