-
Notifications
You must be signed in to change notification settings - Fork 36
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
Support for fixed points of n-th order in fixedpoints
#326
Conversation
src/stability/fixedpoints.jl
Outdated
for i in 2:o | ||
trajectory[i] = ds.f(trajectory[i-1], p, 0.0) # trajectory from DynamicalSystemsBase won't work because of IntervalRootFinding | ||
end | ||
reverse!(trajectory) | ||
jacobians = map(x -> Jf(x, p, 0.0), trajectory) |
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.
I am not sure I understand the necessity of this trajectory
. Why make a trajectory
container and loop over, instead of initializing a Jacobian matrix and multiplying in the loop? J *= Jf(x, p, 0)
and x = ds.f(x, p, 0)
would do the trick in the loop.
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.
Good idea. I simplified the function.
Codecov ReportAttention: Patch coverage is
Additional details and impacted files@@ Coverage Diff @@
## main #326 +/- ##
==========================================
+ Coverage 54.86% 56.05% +1.18%
==========================================
Files 22 22
Lines 904 917 +13
==========================================
+ Hits 496 514 +18
+ Misses 408 403 -5 ☔ View full report in Codecov by Sentry. |
I have added support for computing fixed points of n-th iteration of a discrete dynamical system.