Function constantCuspToolPath

source/MRVoxels/MRToolPath.cpp:1022–1357

}

Expected<ToolPathResult> constantCuspToolPath( const MeshPart& mp, const ConstantCuspParams& params )
{
    ToolPathResult  res;

    if ( !params.offsetMesh )
    {
        auto preprocessedMesh = preprocessMesh( mp.mesh, params, false );
        if ( !preprocessedMesh )
            return unexpected( preprocessedMesh.error() );

        res.modifiedMesh = std::move( *preprocessedMesh );
    }
    else
    {
        res.modifiedMesh = params.offsetMesh->mesh;
    }
    const auto box = res.modifiedMesh.computeBoundingBox();

    const Vector3f normal = Vector3f::plusZ();
    float minZ = box.min.z + params.millRadius;
    float safeZ = std::max( box.max.z + params.millRadius, params.safeZ );

    const auto undercutPlane = MR::Plane3f::fromDirAndPt( normal, { 0.0f, 0.0f, minZ } );
    
    //compute the lowest contour that might be processed
    const auto undercutSection = extractPlaneSections( res.modifiedMesh, undercutPlane ).front();
    Polyline3 undercutPolyline;
    undercutPolyline.addFromSurfacePath( res.modifiedMesh, undercutSection );
    const auto undercutContour = undercutPolyline.contours().front();

    // if there are multiple independent zones selected we need to process them separately
    const auto processZone = [&] ( const std::vector<SurfacePath>& startSurfacePaths, Vector3f lastPoint, ProgressCallback cb ) -> std::string
    {
        ExtractIsolinesParams extractionParams
        {
            .startSurfacePaths = startSurfacePaths,
            .startContours = params.startContours,
            .startVertices = params.startVertices,
            .sectionStep = params.sectionStep,
            .bypassDir = params.bypassDir,
            .cb = subprogress( cb, 0.0f, 0.4f )
        };
        //compute isolines based on the start point or the bounding contour
        auto extractRes = extractAllIsolines( res.modifiedMesh, extractionParams );
        if ( !extractRes.has_value() )
            return extractRes.error();

        auto& extract = *extractRes;
        res.modifiedMesh = extract.meshAfterCut;
        
        if ( !extract.old2NewMap.empty() )
        {
            const auto oldRegion = res.modifiedRegion;
            res.modifiedRegion = extract.region;

            BitSetParallelFor( oldRegion, [&] ( FaceId f )
            {
                for ( auto& newFace : extract.old2NewMap[f] )