为什么VisualC++自动矢量器不将这个简单的循环矢量化？

Question

我不明白为什么VisualC++不能自动矢量化这个循环.有什么想法吗？

我得到：

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来修复这个问题。

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

如果有人知道，请张贴！