c++ - 多精度初等函数实现

Question

我正在尝试寻找 sin 函数的多精度版本的实现。我正在尝试通读 mpfr 库(C/C++ 源代码)，但找不到与实际实现相关的任何代码。

我基本上想知道的是在这种情况下(多精度)如何执行参数缩减，可能会存储一个巨大的 LUT 来表示 pi 的倒数，这将允许处理所有可能的情况(我指的是尾数大小)。

您能指出要查看的具体来源吗？代码阅读起来相当复杂。

Answer 1

范围缩减不是通过巨大的 LUT 实现的:这是不可能的，因为指数范围太大。所做的只是足够准确地计算 pi； pi 的值被缓存以供下次调用三角函数时使用，并且每次需要更高的精度时都会重新计算 pi。

您可以通过查看 MPFR 的源代码和/或使用 LaTeX 从 doc 目录构建 algorithms.pdf 来获取信息。但此文档可能并不总是最新的，并且可能会丢失一些信息。