43 dpo.dofs_per_object_exclusive = {{1}, {0}, {0, 0, 0, 0, 0}, {0}};
44 dpo.dofs_per_object_inclusive = {{1}, {2}, {3, 3, 4, 4, 4}, {6}};
45 dpo.object_index = {{}, {6}, {6}, {6}};
46 dpo.first_object_index_on_face = {{}, {3, 3, 4, 4, 4}, {3, 3, 4, 4, 4}};
50 dpo.dofs_per_object_exclusive = {{1}, {1}, {0, 0, 1, 1, 1}, {0}};
51 dpo.dofs_per_object_inclusive = {{1}, {3}, {6, 6, 9, 9, 9}, {18}};
52 dpo.object_index = {{}, {6}, {15, 15, 15, 16, 17}, {18}};
53 dpo.first_object_index_on_face = {{}, {3, 3, 4, 4, 4}, {6, 6, 8, 8, 8}};
69 unsigned int n_dofs = 0;
84template <
int dim,
int spacedim>
86 const unsigned int degree,
119 fe_triangle.get_unit_support_points()[
pair[0]][0],
129 for (
const auto &p :
fe_triangle.get_unit_support_points())
133 for (
const auto &p :
fe_quad.get_unit_support_points())
141template <
int dim,
int spacedim>
149 this->get_unit_support_points().size());
153 for (
unsigned int i = 0; i < this->dofs_per_cell; ++i)
163template <
int dim,
int spacedim>
172template <
int dim,
int spacedim>
173std::unique_ptr<FiniteElement<dim, spacedim>>
176 return std::make_unique<FE_WedgeP<dim, spacedim>>(*this);
181template <
int dim,
int spacedim>
187 << this->degree <<
")";
194template <
int dim,
int spacedim>
198 const unsigned int codim)
const
250 if (fe_nothing->is_dominating())
265template <
int dim,
int spacedim>
266std::vector<std::pair<unsigned int, unsigned int>>
281template <
int dim,
int spacedim>
282std::vector<std::pair<unsigned int, unsigned int>>
292 std::vector<std::pair<unsigned int, unsigned int>>
result;
294 for (
unsigned int i = 0; i < this->degree - 1; ++i)
295 result.emplace_back(i, i);
302template <
int dim,
int spacedim>
303std::vector<std::pair<unsigned int, unsigned int>>
306 const unsigned int face_no)
const
323 std::vector<std::pair<unsigned int, unsigned int>>
result;
325 for (
unsigned int i = 0; i < this->n_dofs_per_quad(
face_no); ++i)
326 result.emplace_back(i, i);
333template <
int dim,
int spacedim>
342template <
int dim,
int spacedim>
343std::unique_ptr<FiniteElement<dim, spacedim>>
346 return std::make_unique<FE_WedgeDGP<dim, spacedim>>(*this);
351template <
int dim,
int spacedim>
357 << this->degree <<
")";
363#include "fe/fe_wedge_p.inst"
std::unique_ptr< FiniteElement< dim, spacedim > > clone() const override
std::string get_name() const override
FE_WedgeDGP(const unsigned int degree)
std::vector< std::pair< unsigned int, unsigned int > > hp_quad_dof_identities(const FiniteElement< dim, spacedim > &fe_other, const unsigned int face_no=0) const override
FE_WedgeP(const unsigned int degree)
std::vector< std::pair< unsigned int, unsigned int > > hp_vertex_dof_identities(const FiniteElement< dim, spacedim > &fe_other) const override
std::vector< std::pair< unsigned int, unsigned int > > hp_line_dof_identities(const FiniteElement< dim, spacedim > &fe_other) const override
FiniteElementDomination::Domination compare_for_domination(const FiniteElement< dim, spacedim > &fe_other, const unsigned int codim) const override
std::unique_ptr< FiniteElement< dim, spacedim > > clone() const override
std::string get_name() const override
FE_WedgePoly(const unsigned int degree, const internal::GenericDoFsPerObject &dpos, const typename FiniteElementData< dim >::Conformity conformity)
virtual void convert_generalized_support_point_values_to_dof_values(const std::vector< Vector< double > > &support_point_values, std::vector< double > &nodal_values) const override
const unsigned int degree
unsigned int n_dofs_per_cell() const
ReferenceCell reference_cell() const
std::vector< std::vector< Point< dim - 1 > > > unit_face_support_points
std::vector< Point< dim > > unit_support_points
unsigned int n_faces() const
ReferenceCell face_reference_cell(const unsigned int face_no) const
std_cxx20::ranges::iota_view< unsigned int, unsigned int > face_indices() const
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
#define DEAL_II_ASSERT_UNREACHABLE()
#define DEAL_II_NOT_IMPLEMENTED()
static ::ExceptionBase & ExcNotImplemented()
#define Assert(cond, exc)
static ::ExceptionBase & ExcImpossibleInDim(int arg1)
#define AssertDimension(dim1, dim2)
#define AssertIndexRange(index, range)
@ either_element_can_dominate
@ other_element_dominates
@ neither_element_dominates
constexpr const ReferenceCell Quadrilateral
constexpr const ReferenceCell Wedge
constexpr const ReferenceCell Triangle
std::string dim_string(const int dim, const int spacedim)
constexpr ::ndarray< unsigned int, 6, 2 > wedge_table_1
constexpr ::ndarray< unsigned int, 18, 2 > wedge_table_2
internal::GenericDoFsPerObject expand(const unsigned int dim, const std::vector< unsigned int > &dofs_per_object, const ReferenceCell reference_cell)