我有一个矩阵乘法代码，看起来像这样：

for(i = 0; i < dimension; i++)
    for(j = 0; j < dimension; j++)
        for(k = 0; k < dimension; k++)
            C[dimension*i+j] += A[dimension*i+k] * B[dimension*k+j];

这里，矩阵的大小由维度表示。
现在，如果矩阵的大小是2000，运行这段代码需要147秒，而如果矩阵的大小是2048，则需要447秒。所以虽然不同。的乘法是(2048 * 2048 * 2048)/(2000 * 2000 * 2000)= 1.073，时间差是447/147 = 3。有人可以解释为什么会发生这种情况吗？我预计它会线性扩展，这不会发生。我不是想做最快的矩阵乘法代码，只是试图理解为什么会发生。

规格：AMD Opteron双核节点(2.2GHz)，2G RAM，gcc v 4.5.0

程序编译为gcc -O3 simple.c

我也在英特尔的icc编译器上运行，并看到类似的结果。

编辑：

正如在意见/答案建议，我运行代码与维度= 2060，它需要145秒。

完整的程序：

#include <stdlib.h>
#include <stdio.h>
#include <sys/time.h>

/* change dimension size as needed */
const int dimension = 2048;
struct timeval tv; 

double timestamp()
{
        double t;
        gettimeofday(&tv, NULL);
        t = tv.tv_sec + (tv.tv_usec/1000000.0);
        return t;
}

int main(int argc, char *argv[])
{
        int i, j, k;
        double *A, *B, *C, start, end;

        A = (double*)malloc(dimension*dimension*sizeof(double));
        B = (double*)malloc(dimension*dimension*sizeof(double));
        C = (double*)malloc(dimension*dimension*sizeof(double));

        srand(292);

        for(i = 0; i < dimension; i++)
                for(j = 0; j < dimension; j++)
                {   
                        A[dimension*i+j] = (rand()/(RAND_MAX + 1.0));
                        B[dimension*i+j] = (rand()/(RAND_MAX + 1.0));
                        C[dimension*i+j] = 0.0;
                }   

        start = timestamp();
        for(i = 0; i < dimension; i++)
                for(j = 0; j < dimension; j++)
                        for(k = 0; k < dimension; k++)
                                C[dimension*i+j] += A[dimension*i+k] *
                                        B[dimension*k+j];

        end = timestamp();
        printf("\nsecs:%f\n", end-start);

        free(A);
        free(B);
        free(C);

        return 0;
}