缓存阻塞比原始代码最差吗？

Question

我正在研究如何改善矩阵操作（使用双重类型），并且正在尝试一些技术，例如缓存和循环展开。第二个确实取得了成功，但是我无法使用阻止来提高自己的性能。我不是因为我做错了什么，或者是由于阻塞而引起的，在这种情况下根本没有用。

没有操作的原始代码是：

    for (int i=0; i<N; i++){
        for (int j=0; j<N; j++){
            d[i][j] = 0.0;
        }
    }

    for (int i=0; i<N; i++){
        for (int j=0; j<N; j++){
            for (int k=0; k<K_MAX; k++){
                d[i][j] += 2 * a[i][k] * (b[k][j] - c[k]);
            }   
        }
    }

k_max始终为8，n从250 500 750 1000 1500 2000 2000 2550 3000中获取值

，而我试图处理的封锁是：

    for (int i= 0 ; i<N; i+=block_size){
        for (int j=0; j<N; j+=block_size){      
            for (int ii=i; ii<min (i+block_size, N); ii++){
                for (int jj=j; jj<min(j+block_size, N); jj++){      
                    d[ii][jj] = 0.0;
                    for (int k = 0; k<K_MAX; k++){
                        d[ii][jj] += 2 * a[ii][k] * ( b[k][jj]- c[k]);
                    }
                }
            }
        }
    }

我可能为block_size选择一个不好的值，因为我做了不了解如何选择一个不错的一个，但是我尝试了N的所有分隔线来选择一个块大小，从1到N。而且，我尝试使用适合缓存线上的元素数量的倍数（8倍）就像8、64、128、256和512（我知道n并不总是该值的倍数，有必要处理块无法达到的元素，我尝试并做得很好，因为我已经获得了正确的输出），但表现没有改善。我还尝试使用所有n个块大小值使用相同的块大小值，但是您可以猜到，什么都没有实现。

我的处理器是Intel Core i7-10870h。

先感谢您