Full argument reduction for the elliptic integrals #2368
Merged
+2,440
−138
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.
Amusingly, I think that this uncovered a bug in Mathematica 11.0. You'd expect EllipticPi[3/2, N[3 Pi,200], 9/10] and N[EllipticPi[3/2, 3 Pi, 9/10], 200] to return roughly the same value, but they don't. I believe that the latter is correct.
Benchmarks before:
after: