MRMeshMath.h

Go to the documentation of this file.

#pragma once
 
#include "MRVector.h"
#include "MRMeshTopology.h"
#include "MRLineSegm.h"
#include "MRPointOnFace.h"
 
namespace MR
{
 
[[nodiscard]] inline Vector3f orgPnt( const MeshTopology & topology, const VertCoords & points, EdgeId e )
{
    return points[ topology.org( e ) ];
}
 
[[nodiscard]] inline Vector3f destPnt( const MeshTopology & topology, const VertCoords & points, EdgeId e )
{
    return points[ topology.dest( e ) ];
}
 
[[nodiscard]] inline Vector3f edgeVector( const MeshTopology & topology, const VertCoords & points, EdgeId e )
{
    return destPnt( topology, points, e ) - orgPnt( topology, points, e );
}
 
[[nodiscard]] inline LineSegm3f edgeSegment( const MeshTopology & topology, const VertCoords & points, EdgeId e )
{
    return { orgPnt( topology, points, e ), destPnt( topology, points, e ) };
}
 
[[nodiscard]] inline Vector3f edgePoint( const MeshTopology & topology, const VertCoords & points, EdgeId e, float f )
{
    return f * destPnt( topology, points, e ) + ( 1 - f ) * orgPnt( topology, points, e );
}
 
[[nodiscard]] inline Vector3f edgePoint( const MeshTopology & topology, const VertCoords & points, const MeshEdgePoint & ep )
{
    return edgePoint( topology, points, ep.e, ep.a );
}
 
[[nodiscard]] inline Vector3f edgeCenter( const MeshTopology & topology, const VertCoords & points, UndirectedEdgeId e )
{
    return edgePoint( topology, points, e, 0.5f );
}
 
MRMESH_API void getLeftTriPoints( const MeshTopology & topology, const VertCoords & points, EdgeId e, Vector3f & v0, Vector3f & v1, Vector3f & v2 );
 
inline void getLeftTriPoints( const MeshTopology & topology, const VertCoords & points, EdgeId e, Vector3f (&v)[3] )
{
    getLeftTriPoints( topology, points, e, v[0], v[1], v[2] );
}
 
[[nodiscard]] inline Triangle3f getLeftTriPoints( const MeshTopology & topology, const VertCoords & points, EdgeId e )
{
    Triangle3f res;
    getLeftTriPoints( topology, points, e, res[0], res[1], res[2] );
    return res;
}
 
inline void getTriPoints( const MeshTopology & topology, const VertCoords & points, FaceId f, Vector3f & v0, Vector3f & v1, Vector3f & v2 )
{
    getLeftTriPoints( topology, points, topology.edgeWithLeft( f ), v0, v1, v2 );
}
 
inline void getTriPoints( const MeshTopology & topology, const VertCoords & points, FaceId f, Vector3f (&v)[3] )
{
    getTriPoints( topology, points, f, v[0], v[1], v[2] );
}
 
[[nodiscard]] inline Triangle3f getTriPoints( const MeshTopology & topology, const VertCoords & points, FaceId f )
{
    Triangle3f res;
    getTriPoints( topology, points, f, res[0], res[1], res[2] );
    return res;
}
 
[[nodiscard]] MRMESH_API Vector3f triPoint( const MeshTopology & topology, const VertCoords & points, const MeshTriPoint & p );
 
[[nodiscard]] MRMESH_API Vector3f triCenter( const MeshTopology & topology, const VertCoords & points, FaceId f );
 
[[nodiscard]] MRMESH_API float triangleAspectRatio( const MeshTopology & topology, const VertCoords & points, FaceId f );
    
[[nodiscard]] MRMESH_API float circumcircleDiameterSq( const MeshTopology & topology, const VertCoords & points, FaceId f );
 
[[nodiscard]] MRMESH_API float circumcircleDiameter( const MeshTopology & topology, const VertCoords & points, FaceId f );
 
[[nodiscard]] MRMESH_API MeshTriPoint toTriPoint( const MeshTopology & topology, const VertCoords & points, FaceId f, const Vector3f & p );
 
[[nodiscard]] MRMESH_API MeshTriPoint toTriPoint( const MeshTopology & topology, const VertCoords & points, const PointOnFace& p );
 
[[nodiscard]] MRMESH_API MeshEdgePoint toEdgePoint( const MeshTopology & topology, const VertCoords & points, EdgeId e, const Vector3f & p );
 
[[nodiscard]] MRMESH_API VertId getClosestVertex( const MeshTopology & topology, const VertCoords & points, const PointOnFace & p );
 
[[nodiscard]] inline VertId getClosestVertex( const MeshTopology & topology, const VertCoords & points, const MeshTriPoint & p )
{
    return getClosestVertex( topology, points, PointOnFace{ topology.left( p.e ), triPoint( topology, points, p ) } );
}
 
[[nodiscard]] MRMESH_API UndirectedEdgeId getClosestEdge( const MeshTopology & topology, const VertCoords & points, const PointOnFace & p );
 
[[nodiscard]] inline UndirectedEdgeId getClosestEdge( const MeshTopology & topology, const VertCoords & points, const MeshTriPoint & p )
{
    return getClosestEdge( topology, points, PointOnFace{ topology.left( p.e ), triPoint( topology, points, p ) } );
}
 
[[nodiscard]] inline float edgeLength( const MeshTopology & topology, const VertCoords & points, UndirectedEdgeId e )
{
    return edgeVector( topology, points, e ).length();
}
 
[[nodiscard]] MRMESH_API UndirectedEdgeScalars edgeLengths( const MeshTopology & topology, const VertCoords & points );
 
[[nodiscard]] inline float edgeLengthSq( const MeshTopology & topology, const VertCoords & points, UndirectedEdgeId e )
{
    return edgeVector( topology, points, e ).lengthSq();
}
 
[[nodiscard]] MRMESH_API Vector3f leftDirDblArea( const MeshTopology & topology, const VertCoords & points, EdgeId e );
 
[[nodiscard]] inline Vector3f dirDblArea( const MeshTopology & topology, const VertCoords & points, FaceId f )
{
    return leftDirDblArea( topology, points, topology.edgeWithLeft( f ) );
}
 
[[nodiscard]] MRMESH_API Vector<Vector3f, VertId> dirDblAreas( const MeshTopology & topology, const VertCoords & points, const VertBitSet * region = nullptr );
 
[[nodiscard]] inline float dblArea( const MeshTopology & topology, const VertCoords & points, FaceId f )
{
    return dirDblArea( topology, points, f ).length();
}
 
[[nodiscard]] inline float area( const MeshTopology & topology, const VertCoords & points, FaceId f )
{
    return 0.5f * dblArea( topology, points, f );
}
 
[[nodiscard]] MRMESH_API double area( const MeshTopology & topology, const VertCoords & points, const FaceBitSet & fs );
 
[[nodiscard]] inline double area( const MeshTopology & topology, const VertCoords & points, const FaceBitSet * fs = nullptr )
{
    return area( topology, points, topology.getFaceIds( fs ) );
}
 
[[nodiscard]] MRMESH_API Vector3d dirArea( const MeshTopology & topology, const VertCoords & points, const FaceBitSet & fs );
 
[[nodiscard]] inline Vector3d dirArea( const MeshTopology & topology, const VertCoords & points, const FaceBitSet * fs = nullptr )
{
    return dirArea( topology, points, topology.getFaceIds( fs ) );
}
 
[[nodiscard]] MRMESH_API double projArea( const MeshTopology & topology, const VertCoords & points, const Vector3f & dir, const FaceBitSet & fs );
 
[[nodiscard]] inline double projArea( const MeshTopology & topology, const VertCoords & points, const Vector3f & dir, const FaceBitSet * fs = nullptr )
{
    return projArea( topology, points, dir, topology.getFaceIds( fs ) );
}
 
[[nodiscard]] MRMESH_API double volume( const MeshTopology & topology, const VertCoords & points, const FaceBitSet* region = nullptr );
 
[[nodiscard]] MRMESH_API double holePerimiter( const MeshTopology & topology, const VertCoords & points, EdgeId e );
 
[[nodiscard]] MRMESH_API Vector3d holeDirArea( const MeshTopology & topology, const VertCoords & points, EdgeId e );
 
[[nodiscard]] MRMESH_API Vector3f leftTangent( const MeshTopology & topology, const VertCoords & points, EdgeId e );
 
[[nodiscard]] inline Vector3f leftNormal( const MeshTopology & topology, const VertCoords & points, EdgeId e )
{
    return leftDirDblArea( topology, points, e ).normalized();
}
 
[[nodiscard]] inline Vector3f normal( const MeshTopology & topology, const VertCoords & points, FaceId f )
{
    return dirDblArea( topology, points, f ).normalized();
}
 
[[nodiscard]] MRMESH_API Vector3f dirDblArea( const MeshTopology & topology, const VertCoords & points, VertId v );
 
[[nodiscard]] inline float dblArea( const MeshTopology & topology, const VertCoords & points, VertId v )
{
    return dirDblArea( topology, points, v ).length();
}
 
[[nodiscard]] inline Vector3f normal( const MeshTopology & topology, const VertCoords & points, VertId v )
{
    return dirDblArea( topology, points, v ).normalized();
}
 
[[nodiscard]] MRMESH_API Vector3f normal( const MeshTopology & topology, const VertCoords & points, const MeshTriPoint & p );
 
[[nodiscard]] MRMESH_API Vector3f pseudonormal( const MeshTopology & topology, const VertCoords & points, VertId v, const FaceBitSet * region = nullptr );
 
[[nodiscard]] MRMESH_API Vector3f pseudonormal( const MeshTopology & topology, const VertCoords & points, UndirectedEdgeId e, const FaceBitSet * region = nullptr );
 
[[nodiscard]] MRMESH_API Vector3f pseudonormal( const MeshTopology & topology, const VertCoords & points, const MeshTriPoint & p, const FaceBitSet * region = nullptr );
 
[[nodiscard]] MRMESH_API float sumAngles( const MeshTopology & topology, const VertCoords & points, VertId v, bool * outBoundaryVert = nullptr );
 
[[nodiscard]] MRMESH_API Expected<VertBitSet> findSpikeVertices( const MeshTopology & topology, const VertCoords & points, float minSumAngle, const VertBitSet* region = nullptr, const ProgressCallback& cb = {} );
 
[[nodiscard]] MRMESH_API float dihedralAngleSin( const MeshTopology & topology, const VertCoords & points, UndirectedEdgeId e );
 
[[nodiscard]] MRMESH_API float dihedralAngleCos( const MeshTopology & topology, const VertCoords & points, UndirectedEdgeId e );
 
[[nodiscard]] MRMESH_API float dihedralAngle( const MeshTopology & topology, const VertCoords & points, UndirectedEdgeId e );
 
[[nodiscard]] MRMESH_API float discreteMeanCurvature( const MeshTopology & topology, const VertCoords & points, VertId v );
 
[[nodiscard]] MRMESH_API float discreteMeanCurvature( const MeshTopology & topology, const VertCoords & points, UndirectedEdgeId e );
 
[[nodiscard]] inline float discreteGaussianCurvature( const MeshTopology & topology, const VertCoords & points, VertId v, bool * outBoundaryVert = nullptr )
{
    return 2 * PI_F - sumAngles( topology, points, v, outBoundaryVert );
}
 
[[nodiscard]] MRMESH_API UndirectedEdgeBitSet findCreaseEdges( const MeshTopology & topology, const VertCoords & points, float angleFromPlanar );
 
[[nodiscard]] MRMESH_API float leftCotan( const MeshTopology & topology, const VertCoords & points, EdgeId e );
 
[[nodiscard]] inline float cotan( const MeshTopology & topology, const VertCoords & points, UndirectedEdgeId ue )
{
    EdgeId e{ ue };
    return leftCotan( topology, points, e ) + leftCotan( topology, points, e.sym() );
}
 
[[nodiscard]] MRMESH_API QuadraticForm3f quadraticForm( const MeshTopology & topology, const VertCoords & points, VertId v, bool angleWeigted,
    const FaceBitSet * region = nullptr, const UndirectedEdgeBitSet * creases = nullptr );
 
[[nodiscard]] MRMESH_API Box3f computeBoundingBox( const MeshTopology & topology, const VertCoords & points, const FaceBitSet* region, const AffineXf3f* toWorld = nullptr );
 
[[nodiscard]] MRMESH_API float averageEdgeLength( const MeshTopology & topology, const VertCoords & points );
 
[[nodiscard]] MRMESH_API Vector3f findCenterFromPoints( const MeshTopology & topology, const VertCoords & points );
 
[[nodiscard]] MRMESH_API Vector3f findCenterFromFaces( const MeshTopology & topology, const VertCoords & points );
 
[[nodiscard]] MRMESH_API Vector3f findCenterFromBBox( const MeshTopology & topology, const VertCoords & points );
 
} //namespace MR