Почему параллельная версия работает за то же время или больше чем последовательная версия алгоритма? Что я сделал не так?
LU разложение с помощью OpenMP:
void lup_od_omp(double* a, int n){
int i,j,k;
for(k = 0; k < n - 1; ++k)
{
#pragma omp parallel for shared(a,n,k) private(i,j)
for(i = k + 1; i < n; i++)
{
a[i*n + k] /= a[k*n + k];
for(j = k + 1; j < n; j++)
{
a[i*n + j] -= a[i*n + k]*a[k*n + j];
}
}
}
}
P.S. Матрица хранится в одномерном массиве. Та же самая ситуация с временем выполнения (параллельно >= последовательно) и с матрицей, хранящейся в двумерном массиве.