Skip to content

Commit

Permalink
Update sde_nondiag_stratonovich.jl
Browse files Browse the repository at this point in the history
  • Loading branch information
ChrisRackauckas authored Jan 10, 2023
1 parent d3e0de2 commit 34c1809
Showing 1 changed file with 5 additions and 6 deletions.
11 changes: 5 additions & 6 deletions test/sde_nondiag_stratonovich.jl
Original file line number Diff line number Diff line change
Expand Up @@ -161,19 +161,18 @@ end
dt = dtnd, adaptive = false,
sensealg = BacksolveAdjoint(autojacvec = ZygoteVJP()))

@test all(isapprox(res_sde_u0a, res_sde_u0, rtol = 1e-6),isapprox(res_sde_pa, res_sde_p, rtol = 1e-6))
@test all((isapprox(res_sde_u0a, res_sde_u0, rtol = 1e-6),isapprox(res_sde_pa, res_sde_p, rtol = 1e-6)))
end

@info res_sde_pa

@test_broken begin
res_sde_u0a, res_sde_pa = adjoint_sensitivities(sol, EulerHeun(), t = Array(t),
res_sde_u0a, res_sde_pa = adjoint_sensitivities(sol, EulerHeun(), t = Array(t),
dgdu_discrete = dg!,
dt = dtnd, adaptive = false,
sensealg = BacksolveAdjoint(autojacvec = false))

@test all(isapprox(res_sde_u0a, res_sde_u0, rtol = 1e-6),isapprox(res_sde_pa, res_sde_p, rtol = 1e-6))
end
@test isapprox(res_sde_u0a, res_sde_u0, rtol = 1e-6)
@test isapprox(res_sde_pa, res_sde_p, rtol = 1e-6)

@info res_sde_pa

Expand All @@ -195,7 +194,7 @@ end
dt = dtnd, adaptive = false,
sensealg = InterpolatingAdjoint(autojacvec = ZygoteVJP()))

@test all(isapprox(res_sde_u0a, res_sde_u0, rtol = 1e-5),isapprox(res_sde_pa, res_sde_p, rtol = 1e-4))
@test all((isapprox(res_sde_u0a, res_sde_u0, rtol = 1e-5),isapprox(res_sde_pa, res_sde_p, rtol = 1e-4)))
end

@info res_sde_pa
Expand Down

0 comments on commit 34c1809

Please sign in to comment.