《深入理解计算机系统第三版》— 第六章 — cache lab

这章的实验还挺有意思，简单记录下做法。

一、partA

partA 相对来说比较简单，需要注意的就是 LRU 的使用，以及对于同一个 set 中的不同行，它们并不是存放连续的地址。以 s = 1，E = 2，b = 1 为例：第一个 set 的第一行存放的是 0x2 0xa ···，第二行存放的是 0x6 0xd ···。

二、partB 32*32

像手册中提到的，这里的关键在于对角线的冲突。

首先，我们不难发现，由于缓存中共有 32*32 = 1024 字节，因此在矩阵中，每 1024 / 4 / 32 = 8 行就会引发冲突。因此，我们最好选择 blocksize = 8*8。然而，这样会导致对角线上的 a 和 b 的 block 同时操作，引发地址冲突，即 A[0][0]-A[7][7] 和 B[0][0]-B[7][7] ··· 是冲突的。经过观察，可以看到，当我们加载 A[0][0] 后，由于需要写入 B[0][0]，引发了冲突（其后的冲突都类似）。这导致后续对A[0][1]的读取产生 miss。因此，我们希望能够使得通过改变 A[0][0]-A[0][7] 只出现 A[0][0] 一次 miss 来改善程序。最终我们选择使用一个临时变量来保存对角线元素。行读取结束后，再将该元素写回 B 中。

int i, j, ii, jj, tmp, diagonal;
    int msize = 8;
    int nsize = 8;

    for (ii = 0; ii < N; ii += nsize)
    {
   
        for (jj = 0; jj < M; jj += msize)
        {
   
            if (ii != jj)
            {
   
                for (i = ii; i < ii + nsize; i++)
                {
   
                    for (j = jj; j < jj + msize; j++)
                    {
   

                        tmp = A[i][j];
                        B[j][i] = tmp;
                    }
                }

                continue;
            }

            for (i = ii; i < ii + nsize; i++)
            {
   
                diagonal = A[i][i];
                for (j = jj; j < jj + msize; j++)
                {
   
                    if (i == j)
                    {
   
                        continue;
                    }

                    tmp = A[i][j];
                    B[j][i] = tmp;
                }
                B[i][i] = diagonal;
            }
        }
    }

通过 ii != jj 的判断来分辨是否为对角线 block。

三、partB 64*64

这部分是这个实验最有趣的地方。类似地，我们不难分析出来，在每个矩阵中每4行就会引发一次冲突。因此这里我们首先想到的就是使用 4*4 矩阵。然而，由于每个 set 中的 block 能保存 32 字节，因此我们需要分析如何能够最大化的利用剩余在缓存中的字节。而如果想要充分利用所有字节，就肯定需要考虑 8*8 矩阵。以下分两种情况讨论：

1、位于对角线上的 8*8 矩阵

对于此类矩阵，A 矩阵和 B 矩阵加载到相同的缓冲区中，会相互冲突。其中 A 矩阵的1-4行，A 矩阵的5-8行，B 矩阵的1-4行，B 矩阵的5-8行都会映射到相同的 block 中。假设这些 block 的序号分别为1，9，17，25。我们可以把每个 8*8 矩阵划分为4个 4*4 矩阵：

首先考虑区域1。如果我们按照如下图所示的方式，对矩阵 A 中的1号子矩阵从左至右，从上至下使用前面的对角线优化方式进行访问：

那么访问 A 矩阵中的第一行时，其命中缺失（忽略冲突）情况应该为：

A[0][0] m
A[0][1] h
B[1][0] m
A[0][2] h
B