edge_insertion.hpp

Go to the documentation of this file.

#pragma once
 
#include <wmtk/EdgeMesh.hpp>
#include <wmtk/TriMesh.hpp>
 
#include <wmtk/Types.hpp>
#include <wmtk/utils/Rational.hpp>
 
 
namespace wmtk::components::internal {
 
struct bbox
{
    Rational x_min, x_max, y_min, y_max;
 
    bbox(const Vector2r& p0, const Vector2r& p1)
    {
        if (p0[0] < p1[0]) {
            x_min = p0[0];
            x_max = p1[0];
        } else {
            x_min = p1[0];
            x_max = p0[0];
        }
 
        if (p0[1] < p1[1]) {
            y_min = p0[1];
            y_max = p1[1];
        } else {
            y_min = p1[1];
            y_max = p0[1];
        }
    }
 
    bbox(const Vector2r& p0, const Vector2r& p1, const Vector2r& p2)
    {
        if (p0[0] < p1[0]) {
            x_min = p0[0];
            x_max = p1[0];
        } else {
            x_min = p1[0];
            x_max = p0[0];
        }
 
        if (p0[1] < p1[1]) {
            y_min = p0[1];
            y_max = p1[1];
        } else {
            y_min = p1[1];
            y_max = p0[1];
        }
 
        x_min = (x_min > p2[0]) ? p2[0] : x_min;
        x_max = (x_max < p2[0]) ? p2[0] : x_max;
        y_min = (y_min > p2[1]) ? p2[1] : y_min;
        y_max = (y_max < p2[1]) ? p2[1] : y_max;
    }
};
class Segment
{
public:
    const Vector2r p0, p1;
    const int64_t idx0, idx1;
    Vector2r bbox_min, bbox_max;
    bool deprecated = false;
    bool is_on_input = false;
 
    // std::vector<std::pair<int64_t, Rational>> points_on_segment;
    std::vector<std::pair<int64_t, Vector2r>> points_on_segment;
 
    Segment(const Vector2r& _p0, const Vector2r& _p1, const int64_t _idx0, const int64_t _idx1)
        : p0(_p0)
        , p1(_p1)
        , idx0(_idx0)
        , idx1(_idx1)
    {
        Rational x_min, x_max, y_min, y_max;
 
        if (p0[0] < p1[0]) {
            x_min = p0[0];
            x_max = p1[0];
        } else {
            x_min = p1[0];
            x_max = p0[0];
        }
 
        if (p0[1] < p1[1]) {
            y_min = p0[1];
            y_max = p1[1];
        } else {
            y_min = p1[1];
            y_max = p0[1];
        }
 
        bbox_min = Vector2r(x_min, y_min);
        bbox_max = Vector2r(x_max, y_max);
 
        points_on_segment.emplace_back(idx0, p0);
        points_on_segment.emplace_back(idx1, p1);
    }
};
 
int is_point_inside_triangle(
    const wmtk::Vector2r& P,
    const wmtk::Vector2r& A,
    const wmtk::Vector2r& B,
    const wmtk::Vector2r& C);
 
 
bool segment_segment_inter(
    const Vector2r& s0,
    const Vector2r& e0,
    const Vector2r& s1,
    const Vector2r& e1,
    Vector2r& res);
 
void edge_insertion(
    TriMesh& _trimesh,
    EdgeMesh& edgemesh,
    std::vector<Vector2r>& v_final,
    std::vector<std::array<int64_t, 3>>& FV_new,
    std::vector<std::array<int, 3>>& local_e_on_input);
 
} // namespace wmtk::components::internal