-
-
Notifications
You must be signed in to change notification settings - Fork 210
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
fix out of place WOperator
#2291
fix out of place WOperator
#2291
Conversation
9a99596
to
a5f0be8
Compare
0239364
to
0772416
Compare
SciML/StochasticDiffEq.jl#574 fixes the StochasticDiffEq test failure (which was caused by testing that you can call calc_W! on an out of place ODE which you just can't do). The other test failures look unrelated. I think this is finally ready to merge. |
@@ -908,42 +873,36 @@ function build_J_W(alg, u, uprev, p, t, dt, f::F, ::Type{uEltypeNoUnits}, | |||
else | |||
deepcopy(f.jac_prototype) | |||
end | |||
__f = if IIP | |||
(du, u, p, t) -> _f(du, u, p, t) | |||
W = if J isa StaticMatrix && alg isa OrdinaryDiffEqRosenbrockAdaptiveAlgorithm |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
why is this Rosenbrock specific?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
not sure. I was vaguely following how this worked previously
#2252 made it so there were some cases in which we weren't updating J correctly leading to convergence failures. Debugging this was hell, and all of this code needs to be ripped out in the nearish future and replaced with DifferentiationInterface and a sane layout. While I was at it, I also fixed a couple places where the stats were completely broken for out of place. Tests incoming.