PolygonalBasis2d.hpp

Go to the documentation of this file.

#pragma once
 
#include <polyfem/mesh/mesh2D/CMesh2D.hpp>
#include <polyfem/mesh/mesh2D/NCMesh2D.hpp>
#include <polyfem/assembler/Assembler.hpp>
#include <polyfem/basis/ElementBases.hpp>
#include <polyfem/assembler/ElementAssemblyValues.hpp>
#include <polyfem/basis/InterfaceData.hpp>
 
#include <Eigen/Dense>
#include <vector>
#include <map>
 
namespace polyfem
{
    namespace basis
    {
        class PolygonalBasis2d
        {
        public:
            // Compute the integral constraints for each basis of the mesh. This step is PDE-dependent.
            //
            // @param[in]  mesh             Input surface mesh
            // @param[in]  n_bases          Number of bases/nodes in the mesh
            // @param[in]  values           Shape functions for the PDE field evaluated over each
            //                              elements
            // @param[in]  gvalues          Shape functions for the geometric mapping evaluated over
            //                              each elements
            // @param[out] basis_integrals  Integral constraints over polygon that each basis of the
            //                              mesh must verify
            //
            static void compute_integral_constraints(
                const assembler::LinearAssembler &assembler,
                const mesh::Mesh2D &mesh,
                const int n_bases,
                const std::vector<ElementBases> &bases,
                const std::vector<ElementBases> &gbases,
                Eigen::MatrixXd &basis_integrals);
 
            static int build_bases(
                const assembler::LinearAssembler &assembler,
                const int n_samples_per_edge,
                const mesh::Mesh2D &mesh,
                const int n_bases,
                const int quadrature_order,
                const int mass_quadrature_order,
                const int integral_constraints,
                std::vector<ElementBases> &bases,
                const std::vector<ElementBases> &gbases,
                const std::map<int, InterfaceData> &poly_edge_to_data,
                std::map<int, Eigen::MatrixXd> &mapped_boundary);
        };
    } // namespace basis
} // namespace polyfem