From 1a8c408e646585e8614d1dc998bc5f808959a030 Mon Sep 17 00:00:00 2001 From: Christopher Rackauckas Date: Thu, 6 Jun 2024 14:13:50 -0400 Subject: [PATCH] a few fixes --- src/gauss_adjoint.jl | 1 - test/sde_nondiag_stratonovich.jl | 2 +- 2 files changed, 1 insertion(+), 2 deletions(-) diff --git a/src/gauss_adjoint.jl b/src/gauss_adjoint.jl index fbbc9fe88..8beb698ae 100644 --- a/src/gauss_adjoint.jl +++ b/src/gauss_adjoint.jl @@ -424,7 +424,6 @@ function GaussIntegrand(sol, sensealg, checkpoints, dgdp = nothing) else pf = DiffEqBase.ParamJacobianWrapper(unwrappedf, tspan[1], y) pJ = similar(u0, length(u0), numparams) - pJ .= 0 paramjac_config = build_param_jac_config(sensealg, pf, y, p) end diff --git a/test/sde_nondiag_stratonovich.jl b/test/sde_nondiag_stratonovich.jl index 0e06eb126..ec89c2372 100644 --- a/test/sde_nondiag_stratonovich.jl +++ b/test/sde_nondiag_stratonovich.jl @@ -677,7 +677,7 @@ end _dp1 = compute_dp(p, prob!, ForwardDiffSensitivity()) _dp2 = compute_dp(p, prob!, BacksolveAdjoint(autojacvec = ReverseDiffVJP())) _dp3 = compute_dp(p, prob!, InterpolatingAdjoint(autojacvec = ReverseDiffVJP())) - @test_broken !any(isnan, compute_dp(p, prob!, InterpolatingAdjoint())) + @test !any(isnan, compute_dp(p, prob!, InterpolatingAdjoint())) @test dp1≈_dp1 rtol=1e-8 @test dp2≈_dp2 rtol=1e-8