小编Sco*_*ttD的帖子

sinl当参数接近pi的非零倍数时，为什么给出不正确的结果？sinl当参数很大时，为什么会给出错误的结果？以下代码说明了这一点。

注意，用于初始化变量pi的数字与任何64位long double值都不完全匹配。编译器选择最接近的值3.14159265358979323851280895940618620443274267017841339111328125。可以使用libquadmath，gnu MPFR lib或在线计算器（例如http://www.ttmath.org/online_calculator）找到期望的正弦值。

#include <stdio.h>
#include <math.h>

int main (int argc, char *argv [])
    {
    volatile long double pi = 3.14159265358979323846L;
    volatile long double big = 9223372035086174241L;
    volatile long double expected1 = -5.0165576126683320235E-20L;
    volatile long double expected2 = -4.2053336735954077951E-10L;
    double result;
    double ex1 = expected1, ex2 = expected2;

    result = sinl (pi);
    printf("expected: %g, \nreturned: %g\n\n", ex1, result);
    result = sinl (big);
    printf("expected: %g, \nreturned: %g\n\n", ex2, result);
    return 0;
    } …