复现std::complex存疑：无穷大除法运算结果为何不一致？

问题场景 #

正在为HLSL库实现std::complex，基础算术运算正常，但处理无穷大时出现异常。以下是Godbolt上的测试代码：

标准库/=实现测试 #

std::complex<float> a(-3.0, 4.0);
std::complex<float> inf(1.0f / 0.0f, 0.0f);
a /= inf;
std::cout << a << std::endl;

输出：(-0,0)，符合预期。为复刻逻辑，复制了GCC版本的operator/=源码。

手动实现测试 #

std::complex<float> a(-3.0, 4.0);
std::complex<float> inf(1.0f / 0.0f, 0.0f);
float r = a.real() * inf.real() + a.imag() * inf.imag();
float n = std::norm(inf);
float imag = (a.imag() * inf.real() - a.real() * inf.imag()) / n;
float real = r / n;
std::cout << std::complex<float>(real, imag) << std::endl;

输出：(-nan, -nan)，与标准库结果不符。

问题原因 #

手动实现直接套用了复数除法的通用公式，但忽略了标准库对**特殊数值（无穷大、NaN）**的分支处理：

1. 当除数是实部无穷大、虚部为0的复数时，std::norm(inf)的结果是inf * inf = inf
2. 分子r = -3.0 * inf + 4.0 * 0 = -inf，-inf / inf的结果为NaN；虚部计算同理，最终输出NaN

而GCC的std::complex::operator/=会先做特殊数值检测：
当识别到除数是纯实无穷大时，直接将被除数的实部、虚部分别除以除数的实部，即-3.0/inf和4.0/inf，这两个运算的结果分别是-0和0，与预期一致。

解决思路 #

你复制的源码大概率只包含通用情况的实现，缺少特殊数值的前置检查。在HLSL中实现时需要补充这些分支：

• 优先判断除数是否为纯实无穷大/纯虚无穷大，直接执行对应除法
• 补充被除数为零、除数为NaN等其他特殊场景的处理，确保数值行为符合C++标准

内容的提问来源于stack exchange，提问作者Francisco José Letterio