多精度初等函数实现

Question

我试图找到sin函数的多精度版本的实现。我想通读mpfr库（C/C++源代码），但我找不到与实际实现相关的任何代码。

我基本上想知道的是如何在这种情况下执行参数缩减（多精度），这可能是一个巨大的LUT被存储来代表pi的倒数，这将允许处理所有可能的情况与案件我的意思是尾数大小）。

你能指出我看什么具体的来源？代码阅读起来相当复杂。

Answer 1

范围缩减不是用一个巨大的LUT实现的：这是不可能的，因为指数范围太大。所做的只是以足够的精度计算pi; pi的值被缓存以用于下一次对三角函数的调用，并且每次需要更多精度时重新计算pi。

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