Skip to content

Commit

Permalink
Merge pull request #577 from SciML/os/fix-calc_W!-rosenbrock-transfor…
Browse files Browse the repository at this point in the history
…m-test

fix calc_W! test with Rosenbrock transform
  • Loading branch information
ChrisRackauckas authored Aug 2, 2024
2 parents 8c5826b + ef4a16c commit 85c4636
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions test/utility_tests.jl
Original file line number Diff line number Diff line change
Expand Up @@ -17,7 +17,7 @@ using StochasticDiffEq.SciMLOperators: MatrixOperator
jac=(u,p,t) -> A)
prob = SDEProblem(fun, u0, tspan)
integrator = init(prob, ImplicitEM(theta=1); adaptive=false, dt=dt)
W = calc_W(integrator, integrator.cache.nlsolver, dtgamma, true)
W = calc_W(integrator, integrator.cache.nlsolver, dtgamma, #=repeat_step=#false, #=W_transform=#true)
@test convert(AbstractMatrix, W) concrete_W
@test W \ u0 concrete_W \ u0

Expand All @@ -29,7 +29,7 @@ using StochasticDiffEq.SciMLOperators: MatrixOperator
prob = SDEProblem(fun, u0, tspan)
integrator = init(prob, ImplicitEM(theta=1); adaptive=false, dt=dt)
W = integrator.cache.nlsolver.cache.W
calc_W!(W, integrator, integrator.cache.nlsolver, integrator.cache, dtgamma, true)
calc_W!(W, integrator, integrator.cache.nlsolver, integrator.cache, dtgamma, #=repeat_step=#false, #=W_transform=#true)

# Did not update because it's an array operator
# We don't want to build Jacobians when we have operators!
Expand Down

0 comments on commit 85c4636

Please sign in to comment.