MRFastWindingNumber.h

Go to the documentation of this file.

#pragma once
 
#include "MRProgressCallback.h"
#include "MRExpected.h"
#include "MRId.h"
 
namespace MR
{
 

class IFastWindingNumber
{
public:
    IFastWindingNumber() = default;
    IFastWindingNumber( const IFastWindingNumber & ) = default;
    IFastWindingNumber( IFastWindingNumber && ) noexcept = default;
 
    virtual ~IFastWindingNumber() = default;

    virtual Expected<void> calcFromVector( std::vector<float>& res, const std::vector<Vector3f>& points, float beta, FaceId skipFace = {}, const ProgressCallback& cb = {} ) = 0;

    virtual Expected<void> calcSelfIntersections( FaceBitSet& res, float beta, const ProgressCallback& cb = {} ) = 0;

    virtual Expected<void> calcFromGrid( std::vector<float>& res, const Vector3i& dims, const AffineXf3f& gridToMeshXf, float beta, const ProgressCallback& cb = {} ) = 0;

    virtual Expected<void> calcFromGridWithDistances( std::vector<float>& res, const Vector3i& dims, const AffineXf3f& gridToMeshXf, const DistanceToMeshOptions& options, const ProgressCallback& cb ) = 0;
};

class MRMESH_CLASS FastWindingNumber : public IFastWindingNumber
{
public:
    [[nodiscard]] MRMESH_API FastWindingNumber( const Mesh & mesh );

    MRMESH_API Expected<void> calcFromVector( std::vector<float>& res, const std::vector<Vector3f>& points, float beta, FaceId skipFace, const ProgressCallback& cb ) override;
    MRMESH_API Expected<void> calcSelfIntersections( FaceBitSet& res, float beta, const ProgressCallback& cb ) override;
    MRMESH_API Expected<void> calcFromGrid( std::vector<float>& res, const Vector3i& dims, const AffineXf3f& gridToMeshXf, float beta, const ProgressCallback& cb ) override;
    MRMESH_API float calcWithDistances( const Vector3f& p, const DistanceToMeshOptions& options );
    MRMESH_API Expected<void> calcFromGridWithDistances( std::vector<float>& res, const Vector3i& dims, const AffineXf3f& gridToMeshXf, const DistanceToMeshOptions& options, const ProgressCallback& cb ) override;
 
private:
    [[nodiscard]] float calc_( const Vector3f & q, float beta, FaceId skipFace = {} ) const;
 
    const Mesh & mesh_;
    const AABBTree & tree_;
    const Dipoles & dipoles_;
};

class IFastWindingNumberByParts
{
public:
    IFastWindingNumberByParts() = default;
    IFastWindingNumberByParts( const IFastWindingNumberByParts& ) = default;
    IFastWindingNumberByParts( IFastWindingNumberByParts&& ) noexcept = default;
 
    virtual ~IFastWindingNumberByParts() = default;

    using GridByPartsFunc = std::function<Expected<void> ( std::vector<float>&& data, const Vector3i& dims, int zOffset )>;

    virtual Expected<void> calcFromGridByParts( GridByPartsFunc resFunc, const Vector3i& dims,
        const AffineXf3f& gridToMeshXf, float beta, int layerOverlap, const ProgressCallback& cb ) = 0;

    virtual Expected<void> calcFromGridWithDistancesByParts( GridByPartsFunc resFunc, const Vector3i& dims,
        const AffineXf3f& gridToMeshXf, const DistanceToMeshOptions& options, int layerOverlap, const ProgressCallback& cb ) = 0;
};
 
}