使用std::complex<_Float16>乘法时__constexpr_copysign调用歧义问题求助

问题描述 #

使用std::complex<_Float16>类型执行乘法运算时，编译器报错：call to '__constexpr_copysign' is ambiguous。错误触发于实例化std::operator*<_Float16>模板阶段，编译器列出float、double、long double版本的__constexpr_copysign候选函数，无法确定匹配的重载。

可复现代码：

#include <complex>

int main(){
    std::complex<_Float16> a(10.0, 11.2);
    std::complex<_Float16> b(3.0, 4.0);
    std::complex<_Float16> c = a * b;
    return 0;
}

解决方案 #

1. 补充_Float16版本的__constexpr_copysign重载 #

在代码中为_Float16类型添加__constexpr_copysign的重载，让编译器能匹配到精确类型的函数：

#include <complex>
#include <cmath>

// 为_Float16补充__constexpr_copysign重载
namespace std {
    constexpr _Float16 __constexpr_copysign(_Float16 magnitude, _Float16 sign) noexcept {
        // 使用GCC内置的_Float16版copysign实现，其他编译器需替换对应内置函数
        return __builtin_copysignf16(magnitude, sign);
    }
}

int main(){
    std::complex<_Float16> a(10.0, 11.2);
    std::complex<_Float16> b(3.0, 4.0);
    std::complex<_Float16> c = a * b;
    return 0;
}

2. 手动实现std::complex<_Float16>的乘法运算符 #

绕过标准库的乘法实现，手动编写符合复数乘法规则的运算符重载：

#include <complex>

// 手动实现_Float16复数的乘法
std::complex<_Float16> operator*(const std::complex<_Float16>& lhs, const std::complex<_Float16>& rhs) {
    const _Float16 real_part = lhs.real() * rhs.real() - lhs.imag() * rhs.imag();
    const _Float16 imag_part = lhs.real() * rhs.imag() + lhs.imag() * rhs.real();
    return {real_part, imag_part};
}

int main(){
    std::complex<_Float16> a(10.0, 11.2);
    std::complex<_Float16> b(3.0, 4.0);
    std::complex<_Float16> c = a * b;
    return 0;
}

3. 升级编译器版本 #

该问题多源于旧版本编译器对_Float16类型的std::complex支持不完善。升级到GCC 12+、Clang 15+等较新版本，编译器已完善半精度浮点类型的标准库函数支持，可直接解决歧义问题。

问题根源 #

标准库中std::complex的乘法实现依赖于__constexpr_copysign这类内部函数，但旧版本编译器未提供_Float16类型的重载实例。当模板实例化时，编译器只能在float、double、long double版本中选择，因类型转换存在歧义而报错。

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