c++ - c=c+a*b 的 OpenMP 4 simd 矢量化

Question

我不知道 OpenMP 4 是否支持这个 for 循环。使用和不使用 pragma 的速度是一样的。

#pragma omp  for simd
for (size_t i = 0; i < col; i++)
{
    C[i] += A[i]* B[i];
}

Answer 1

编译指示无效的原因(我猜)是双重的:

• 代码在没有simd 指令的情况下已经向量化；和
• 代码无论如何都是内存限制的，因此添加更多线程来计算它不会产生太大差异，除非它可以让您访问更多内存带宽。参见 this excellent answer了解更多详情。