C ++：正弦函数的工作

Question

我试图用C ++绘制一条正弦曲线并遇到一些有趣的东西。

我有一个函数，它返回一个以度为单位的值的正弦值

double sind(double a)
{
    return sin(a*3.14159/180);
}

现在在主要功能

sind(18)==sind(18)?cout<<1:cout<<0;

我写了上面的代码。结果似乎是假的，它在控制台上打印0。 但据我说，罪（18）和罪（18）是平等的。那么计算机的想法是怎么回事？

另外，如果我想检查两个正弦值的相等性，我该怎么做呢？

Answer 1

在PC上，至少从旧版开始，浮点值是用80位计算的，但是对于主存储器则向下舍入到64位。当编译器认识到它可以重用80位结果进行80位比较时，你会得到令人困惑的结果。是的，它是神圣标准所允许的。

---

顺便说一句，void main无效。这意味着代码具有未定义的行为，原则上可以做任何事情，包括什么都不做，或者你期望什么。在实践中它并没有那么糟糕，它只是使代码不可移植，但仍然不能这样做：添加一个字符以使代码不可移植是愚蠢的，所以写int main。

Answer 2

  

另外，如果我想检查两个正弦值的相等性，我该怎么做呢？

有一个非常好的浮点示例＆＃34; almost_equal＆＃34;函数here

按要求，从原始来源复制的代码：

#include <cmath>
#include <limits>
#include <iomanip>
#include <iostream>
#include <type_traits>
#include <algorithm>

template<class T>
typename std::enable_if<!std::numeric_limits<T>::is_integer, bool>::type
    almost_equal(T x, T y, int ulp)
{
    // the machine epsilon has to be scaled to the magnitude of the values used
    // and multiplied by the desired precision in ULPs (units in the last place)
    return std::abs(x-y) < std::numeric_limits<T>::epsilon() * std::abs(x+y) * ulp
    // unless the result is subnormal
           || std::abs(x-y) < std::numeric_limits<T>::min();
}
int main()
{
    double d1 = 0.2;
    double d2 = 1 / std::sqrt(5) / std::sqrt(5);

    if(d1 == d2)
            std::cout << "d1 == d2\n";
    else
            std::cout << "d1 != d2\n";

    if(almost_equal(d1, d2, 2))
            std::cout << "d1 almost equals d2\n";
    else
            std::cout << "d1 does not almost equal d2\n";
}