为什么Visual C ++自动矢量化器不能矢量化这个简单的循环？

Question

我无法弄清楚为什么Visual C ++无法自动向量化这个循环...任何想法？

我明白了：

testvec.cpp:12: info C5002: loop not vectorized due to reason '1200'

原因代码1200是：

  

循环包含阻止矢量化的循环携带数据依赖性。循环的不同迭代相互干扰，使得对循环进行矢量化会产生错误的答案，并且自动矢量化器无法证明自己没有这样的数据依赖性。

但为什么？

#include <stdlib.h>

int main(int argc, char *argv[])
{
    int const n = argc;
    double
        *const p1 = (double *)malloc(n * n * sizeof(*p1)),
        *const p2 = (double *)malloc(n * n * sizeof(*p2));
    for (int j = 0; j < n; ++j)
    {
        double const sj = p1[n * j];
        for (int i = 0; i < n; ++i)
        {
            double const sum = p1[i] + sj, old = p1[i + n * j];
            p2[i + n * j] = sum < old ? sum : old;
        }
    }
}

Answer 1

我终于找到了如何修复它...似乎n * j中的乘法是罪魁祸首。
将其作为int nj = n * j;提取到外部并在内部循环中使用nj来修复问题。

我仍然不知道为什么会发生这种情况。

如果有人知道，请发布！