LagrangeBasis3d.hpp

Go to the documentation of this file.

#pragma once
 
#include <polyfem/basis/ElementBases.hpp>
#include <polyfem/mesh/mesh3D/CMesh3D.hpp>
#include <polyfem/mesh/mesh3D/NCMesh3D.hpp>
#include <polyfem/mesh/LocalBoundary.hpp>
#include <polyfem/basis/InterfaceData.hpp>
#include <polyfem/mesh/mesh3D/Navigation3D.hpp>
#include <polyfem/mesh/MeshNodes.hpp>
 
#include <Eigen/Dense>
#include <vector>
 
namespace polyfem
{
    namespace basis
    {
        class LagrangeBasis3d
        {
        public:
            static int build_bases(
                const mesh::Mesh3D &mesh,
                const std::string &assembler,
                const int quadrature_order,
                const int mass_quadrature_order,
                const int discr_orderp,
                const int discr_orderq,
                const bool bernstein,
                const bool serendipity,
                const bool has_polys,
                const bool is_geom_bases,
                const bool use_corner_quadrature,
                std::vector<ElementBases> &bases,
                std::vector<mesh::LocalBoundary> &local_boundary,
                std::map<int, InterfaceData> &poly_face_to_data,
                std::shared_ptr<mesh::MeshNodes> &mesh_nodes);
 
            static int build_bases(
                const mesh::Mesh3D &mesh,
                const std::string &assembler,
                const int quadrature_order,
                const int mass_quadrature_order,
                const Eigen::VectorXi &discr_orderp,
                const Eigen::VectorXi &discr_orderq,
                const bool bernstein,
                const bool serendipity,
                const bool has_polys,
                const bool is_geom_bases,
                const bool use_corner_quadrature,
                std::vector<ElementBases> &bases,
                std::vector<mesh::LocalBoundary> &local_boundary,
                std::map<int, InterfaceData> &poly_face_to_data,
                std::shared_ptr<mesh::MeshNodes> &mesh_nodes);
 
            // return the local faces nodes for a tet or a hex of order p, index points to a face
            static Eigen::VectorXi tet_face_local_nodes(const int p, const mesh::Mesh3D &mesh, mesh::Navigation3D::Index index);
            static Eigen::VectorXi hex_face_local_nodes(const bool serendipity, const int q, const mesh::Mesh3D &mesh, mesh::Navigation3D::Index index);
            static Eigen::VectorXi prism_face_local_nodes(const int p, const int q, const mesh::Mesh3D &mesh, mesh::Navigation3D::Index index);
 
        private:
            static Eigen::MatrixXd linear_hex_face_local_nodes_coordinates(const mesh::Mesh3D &mesh, mesh::Navigation3D::Index index);
 
            static Eigen::RowVector3d quadr_hex_local_node_coordinates(int local_index);
        };
    } // namespace basis
} // namespace polyfem