c++ - 复数的平方根 - g++ vs icpc

Question

一个相当基本的问题，但我在这里有点迷路。我正在开发一个处理复数的库。到目前为止，我在 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 结果。

Answer 1

sqrt 通常是非负平方根，也就是数学中用“平方根”符号表示的平方根。

根据(草案)标准(26.4.8，复先验数)，

If the argument is a negative real number, the value returned lies on the positive imaginary axis.

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