CDT.cpp

Go to the documentation of this file.

#include "CDT.hpp"
 
#include <wmtk/Mesh.hpp>
#include <wmtk/TetMesh.hpp>
#include <wmtk/TriMesh.hpp>
 
// #include <wmtk/components/multimesh_from_tag/internal/MultiMeshFromTag.hpp>
#include <wmtk/operations/attribute_update/AttributeTransferStrategy.hpp>
 
#include "internal/CDT.hpp"
// #include "internal/CDTOptions.hpp"
 
namespace wmtk {
namespace components {
 
using namespace internal;
 
std::shared_ptr<Mesh> CDT(const TriMesh& trimesh, const bool inner_only, const bool rational_output)
{
    constexpr static PrimitiveType PV = PrimitiveType::Vertex;
    constexpr static PrimitiveType PE = PrimitiveType::Edge;
    constexpr static PrimitiveType PF = PrimitiveType::Triangle;
    constexpr static PrimitiveType PT = PrimitiveType::Tetrahedron;
 
    std::vector<std::array<bool, 4>> local_f_on_input;
 
    wmtk::logger().info("start CDT ...");
 
    std::shared_ptr<TetMesh> tm =
        CDT_internal(trimesh, local_f_on_input, inner_only, rational_output);
 
    wmtk::logger().info("finished CDT");
 
    // auto surface_handle = tm->register_attribute<int64_t>("surface", PrimitiveType::Triangle, 1);
    // auto surface_accessor = tm->create_accessor<int64_t>(surface_handle);
 
    // const auto& tets = tm->get_all(PrimitiveType::Tetrahedron);
 
    // for (int64_t i = 0; i < tets.size(); ++i) {
    //     const auto& t = tets[i];
    //     std::array<Tuple, 4> fs = {
    //         {tm->Mesh::switch_tuples(t, {PV, PE, PF}),
    //          tm->Mesh::switch_tuples(t, {PE, PF}),
    //          t,
    //          tm->Mesh::switch_tuples(t, {PF})}};
 
    //     for (int64_t k = 0; k < 4; ++k) {
    //         if (local_f_on_input[i][k]) {
    //             surface_accessor.scalar_attribute(fs[k]) = 1;
    //         } else {
    //             surface_accessor.scalar_attribute(fs[k]) = 0;
    //         }
    //     }
    // }
 
    // std::shared_ptr<Mesh> surface_mesh;
    // internal::MultiMeshFromTag mmft(*tm, surface_handle, 1);
    // mmft.compute_substructure_mesh();
    // surface_mesh = tm->get_child_meshes().back();
    // mmft.remove_soup();
 
    // wmtk::logger().info("registered surface child mesh to tetmesh {}", options.output);
 
    // // propagate position to all child meshes
    // auto pt_attribute = tm->get_attribute_handle<Rational>("vertices", PrimitiveType::Vertex);
 
    // for (auto child : tm->get_child_meshes()) {
    //     auto child_position_handle = child->register_attribute<Rational>(
    //         "vertices",
    //         PrimitiveType::Vertex,
    //         tm->get_attribute_dimension(pt_attribute.as<Rational>()));
 
    //     auto propagate_to_child_position =
    //         [](const Eigen::MatrixX<Rational>& P) -> Eigen::VectorX<Rational> { return P; };
    //     auto update_child_positon =
    //         std::make_shared<wmtk::operations::SingleAttributeTransferStrategy<Rational,
    //         Rational>>(
    //             child_position_handle,
    //             pt_attribute,
    //             propagate_to_child_position);
    //     update_child_positon->run_on_all();
    // }
 
    return tm;
}
 
} // namespace components
} // namespace wmtk