GU_WindingNumber.h

Go to the documentation of this file.

#ifndef __GU_WindingNumber_h__
#define __GU_WindingNumber_h__

#include "GU_API.h"

#include <GEO/GEO_Detail.h>
#include <GA/GA_SplittableRange.h>
#include <UT/UT_SolidAngle.h>


class GU_API GU_WindingNumber3DExact
{
public:
    /// Initializes WindingNumber3DExact class
    /// 
    /// NOTE: The class keeps non-owning pointer `mesh_geo` to the geometry.
    ///       The pointer has to stay valid as long as you want to call `eval`
    void init(const GEO_Detail *const mesh_geo,
              const GA_PrimitiveGroup *const mesh_prim_group);

    void clear();

    double eval(const UT_Vector3 &query_point,
                const bool as_solid_angle = false,
                const bool negate = false) const;

    void eval(const GEO_Detail *const query_points,
              const GA_SplittableRange &point_range,
              const GA_RWHandleF &winding_number_attrib,
              const bool as_solid_angle = false,
              const bool negate = false) const;

private:
    const GEO_Detail * myMeshGeo = nullptr;
    GA_OffsetList myPrimoffs{};
};

class GU_API GU_WindingNumber3DApprox
{
public:
    void init(const GA_Detail &mesh_geo,
              const GA_PrimitiveGroup *const mesh_prim_group,
              const int approx_order);

    void clear();

    bool isClear() const { return mySolidAngleTree.isClear(); }    
    

    double eval(const UT_Vector3 &query_point,
                const double accuracy_scale,
                const bool as_solid_angle = false,
                const bool negate = false) const;

    void eval(const GEO_Detail *const query_points,
              const GA_SplittableRange &point_range,
              const GA_RWHandleF &winding_number_attrib,
              const double accuracy_scale,            
              const bool as_solid_angle = false,
              const bool negate = false) const;

private:
    UT_SolidAngle<float,float> mySolidAngleTree{};
    UT_Array<int> myTrianglePoints{};
    UT_Array<UT_Vector3> myPositions3D{};
};


class GU_API GU_WindingNumber2DExact
{
public:
    /// Initializes WindingNumber2DExact class
    /// 
    /// NOTE: The class keeps non-owning pointer `mesh_geo` to the geometry.
    ///       The pointer has to stay valid as long as you want to call `eval`
    void init(const GEO_Detail *const mesh_geo,
              const GA_PrimitiveGroup *const mesh_prim_group,
              const int axis0, const int axis1);

    void clear();

    double eval(const UT_Vector2 &query_point,
                const bool as_angle = false,
                const bool negate = false) const;

    void eval(const GEO_Detail *const query_points,
              const GA_SplittableRange &point_range,
              const GA_RWHandleF &winding_number_attrib,
              const bool as_angle = false,
              const bool negate = false) const;

private:
    const GEO_Detail * myMeshGeo = nullptr;
    GA_OffsetList myPrimoffs{};
    int myAxis0, myAxis1;
};


class GU_API GU_WindingNumber2DApprox
{
public:
    void init(const GA_Detail &mesh_geo,
              const GA_PrimitiveGroup *const prim_group,
              const int approx_order,
              const int axis0, const int axis1);

    void clear();

    bool isClear() const { return mySubtendedAngleTree.isClear(); }    

    double eval(const UT_Vector2 &query_point,
                const double accuracy_scale,
                const bool as_angle = false,
                const bool negate = false) const;

    void eval(const GEO_Detail *const query_points,
              const GA_SplittableRange &point_range,
              const GA_RWHandleF &winding_number_attrib,
              const double accuracy_scale,            
              const bool as_angle = false,
              const bool negate = false) const;

private:
    UT_SubtendedAngle<float,float> mySubtendedAngleTree{};
    UT_Array<int> mySegmentPoints{};
    UT_Array<UT_Vector2> myPositions2D{};
    int myAxis0, myAxis1;
};

#endif