векторизация умножения матриц

[alek3s]

имеется две матрицы 1000 X 1000 которые нужно стандартно перемножить. элементы этих матриц - матрицы типа double 5X5.

вот какой получился алгоритм перемножения

for(int k = 0;k < N;k++)
for(int g2 = 0;g2 < (N*N);g2+= N)
for(int i = 0;i < N;i++)
//теперь работаем с внутренней матрицей
for(int j=0;j<M;j++)
for(int h1 = 0;h1 < M;h1++)
for(int g3 = 0;g3 < M;g3++)
{
rez_matrix[g2+k].MM[g3][j] += first_matrix[i+g2].MM[g3][h1] * second_matrix[i*N+k].MM[h1][j];
}

вот туточки у меня вопросик.
данный алгоритм векторизуется на ура в Intel C++ без всяких прагм и подсказок. но почемуто время выполнения векторизированного цикла в 2 раза больше чем того-же цикла только без векторизации.
может у кого мысли есть?

[alek3s]

неужели никто не знает про векторизацию?