循环嵌套之美

for循环的嵌套顺序不一样，对于代码的运行效率完全不一样。如下面的代码：

void matrixMultiply(int **a, int **b, int **c, int n)
{
    for (int i =0; i < n; i++)
        for (int j = 0; j < n; j++)
            for (int k = 0; k < n; k++)
                c[i][j] += a[i][k] * b[k][j];
}

如果我们将循环嵌套从ijk，换成ikj，则运算效率大大提高。原因在于现代计算机结果都会有一级，二级等缓存机理。如下图所示：

                                                        图1. 计算机缓存图

由于计算器数据从L1到CPU的速度是 L2到CPU 速度的好多倍，而主存到CPU的速度更慢。因此计算机采用了一种策略，把所需要计算的数据预载到缓存中，当出现一个缓存未命中的时候，把需要的数据和相邻字节中的数据装入到缓存中，以此来加速运行速度。而题目中的同行元素也就是一维数组的元素会被预载到缓存中，当数据量大的时候，非同行的数据并不一定会被预载到缓存中，因此ikj运算时，会将a数组和b数组中的同行元素预载到缓存中。这样每次在循环j的数据的时候都是已经缓存过的，运算速度提升很多倍。以下是运行速度比较：