如何在-O3优化下保留Split函数算术逻辑并实现内联

问题背景 #

需要实现某论文第5页的Split函数，代码如下：

void Split(double a, double *x, double *y) {
   double z = a * ((1 << r) + 1);
   *x = z - (z - a);
   *y = a - *x;
}

启用-O3优化后，GCC内联该函数时会将z - (z - a)直接简化为a，导致原本的浮点数拆分逻辑失效。给z添加volatile修饰虽能阻止优化，但会破坏向量化效果，需寻找替代内联汇编的编译器提示方案。

可行方案 #

1. 针对性禁用死代码消除优化 #

给Split函数添加GCC扩展属性，禁用会消除冗余计算的tree-dce优化，同时保留-O3的其他优化项：

void __attribute__((optimize("no-tree-dce"))) Split(double a, double *x, double *y) {
   double z = a * ((1 << r) + 1);
   *x = z - (z - a);
   *y = a - *x;
}

注意：optimize属性是GCC专属扩展，不同版本的支持细节可能有差异，需实际测试验证是否仅阻止目标优化且不影响向量化。

2. 用__builtin_assume制造伪依赖 #

利用GCC内置函数向编译器传递“z的地址可能被外部访问”的提示，让编译器保留z的计算及后续算术步骤：

void Split(double a, double *x, double *y) {
   double z = a * ((1 << r) + 1);
   __builtin_assume(&z != NULL); // 提示编译器z的地址可能被使用，避免优化掉相关计算
   *x = z - (z - a);
   *y = a - *x;
}

该方案不会生成额外指令，也不会标记变量为volatile，因此对向量化没有负面影响。

3. 精细使用volatile限定临时变量 #

仅在关键计算步骤使用volatile修饰临时变量，减少对整体向量化的影响：

void Split(double a, double *x, double *y) {
   double z = a * ((1 << r) + 1);
   volatile double temp = z - a;
   *x = z - temp;
   *y = a - *x;
}

这种方式只保护z - a的计算不被优化，z的计算仍可参与向量化，比直接给z加volatile的影响更小。

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