Я новичок в OpenMP, и я пытаюсь парализовать следующий код с помощью OpenMP:
#pragma omp parallel for
for(int k=0;k<m;k++)
{
for(int j=n-1;j>=0;j--)
{
outX[k+j*m] = inB2[j+n * k] / inA2[j*n + j];
for(int i=0;i<j;i++)
{
inB2[k*n+i] -= inA2[i+n * j] * outX[k + m*j];
}
}
}
Параллелизировать внешний цикл довольно просто, но чтобы оптимизировать его, я хотел также парализовать самый внутренний цикл (тот, который повторяется по i). Но когда я пытаюсь сделать это вот так:
#pragma omp parallel for
for(int i=0;i<j;i++)
{
inB2[k*n+i] -= inA2[i+n * j] * outX[k + m*j];
}
компилятор не векторизует внутренний цикл («версия цикла для векторизации из-за возможного сглаживания»), что замедляет его работу. Я скомпилировал его с помощью gcc -ffast-math -std=c++11 -fopenmp -O3 -msse2 -funroll-loops -g -fopt-info-vec prog.cpp
Спасибо за любой совет!
РЕДАКТИРОВАТЬ: я использую ключевое слово __restrict для массивов.
EDIT2: Интересно то, что когда я сохраняю только прагму во внутреннем цикле и удаляю ее из внешнего, gcc векторизует ее. Так что проблема возникает только тогда, когда я пытаюсь парализовать оба цикла.
EDIT3: компилятор векторизует цикл, когда я использую #pragma omp parallel для simd. Но это все равно медленнее, чем без распараллеливания внутреннего цикла.
