Function getBooleanPoints

source/MRMesh/MRMeshBoolean.cpp:672–808

}

Expected<BooleanResultPoints> getBooleanPoints( const Mesh& meshA, const Mesh& meshB,
    BooleanOperation operation, const AffineXf3f* rigidB2A )
{
    MR_TIMER;

    BooleanResultPoints result;
    result.meshAVerts.resize( meshA.topology.lastValidVert() + 1 );
    result.meshBVerts.resize( meshB.topology.lastValidVert() + 1 );

    const auto converters = getVectorConverters( meshA, meshB, rigidB2A );
    const auto intersections = findCollidingEdgeTrisPrecise( meshA, meshB, converters.toInt, rigidB2A );
    result.intersectionPoints.reserve( intersections.size() );

    FaceBitSet collFacesA, collFacesB;
    VertBitSet destVertsA, destVertsB, orgVertsA, orgVertsB;
    collFacesA.resize( meshA.topology.lastValidFace() + 1 );
    collFacesB.resize( meshB.topology.lastValidFace() + 1 );
    orgVertsA.resize( meshA.topology.lastValidVert() + 1 );
    orgVertsB.resize( meshB.topology.lastValidVert() + 1 );
    destVertsA.resize( meshA.topology.lastValidVert() + 1 );
    destVertsB.resize( meshB.topology.lastValidVert() + 1 );
    for ( const auto& et : intersections )
    {
        if ( et.isEdgeATriB() )
        {
            gatherEdgeInfo( meshA.topology, et.edge, collFacesA, orgVertsA, destVertsA );
            collFacesB.set( et.tri() );

            const auto isect = findEdgeTriIntersectionPoint( meshA, et.edge, meshB, et.tri(), converters, rigidB2A, EdgeTriComponent::Tri );
            result.intersectionPoints.emplace_back( isect );
        }
        else
        {
            gatherEdgeInfo( meshB.topology, et.edge, collFacesB, orgVertsB, destVertsB );
            collFacesA.set( et.tri() );

            const auto isect = findEdgeTriIntersectionPoint( meshB, et.edge, meshA, et.tri(), converters, rigidB2A, EdgeTriComponent::Edge );
            result.intersectionPoints.emplace_back( isect );
        }
    }

    if ( orgVertsA.intersects( destVertsA ) || orgVertsB.intersects( destVertsB ) )
    {
        // in this case we are not able to detect inside outside correctly
        BooleanResultMapper mapper;
        auto boolRes = MR::boolean( meshA, meshB, operation, rigidB2A, &mapper );
        if ( !boolRes.valid() )
            return unexpected( boolRes.errorString );

        if ( !mapper.maps[int( BooleanResultMapper::MapObject::A )].old2newVerts.empty() )
            for ( auto v : meshA.topology.getValidVerts() )
            {
                auto vn = mapper.maps[int( BooleanResultMapper::MapObject::A )].old2newVerts[v];
                if ( vn.valid() )
                    result.meshAVerts.set( v );
            }

        if ( !mapper.maps[int( BooleanResultMapper::MapObject::B )].old2newVerts.empty() )