Method calcSelfIntersections

source/MRCuda/MRCudaFastWindingNumber.cpp:129–173

}

Expected<void> FastWindingNumber::calcSelfIntersections( FaceBitSet& res, float beta, const ProgressCallback& cb )
{
    MR_TIMER;
    return prepareData_( subprogress( cb, 0.0, 0.5f ) ).and_then( [&]() -> Expected<void>
    {
        const auto totalSize = mesh_.topology.faceSize();
        const auto bufferSize = maxBufferSize( getCudaSafeMemoryLimit(), totalSize, sizeof( float ) );

        DynamicArrayF cudaResult;
        CUDA_LOGE_RETURN_UNEXPECTED( cudaResult.resize( bufferSize ) );
        if ( !reportProgress( cb, 0.6f ) )
            return unexpectedOperationCanceled();

        std::vector<float> wns;
        wns.resize( totalSize );

        const auto cb1 = subprogress( cb, 0.60f, 0.90f );
        const auto iterCount = chunkCount( totalSize, bufferSize );
        size_t iterIndex = 0;

        for ( const auto [offset, size] : splitByChunks( totalSize, bufferSize ) )
        {
            const auto cb2 = subprogress( cb1, iterIndex++, iterCount );

            fastWindingNumberFromMesh( data_->toData(), cudaResult.data(), beta, size, offset );
            CUDA_LOGE_RETURN_UNEXPECTED( cudaGetLastError() );
            if ( !reportProgress( cb2, 0.33f ) )
                return unexpectedOperationCanceled();

            CUDA_LOGE_RETURN_UNEXPECTED( cudaResult.copyTo( wns.data() + offset, size ) );
            if ( !reportProgress( cb2, 1.00f ) )
                return unexpectedOperationCanceled();
        }

        res.resize( totalSize );
        if ( !BitSetParallelForAll( res, [&] (FaceId f)
        {
            if ( wns[f] < 0 || wns[f] > 1 )
                res.set( f );
        }, subprogress( cb, 0.9f, 1.0f ) ) )
            return unexpectedOperationCanceled();

        return {};
    } );
}

Expected<void> FastWindingNumber::calcFromGrid( std::vector<float>& res, const Vector3i& dims, const AffineXf3f& gridToMeshXf, float beta, const ProgressCallback& cb )
{