UpsampleMesh.hpp

Go to the documentation of this file.

#pragma once
 
#include <polyfem/mesh/Mesh.hpp>
#include <polyfem/mesh/LocalBoundary.hpp>
#include <polyfem/basis/ElementBases.hpp>
#include <polyfem/utils/Types.hpp>
 
#include <Eigen/Core>
 
namespace polyfem::mesh
{
    void stitch_mesh(
        const Eigen::MatrixXd &V,
        const Eigen::MatrixXi &F,
        Eigen::MatrixXd &V_out,
        Eigen::MatrixXi &F_out,
        const double epsilon = 1e-5);
 
    void stitch_mesh(
        const Eigen::MatrixXd &V,
        const Eigen::MatrixXi &F,
        const std::vector<Eigen::Triplet<double>> &W,
        Eigen::MatrixXd &V_out,
        Eigen::MatrixXi &F_out,
        std::vector<Eigen::Triplet<double>> &W_out,
        const double epsilon = 1e-5);
 
    double max_edge_length(const Eigen::MatrixXd &V, const Eigen::MatrixXi &F);
 
    // Regular tessellation
 
    void regular_grid_triangle_barycentric_coordinates(
        const int n, Eigen::MatrixXd &V, Eigen::MatrixXi &F);
 
    void regular_grid_tessellation(
        const Eigen::MatrixXd &V,
        const Eigen::MatrixXi &F,
        const double max_edge_length,
        Eigen::MatrixXd &V_out,
        Eigen::MatrixXi &F_out);
 
    // Irregular tessellation
 
    Eigen::MatrixXd
    refine_edge(const VectorNd &a, const VectorNd &b, const double max_edge_length);
 
    void refine_triangle_edges(
        const VectorNd &a,
        const VectorNd &b,
        const VectorNd &c,
        const double max_edge_len,
        Eigen::MatrixXd &V,
        Eigen::MatrixXi &E);
 
    void irregular_triangle_barycentric_coordinates(
        const Eigen::Vector3d &a,
        const Eigen::Vector3d &b,
        const Eigen::Vector3d &c,
        const double max_edge_length,
        Eigen::MatrixXd &UV,
        Eigen::MatrixXi &F);
 
    void irregular_triangle(
        const Eigen::Vector3d &a,
        const Eigen::Vector3d &b,
        const Eigen::Vector3d &c,
        const double max_edge_length,
        Eigen::MatrixXd &V,
        Eigen::MatrixXi &F);
 
    void irregular_tessellation(
        const Eigen::MatrixXd &V,
        const Eigen::MatrixXi &F,
        const double max_edge_length,
        Eigen::MatrixXd &V_out,
        Eigen::MatrixXi &F_out);
} // namespace polyfem::mesh