OpenMP parallelization

我正在写一个科学用途的c++程序。程序运行良好，返回良好的结果，所以我决定使用OpenMP提高其性能。我要优化的循环如下:

        //== #pragma omp parallel for private(i,j)
        for (k=0; k < number; k++)
        {
        for (i=0; i < L; i++)
        {       
                for (j=0; j < L; j++)
                {       
                        red[i][j] = UNDEFINED;
                }
        }

        Point inicial =  {L/2, L/2, OCCUPIED};
        red[L/2][L/2] = OCCUPIED;
        addToList(inicial, red, list, L,f); 
        oc.push_back(inicial); 
        while (list.size() > 0 && L > 0)
        {       
                punto = selectPoint(red, list, generator, prob, p); 
                if (punto.state == OCCUPIED)
                {       
                        addToList(punto, red, list, L,f);
                        oc.push_back(punto);
                }    
                else
                {
                        out.push_back(punto);
                }

        }
        L = auxL; 
        oc.clear();
        out.clear();
        list.clear();
        }
 f = f*1.0/(number*1.0);
        if (f > 0.5)
        {
                inta = inta;
                intb = p;
                p = (inta + intb) / 2.0;
        }
        else if (f < 0.5)
        {
                intb = intb;
                inta = p;
                p = (inta + intb) / 2.0;
        }
        cout << p << endl;

        }

我对OpenMP的尝试如上所述。正如您所看到的，我已经将i和j声明为私有，因为它们是在并行部分之前声明的。我也试图使L私人，没有结果。到处都是分段错误和坏指针。我认为问题在于while循环嵌套在里面。我的问题是:在这种情况下omp parallel for是否正确?还是我应该只优化那个while循环?std::vector是否干扰OpenMP?

注意:list, oc和out是std::vector<Point>, Point是一个简单的结构体，有三个int属性。addToList是一个内部没有循环的函数