LocalOperations.h

// This file is part of TetWild, a software for generating tetrahedral meshes.
//
// Copyright (C) 2018 Yixin Hu <yixin.hu@nyu.edu>
//
// This Source Code Form is subject to the terms of the Mozilla Public License
// v. 2.0. If a copy of the MPL was not distributed with this file, You can
// obtain one at http://mozilla.org/MPL/2.0/.
//
// Created by Yixin Hu on 5/6/17.
//
 
#ifndef NEW_GTET_LOCALOPERATIONS_H
#define NEW_GTET_LOCALOPERATIONS_H
 
#include <igl/Timer.h>
#include <igl/grad.h>
#include <wmtk/Types.hpp>
#include <wmtk/envelope/Envelope.hpp>
#include "ForwardDecls.h"
#include "TetmeshElements.h"
 
namespace wmtk::components::tetwild::orig {
 
enum class EnvelopSide { OUTSIDE = 0, INSIDE = 1, UNCERTAIN = 2 };
 
class LocalOperations
{
public:
    const Args& args;
    State& state;
 
    std::vector<TetVertex>& tet_vertices;
    std::vector<std::array<int, 4>>& tets;
    std::vector<std::array<int, 4>>& is_surface_fs;
    std::vector<bool>& v_is_removed;
    std::vector<bool>& t_is_removed;
    std::vector<TetQuality>& tet_qualities;
 
    const SampleEnvelope& geo_sf_tree;
    const SampleEnvelope& geo_b_tree;
 
    int counter = 0;
    int suc_counter = 0;
 
    LocalOperations(
        std::vector<TetVertex>& t_vs,
        std::vector<std::array<int, 4>>& ts,
        std::vector<std::array<int, 4>>& is_sf_fs,
        std::vector<bool>& v_is_rm,
        std::vector<bool>& t_is_rm,
        std::vector<TetQuality>& tet_qs,
        const SampleEnvelope& geo_tree,
        const SampleEnvelope& b_t,
        const Args& ar,
        State& st)
        : tet_vertices(t_vs)
        , tets(ts)
        , is_surface_fs(is_sf_fs)
        , v_is_removed(v_is_rm)
        , t_is_removed(t_is_rm)
        , tet_qualities(tet_qs)
        , geo_sf_tree(geo_tree)
        , geo_b_tree(b_t)
        , args(ar)
        , state(st)
    {}
 
    void outputInfo(int op_type, double time, bool is_log = true);
 
    void calTetQualities(
        const std::vector<std::array<int, 4>>& new_tets,
        std::vector<TetQuality>& tet_qs);
 
    double calEdgeLength(int v1_id, int v2_id);
    void calTetQuality_AMIPS(const std::array<int, 4>& tet, TetQuality& t_quality);
 
    bool isFlip(const std::vector<std::array<int, 4>>& new_tets);
    bool isTetFlip(const std::array<int, 4>& t);
    bool isTetFlip(int t_id);
 
    void getAvgMaxEnergy(double& avg_tq, double& max_tq);
    double getMaxEnergy();
 
    void getCheckQuality(const std::vector<TetQuality>& tet_qs, TetQuality& tq);
    void getCheckQuality(const std::vector<int>& t_ids, TetQuality& tq);
 
    bool isEdgeOnSurface(int v1_id, int v2_id);
    bool isEdgeOnBbox(int v1_id, int v2_id);
    bool isEdgeOnSurface(int v1_id, int v2_id, const std::vector<int>& t_ids);
    bool isEdgeOnBbox(int v1_id, int v2_id, const std::vector<int>& t_ids);
    bool isEdgeOnBoundary(int v1_id, int v2_id);
 
    bool isFaceOutEnvelop(const std::array<Vector3d, 3>& tri);
    bool isPointOutEnvelop(const Vector3d& p);
    bool isPointOutBoundaryEnvelop(const Vector3d& p);
 
    void getFaceConnTets(int v1_id, int v2_id, int v3_id, std::vector<int>& t_ids);
    bool isIsolated(int v_id);
    bool isBoundaryPoint(int v_id);
 
    igl::Timer igl_timer0;
    int id_sampling = 0;
    int id_aabb = 1;
 
    int mid_id = 0;
 
    bool isLocked_ui(const std::array<int, 2>& e);
    bool isTetLocked_ui(int tid);
};
 
} // namespace wmtk::components::tetwild::orig
 
#endif // NEW_GTET_LOCALOPERATIONS_H