为什么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代替修复问题。

我仍然不知道为什么会这样。

如果有人知道，请贴出来！