Method run

source/MRMesh/MRPolylineDecimate.cpp:295–352

template<typename V>
DecimatePolylineResult PolylineDecimator<V>::run()
{
    MR_TIMER;

    initializeQueue_();

    res_.errorIntroduced = settings_.maxError;
    while ( !queue_.empty() )
    {
        auto topQE = queue_.top();
        assert( presentInQueue_.test( topQE.uedgeId ) );
        queue_.pop();
        if( res_.vertsDeleted >= settings_.maxDeletedVertices )
        {
            res_.errorIntroduced = std::sqrt( topQE.c );
            break;
        }

        if ( polyline_.topology.isLoneEdge( topQE.uedgeId ) )
        {
            // edge has been deleted by this moment
            presentInQueue_.reset( topQE.uedgeId );
            continue;
        }

        QuadraticForm<V> collapseForm;
        V collapsePos;
        auto qe = computeQueueElement_( topQE.uedgeId, &collapseForm, &collapsePos );
        if ( !qe )
        {
            presentInQueue_.reset( topQE.uedgeId );
            continue;
        }

        if ( qe->c > topQE.c )
        {
            queue_.push( *qe );
            continue;
        }

        presentInQueue_.reset( topQE.uedgeId );
        VertId collapseVert = collapse_( topQE.uedgeId, collapsePos );
        if ( !collapseVert )
            continue;

        vertForms_[collapseVert] = collapseForm;

        EdgeId e = polyline_.topology.edgeWithOrg( collapseVert );
        addInQueueIfMissing_( e.undirected() );
        EdgeId eNext = polyline_.topology.next( e );
        if ( e != eNext )
            addInQueueIfMissing_( eNext.undirected() );
    }

    if ( settings_.vertForms )
        *settings_.vertForms = std::move( vertForms_ );
    return res_;
}