从矩阵乘法的不同计算方式来看局部性原理

今天碰到的关于矩阵乘法不同情况下运算速度的问题，隐约记得是因为缓存的问题，后来突然想起来CSAPP那本书上讲过这个东西的，就是通过矩阵乘法三重循环的不同顺序来讲的局部性原理的，所以翻过来又看了一下。

两个矩阵A,B相乘得到C【为了方便起见，把它们都看成n*n的方阵】经典的做法就是用三重循环来实现，但是具体这三重循环如何进行排列，就非常的有讲究。
假设n是一个非常大的数,也就意味着如果跨行的话,必然不会缓存命中,假设只有一个告诉缓存,其块大小为32字节,也就说一个块只能存4个double

经典做法ijk和jik

for(i=0;i<n;++i)
    for(j=0;j<n;++j)
    {
        double sum=0;
        for(k=