MMGRemesh.hpp

Go to the documentation of this file.

// Original source from cellogram (https://github.com/cellogram/cellogram/blob/master/src/cellogram/remesh_adaptive.h)
// Authors: Tobias Lendenmann, Teseo Schneider, Jérémie Dumas, Marco Tarini
// License: MIT (https://github.com/cellogram/cellogram/blob/master/LICENSE)
 
#pragma once
 
#ifdef POLYFEM_WITH_MMG
 
#include <Eigen/Dense>
 
#include <vector>
 
namespace polyfem::mesh
{
 
    // See MmgTools documentation for interpreation
    // https://www.mmgtools.org/mmg-remesher-try-mmg/mmg-remesher-options
    //
    struct MmgOptions
    {
        /* Remeshing */
        bool angle_detection = true;
        double angle_value = 45.;
        double hausd = 0.01;
        double hsiz = 0.0; /* using hmin and hmax if set to 0 */
        double hmin = 0.01;
        double hmax = 2.;
        double hgrad = 1.105171;
        bool enable_anisotropy = false;
        bool optim = false;
        bool optimLES = false;
        bool opnbdy = false;
        bool noinsert = false;
        bool noswap = false;
        bool nomove = false;
        bool nosurf = false;
        std::string metric_attribute = "no_metric";
        /* Level set extraction */
        bool level_set = false;
        std::string ls_attribute = "no_ls";
        double ls_value = 0.;
    };
 
    bool remesh_2d(const Eigen::MatrixXd &V, const Eigen::MatrixXi &F,
                   Eigen::MatrixXd &OV, Eigen::MatrixXi &OF,
                   MmgOptions opt = MmgOptions(),
                   const std::vector<int> *pinned_vertices = nullptr);
 
    void remesh_adaptive_2d(const Eigen::MatrixXd &V, const Eigen::MatrixXi &F, const Eigen::VectorXd &S,
                            Eigen::MatrixXd &OV, Eigen::MatrixXi &OF, MmgOptions opt = MmgOptions());
 
    bool remesh_3d(const Eigen::MatrixXd &V, const Eigen::MatrixXi &T,
                   Eigen::MatrixXd &OV, Eigen::MatrixXi &OF, Eigen::MatrixXi &OT,
                   MmgOptions opt = MmgOptions(),
                   const std::vector<int> *pinned_vertices = nullptr);
 
    void remesh_adaptive_3d(const Eigen::MatrixXd &V, const Eigen::MatrixXi &T, const Eigen::VectorXd &S,
                            Eigen::MatrixXd &OV, Eigen::MatrixXi &OF, Eigen::MatrixXi &OT, MmgOptions opt = MmgOptions());
 
} // namespace polyfem::mesh
 
#endif