Function selfBoolean

source/MRMesh/MRMeshBoolean.cpp:254–386

}

Expected<MR::Mesh> selfBoolean( const Mesh& inMesh )
{
    auto box = Box3d( inMesh.computeBoundingBox() );
    CoordinateConverters converters{ .toInt = getToIntConverter( box ),.toFloat = getToFloatConverter( box ) };
    std::vector<EdgeTri> intersections;
    ContinuousContours contours;

    Mesh mesh = inMesh;

    int iters = 0;
    const int cMaxFixLoneIterations = 3; // lets not do it many times
    for ( ;; ++iters )
    {
        // find intersections
        intersections = findSelfCollidingEdgeTrisPrecise( mesh, converters.toInt );
        // order intersections
        contours = orderSelfIntersectionContours( mesh.topology, intersections );

        // find lone
        auto loneProcessing = detail::subdivideSelfLone( mesh, converters, contours );
        if ( loneProcessing == detail::LoneProccessingState::None )
            break;
        else if ( loneProcessing == detail::LoneProccessingState::DegenerateOnly || iters == cMaxFixLoneIterations )
        {
            // in some rare cases there are lone contours with zero area that cannot be resolved
            // they lead to infinite loop, so just try to remove them
            removeLoneContours( contours, true );
            break;
        }
    }

    auto holePairs = detail::findSelfContoursMapping( contours );
    auto meshCpy = mesh; // for now lets do copy each time, TODO: do copy only if needed
    auto intContours = getOneMeshSelfIntersectionContours( mesh, contours, converters );

    // update non-closed
    ParallelFor( holePairs, [&] ( size_t hpId )
    {
        auto [f, s, r] = holePairs[hpId];
        auto& cf = contours[f];
        if ( isClosed( cf ) )
            return;

        if ( cf.capacity() < 2 * cf.size() + 3 )
            cf.reserve( 2 * cf.size() + 3 );
        auto& icf = intContours[f].intersections;
        if ( icf.capacity() < 2 * icf.size() + 3 )
            icf.reserve( 2 * icf.size() + 3 );

        auto& cs = contours[s];
        auto& ics = intContours[s].intersections;

        auto prev = mesh.topology.prev( cf.front().edge );
        auto next = mesh.topology.next( cf.back().edge );
        auto svId = mesh.topology.dest( prev );
        auto fvId = mesh.topology.dest( next );

        auto triFVerts = mesh.topology.getTriVerts( cf.front().tri() );