КЭШ-промахи при умножении квадратных матриц (Си)

[Stasiksis]

Задача состоит в том, чтобы составить аналитически функцию КЭШ-промахов, которая будет зависеть от размера L2 КЭША (у меня он 2 Mb) и от размера матриц. Есть программа для умножения матриц, написанная на СИ:

#include <stdio.h>
#include <stdlib.h>
#include <stdint.h>
#include <locale.h>
#include <string.h>
#include <time.h>

typedef int32_t T;
	T A[M][N], B[N][H], C[M][H];
int mul () {
	int i, j, k;
	for (i = 0; i < M; i++)
		for (j = 0; j < H; j++) 
			for (k = 0; k < N; k++)
				C[i][j] += A[i][k] * B[k][j];
}
int main() {
	static struct timespec t1, t2;
	int i, j, k, n;
	T temp;
	memset( C, 0, sizeof(C) );
	for (i = 0; i < M; i++ )
		for (j = 0; j < N; j++) 
			A[i][j] = 1;
	for (i = 0; i < N; i++ )
		for (j = 0; j < H; j++) 
			B[i][j] = 1;

	setlocale(LC_ALL, "");
	clock_gettime(CLOCK_PROCESS_CPUTIME_ID, &t1 );
	
	mul();

	clock_gettime(CLOCK_PROCESS_CPUTIME_ID, &t2);
	
	printf( "%.05g\n",((t2.tv_sec+t2.tv_nsec*1.e-9)-(t1.tv_sec+t1.tv_nsec*1.e-9)));
	return 0;
}

Очевидно, что если размер матрицы меньше, чем

3 * N^2 * 8 sizeof(int) <= 2 Mb

, где N - порядок матрицы, то получим примерно, что при порядке матриц, меньше 418 КЭШ-промахов не будет, ну это только начало, кто-нибудь может помочь и подсказать, как для данного случая перемножения матриц определить количество КЭШ-промахов при большем порядке матриц?

[Stasiksis]

Неужели никто не в силах помочь?