9

我有一些从 Java 移植到 C++ 的代码

// since this point is a vector from (0,0,0), we can just take the
// dot product and compare
double r = point.dot(normal);
return (r>=0.0);

但在 C++r中，可以是+0.0or -0.0，当r等于-0.0时，检查失败。

我试图在下面的代码中调整负零，但它从未到达 DEBUG("Negative zero") 行。但是r2打印出来等于+0.0.

// since this point is a vector from (0,0,0), we can just take the
// dot product and compare
double r = point.dot(normal);
if (std::signbit(r)){
    double r2 = r*-1;
    DEBUG("r=%f r=%f", r,r2);
    if (r2==0.0) {
        DEBUG("Negative zero");
        r = 0.0; //Handle negative zero
    }
}
return (r>=0.0);

有什么建议吗？

测试代码：

DEBUG("point=%s", point.toString().c_str());
DEBUG("normal=%s", normal->toString().c_str());
double r = point.dot(normal);
DEBUG("r=%f", r);
bool b = (r>=0.0);
DEBUG("b=%u", b);

测试结果：

DEBUG - point=Vector3D[ x=1,y=0,z=0 ]
DEBUG - normal=Vector3D[ x=0,y=-0.0348995,z=0.0348782 ]
DEBUG - r=0.000000
DEBUG - b=1
DEBUG - point=Vector3D[ x=1,y=0,z=0 ]
DEBUG - normal=Vector3D[ x=-2.78269e-07,y=0.0174577,z=-0.0174391 ]
DEBUG - r=-0.000000
DEBUG - b=0

海合会：

Target: x86_64-linux-gnu
--enable-languages=c,c++,fortran,objc,obj-c++ 
--prefix=/usr 
--program-suffix=-4.6 
--enable-shared 
--enable-linker-build-id 
--with-system-zlib 
--libexecdir=/usr/lib 
--without-included-gettext 
--enable-threads=posix 
--with-gxx-include-dir=/usr/include/c++/4.6 
--libdir=/usr/lib 
--enable-nls 
--with-sysroot=/ 
--enable-clocale=gnu 
--enable-libstdcxx-debug 
--enable-libstdcxx-time=yes 
--enable-gnu-unique-object 
--enable-plugin 
--enable-objc-gc 
--disable-werror 
--with-arch-32=i686 
--with-tune=generic 
--enable-checking=release 
--build=x86_64-linux-gnu 
--host=x86_64-linux-gnu --target=x86_64-linux-gnu
Thread model: posix
gcc version 4.6.3 (Ubuntu/Linaro 4.6.3-1ubuntu5)

旗帜：

CXXFLAGS += -g -Wall -fPIC

回答：

我已经使用@amit 的答案来执行以下操作。

return (r>=(0.0-std::numeric_limits<double>::epsilon()));

这似乎有效。

4

6 回答 6

10

好吧，使用doubles 时的一般建议是记住它们并不准确。因此，如果平等很重要 - 通常建议使用一些容差因子。

在你的情况下：

if (|r - 0.0| >= EPSILON)

EPSILON您的公差因子在哪里，如果r不是 0.0，将产生 true，至少有EPSILON间隔。

于 2012-12-07T17:12:09.200 回答
10

在某些较旧的系统（即 IEE754 之前）上，您可能会发现针对 0 的相等检查因负 0 失败：

if (a == 0.0) // when a==-0.0, fails

您可以通过在比较之前将 0.0 添加到值来解决此问题：

if ((a+0.0) == 0.0) // when a == -0.0, succeeds

但是，我要提醒的是，真正需要这样做的硬件/软件组合是非常不寻常的。我最后一次必须这样做是在 Control Data 大型机上。即使在那里，它也只是在一些不寻常的情况下出现：Fortran 编译器允许生成负零，并且知道在比较中补偿它们。作为计算的一部分，Pascal 编译器生成了将负零转换为正常零的代码。

因此，如果您在 Fortran 中编写了一个例程并从 Pascal 中调用它，您可能会遇到这个问题，并通过在进行比较之前添加 0.0 来防止它如上所述。

不过，我很有可能您的问题并不是真正源于与负零的比较。我所知道的所有相当现代的硬件都会完全自动处理这个问题，所以软件根本不需要考虑它。

于 2012-12-07T17:50:54.497 回答
5

OPstd::signbit在他的问题中提到，但它应该在一个可以快速找到的答案中，所以这里是：

这从 C++11 开始就可以区分-0.0和+0.0：

#include <cmath>
std::signbit(x) // true iif sign bit of x is set i.e. x is negative

cppreference 上的 std::signbit

于 2020-08-17T13:40:18.560 回答
3

大概你的意思是if (r2==-0.0). 尽管如此，负 0 和正 0 将比较相等。出于所有意图和目的，两者之间没有区别。对于负 0，您可能不需要特殊情况。对于负 0或正 0，您的比较r >= 0 应该为真。

于 2012-12-07T17:13:50.550 回答
1

考虑：

#include <stdio.h>
#include <iostream>
using namespace std;

int main()
{
    double const x = -2.78269e-07;

    printf( "printf: x=%f\n", x );
    cout << "cout: x=" << x << endl;
}

结果（使用 Visual C++ 11.0）：

[D:\开发\测试]
> cl foo.cpp
foo.cpp

[D:\开发\测试]
>富
printf: x=-0.000000
cout: x=-2.78269e-007

[D:\开发\测试]
> _

这似乎与问题中的神秘结果非常相似。

我认为它像问题的结果一样嘎嘎作响，看起来像问题的结果，像问题的结果一样摇摆不定。

所以，我相信未显示的计算代码产生了 value 2.78269e-007。

因此，总而言之，显然只有所描述的行为，即“负零”比较不等于零，这是非标准的。实际上，显然没有负零，只有一个非常小的负值。对于给定的输出格式，它显示为前面带减号的全零数字。

于 2012-12-07T17:56:13.243 回答
0

这只是一个很小的负数，控制台无法为您打印它
您可以尝试根据您的精度检查它，以便您可以用纯零替换它。

std::cout << std::setprecision(7) << (abs(value) < 0.0000005f ? 0 : value);

请注意我是如何将 7 个位置添加到我指定的浮动精度中的std::setprecision()。

这取决于您想要打印浮点数/双精度的准确度。

于 2018-07-19T09:10:03.323 回答