複数精度の基本関数の実装

Question

sin関数の倍精度バージョンの実装を探しています。 mpfrライブラリ（C/C++ソースコード）を読み込もうとしていますが、実際の実装に関連するコードは見つかりません。

このような場合に引数の削減がどのように実行されるかは、複数の精度で表現されています。piの逆数を表す巨大なLUTが格納されている可能性があります。場合によっては仮数サイズを意味する）。

具体的な情報源を教えてください。このコードは読み込みが非常に複雑です。

Answer 1

指数範囲が大きすぎるため、巨大なLUTでは範囲の縮小が実装されていません。これは不可能です。実行されるのは、十分な精度でpiを計算することです。次の三角関数の呼び出しのためにpiの値がキャッシュされ、さらに精度が必要になるたびにpiが再計算されます。

MPFRのソースを参照して、および/またはalgorithms.pdfをLaTeXのdocディレクトリからビルドして情報を入手できます。しかし、この文書は常に最新のものではなく、一部の情報が欠落している可能性があります。