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的数据的时候都是已经缓存过的,运算速度提升很多倍。以下是运行速度比较: