Возникла следующая проблема: имеется вот такой участок кода
start = clock();
LLTWithBlock(n, 4, a);
end = clock();
duration = (double)(end - start);
std::cout << "Alltime: " << duration << std::endl;
Время выполнения 8
секунд.
После чего заменяю константу "4" на переменную "r" (тоже равную 4, считывается из файла)
LLTWithBlock(n, r, a);
И новое время выполнения становится 22
секунды. Я отдаленно понимаю, что в ассемблерном коде будут создаваться разные экземпляры функции, но откуда такая разница во времени выполнения?
P.S. Время замерялось не один раз, данная ситуация повторяется при каждом запуске программы.
upd: меняется конкретно время выполнения следующей функции (вызываемой внутри)
void A_LLT(unsigned n, unsigned r, unsigned numBlock, double *b, double *bb, double *a){
unsigned adr;
unsigned cnt = r*(numBlock + 1);//номер 1 строки в текущем блоке(прямоугольном) в большой матрице
unsigned nb = (n - cnt) / r, end;//n/r число целых блоков размера р в матрице.
//Отнимаем 1(мы его раскладываем) и нумблок(число уже разложенных блоков)
for (unsigned ib = 0; ib < nb; ib++)//ход по квадратным блокам
{
end = ib * r + r;
for (unsigned i = ib * r; i < end; i++)//по строкам в блоке
{
adr = (cnt + i) * (cnt + i + 1) / 2 + cnt;
for (unsigned j = 0; j <= i; j++, adr++)//по строкам
{
for (unsigned k = 0; k < r; k += 4)
{
a[adr] -= b[i*r + k] * b[j*r + k];
a[adr] -= b[i*r + k + 1] * b[j*r + k + 1];
a[adr] -= b[i*r + k + 2] * b[j*r + k + 2];
a[adr] -= b[i*r + k + 3] * b[j*r + k + 3];
}
}
}
}
for (unsigned i = nb * r, end = i + n % r; i < end; i++)//по строкам в блоке n - r - numBlock * r
{
adr = (cnt + i) * (cnt + i + 1) / 2 + cnt;
for (unsigned j = 0; j <= i; j++, adr++)//по строкам
{
for (unsigned k = 0; k < r; k += 4)
{
a[adr] -= b[i*r + k] * b[j*r + k];
a[adr] -= b[i*r + k + 1] * b[j*r + k + 1];
a[adr] -= b[i*r + k + 2] * b[j*r + k + 2];
a[adr] -= b[i*r + k + 3] * b[j*r + k + 3];
}
}
}}
r
после строчкиstart = ...
?r
точно равно 4? отлладочный вывод делали?