Complete loss of precision for trig reduction in non-default rounding modes
まだ誰も着手していません。
評価
調査の方向性
まず rem_pio2f の中程度の引数のパスと rem_pio2 のフォールバックから始め、次に round-to-zero と round-down で sinf(-0x1.4665d2p+27) を再現します。結果を記載された正確な値と比較し、各丸めモードが三角関数のリダクションにどのように影響するかを追跡します。round-to-nearest 以外でもリダクションが正確な結果を生成し、報告されたケースのリグレッションカバレッジがあることが完了の条件です。
索引モデルが issue の本文から書いたものです。
説明
The trig-reduction implementations assume that the rounding mode is round-to-nearest, and do not function correctly if it isn't. Example: sinf(-0x1.4665d2p+27). The mathematically precise result is approximately -0x1.927bcc77af475p-25, but in round-toward-zero or round-down, openlibm computes a result of -0x1.a2c7eep-17.
This behavior is especially bad for the "medium argument" path in rem_pio2f. When the argument is so large that you punt to rem_pio2, double-precision has enough extra bits to still deliver relatively good results.
- 主要言語
- C
- スター
- 637
- フォーク
- 165
- PR マージ指標
- 30日以内にマージされた PR はありません
環境構築
このプロジェクトには開発コンテナ、Dockerfile、コントリビューションガイドがありません。まず README を読み、一般的な手順ははじめてのコントリビューションガイドを参照してください。
はじめの一歩
- issue を最後まで読み、次にプロジェクトのコントリビューションガイドを読みます。
- 着手することを issue にコメントします — 二人が同じ作業をするのを防げます。
- リポジトリをフォークし、ブランチを切って変更します。
- issue 番号を参照したプルリクエストを送ります。
JuliaMath/openlibm のほかの issue
-
nextafter and nexttoward return the finite endpoint in directed rounding modes対応中かも @vera-rubin が 4 日前に担当しました。 オープン
難易度 3/5 1〜2日 初心者へのやさしさ 68/100
-
難易度 4/5 3〜5日 初心者へのやさしさ 38/100
-
難易度 5/5 1週間以上 初心者へのやさしさ 25/100
-
難易度 4/5 3〜5日 初心者へのやさしさ 35/100
-
Non-standard behavior of tgamma()対応中かも @ViralBShah が 117 日前に担当しました。 オープン
難易度 3/5 1〜2日 初心者へのやさしさ 45/100
JuliaMath/openlibm の issue をすべて見る
似ている issue
-
難易度 2/5 1〜3時間 初心者へのやさしさ 70/100
OpenPrinting/cups#1772 ·
メンテナーはふだん 1 日以内に返信
-
難易度 2/5 1〜3時間 初心者へのやさしさ 88/100
snowballstem/snowball#311 ·
-
tech-debt
難易度 2/5 1〜3時間 初心者へのやさしさ 78/100
メンテナーはふだん 1 日以内に返信
-
難易度 2/5 1〜3時間 初心者へのやさしさ 68/100
メンテナーはふだん 1 日以内に返信
-
難易度 2/5 1〜3時間 初心者へのやさしさ 70/100