-
-
Notifications
You must be signed in to change notification settings - Fork 214
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
Don't setup a JacVec
if concrete_jac==true
, and use non SciMLOperator
W_prototype
#2491
base: master
Are you sure you want to change the base?
Don't setup a JacVec
if concrete_jac==true
, and use non SciMLOperator
W_prototype
#2491
Conversation
JacVec
if concrete=jac==true
, and use non SciMLOperator
W_prototype
JacVec
if concrete=jac==true
, and use non SciMLOperator
W_prototype
JacVec
if concrete_jac==true
, and use non SciMLOperator
W_prototype
jacvec = JacVec(__f, copy(u), p, t; | ||
autodiff = alg_autodiff(alg), tag = OrdinaryDiffEqTag()) | ||
WOperator{IIP}(f.mass_matrix, dt, J, u, jacvec) | ||
WOperator{IIP}(f.mass_matrix, dt, J, u, nothing) |
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.
how would this work?
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.
by using the Jacobian. if you scroll up to the w_operator ops, they all check for whether ja cacvec exists
@@ -391,7 +391,7 @@ function do_newJW(integrator, alg, nlsolver, repeat_step)::NTuple{2, Bool} | |||
# TODO: add `isJcurrent` support for Rosenbrock solvers | |||
if !isnewton(nlsolver) | |||
isfreshJ = !(integrator.alg isa CompositeAlgorithm) && | |||
(integrator.iter > 1 && errorfail && !integrator.u_modified) |
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.
how would this not be an issue to drop?
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.
if iter<=1 we already returned on the first line (I changed that a few months ago)
Rebase? |
6a2cb4d
to
01b0c61
Compare
01b0c61
to
851e3aa
Compare
rebased now that CI ismore stable |
I'm not sure if this is a good idea yet, but I want to run it past CI.