Параллельная программа на OpenMP

[Dmitriy40]

Mihaylo
С вашей логикой должно быть очевидно 2n потоков, может быть даже плюс ещё несколько.
Так что да, два потока совершенно не очевидно, photon абсолютно прав (да ещё о многом важном и не упомянул).

[Mihaylo]

Два цикла с сугубо последовательными вычислениями. Последовательные вычисления плохо параллелятся. Требуются какие-то нетривиальные преобразования, чтобы эти циклы эффективно разбить по потокам.

[photon]

Mihaylo в сообщении #1434707 писал(а):

Последовательные вычисления плохо параллелятся.

Это верно, в частности, параллелить строки типа

P[i] = c / (-a * P[i - 1] + b);

как это сделано у автора нельзя - результат работы программы станет непредсказуем, но это не значит, что нельзя разбить на отдельные потоки внутренние циклы типа

Используется синтаксис C++

 for (i = n - 2; i >= 1; i--)
{
    C1[i][j] = P[i] * C1[i + 1][j] + Q[i];
}

[Mihaylo]

Да, преобразование только надо сделать. OpenMP сможет сам справиться?

[photon]

А, я соврал, тут же снова есть зависимость C1[i] от С1[i+1]...

[Dmitriy40]

А строки типа P[i] = c / (-a * P[i - 1] + b); параллелятся созданием копии массивов для правой части выражения, можно и лучше даже общей для всех потоков (многопоточный readonly доступ не конфликтует). Копию кстати тоже можно создавать параллельно.

(Оффтоп)

При дальнейшей оптимизации по памяти и быстродействию дублировать можно только пересекающиеся области массивов в точках разрывов между потоками, но это уже к учебной задаче по OpenMP не относится.