为何C语言中carg(-0)返回值为-M_PI？

在C语言中使用double complex类型时，发现将带有负0虚部的实数传入carg()函数时，返回值为-M_PI而非0或M_PI，这是为什么？该行为影响了newlib库中的cpow函数，导致出现虚部符号错误。

最小复现代码 #

#include <stdio.h>
#include <complex.h>
#include <math.h>

int main(void) {

    double complex a = 2;

    printf("creal(a): %f\n", creal(a));
    printf("cimag(a): %f\n", cimag(a));
    printf("carg(a): %f\n", carg(a));

    printf("\n");

    printf("creal(-a): %f\n", creal(-a));
    printf("cimag(-a): %f\n", cimag(-a));
    printf("carg(-a): %f\n", carg(-a));

    printf("\n");
    a = 0;
    printf("creal(-a): %f\n", creal(-a));
    printf("cimag(-a): %f\n", cimag(-a));
    printf("carg(-a): %f\n", carg(-a));

}

运行结果（gcc 14.2.0编译） #

creal(a): 2.000000
cimag(a): 0.000000
carg(a): 0.000000

creal(-a): -2.000000
cimag(-a): -0.000000
carg(-a): -3.141593

creal(-a): -0.000000
cimag(-a): -0.000000
carg(-a): -3.141593

原因解释 #

1. IEEE 754浮点数的正负0特性：IEEE 754标准明确区分正0和负0，负0是合法的浮点数表示，对正0取负会得到负0。代码中-a的虚部是-cimag(a)，原a的虚部为正0，取负后就变成了负0。

2. carg()函数的标准行为：C语言标准（如C11及后续版本）和IEEE 754规定，carg(z)返回复数z的辐角，范围在[-π, π]之间：

   • 当虚部y为正（包括正0）时，辐角符号为正；当y为负（包括负0）时，辐角符号为负。
   • 即使实部x为负数，只要虚部是负0，carg()会返回-π而非π——负0表示虚部从下半平面趋近于0，对应复平面中第四象限方向趋近负实轴，因此辐角取-π；若虚部是正0，则是从第二象限趋近负实轴，辐角取π。

3. 特殊情况：-0 + (-0)i：当复数为-0 + (-0)i时，carg()返回-π，这同样符合标准——虚部的负0决定了辐角的符号，即使实部也是负0。

这种行为是标准规定的，并非实现bug，它保证了复数运算的连续性和符号一致性，尤其是在涉及分支切割（比如幂函数cpow）的场景中，符号的正确性会直接影响计算结果。

内容的提问来源于stack exchange，提问作者Leo