Function reduceTotalAngle

source/MRMesh/MRMeshTotalAngle.cpp:134–205

} //anonymous namespace

int reduceTotalAngle( MeshTopology& topology, const VertCoords& points, int numIters, const ReduceTotalAngleParams& region, const ProgressCallback& progressCallback )
{
    if ( numIters <= 0 )
        return 0;
    MR_TIMER;

    UndirectedEdgeBitSet flipCandidates( topology.undirectedEdgeSize() );
    UndirectedEdgeBitSet nextFlipCandidates( topology.undirectedEdgeSize(), true );

    int flipsDone = 0;
    for ( int iter = 0; iter < numIters; ++iter )
    {
        if ( progressCallback && !progressCallback( float( iter ) / numIters ) )
            return flipsDone;

        flipCandidates.reset();
        BitSetParallelFor( nextFlipCandidates, [&] ( UndirectedEdgeId e )
        {
            if ( totalAngleIncreaseOnFlip( topology, points, e, region ) < 1 )
                flipCandidates.set( e );
        } );
        nextFlipCandidates.reset();
        int flipsDoneBeforeThisIter = flipsDone;
        for ( UndirectedEdgeId ue : flipCandidates )
        {
            const EdgeId e = ue;
            if ( totalAngleIncreaseOnFlip( topology, points, e, region ) >= 1 )
                continue;

            ++flipsDone;
            topology.flipEdge( e );

            if ( iter + 1 >= numIters )
                continue;

            const auto e0 = topology.prev( e );
            nextFlipCandidates.set( e0 );
            if ( topology.right( e0 ) )
            {
                nextFlipCandidates.set( topology.prev( e0 ) );
                nextFlipCandidates.set( topology.next( e0.sym() ) );
            }

            const auto e1 = topology.next( e.sym() );
            nextFlipCandidates.set( e1 );
            if ( topology.left( e1 ) )
            {
                nextFlipCandidates.set( topology.next( e1 ) );
                nextFlipCandidates.set( topology.prev( e1.sym() ) );
            }

            const auto e2 = topology.prev( e.sym() );
            nextFlipCandidates.set( e2 );
            if ( topology.right( e2 ) )
            {
                nextFlipCandidates.set( topology.prev( e2 ) );
                nextFlipCandidates.set( topology.next( e2.sym() ) );
            }