17 #include <deal.II/base/std_cxx14/memory.h> 19 #include <deal.II/fe/fe_dgp.h> 20 #include <deal.II/fe/fe_tools.h> 25 DEAL_II_NAMESPACE_OPEN
27 template <
int dim,
int spacedim>
31 Polynomials::Legendre::generate_complete_basis(degree)),
41 std::vector<bool>(1, true)))
56 template <
int dim,
int spacedim>
64 std::ostringstream namebuf;
73 template <
int dim,
int spacedim>
74 std::unique_ptr<FiniteElement<dim, spacedim>>
77 return std_cxx14::make_unique<FE_DGP<dim, spacedim>>(*this);
87 template <
int dim,
int spacedim>
88 std::vector<unsigned int>
91 std::vector<unsigned int> dpo(dim + 1, 0U);
93 for (
unsigned int i = 1; i < dim; ++i)
95 dpo[dim] *= deg + 1 + i;
103 template <
int dim,
int spacedim>
113 (void)interpolation_matrix;
117 (dynamic_cast<const FEDGP *>(&x_source_fe) !=
nullptr),
118 typename FE::ExcInterpolationNotImplemented());
120 Assert(interpolation_matrix.
m() == 0,
122 Assert(interpolation_matrix.
n() == 0,
128 template <
int dim,
int spacedim>
139 (void)interpolation_matrix;
143 (dynamic_cast<const FEDGP *>(&x_source_fe) !=
nullptr),
144 typename FE::ExcInterpolationNotImplemented());
146 Assert(interpolation_matrix.
m() == 0,
148 Assert(interpolation_matrix.
n() == 0,
154 template <
int dim,
int spacedim>
163 template <
int dim,
int spacedim>
164 std::vector<std::pair<unsigned int, unsigned int>>
170 return std::vector<std::pair<unsigned int, unsigned int>>();
174 return std::vector<std::pair<unsigned int, unsigned int>>();
180 template <
int dim,
int spacedim>
181 std::vector<std::pair<unsigned int, unsigned int>>
187 return std::vector<std::pair<unsigned int, unsigned int>>();
191 return std::vector<std::pair<unsigned int, unsigned int>>();
197 template <
int dim,
int spacedim>
198 std::vector<std::pair<unsigned int, unsigned int>>
204 return std::vector<std::pair<unsigned int, unsigned int>>();
208 return std::vector<std::pair<unsigned int, unsigned int>>();
214 template <
int dim,
int spacedim>
230 template <
int dim,
int spacedim>
233 const unsigned int)
const 241 template <
int dim,
int spacedim>
242 std::pair<Table<2, bool>, std::vector<unsigned int>>
246 constant_modes(0, 0) =
true;
247 return std::pair<Table<2, bool>, std::vector<unsigned int>>(
248 constant_modes, std::vector<unsigned int>(1, 0));
253 template <
int dim,
int spacedim>
264 #include "fe_dgp.inst" 267 DEAL_II_NAMESPACE_CLOSE
virtual std::size_t memory_consumption() const override
std::vector< std::vector< FullMatrix< double > > > restriction
virtual std::vector< std::pair< unsigned int, unsigned int > > hp_vertex_dof_identities(const FiniteElement< dim, spacedim > &fe_other) const override
const unsigned int degree
FE_DGP(const unsigned int p)
#define AssertThrow(cond, exc)
virtual std::string get_name() const override
std::vector< std::vector< FullMatrix< double > > > prolongation
#define Assert(cond, exc)
static::ExceptionBase & ExcDimensionMismatch(std::size_t arg1, std::size_t arg2)
static std::vector< unsigned int > get_dpo_vector(const unsigned int degree)
virtual FiniteElementDomination::Domination compare_for_face_domination(const FiniteElement< dim, spacedim > &fe_other) const override
virtual std::string get_name() const =0
virtual bool has_support_on_face(const unsigned int shape_index, const unsigned int face_index) const override
std::string dim_string(const int dim, const int spacedim)
const unsigned int dofs_per_cell
virtual std::pair< Table< 2, bool >, std::vector< unsigned int > > get_constant_modes() const override
void reinit_restriction_and_prolongation_matrices(const bool isotropic_restriction_only=false, const bool isotropic_prolongation_only=false)
static::ExceptionBase & ExcNotImplemented()
virtual void get_subface_interpolation_matrix(const FiniteElement< dim, spacedim > &source, const unsigned int subface, FullMatrix< double > &matrix) const override
virtual void get_face_interpolation_matrix(const FiniteElement< dim, spacedim > &source, FullMatrix< double > &matrix) const override
virtual bool hp_constraints_are_implemented() const override
virtual std::vector< std::pair< unsigned int, unsigned int > > hp_line_dof_identities(const FiniteElement< dim, spacedim > &fe_other) const override
virtual std::vector< std::pair< unsigned int, unsigned int > > hp_quad_dof_identities(const FiniteElement< dim, spacedim > &fe_other) const override
virtual std::unique_ptr< FiniteElement< dim, spacedim > > clone() const override