305
305
@variables A (t) B (t) C (t)
306
306
@parameters k
307
307
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]))
309
309
oprob = ODEProblem (js, [A => 0 , C => 0 ], (0.0 , 10.0 ), [k => 1.0 ])
310
310
jprob = JumpProblem (js, oprob, Direct (); rng)
311
311
sol = solve (jprob, Tsit5 ())
329
329
nothing
330
330
end
331
331
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 ,))
333
333
jprob2 = JumpProblem (oprob2, Direct (), vrj2; rng)
334
334
cmean2 = zeros (11 )
335
335
for n in 1 : N
338
338
end
339
339
cmean2 ./= N
340
340
341
- @test all ( abs .(cmean .- cmean2) .<= .05 .* cmean)
342
- end
341
+ @test all (abs .(cmean .- cmean2) .<= 0 .05 .* cmean)
342
+ end
0 commit comments