如何为__m128向量实现fabs()函数？

不需要通过与1.0f/-1.0f相乘来实现，位掩码操作是更高效且原生支持的方案，具体实现如下：

• 核心原理：单精度浮点数（float）的符号位是最高位（第31位），将该位清零即可得到对应值的绝对值。对于__m128向量，只需构造一个掩码向量（每个float元素的符号位为0，其余位为1），通过位与操作清零原向量的符号位。

• 代码实现：

#include <xmmintrin.h>

__m128 sse_fabs_ps(__m128 vec) {
    // 构造掩码：0x7FFFFFFF 是单精度浮点数的无符号掩码，复制到__m128的4个元素中
    const __m128 mask = _mm_set1_ps(*(const float*)&(const unsigned int){0x7FFFFFFF});
    // 位与操作清零所有元素的符号位
    return _mm_and_ps(vec, mask);
}

• 指令集说明：_mm_and_ps对应SSE1指令集中的andps指令，属于基础SSE指令，完全适配__m128类型，不需要依赖更高版本的AVX指令集。

• 为什么不推荐相乘方案：相乘方案需要额外的符号判断、条件选择操作，会增加指令延迟和周期数；而位与操作是单周期指令，执行效率远高于相乘方案。

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