Find boundary intersections about meshlib HOT 2 CLOSED

MeshEdgePoint a,b; // input data
Mesh mesh; // input mesh
UndirectedEdgeBitSet selectedEdges; // input edges to find intersections on

//define plane
auto plane = Plane3f::fromDirAndPt( cross( mesh.edgePoint(b) - mesh.edgePoint(a), mesh.normal( MeshTriPoint(a) ) ), mesh.edgePoint(a) );

// find intersections
std::vector<MeshEdgePoint> intersections;
for ( auto ue : selectedEdges )
{
    float orgD = plane.distance( mesh.orgPt( EdgeId( ue ) ) );
    float destD = plane.distance( mesh.destPt( EdgeId( ue ) ) );
    if ( orgD * destD > 0.0f ) // same sign - no intersection
        continue;
    intersections.push_back( MeshEdgePoint( EdgeId( ue ), std::abs(orgD)/(std::abs(orgD)+std::abs(destD)) ) );
}