OpenMP. Как распараллелить цикл?

Благодарю за помощь.

OpenMP. Как распараллелить цикл?

https://pastebin.com/XE4UbeB8
Матрица у меня хранится в файле, в формате mtx (Matrix Market), скачала тут https://sparse.tamu.edu/VDOL/freeFlyingRobot_7

OpenMP. Как распараллелить цикл?

Сделала вот так, но теперь вернулась проблема с неверными значениями...

#pragma omp parallel private(i, j)
{
    #pragma omp for 
    for (i = 0; i < N; i++)
    {
        for (j = RowIndex[i]; j < RowIndex[i + 1]; j++)
            result[i] += val[j] * vector[col[j]];    
    }
}

#pragma omp parallel private(i, j)
{
    #pragma omp for 
    for (i = 0; i < N; i++)
    {
        if (i != col[j])
        {
            result[col[j]] += val[j] * vector[i];
        }    
    }
}

OpenMP. Как распараллелить цикл?

Roman, Спасибо, так работает. А возможно как-то избавиться от atomic? На сколько я поняла, его использование сильно замедляет программу.

OpenMP. Как распараллелить цикл?

Roman,

#pragma omp parallel for private (i, j) 
        for (i = 0; i < N; i++)
        {
            for ( j=RowIndex[i]; j<RowIndex[i + 1]; j++)
            {
                if (i == col[j])
                    result[i] += val[j] * vector[col[j]];
                else
                {
                    result[i] += val[j] * vector[col[j]];
                    #pragma omp atomic
                    result[col[j]] += val[j] * vector[i];
                }
}

OpenMP. Как распараллелить цикл?

Спасибо, но это я уже читала. Смогла распараллелить, но значения правильные получаю через раз. Поэтому помощи и прошу.