MRPrecisePredicates3.h

Go to the documentation of this file.

#pragma once
 
#include "MRVector3.h"
#include "MRId.h"
#include <array>
#include <optional>
 
namespace MR
{
 


MRMESH_API bool orient3d( const Vector3i & a, const Vector3i & b, const Vector3i & c );

inline bool orient3d( const Vector3i & a, const Vector3i & b, const Vector3i & c, const Vector3i & d )
    { return orient3d( a - d, b - d, c - d ); }
 
struct PreciseVertCoords
{
    VertId id;   
    Vector3i pt; 
};

MRMESH_API bool orient3d( const std::array<PreciseVertCoords, 4> & vs );
MRMESH_API bool orient3d( const PreciseVertCoords* vs );
 
struct TriangleSegmentIntersectResult
{
    bool doIntersect = false;    
    bool dIsLeftFromABC = false; 
 
    explicit operator bool() const { return doIntersect; }
};

[[nodiscard]] MRMESH_API TriangleSegmentIntersectResult doTriangleSegmentIntersect(
    const std::array<PreciseVertCoords, 5> & vs );

[[nodiscard]] MRMESH_API bool segmentIntersectionOrder( const std::array<PreciseVertCoords, 8> & vs );

struct ConvertToIntVector
{
    Vector3d center;
    double invRange = 1;

    [[nodiscard]] int scaleOnly( double v ) const
    {
        return (int)std::round( v * invRange );
    }

    [[nodiscard]] Vector3i operator()( const Vector3d& p ) const
    {
        const auto d = p - center;
        return { scaleOnly( d.x ), scaleOnly( d.y ), scaleOnly( d.z ) };
    }

    [[nodiscard]] Vector3i operator()( const Vector3f& p ) const
    {
        return operator()( Vector3d{ p } );
    }
};

struct ConvertToFloatVector
{
    double range = 1;
    Vector3d center;

    [[nodiscard]] Vector3d convert( const Vector3i& v ) const
    {
        return Vector3d{ v } * range + center;
    }

    [[nodiscard]] Vector3f operator()( const Vector3i& v ) const
    {
        return Vector3f( convert( v ) );
    }
};

struct CoordinateConverters
{
    ConvertToIntVector toInt;
    ConvertToFloatVector toFloat;
};

MRMESH_API ConvertToIntVector getToIntConverter( const Box3d& box );
MRMESH_API ConvertToFloatVector getToFloatConverter( const Box3d& box );

[[nodiscard]] MRMESH_API Vector<Vector3i, VertId> computeIntCoords( const ConvertToIntVector& conv,
    const VertCoords& points, const VertBitSet* valid = nullptr );

[[nodiscard]] MRMESH_API VertCoords computeFloatCoords( const ConvertToFloatVector& conv,
    const Vector<Vector3i, VertId>& intCoords, const VertBitSet* valid = nullptr );

[[nodiscard]] MRMESH_API std::optional<Vector3i> findTwoSegmentsIntersection( const Vector3i& ai, const Vector3i& bi, const Vector3i& ci, const Vector3i& di );

[[nodiscard]] MRMESH_API Vector3f findTriangleSegmentIntersectionPrecise( 
    const Vector3f& a, const Vector3f& b, const Vector3f& c,
    const Vector3f& d, const Vector3f& e, 
    CoordinateConverters converters );

 
}