Fortran中Sin与Cos函数精度差异：为何Sin(π/2)为1而Cos非0？

先看一下用户提供的代码和运行结果：

PROGRAM TRIG_TEST
IMPLICIT NONE
DOUBLE PRECISION, PARAMETER :: PI=4.D0*DATAN(1.0D)
print *, sin(PI/2.0), cos(PI/2.0)
END PROGRAM TRIG_TEST

用gfortran编译运行后输出为：1.0000000000000000 6.1232339957367660E-017

这个问题问得特别到位，其实背后是浮点数表示限制和数学库的特殊优化逻辑共同作用的结果，我来一步步拆解：

首先：你的PI不是精确的数学π #

你用4.D0*DATAN(1.0D)定义的PI，是双精度浮点数能表示的最接近真实π的数值，但它和数学上的π之间依然存在一个极小的误差（大约是1.2e-16）。也就是说，PI/2.0也不是精确的π/2，而是一个非常接近它的浮点数，我们暂时叫它x。

为什么sin(x)返回精确的1.0？ #

绝大多数Fortran编译器（包括gfortran）依赖的数学库（比如glibc的数学库）会对特殊三角函数值做针对性优化：

• 当检测到输入值x落在π/2的一个极小误差范围内时，库函数会直接返回精确的1.0——毕竟这是数学上sin(π/2)的精确结果，这么做既符合用户的直觉，也不会引入额外的计算误差。
• 简单说，库知道你想算的是sin(π/2)，哪怕输入有一点点偏差，也直接给你正确的精确值。

为什么cos(x)返回的是6.12e-17而不是0？ #

这就要结合三角函数的特性和浮点数运算的敏感性来看：

1. cos在π/2附近的敏感性：数学上，当ε极小时，cos(π/2 + ε) ≈ -ε。你的x = PI/2.0比真实的π/2大了约6e-17（因为PI和真实π的差是1.2e-16，除以2就是这个数），所以cos(x)的结果自然就是这个微小偏差的近似值，也就是你看到的6.1232339957367660E-017。
2. 库的优化策略差异：数学库没有把这个微小值近似为0，是因为：
   • 这个结果本身就是浮点数运算的真实输出，反映了输入和精确π/2的偏差；
   • 如果强行返回0，反而会掩盖浮点数近似带来的误差，对于一些需要高精度的场景可能会有问题；
   • 相比sin(π/2)的结果是明确的1，cos(π/2)的0在浮点数中是一个“可表示但需要反映输入误差”的值，所以库选择保留这个微小的结果。

举个直观的例子：你用一把精度极高的尺子量π/2，虽然量出来的结果和真实值只差一根头发丝的万分之一，但cos对这个偏差的反应就是这么灵敏，直接把这个微小偏差放大成了输出里的那个小数值。

内容的提问来源于stack exchange，提问作者user1887919