Function findMaxDistanceSqOneWay

source/MRMesh/MRPointCloudDistance.cpp:11–38

{

float findMaxDistanceSqOneWay( const PointCloud& a, const PointCloud& b, const AffineXf3f* rigidB2A, float maxDistanceSq )
{
    MR_TIMER;

    return tbb::parallel_reduce
    (
        tbb::blocked_range( 0_v, b.validPoints.endId() ),
        0.0f,
        [&] ( const tbb::blocked_range<VertId>& range, float init )
        {
            for ( VertId i = range.begin(); i < range.end(); ++i )
            {
                if ( !b.validPoints.test( i ) )
                    continue;

                auto distSq = findProjectionOnPoints( rigidB2A ? (*rigidB2A)( b.points[i] ) : b.points[i], PointCloudPart{ a }, maxDistanceSq ).distSq;
                if ( distSq > init )
                    init = distSq;
            }

            return  init;
        },
        [] ( float a, float b ) -> float
        {
            return a > b ? a : b;
        }
    );
}

float findMaxDistanceSq( const PointCloud& a, const PointCloud& b, const AffineXf3f* rigidB2A, float maxDistanceSq )
{