Add cutMeshByPolyline function. (#4227)

* Add `cutMeshByPolyline` function.

* Remove unused parameters.

* PR fixes.

* Grammar.

* Make parallel and check cut.

* PR fixes.

Author: andrewerf

source/MRMesh/MRContoursCut.cpp

@@ -1,6 +1,7 @@
 #include "MRContoursCut.h"
 #include "MRAffineXf3.h"
 #include "MRMesh.h"
+#include "MRPolyline.h"
 #include "MRTriangleIntersection.h"
 #include "MRMeshTopology.h"
 #include "MRMeshDelone.h"
@@ -2190,6 +2191,37 @@ CutMeshResult cutMesh( Mesh& mesh, const OneMeshContours& contours, const CutMes
     return res;
 }
 
+Expected<FaceBitSet> cutMeshByContour( Mesh& mesh, const Contour3f& contour, const AffineXf3f& xf )
+{
+    std::vector<MeshTriPoint> surfaceLine( contour.size() );
+    tbb::task_group_context ctx;
+    bool ok = true;
+    ParallelFor( (size_t)0, contour.size(), [&] ( size_t i )
+    {
+        PointOnFace projPt;
+        if ( !mesh.projectPoint( xf( contour[i] ), projPt ) )
+        {
+            if ( ctx.cancel_group_execution() )
+                ok = false;
+            return;
+        }
+        surfaceLine[i] = mesh.toTriPoint( projPt );
+    } );
+    if ( !ok )
+        return unexpected( "Cannot project point to mesh" );
+
+    auto meshContour = convertMeshTriPointsToMeshContour( mesh, surfaceLine );
+    if ( !meshContour )
+        return unexpected( "Cannot convert tri points to mesh contour: " + meshContour.error() );
+
+    auto cutRes = cutMesh( mesh, { *meshContour } );
+    if ( !cutRes.fbsWithCountourIntersections.none() )
+        return unexpected( "Cannot cut mesh because of contour self intersections" );
+    auto sideFbv = fillContourLeft( mesh.topology, cutRes.resultCut );
+    return sideFbv;
+}
+
+
 TEST( MRMesh, BooleanIntersectionsSort )
 {
     Mesh meshA;

source/MRMesh/MRContoursCut.h

@@ -1,6 +1,7 @@
 #pragma once
 
 #include "MRVector3.h"
+#include "MRAffineXf3.h"
 #include "MRId.h"
 #include "MRBitSet.h"
 #include "MRIntersectionContour.h"
@@ -233,4 +234,11 @@ struct CutMeshResult
   */
 MRMESH_API CutMeshResult cutMesh( Mesh& mesh, const OneMeshContours& contours, const CutMeshParameters& params = {} );
 
+
+/// Cuts \p mesh by \p contour by projecting all the points
+/// \param xf transformation from the CSYS of \p contour to the CSYS of \p mesh
+/// \note \p mesh is modified, see \ref cutMesh for info
+/// \return Faces to the left of the polyline
+MRMESH_API Expected<FaceBitSet> cutMeshByContour( Mesh& mesh, const Contour3f& contour, const AffineXf3f& xf = {} );
+
 } //namespace MR