使用 GCC vector 内在函数计算 x^y

Question

假设我有一个 2 元素 vector 定义如下(使用打包 vector 的 GCC 语法)

// packed vector of 2-elements
typedef double v2d __attribute__((vector_size(sizeof(double)*2)));

v2d x = ...;
double y = ...;

x[0] = pow(x[0], y)
x[1] = pow(x[1], y)

我想知道是否有更快的方法来使用 vector 运算来进行两次幂计算。该体系结构是 x86-64 上的 GCC，平台特定代码是可以的。

Answer 1

是的，如果您没有特殊情况(负数、0、1、NaN 等...)，这应该是可能的，这样代码路径是线性的。

Here是 IEEE754 double 的 pow 函数的通用代码，它没有循环结构，因此如果您充实所有特殊情况，向量化似乎很简单。玩得开心。