Как правильно распараллелить вложенные циклы for

Я работа с OpenMP для распараллеливания скалярного вложенного цикла for:

double P[N][N];
double x=0.0,y=0.0;

for (int i=0; i<N; i++)
{
    for (int j=0; j<N; j++)
    {
        P[i][j]=someLongFunction(x,y);
        y+=1;
    }
    x+=1;
}

В этом цикле важно то, что матрица P должна быть одинаковой как в скалярной, так и в параллельной версиях:

Все мои возможные испытания не увенчались успехом ...