防止C++矢量化的数学函数
防止C++矢量化的数学函数
提问于 2020-09-17 07:38:07
我有如下所示的代码，并希望将其矢量化。我用一个情报顾问分析了整个代码。它说我不能把它矢量化，因为它里面有数学函数。当然，它指出了循环中使用的sin和cos函数。
如何在不使用英特尔短向量数学库的情况下矢量化这个循环？
代码：
for (size_t j = 0; j < NA; ++j) {
esf = sfs[j];
x = p_data[3 * j];
y = p_data[3 * j + 1];
z = p_data[3 * j + 2];
p = x * qx + y * qy + z * qz;
Ar += esf * cos(p);
Ai += esf * sin(p);
}回答 1
Stack Overflow用户
回答已采纳
发布于 2020-09-17 07:58:20
，它说我不能把它矢量化，因为里面有数学函数。
实际上，阻止矢量化的是Ar +=和Ai +=术语，因为这意味着j = 2的输出将取决于j = 1的输出。如果Ar和Ai只是输出，则可以将它们改为数组，并在运行循环后对它们进行求和。
// init to 0
double Ar_elem[NA] = {0.0};
double Ai_elem[NA] = {0.0};
for (size_t j = 0; j < NA; ++j) {
esf = sfs[j];
x = p_data[3 * j];
y = p_data[3 * j + 1];
z = p_data[3 * j + 2];
p = x * qx + y * qy + z * qz;
Ar_elem[j] = esf * cos(p);
Ai_elem[j] = esf * sin(p);
}
//sum
double Ar = std::accumulate(begin(Ar_elem), end(Ar_elem), 0, plus<double>());
double Ai = std::accumulate(begin(Ai_elem), end(Ai_elem), 0, plus<double>());页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/63933350
复制相关文章
相似问题
腾讯云开发者
