MinGW GCC中atan2/cos等函数用x87库，如何切换至AVX2/SSE2？

问题背景 #

将Intel 2023及MSVC Visual C++ 2022下的数值代码移植到Windows平台MinGW GCC 13.1.0时，发现编译后atan2、cos、exp、sin、sqrt等函数会调用x87 80位浮点库实现，单调用耗时约100周期，远慢于SSE/AVX2版本的20-30周期；但tan、atan却正常使用AVX2指令。已尝试-O3、-Ofast、-march=native、-mavx2、-mfpmath=sse等编译选项，问题仍存在。当前临时方案是集成BSD三角函数库，但不够优雅，希望找到强制使用AVX2/SSE2浮点库或重新编译默认数学库的方法。

强制使用AVX2/SSE2数学函数的直接方法 #

1. 优化编译参数组合 #

在编译时补充更明确的ABI参数，确保GCC优先选择SIMD版本的数学库：

-O3 -march=native -mavx2 -mfpmath=sse -mabi=ms

-mabi=ms强制使用Windows系统调用约定，避免因ABI不匹配导致GCC fallback到x87实现。

2. 显式调用SIMD专用函数 #

MinGW的数学库中部分函数提供了SSE/AVX2专用变体，可直接在代码中调用：

• SSE2版本：单精度用atan2f、cosf、expf、sinf、sqrtf；双精度对应atan2l、cosl等变体
• AVX2版本：包含头文件<immintrin.h>后，使用_mm256_atan2_pd、_mm256_cos_pd等内置函数

注意：直接调用内置函数需手动处理向量打包/解包逻辑，适合性能敏感的热点代码段。

3. 替换第三方优化数学库 #

替换MinGW默认的libm库为SIMD优化的第三方库：

• Intel SVML：通过编译选项-msvml启用，GCC会链接Intel短向量数学库，自动优先使用SSE/AVX2实现
• AMD LibM：针对AMD架构优化的数学库，编译时指定链接该库即可

重新编译MinGW默认数学库的步骤 #

如果需要自定义默认数学库实现，可按以下操作：

1. 获取MinGW-w64的libm源码（通常包含在mingw-w64-src包中）
2. 修改编译配置，添加-mavx2 -mfpmath=sse等参数，强制生成SIMD优化代码
3. 编译生成新的libm.a/libm.dll.a库文件
4. 项目编译时，指定链接路径优先使用自定义编译的libm库

注意：重新编译库需熟悉MinGW编译链，且要保证新库与系统ABI兼容，避免运行时错误。

额外注意事项 #

• 检查代码中是否存在long double这类80位浮点类型，该类型会强制GCC使用x87实现，需统一替换为double或float
• -Ofast已包含-ffast-math，该选项允许GCC进行激进浮点优化，包括优先选择SIMD数学函数实现

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