DelaunayMesh.h

#pragma once

#include <Euclid/MeshUtil/MeshDefs.h>

namespace Euclid
{
template<typename Mesh>
bool is_delaunay(const Mesh& mesh, double eps = 1e-14);

template<typename Mesh>
bool is_delaunay(const Mesh& mesh, edge_t<Mesh> e, double eps = 1e-14);

enum class RemeshDelaunayScheme
{
    SimpleFlip,

    GeometryPreserving,

    FeaturePreserving,
};

template<typename Mesh>
struct SplitSite
{
    edge_t<Mesh> epq;
    vertex_t<Mesh> vs;
    halfedge_t<Mesh> hps;
    halfedge_t<Mesh> hqs;
    halfedge_t<Mesh> hus;
    halfedge_t<Mesh> hvs;
};

template<typename Mesh>
struct RemeshDelaunayVisitor
{
    void on_started(Mesh&) {}

    void on_finished(Mesh&) {}

    void on_flipping(Mesh&, const edge_t<Mesh>&) {}

    void on_flipped(Mesh&, const edge_t<Mesh>&) {}

    void on_nonflippable(Mesh&, const edge_t<Mesh>&) {}

    void on_splitting(Mesh&, const edge_t<Mesh>&) {}

    void on_split(Mesh&, const SplitSite<Mesh>&) {}
};

template<typename Mesh>
void remesh_delaunay(
    Mesh& mesh,
    RemeshDelaunayScheme scheme = RemeshDelaunayScheme::SimpleFlip,
    double dihedral_angle = 10.0);

template<typename Mesh, typename Visitor>
void remesh_delaunay(
    Mesh& mesh,
    Visitor& visitor,
    RemeshDelaunayScheme scheme = RemeshDelaunayScheme::SimpleFlip,
    double dihedral_angle = 10.0);

} // namespace Euclid

#include "src/DelaunayMesh.cpp"