MRDipole.h

Go to the documentation of this file.

#pragma once
 
#include "MRVector3.h"
 
namespace MR
{
 
struct Dipole
{
    Vector3f pos;
    float area = 0;
    Vector3f dirArea;
    float rr = 0; // maximum squared distance from pos to any corner of the bounding box
    [[nodiscard]] bool addIfGoodApprox( const Vector3f& q, float betaSq, float& addTo ) const
    {
        const auto dp = pos - q;
        const auto dd = dp.lengthSq();
        if ( dd <= betaSq * rr )
            return false;
        if ( const auto d = std::sqrt( dd ); d > 0 )
            addTo += dot( dp, dirArea ) / ( d * dd );
        return true;
    }
};
 
static_assert( sizeof( Dipole ) == 8 * sizeof( float ) );
 
MRMESH_API void calcDipoles( Dipoles& dipoles, const AABBTree& tree, const Mesh& mesh );
[[nodiscard]] MRMESH_API Dipoles calcDipoles( const AABBTree& tree, const Mesh& mesh );
 
[[nodiscard]] MRMESH_API float calcFastWindingNumber( const Dipoles& dipoles, const AABBTree& tree, const Mesh& mesh,
    const Vector3f & q, float beta, FaceId skipFace );
 
} //namespace MR