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

Collection of the Lyapunov exponents over time #318

Closed
rseydam opened this issue Oct 31, 2023 · 1 comment
Closed

Collection of the Lyapunov exponents over time #318

rseydam opened this issue Oct 31, 2023 · 1 comment

Comments

@rseydam
Copy link

rseydam commented Oct 31, 2023

Is it possible to collect the computed Lyapunovspectra over time to check on the convergence of the exponents?

@Datseris
Copy link
Member

Closing in favor of the discussion of #317 , as the answer here is similar:

This option is not available out of the box. It will likely not be available out of the box in the near future either. That is because adding auxilary options for every possibility some user will find useful would obfuscate the code and reduce code clarity, making more than 50% of the lines of code auxilary options rather than computing Lyapunov exponents. To add this option you can trivially create a custom function by

  1. You copy paste this function: https://github.com/JuliaDynamics/ChaosTools.jl/blob/main/src/chaosdetection/lyapunovs/lyapunovspectrum.jl#L67-L99 and rename it to mylyapunovspectrum,
  2. You add the exponents_history = [zeros(k)] at the start of the function.
  3. In the core step loop, right after https://github.com/JuliaDynamics/ChaosTools.jl/blob/main/src/chaosdetection/lyapunovs/lyapunovspectrum.jl#L89 you add the line push!(exponents_history, λ ./ (current_time(tands) - t0)).
  4. You change the return statement to return λ, exponent_history.

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