Clang 9下__fp16转float链接失败问题求助

你遇到的这个链接错误，是因为Clang在处理__fp16到float的转换时，需要依赖外部库提供的半精度浮点转换函数。这里有几个实用的解决办法：

1. 添加链接数学库 -lm #

__gnu_h2f_ieee是GNU实现的半精度转单精度浮点的函数，它通常包含在数学库libm中。你只需要在链接命令里加上-lm就能解决这个问题。

修改后的链接命令示例：

&& /usr/bin/clang++-9 -g -fsanitize=address,undefined -fno-omit-frame-pointer -fno-limit-debug-info CMakeFiles/run_model.dir/src/run_model.cc.o -o run_model libfused_transformer.a ../thirdparty/OpenBLAS/libopenblas.a ../thirdparty/icu/icu4c/linux/prebuilt/lib/libicui18n.a ../thirdparty/icu/icu4c/linux/prebuilt/lib/libicuuc.a ../thirdparty/icu/icu4c/linux/prebuilt/lib/libicudata.a -lpthread /usr/lib/llvm-9/lib/libomp.so -lpthread -lm && :

把-lm放在链接命令的末尾即可，确保编译器能找到该库的符号。

2. 添加Clang特定编译选项 -ffp16-format=ieee #

有时候编译器默认可能没有正确识别IEEE格式的半精度浮点，加上-ffp16-format=ieee编译选项，可以让Clang直接生成转换的机器码，无需依赖外部库函数。你可以把这个选项加入编译flags，比如在CMake中设置：

set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -ffp16-format=ieee")

3. 手动实现转换函数（作为备选方案） #

如果上面两个方法都不生效，你可以自己实现一个__gnu_h2f_ieee的替代函数，核心是把binary16的位模式转换成float的位模式，代码示例如下：

float __gnu_h2f_ieee(unsigned short h) {
    unsigned int sign = (h >> 15) & 0x1;
    unsigned int exponent = (h >> 10) & 0x1F;
    unsigned int mantissa = h & 0x3FF;

    unsigned int f_sign = sign << 31;
    unsigned int f_exponent;
    unsigned int f_mantissa;

    if (exponent == 0) {
        if (mantissa == 0) {
            // 处理零值
            return *reinterpret_cast<float*>(&f_sign);
        } else {
            // 处理亚正规数
            int shift = __builtin_clz(mantissa) - 22;
            exponent = 1 - shift;
            mantissa <<= shift;
            f_exponent = (exponent + 127 - 15) << 23;
            f_mantissa = (mantissa & 0x3FF) << 13;
        }
    } else if (exponent == 0x1F) {
        // 处理无穷大或NaN
        f_exponent = 0xFF << 23;
        f_mantissa = mantissa << 13;
    } else {
        // 处理正规数
        f_exponent = (exponent + 127 - 15) << 23;
        f_mantissa = mantissa << 13;
    }

    unsigned int f = f_sign | f_exponent | f_mantissa;
    return *reinterpret_cast<float*>(&f);
}

把这段代码加入你的项目中，链接时就不会出现未定义引用的错误了，这个实现完全遵循IEEE 754 binary16的转换规则，和标准函数行为一致。

优先推荐前两个方法，尤其是添加-lm，这是最直接简单的解决方案。

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