volume_complex.hpp

#pragma once
 
#include <array>
#include <cstdint>
#include <optional>
#include <set>
#include <string>
#include <utility>
#include <vector>
 
namespace wmtk::components::tracking::volume {
 
using Vid = uint32_t;
using Tet = std::array<Vid, 4>;
using Tri = std::array<Vid, 3>;
using Tets = std::vector<Tet>;
using Faces = std::set<Tri>;
 
Tri sorted(Tri t);
Tet sorted(Tet t);
bool contains(const Tet& t, Vid v);
bool contains(const Tri& t, Vid v);
std::array<Tri, 4> faces_of(const Tet& t);
std::set<Tet> keys(const Tets& T);
 
std::optional<Faces> boundary(const Tets& T);
bool is_ball(const Tets& T);
std::optional<std::vector<Vid>> disk_cycle(const Faces& F);
std::optional<Tets> collapse(const Tets& T, Vid a, Vid b);
std::vector<Tets> sides(const Tets& T, Vid c, const Faces& surf);
std::vector<Tets> side_groups(const Tets& T, const Faces& surf);
Faces surf_after(const Faces& surf, Vid i, Vid j);
 
struct PairSpec
{
    std::string kind;
    Tets before, after;
    Faces fixed, cap_before, cap_after;
};
std::optional<PairSpec> make_pair(
    const std::string& kind,
    const Tets& Tb,
    const Tets& Ta,
    const Faces& sb,
    const Faces& sa,
    std::string* why = nullptr);
 
char side_class(const Tets& M, Vid i, Vid j);
std::vector<Vid> chords(const Tets& M, Vid i, Vid j);
 
struct Decomposition
{
    std::vector<PairSpec> pairs;
    std::vector<std::pair<Vid, Vid>> virtuals;
    Tets split, collapsed;
};
std::optional<Decomposition> decompose_collapse(
    const Tets& K,
    Vid i,
    Vid j,
    const Faces& surf,
    Vid first_virtual,
    std::string* why);
 
} // namespace wmtk::components::tracking::volume