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

Enable setting Δt in minimal_period #25

Merged
merged 5 commits into from
Oct 3, 2024
Merged

Enable setting Δt in minimal_period #25

merged 5 commits into from
Oct 3, 2024

Conversation

JonasKoziorek
Copy link
Contributor

Enables user to select keyword argumet Δt, eg. time step between points in po.points trajectory. This refers to comment in #22. If the keyword is not specified, Δt is set to default T/100. This keyword argument is only valid for continuous-time systems.

@codecov-commenter
Copy link

Welcome to Codecov 🎉

Once you merge this PR into your default branch, you're all set! Codecov will compare coverage reports and display results in all future pull requests.

Thanks for integrating Codecov - We've got you covered ☂️

Copy link
Member

@Datseris Datseris left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

That was not really my comment. My comment was that minimal_period should not care about Dt. It should re-use the same points of the PO it took as input. The points on the PO remain the same if you have 5 times the period or 1. So my suggestion is that minimal_period should never call complete_orbit. Rather it should just remake the same periodic orbit with the same po.points and just the reduced period.

@JonasKoziorek
Copy link
Contributor Author

JonasKoziorek commented Oct 3, 2024

The points on the PO remain the same if you have 5 times the period or 1.

The geometric curve in the state space remains (in theory) the same. In practice if we compute 50 * minimal_T trajectory and the PO is unstable, the ODE solver will diverge from the correct trajectory. Furthermore, N * minimal_T might have duplicate points compared to minimal_T.

@Datseris
Copy link
Member

Datseris commented Oct 3, 2024

Okay. So, if I understand correctly, we, the codebase never call minimal_period. This is aconvenience function we offer to the user. Hence, it is fine if the orbit is re-completed. Correct?

@JonasKoziorek
Copy link
Contributor Author

JonasKoziorek commented Oct 3, 2024

We don't call minimal_period by default anywhere. For now it is a convenience function for the user. In the future I think there might be a situation where minimal_period is called inside of a detection algorithm. Function minimal_period is not something that would be called in a loop or so, hence I think it's fine to recompute the orbit.

@Datseris Datseris merged commit 77b065d into main Oct 3, 2024
2 checks passed
@Datseris Datseris deleted the minperiod_kwarg branch October 3, 2024 18:16
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