The derivative of an answer at a large hyperbolic argument evaluates to NaN or loses its imaginary part
メンテナーはふだん 1 日以内に返信
まだ誰も着手していません。
評価
調査の方向性
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.
索引モデルが issue の本文から書いたものです。
説明
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
- 主要言語
- C#
- スター
- 831
- フォーク
- 79
- 平均マージ
- 1時間 27分
- マージ済み PR（30日）
- 522
環境構築
- Dockerfile・Docker Compose ファイルなし
- プルリクエストのテンプレートなし
- コントリビューションガイドを読む
はじめの一歩
- issue を最後まで読み、次にプロジェクトのコントリビューションガイドを読みます。
- 着手することを issue にコメントします — 二人が同じ作業をするのを防げます。
- リポジトリをフォークし、ブランチを切って変更します。
- issue 番号を参照したプルリクエストを送ります。
ASC-Community/AngouriMath のほかの issue
-
難易度 4/5 3〜5日 初心者へのやさしさ 42/100
ASC-Community/AngouriMath#2081 · コメント 1 件 ·
メンテナーはふだん 1 日以内に返信
-
難易度 5/5 1週間以上 初心者へのやさしさ 35/100
ASC-Community/AngouriMath#2075 · コメント 1 件 ·
メンテナーはふだん 1 日以内に返信
-
難易度 5/5 1週間以上 初心者へのやさしさ 30/100
ASC-Community/AngouriMath#1807 · コメント 2 件 ·
メンテナーはふだん 1 日以内に返信
-
難易度 5/5 1週間以上 初心者へのやさしさ 25/100
ASC-Community/AngouriMath#1692 ·
メンテナーはふだん 1 日以内に返信
-
難易度 5/5 1週間以上 初心者へのやさしさ 30/100
ASC-Community/AngouriMath#1690 ·
メンテナーはふだん 1 日以内に返信
ASC-Community/AngouriMath の issue をすべて見る
似ている issue
-
devtools
難易度 2/5 1〜3時間 初心者へのやさしさ 82/100
メンテナーはふだん 1 日以内に返信
-
clawsweeper:fix-shape-clear clawsweeper:needs-maintainer-review clawsweeper:no-new-fix-pr clawsweeper:not-repro-on-main impact:ux-release-blocker issue-rating: 🦪 silver shellfish P0
難易度 2/5 1〜3時間 初心者へのやさしさ 72/100
openclaw/openclaw-windows-node#1843 · コメント 1 件 · リアクション 1 件 ·
メンテナーはふだん 1 日以内に返信
-
難易度 2/5 1〜3時間 初心者へのやさしさ 66/100
-
kind:ux size:S type:bug
難易度 2/5 1〜3時間 初心者へのやさしさ 78/100
elsa-workflows/elsa-foundation#3256 ·
メンテナーはふだん 1 日以内に返信
-
[Bug]: Revit selector reads process ids 2000-2099 as years (flaky AgentClientTests.AYearOrAProcessIdSelectsTheHost)対応中かも @Wurschdhaud が今日担当しました。 オープン
難易度 2/5 1〜3時間 初心者へのやさしさ 83/100
pyrevitlabs/pyRevit#3817 ·
メンテナーはふだん 1 日以内に返信