Adaptive implicit solvers and sparse Sherman-Morrison iterator#281
Open
psv4 wants to merge 1 commit into
Open
Conversation
commit 9b4c83d Author: psv4 <44118604+psv4@users.noreply.github.com> Date: Sat Jun 6 08:01:52 2026 -0400 Adjust doc placement commit 1a06502 Author: psv4 <44118604+psv4@users.noreply.github.com> Date: Fri Jun 5 21:57:59 2026 -0400 Updated documentation commit 0b9bf64 Author: psv4 <44118604+psv4@users.noreply.github.com> Date: Fri Jun 5 18:17:52 2026 -0400 Added Kvaerno4 and Kvaerno5 commit 6808945 Author: psv4 <44118604+psv4@users.noreply.github.com> Date: Fri Jun 5 12:43:44 2026 -0400 Add Sherman Morrison fixes to fixed FIRK/DIRK methods commit b5c1849 Author: psv4 <44118604+psv4@users.noreply.github.com> Date: Fri Jun 5 08:58:28 2026 -0400 Working FIRK was fine, DIRK was not recreating the full k matrix. GL4 and GL6 still slow, KV3 much faster commit 92e6674 Author: psv4 <44118604+psv4@users.noreply.github.com> Date: Thu Jun 4 22:17:59 2026 -0400 Update adaptive_implicit.py commit 93bae2d Author: psv4 <44118604+psv4@users.noreply.github.com> Date: Thu Jun 4 22:12:55 2026 -0400 Added Kvaerno3 Still not converging :( commit 5fc2c9e Author: psv4 <44118604+psv4@users.noreply.github.com> Date: Thu Jun 4 19:19:42 2026 -0400 Running but slow commit bbd783e Author: psv4 <44118604+psv4@users.noreply.github.com> Date: Thu Jun 4 14:41:49 2026 -0400 Add GL6 interpolant commit e2ee1d9 Author: psv4 <44118604+psv4@users.noreply.github.com> Date: Thu Jun 4 14:41:11 2026 -0400 Iteration Limit commit 7593a10 Author: psv4 <44118604+psv4@users.noreply.github.com> Date: Thu Jun 4 14:36:39 2026 -0400 Typo for DIRK commit 887936d Author: psv4 <44118604+psv4@users.noreply.github.com> Date: Thu Jun 4 14:34:41 2026 -0400 Complete commit 7d67eb5 Author: psv4 <44118604+psv4@users.noreply.github.com> Date: Thu Jun 4 12:54:03 2026 -0400 Template for new classes and functions
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
1 participant
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
Big Update
Closes #267 and Closes #214
Adds adaptive step implicit solver methods including:
adaptive_gl4adaptive_gl6kvaerno3kvaerno4kvaerno5Only
kvaerno3was added to the test suite as the others take a long time to evaluate.Fixed and adaptive fully implicit Runge-Kutta (FIRK) and diagonally implicit Runge-Kutta (DIRK) now use the good Sherman-Morrison method with a sparse diagonal Jacobian to avoid matrix inversion and scale to large problems.
New solvers have been added to
README.MDandFURTHER_DOCUMENTATION.MD.Side Note:
TestMinMaxStepintests/odeint_tests.pyfails with thescipy_solvernot respecting themin_stepoption in scipy version 1.17.1.