是否可以并行化或展开此循环？

Question

我正在尝试查看是否可以提高C ++中以下循环的性能，C ++使用二维向量（_external和_Table），并且对先前的迭代具有循环依赖关系。 此外，它在最内层的循环中具有计算出的索引访问器，这将使_Table的访问在右侧不连续。

int N = 8000;
int M = 400
int P = 100;
for(int i = 1; i <= N; i++){
    for(int j = 0; j < M; j++){
        for(int k =0; k < P; k++){
            int index = _external.at(j).at(k);
            _Table.at(j).at(i) += _Table.at(index).at(i-1);
        }
    }
}

我该怎么做才能改善这样的循环性能？

Answer 1

在我看来，这些语句的顺序如下：

    int index = _external.at(j).at(k);
    _Table.at(j).at(i) += _Table.at(index).at(i-1);

被执行对于正确性至关重要。 （也就是说，如果i，j，k的迭代顺序发生变化，那么结果将是不同的……并且不正确。）

所以我认为你只留下了微优化，如吊装表达_Table.at(j).at(i)和_external.at(j)从最里面的循环。

考虑一下：

    for(int k =0; k < P; k++){
        int index = _external.at(j).at(k);
        _Table.at(j).at(i) += _Table.at(index).at(i-1);
    }

此循环将数字重复添加到_Table.at(j).at(i) 。 由于（通过检查） _Table.at(index).at(i-1) ） _Table.at(index).at(i-1)必须从表的其他单元格读取（因为i-1与i ），因此可以这样做：

    int temp = 0;
    for(int k =0; k < P; k++){
        int index = _external.at(j).at(k);
        temp += _Table.at(index).at(i-1);
    }
    _Table.at(j).at(i) += temp;

这样可以减少对at的调用次数，也可以稍微提高缓存性能。