ARM Cortex-M4F下避免sqrtf调用双精度函数的替代方案咨询

问题场景 #

给搭载fpv4-sp-d16、使用hard ABI的ARM Cortex-M4F写代码时，发现只要用了sqrtf且没开-ffast-math，ELF文件里就会出现__aeabi_dmul这类双精度函数调用——反汇编显示是sqrtf内部调用了__aeabi_ddiv。不想开-ffast-math（怕带来其他未知影响），得找个纯单精度的开方实现方式。

环境参数 #

• 编译器：ARM-GCC 10 2021.10
• 编译选项：-mcpu=cortex-m4 -mthumb -mfloat-abi=hard -mfpu=fpv4-sp-d16 -Os -fsingle-precision-constant -Wfloat-equal -Wdouble-promotion

踩过的无效尝试 #

试过写个带optimize("fast-math")属性的sqrtf_wrap函数包装sqrtf，结果没用，__aeabi_ddiv照样被调用——后来才明白局部优化属性管不了系统自带的数学函数。

可行解决方案 #

1. 直接调用__ieee754_sqrtf #

直接调用底层的__ieee754_sqrtf函数，编译器会直接生成VSQRT.F32硬件指令，全程单精度计算，完全不会触发双精度函数依赖。

2. 编译时添加-fno-math-errno选项 #

给编译命令加上-fno-math-errno，这个选项会让sqrtf跳过错误码设置逻辑，转而使用纯单精度硬件指令实现，避免调用双精度函数。

最终选择 #

因为局部优化对系统数学函数无效，最终决定直接用__ieee754_sqrtf替换sqrtf。

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