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

Use Krawczyk as default contractor ? #210

Open
Kolaru opened this issue Nov 20, 2024 · 1 comment
Open

Use Krawczyk as default contractor ? #210

Kolaru opened this issue Nov 20, 2024 · 1 comment
Labels

Comments

@Kolaru
Copy link
Collaborator

Kolaru commented Nov 20, 2024

I have been playing around a bit and noticed that the Krawczyk contractor can be much faster than Newton (which is expected as it doesn't use an interval-valued jacobian). For example on the Smiley 5.5 example (4 dimensional):

julia> @time roots(SmileyExample55.f, X ; contractor = Newton) ;
173.267635 seconds (95.31 M allocations: 9.645 GiB, 2.81% gc time, 0.74% compilation time)

julia> @time roots(SmileyExample55.f, X ; contractor = Krawczyk) ;
 16.432576 seconds (8.03 M allocations: 673.482 MiB, 0.51% gc time, 3.66% compilation time)

Therefore, I am tempted to make Krawczyk the default for v1.0.

Is there any reason not to do so ? @dpsanders @lbenet @lucaferranti

@Kolaru Kolaru added the v1.0 label Nov 20, 2024
@dpsanders
Copy link
Member

I can't think of a reason not to do so off the top of my head. This idea had in fact already crossed my mind :) I say go ahead!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

2 participants