00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012 #ifndef __deal2__fe_dgp_h
00013 #define __deal2__fe_dgp_h
00014
00015 #include <deal.II/base/config.h>
00016 #include <deal.II/base/polynomial_space.h>
00017 #include <deal.II/fe/fe_poly.h>
00018
00019 DEAL_II_NAMESPACE_OPEN
00020
00021 template <int dim, int spacedim> class MappingQ;
00022
00025
00107 template <int dim, int spacedim=dim>
00108 class FE_DGP : public FE_Poly<PolynomialSpace<dim>,dim,spacedim>
00109 {
00110 public:
00115 FE_DGP (const unsigned int p);
00116
00126 virtual std::string get_name () const;
00127
00171 virtual
00172 std::vector<std::pair<unsigned int, unsigned int> >
00173 hp_vertex_dof_identities (const FiniteElement<dim,spacedim> &fe_other) const;
00174
00186 virtual
00187 std::vector<std::pair<unsigned int, unsigned int> >
00188 hp_line_dof_identities (const FiniteElement<dim,spacedim> &fe_other) const;
00189
00201 virtual
00202 std::vector<std::pair<unsigned int, unsigned int> >
00203 hp_quad_dof_identities (const FiniteElement<dim,spacedim> &fe_other) const;
00204
00217 virtual bool hp_constraints_are_implemented () const;
00218
00231 virtual
00232 FiniteElementDomination::Domination
00233 compare_for_face_domination (const FiniteElement<dim,spacedim> &fe_other) const;
00234
00260 virtual void
00261 get_face_interpolation_matrix (const FiniteElement<dim,spacedim> &source,
00262 FullMatrix<double> &matrix) const;
00263
00285 virtual void
00286 get_subface_interpolation_matrix (const FiniteElement<dim,spacedim> &source,
00287 const unsigned int subface,
00288 FullMatrix<double> &matrix) const;
00289
00303 virtual bool has_support_on_face (const unsigned int shape_index,
00304 const unsigned int face_index) const;
00305
00317 virtual std::size_t memory_consumption () const;
00318
00319
00330 struct Matrices
00331 {
00336 static const double * const projection_matrices[][GeometryInfo<dim>::max_children_per_cell];
00337
00344 static const unsigned int n_projection_matrices;
00345 };
00346
00347 protected:
00348
00356 virtual FiniteElement<dim,spacedim> *clone() const;
00357
00358 private:
00359
00370 static std::vector<unsigned int> get_dpo_vector (const unsigned int degree);
00371 };
00372
00373
00374 #ifndef DOXYGEN
00375
00376
00377
00378
00379 #ifndef DEAL_II_MEMBER_VAR_SPECIALIZATION_BUG
00380 template <>
00381 const double * const FE_DGP<1>::Matrices::projection_matrices[][GeometryInfo<1>::max_children_per_cell];
00382
00383 template <>
00384 const unsigned int FE_DGP<1>::Matrices::n_projection_matrices;
00385
00386 template <>
00387 const double * const FE_DGP<2>::Matrices::projection_matrices[][GeometryInfo<2>::max_children_per_cell];
00388
00389 template <>
00390 const unsigned int FE_DGP<2>::Matrices::n_projection_matrices;
00391
00392 template <>
00393 const double * const FE_DGP<3>::Matrices::projection_matrices[][GeometryInfo<3>::max_children_per_cell];
00394
00395 template <>
00396 const unsigned int FE_DGP<3>::Matrices::n_projection_matrices;
00397
00398
00399 template <>
00400 const double * const FE_DGP<1,2>::Matrices::projection_matrices[][GeometryInfo<1>::max_children_per_cell];
00401
00402 template <>
00403 const unsigned int FE_DGP<1,2>::Matrices::n_projection_matrices;
00404
00405 template <>
00406 const double * const FE_DGP<2,3>::Matrices::projection_matrices[][GeometryInfo<2>::max_children_per_cell];
00407
00408 template <>
00409 const unsigned int FE_DGP<2,3>::Matrices::n_projection_matrices;
00410
00411 #endif
00412
00413 #endif // DOXYGEN
00414
00415 DEAL_II_NAMESPACE_CLOSE
00416
00417 #endif
00418