我已经在许多正弦/余弦的实现中看到了所谓的扩展模块化精度算法.但是它的用途是什么？例如,在cephes实现中,在缩减到[0,pi/4]范围后,他们正在进行这种模块化精度算法以提高精度.
下面的代码:
z = ((x - y * DP1) - y * DP2) - y * DP3;
其中DP1,DP2和DP3是一些硬编码系数.如何在数学上找到这些系数？我已经理解了大数字"模块化扩展算术"的目的,但这里的确切目的是什么？