CGAL users discussion list

[Andreas Fabri <>]

What kernel do you use?

On 2/22/2020 9:51 PM, Yifang Zhao wrote:

I have a for loop containing approximately 1,000,000 iterations of calling AABBTree.all_intersections() method. Thus, I use #pragma omp parallel for to speed up the program. However, I got an error as "A heap has been corrupted". The cpp file is roughly like the following:

Polyhedron polyhedron;

Facet_tree m_facet_tree(faces(m_polyhedron).first, faces(m_polyhedron).second, m_polyhedon);

#pragma omp parallel for

for (int i = 0; i<100; ++i){

    Plane plane = Plane(1, 0, 0, 0);

    // Compute intersections

    typedef std::vector<Facet_tree::Object_and_primitive_id> Intersections;

    Intersections intersections;

    m_facet_tree.all_intersections(plane, std::back_inserter(intersections));

}

Finally, The program will break here:

    ~Handle()

    {

    if ( PTR && (--PTR->count == 0))

        delete PTR;

    }

What are the possible reasons that cause this problem? Thanks a lot!

Kind Regards,

Yifang Zhao.

 

Sent from Mail for Windows 10

 

-- 
Andreas Fabri, PhD
Chief Officer, GeometryFactory
Editor, The CGAL Project

phone: +33.492.954.912    skype: andreas.fabri