The derivative of an answer at a large hyperbolic argument evaluates to NaN or loses its imaginary part
I maintainer di solito rispondono entro 1 giorno
Nessuno ha ancora preso questa issue.
Valutazione
- Difficoltà
- 4/5
- Tempo stimato
- 3-5 giorni
- Idoneità per principianti
- 43/100
Direzione di ricerca
Start by reproducing the derivative evaluation with the exact parameters and x from the issue, comparing it with the integrand at 50, 300, and 2000 digits. Trace the derivative evaluation path and identify which subexpression causes NaN or loses the imaginary component; the issue does not name a file or test. Done means the derivative agrees with the integrand at the reported large argument and at x = 2 across those precisions.
Scritto dal modello di indicizzazione a partire dal testo della issue.
Descrizione
Found while checking Rubi 6.1.5, row 94 (1/sqrt(a + i a sinh(c + d x))) for 2.6.0. The answer is right. The library's own numerical evaluation of its derivative is not, at large arguments.
With a, c, d = 21/10, 19/10, 17/10 and x = 227283/1250 (so t = c + d x ≈ 311), the true value of the derivative is r(1 − i) with r ≈ 2.0185e-68. That is what mpmath gives at dps 400, and what the library gives when it evaluates the integrand itself.
The derivative of the answer, Differentiate(x) then substituted with exact rationals and evaluated, gives:
| digits | value |
|---|---|
| 50 | NaN |
| 300 | 1.5329…, which is wrong in magnitude too |
| 2000 | 2.0185…e-68: the real part is right and the imaginary part −2.0185…e-68 i is missing |
The relative difference at 2000 digits is exactly −1/2 + i/2. At x = 2 all three precisions agree with the integrand to the last digit.
The derivative holds e^(±t/2) and e^(±t) beside a + i a sinh(t), so large cancellations are expected. Losing a whole component at 2000 digits is not explained by that, though, and NaN at 50 digits should have been "could not evaluate". I have not narrowed it to a single subexpression.
Because of this, the Rubi harness reports this row as wrong, and anything that checks an answer by evaluating would do the same.
🤖 Generated with Claude Code
- Lingua principale
- C#
- Stelle
- 831
- Fork
- 79
- Merge medio
- 1h 24m
- PR unite (30g)
- 540
Preparare l'ambiente
- Nessun Dockerfile né file Docker Compose
- Nessun modello di pull request
- Leggi la guida per i contributori
Come iniziare
- Leggi tutta la issue e poi la guida ai contributi del progetto.
- Commenta sulla issue per dire che te ne occupi tu — evita che due persone facciano lo stesso lavoro.
- Fai un fork del repository e lavora su un branch.
- Apri una pull request che faccia riferimento al numero della issue.
Altre issue di ASC-Community/AngouriMath
-
Opinions wanted
Difficoltà 5/5 Più di una settimana Idoneità per principianti 25/100
ASC-Community/AngouriMath#2120 ·
I maintainer di solito rispondono entro 1 giorno
-
Difficoltà 5/5 Più di una settimana Idoneità per principianti 30/100
ASC-Community/AngouriMath#1807 · 3 commenti ·
I maintainer di solito rispondono entro 1 giorno
-
Difficoltà 5/5 Più di una settimana Idoneità per principianti 25/100
ASC-Community/AngouriMath#1692 ·
I maintainer di solito rispondono entro 1 giorno
-
Difficoltà 5/5 Più di una settimana Idoneità per principianti 30/100
ASC-Community/AngouriMath#1690 ·
I maintainer di solito rispondono entro 1 giorno
-
Difficoltà 5/5 Più di una settimana Idoneità per principianti 25/100
ASC-Community/AngouriMath#1649 · 36 commenti ·
I maintainer di solito rispondono entro 1 giorno
Tutte le issue di ASC-Community/AngouriMath
Issue simili
-
fix(testing): HeadlessTerminal disposes services synchronouslyForse già presa @zidad l’ha presa oggi. Aperta
Difficoltà 2/5 1-3 ore Idoneità per principianti 69/100
RazorConsole/RazorConsole#388 ·
I maintainer di solito rispondono entro 1 giorno
-
A practical example of rulerAperta
Difficoltà 2/5 1-3 ore Idoneità per principianti 65/100
-
Difficoltà 2/5 1-3 ore Idoneità per principianti 72/100
starlinghq/starling#209 ·
I maintainer di solito rispondono entro 1 giorno
-
affected/v17 area/frontend type/bug
Difficoltà 2/5 1-3 ore Idoneità per principianti 63/100
umbraco/Umbraco-CMS#24323 · 1 commento ·
I maintainer di solito rispondono entro 2 giorni
-
Bug
Difficoltà 2/5 1-3 ore Idoneità per principianti 73/100
dotnet/roslynator#1859 ·
I maintainer di solito rispondono entro 1 giorno