GCC启用优化时如何开启FTZ/DAZ？-ffast-math相关疑问

问题背景 #

GCC的-funsafe-math-optimizations是-ffast-math的子选项，原本会开启**FTZ（Flush To Zero，刷新为零）和DAZ（Denormals Are Zero，非正规数视为零）**功能，但启用优化选项（如-O1、-O3、-Ofast）后，这两个功能会被自动禁用。即使在代码中显式调用_MM_SET_DENORMALS_ZERO_MODE(_MM_DENORMALS_ZERO_ON)，优化后该指令也会被忽略，仅当不配合-ffast-math编译时才能生效。

测试代码 #

#include <stdio.h>
#include <pmmintrin.h>

int main(int argc, char** argv) 
{
    //_MM_SET_DENORMALS_ZERO_MODE(_MM_DENORMALS_ZERO_ON);
    float normal_f = 1.18e-38f;
    double normal_d = 2.23e-308;

    normal_f *= 0.1f;
    if (normal_f != 0.0f) 
        printf("FTZ/DAZ disabled for floats (%e)\n", (double) normal_f);
    normal_d *= 0.1;
    if (normal_d != 0.0)
        printf("FTZ/DAZ disabled for doubles (%e)\n", normal_d);

    return 0;
}

编译现象 #

• 仅用gcc foo.c -ffast-math编译：程序无输出，FTZ/DAZ功能正常启用
• 添加优化选项（如gcc foo.c -ffast-math -O3）：程序输出非零值，FTZ/DAZ被禁用

  $ ./a.out 
  FTZ/DAZ disabled for floats (1.180000e-39)
  FTZ/DAZ disabled for doubles (2.230000e-309)
  

解决方案 #

方法1：手动操作MXCSR寄存器 #

高优化级别下GCC会重置MXCSR寄存器到默认状态，可通过直接读写寄存器强制开启FTZ/DAZ，该操作不会被优化器忽略：

#include <stdio.h>
#include <pmmintrin.h>

int main(int argc, char** argv) 
{
    // 读取当前MXCSR寄存器值
    unsigned int mxcsr = __builtin_ia32_stmxcsr();
    // 置位DAZ（0x8000）和FTZ（0x0040）标志位
    mxcsr |= 0x8040;
    // 写入修改后的MXCSR值
    __builtin_ia32_ldmxcsr(mxcsr);

    float normal_f = 1.18e-38f;
    double normal_d = 2.23e-308;

    normal_f *= 0.1f;
    if (normal_f != 0.0f) 
        printf("FTZ/DAZ disabled for floats (%e)\n", (double) normal_f);
    normal_d *= 0.1;
    if (normal_d != 0.0)
        printf("FTZ/DAZ disabled for doubles (%e)\n", normal_d);

    return 0;
}

编译命令：gcc foo.c -ffast-math -O3，此时程序无输出，FTZ/DAZ正常启用。

方法2：使用编译选项强制设置 #

使用GCC的-mdenormal-fp-math选项直接控制DAZ和FTZ行为，该选项优先级高于优化器的默认设置：

gcc foo.c -ffast-math -O3 -mdenormal-fp-math=positive-zero,positive-zero

参数说明：

• 第一个positive-zero：开启DAZ（输入非正规数视为正零）
• 第二个positive-zero：开启FTZ（输出非正规数刷新为正零）

其他可能被优化禁用的-ffast-math特性 #

在高优化级别下，除FTZ/DAZ外，以下-ffast-math的子特性也可能被影响：

• 符号零处理：-fno-signed-zeros（忽略符号零的差异）可能在某些优化路径下被部分还原，导致运算时保留符号零
• 浮点近似运算：-freciprocal-math（使用近似倒数替代除法）在循环或高频运算中，可能被优化器替换为更保守的除法实现
• 精度放宽：-ffinite-math-only（假设无NaN/Inf）可能在优化时被部分限制，尤其是在涉及外部函数调用的场景中

注：以上现象在GCC 6.3和10.2版本中均已观测到

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