Skip to content
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 a typo in Krylov oop tests #1962

Merged
merged 2 commits into from
Jun 11, 2023
Merged

Conversation

gaurav-arya
Copy link
Member

@gaurav-arya gaurav-arya commented Jun 9, 2023

resolves #1956

@gaurav-arya
Copy link
Member Author

@vpuri3 @ChrisRackauckas

@@ -7,7 +7,7 @@ let N = 20
A = diagm(-1 => du, 0 => dd, 1 => du)
_f = (u, p, t) -> A * u - u .^ 3
_f_ip = (du, u, p, t) -> (mul!(du, A, u); du .-= u .^ 3)
_jac = (A, u, p, t) -> A - 3 * diagm(0 => u .^ 2)
_jac = (J, u, p, t) -> A - 3 * diagm(0 => u .^ 2)
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Wait then what is J supposed to be?

Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

It's the old op matrix, which probably shouldn't be provided for OOP update_coefficients since unlike the in place version it doesn't need to be overwritten: SciML/SciMLOperators.jl#202

This PR just gets things to work with the current signature of update_coefficienfs

@ChrisRackauckas ChrisRackauckas merged commit 5de51dd into SciML:master Jun 11, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

Out of place exponential Krylov regression
2 participants