include/deal.II/fe/fe_dgp.h

00001 //---------------------------------------------------------------------------
00002 //    @f$Id: fe_dgp.h 25345 2012-03-31 08:37:04Z bangerth @f$
00003 //
00004 //    Copyright (C) 2002, 2003, 2004, 2005, 2006, 2010, 2012 by the deal.II authors
00005 //
00006 //    This file is subject to QPL and may not be  distributed
00007 //    without copyright and license information. Please refer
00008 //    to the file deal.II/doc/license.html for the  text  and
00009 //    further information on this license.
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 // declaration of explicit specializations of member variables, if the
00378 // compiler allows us to do that (the standard says we must)
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 //codimension 1
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 
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines

deal.II documentation generated on Tue May 22 2012 12:06:07 by doxygen 1.7.3