Repository navigation
Two half-odd powers of secant sums in the half angle's sine - #1828
Merged
Merged
Conversation
Part of #718. Co-Authored-By: Claude Opus 5.5 (1M context) <noreply@anthropic.com> Claude-Session: https://claude.ai/code/session_012sonx8iAspMiwRwokT1Ura
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
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.
Part of #718.
sec(e + f x) sqrt(a + a sec(e + f x))/sqrt(c + d sec(e + f x))was declined, with eight more of Rubi's 4.5.2.1 and 4.5.2.3 that are a half-odd power ofa + a secbeside one ofc + d sec. The half-angle tangent reads a root of the first sum and leaves the second as another root, which it declines. Written in the cosine instead,so the two powers of
cos(y)make a whole one, and inw = sin(y/2), wherecos(y) = 1 - 2w^2anddy = 2 dw/cos(y/2), what is left is rational beside the one root ofc + d - 2c w^2, which is answered at once.SolveTwoHalfOddPowersOfSecantSumsByTheHalfAngleSinetakes the secant's, the cosine's, the minus sign's (incos(y/2)) and the cosecant's (by the complement). The constants are not written: the answer is the integrand times the antiderivative inwover what that differentiates back to, a quotient whose square is one. It is kept only where that square is one at the sampled points.e1a1a3fcsec(x) sqrt(a + a sec(x))/sqrt(c + d sec(x))sec(x) sqrt(c + d sec(x))/sqrt(a + a sec(x))sec(x)/(sqrt(a + a sec(x)) sqrt(c + d sec(x)))sec(x)^2/(sqrt(a + a sec(x)) sqrt(c + d sec(x)))sqrt(c + d sec(x))/sqrt(a + a sec(x))xabove stands fore + f x, which is what was probed. The other four rows answer in 891 to 6,066 characters. The longest answers carry a case for each sign of the root's leading coefficient.Tests:
TwoRootsOfSecantSumsIntegralTest, five rows (shifted argument, minus sign, cosecant), each differentiated back and compared as a complex number at six real points.Measured on every corpus problem with a secant or a cosecant in its integrand, 6,489 of them, at the corpus's 5-second budget, against master
e1a1a3fc:Nine problems are answered here and not on master, five of 4.5.2.1 and four of 4.5.2.3, and none the other way. On the 6,363 both answer the time is 1,686 seconds on master and 1,477 here.
The harness counts no answer wrong in either. The nine problems the builds disagreed on, run again one build at a time: master answers none, this all nine.
The suite passes: 15,176 passed, 13 skipped, none failed. The allocation gate passes: every gated benchmark allocates what the baseline says. The library builds for every target.
🤖 Generated with Claude Code
https://claude.ai/code/session_012sonx8iAspMiwRwokT1Ura