Fix incorrect multiplication #7
Draft
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
Closes #6
The issue is that the early returns in
select_root
can return early for an incorrect root. The multiplication operation leads to the polynomial 1-4x-64x²-64x³-64x³+256x⁴ which has roots at -0.25 and 0.0954.., the second of which is the correct answer. While shrinking the intervals, the lower bound landed on -0.25, which is a root, triggering the early return.This PR fixes the logic in
select_root
which should prevent this problem from occurring for all of the operations, not just multiplication.