1

一个相当基本的问题,但我在这里有点迷失。我正在研究一个使用复数的库。到目前为止,我在 Mac OS 10.8 和 Ubuntu 10.04 下使用了我的实现(在 Man 上使用 llvm-gcc,在 Ubuntu 下使用 gcc)。昨天我在另一台带有icpc(icc/2013.3)的机器上试了一下,复数运算部分失败了。

我用这个简单的例子解决了这个问题:

#include <iostream>
#include <complex>

using namespace std;

int main(int argc, char* argv[])
{
    std::complex<double> tmpcmplxnumber(26.0,0.0);
    cout << "Number: " << tmpcmplxnumber << endl;
    std::complex<double> newcmplxnumber = 1.0 - tmpcmplxnumber;
    cout << "New Number: " << newcmplxnumber << endl;

    cout << "Sqrt of New Number: " << sqrt(newcmplxnumber) << endl;

    return 0;
}

gcc 输出(两台机器):

Number: (26,0)
New Number: (-25,-0)
Sqrt of New Number: (0,-5)

icpc输出:

Number: (26,0)
New Number: (-25,0)
Sqrt of New Number: (0,5)

如果我被正确告知 gcc 实现了 ISO 标准,那么结果应该是正确的(所有符号为零的东西......等等)那么 icpc 实现是“错误的”还是有针对该特定问题的标志?

更新 我知道从数学的角度来看这两个结果都是有效的,但是在使用 icpc 时我也需要 gcc 结果。

4

1 回答 1

1

sqrt通常是非负平方根,在数学中用“平方根”符号表示。

根据(草案)标准(26.4.8,复杂超越),

如果参数是负实数,则返回的值位于正虚轴上。

所以看起来 gcc 在这里有问题。

于 2013-05-24T08:51:02.187 回答