Untitled

Race conditions with OpenMP
testfunc(const double* var, const double* par)
    {
       .......
    }

    #define Nx 10000
    #define Ny 10000
    static double tmp[Ny][Nx];
    double parGlob[2]; //<- Here are they!!!
    #pragma omp threadprivate(parGlob)  // <-Magic directive!!!!

    int main()
    {
 // Not here !!!! ->       double par[2]; // parameters
        double xmin[]={0,0} // limits of 2D integration
        double xmax[]={1,1};// limits of 2D integration

        double val,Tol=1e-7,AbsTol=1e-7;
        int j,NDim=2,NEval=1e5;

    #pragma omp parallel for private(j,val)  // no `i` inside `private` clause
        for (int i=0;i<Nx;i++)
        {
            for (j=0;j<Ny;j++)
            {
                parGlob[0]=i;
                parGlob[1]=j*j;
                adapt_integrate(testfunc,par, NDim, xmin, xmax,
                                NEval, Tol, AbsTol, &val, &err);
                tmp[i][j] = val;
            }
        }
    }