-
Notifications
You must be signed in to change notification settings - Fork 15
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
Newton error propagation based convergence check. #1403
Conversation
Codecov ReportAttention: Patch coverage is
Additional details and impacted files@@ Coverage Diff @@
## master #1403 +/- ##
==========================================
- Coverage 86.56% 86.54% -0.03%
==========================================
Files 182 182
Lines 13549 13559 +10
==========================================
+ Hits 11729 11734 +5
- Misses 1820 1825 +5 ☔ View full report in Codecov by Sentry. |
This comment has been minimized.
This comment has been minimized.
This comment has been minimized.
This comment has been minimized.
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.
Looks sensible.
Instead of `(X - X_old) = dt * dX`. The implemented scaling is the one implemented in NOCMODL. It should prevent the residual from growing linearly with `dt`.
Includes a test that solves the pump equation ~ X + Y <-> Z with extreme compartment sizes and steady state calculation.
Logfiles from GitLab pipeline #228874 (:white_check_mark:) have been uploaded here! Status and direct links: |
Since we have access to the Jacobian, we can use Newton error propagation:
with
sd = eps * X
to estimate the impact of roundoff errors onF
.