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

Support for fixed points of n-th order in fixedpoints #326

Merged
merged 4 commits into from
Mar 18, 2024

Conversation

JonasKoziorek
Copy link
Contributor

I have added support for computing fixed points of n-th iteration of a discrete dynamical system.

src/stability/fixedpoints.jl Outdated Show resolved Hide resolved
Comment on lines 116 to 120
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)
Copy link
Member

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.

Copy link
Contributor Author

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.

test/stability/fixedpoints.jl Outdated Show resolved Hide resolved
@codecov-commenter
Copy link

Codecov Report

Attention: Patch coverage is 92.85714% with 1 lines in your changes are missing coverage. Please review.

Project coverage is 56.05%. Comparing base (928bdb0) to head (2602883).
Report is 2 commits behind head on main.

Files Patch % Lines
src/stability/fixedpoints.jl 92.85% 1 Missing ⚠️
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.
📢 Have feedback on the report? Share it here.

@Datseris Datseris merged commit 80f9b39 into JuliaDynamics:main Mar 18, 2024
1 of 2 checks passed
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.

3 participants