VC ++使用fp:fast导致错误(不仅仅是不准确)结果 - 这是编译器错误吗?

时间:2017-11-23 22:51:47

标签: c++ visual-c++ floating-point visual-studio-2017 compiler-bug

我已经安装了最新的VS2017更新(15.4.4),在编译我们的项目时,单元测试开始失败。 在使用优化(/ O2)和浮点快速模型(/ fp:fast)时,似乎在某些情况下会出现问题。 以前的编译器(VS2017更新15.2)没有发生此问题。

这是一个示例程序:

#include <iostream>

const float FACTOR = 0.01745329251994329576923690768489f;

unsigned long long hoursToMicrosecs(int hours)
{
    return hours * 3600 * 1000000LL;
}

float degToRad(float deg)
{
    return deg * FACTOR;
}

float f(int u1, int u2)
{
    float percent = ((u1 - u2) / 1000.0f) / (hoursToMicrosecs(24) / 1000000.0f);
    return degToRad(percent * 360);
}

int main()
{   
    auto result = f(3600000, 35063681);
    std::cout << result << '\n';
    return (result > -3.0f) ? 0 : -1;
}

result应该是-2.2881,但输出实际上是-394.868,这不仅是不准确的。

如果我执行以下任何操作,它会起作用:

  • 删除优化
  • 更改为fp:exact
  • 返回上一个编译器(15.2)

查看反汇编表明编译器尝试为我们做了一些不错的事情 - 它只是在编译时计算了整个事情并用一个数字替换。

优化的代码只是一行代码:

011F1000  vmovss      xmm0,dword ptr [__real@c3c56f11 (011F2118h)]  

我的问题是:这是一个编译器错误(我应该报告)还是错误使用fp:fast?

1 个答案:

答案 0 :(得分:6)

我在这看到同样的事情。我想我发现了这个问题:hoursToMicroseconds(24)正在泛滥。它不应该,因为表达式应该在与unsigned long long相乘之前转换为1000000LL。但是如果你把结果反馈回unsigned int,你将得到5006540​​80而不是86400000000.所以整个计算最终归结为-394.868。

我肯定会说这是一个编译错误。您似乎可以通过先将unsigned long long结果转换为double来避开它。

修改

你知道还有什么好笑的吗?如果您完成所有功能constexpr,并在result中将constexpr声明为main(),则生成正确的结果。因此,编译器中显然有两个独立的代码路径可以在编译时计算值,其中一个被破坏。