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

Scipy LSODA for stiff ODE #249

Open
janrohleff opened this issue Apr 3, 2024 · 1 comment
Open

Scipy LSODA for stiff ODE #249

janrohleff opened this issue Apr 3, 2024 · 1 comment

Comments

@janrohleff
Copy link

Hey,

thank you for this great package. I have a problem using the scipy LSODA solver together with odeint_adjoint for a nonlinear stiff ODE. The problem is the following:

  1. The forward solve with the scipy solver works well.
  2. When I try to compute the backward pass the solver always get stuck at some specific time-point depending on the tolerance. I don‘t get an error and the file is still running, but it doesn‘t make any progress, i.e. stuck at this t. I think the problem is that you solve in the „backward call“ the state-equation backward in time, which is numericaly not a good idea for stiff ODEs.

But do you have any idea how this could be fixed? Or do you thinking of implementing a Torch version of LSODA, so that one just can use automatic-differentiation and use odeint.

Thank you!

@HoangMinhPhan
Copy link

I also got the same problem for stiff ODE. I hope that the author can introduce new solvers to solve stiff problems, or adapt scipy solvers to be compatible with backward in pytorch.

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

No branches or pull requests

2 participants