Delaunay.hpp

#pragma once
 
#include <array>
#include <cstddef>
#include <vector>
#include <wmtk/Types.hpp>
 
namespace wmtk::delaunay {
using Point3D = std::array<double, 3>;
using Tetrahedron = std::array<size_t, 4>;
using Point2D = std::array<double, 2>;
using Triangle = std::array<size_t, 3>;
 
auto delaunay3D(const std::vector<Point3D>& points)
    -> std::pair<std::vector<Point3D>, std::vector<Tetrahedron>>;
 
 
auto delaunay2D(const std::vector<Point2D>& points)
    -> std::pair<std::vector<Point2D>, std::vector<Triangle>>;
 
void constrained_delaunay2D(
    const MatrixXd& V,
    const MatrixXi& E,
    MatrixXd& V_out,
    MatrixXi& F_out,
    MatrixXi& E_out);
 
} // namespace wmtk::delaunay