PolygonalBasis3d.hpp

Go to the documentation of this file.

#pragma once
 
#include <polyfem/mesh/mesh3D/Mesh3D.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 PolygonalBasis3d
        {
        public:
            // Compute the integral constraints for each basis of the mesh. This step is PDE-dependent.
            //
            // @param[in]  mesh             Input volume 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::Assembler &assembler,
                const mesh::Mesh3D &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::Mesh3D &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_face_to_data,
                std::map<int, std::pair<Eigen::MatrixXd, Eigen::MatrixXi>> &mapped_boundary);
        };
    } // namespace basis
} // namespace polyfem