include/deal.II/fe/fe_q.h

00001 //---------------------------------------------------------------------------
00002 //    @f$Id: fe_q.h 25499 2012-05-06 11:32:37Z kormann @f$
00003 //
00004 //    Copyright (C) 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008, 2009, 2010, 2011, 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_q_h
00013 #define __deal2__fe_q_h
00014 
00015 #include <deal.II/base/config.h>
00016 #include <deal.II/base/tensor_product_polynomials.h>
00017 #include <deal.II/fe/fe_poly.h>
00018 
00019 DEAL_II_NAMESPACE_OPEN
00020 
00021 
00024 
00222 template <int dim, int spacedim=dim>
00223 class FE_Q : public FE_Poly<TensorProductPolynomials<dim>,dim,spacedim>
00224 {
00225   public:
00230     FE_Q (const unsigned int p);
00231 
00252     FE_Q (const Quadrature<1> &points);
00253 
00263     virtual std::string get_name () const;
00264 
00281     virtual void
00282     get_interpolation_matrix (const FiniteElement<dim,spacedim> &source,
00283                               FullMatrix<double>       &matrix) const;
00284 
00285 
00307     virtual void
00308     get_face_interpolation_matrix (const FiniteElement<dim,spacedim> &source,
00309                                    FullMatrix<double>       &matrix) const;
00310 
00332     virtual void
00333     get_subface_interpolation_matrix (const FiniteElement<dim,spacedim> &source,
00334                                       const unsigned int        subface,
00335                                       FullMatrix<double>       &matrix) const;
00336 
00350     virtual bool has_support_on_face (const unsigned int shape_index,
00351                                       const unsigned int face_index) const;
00352 
00371     virtual bool hp_constraints_are_implemented () const;
00372 
00407     virtual
00408     std::vector<std::pair<unsigned int, unsigned int> >
00409     hp_vertex_dof_identities (const FiniteElement<dim,spacedim> &fe_other) const;
00410 
00418     virtual
00419     std::vector<std::pair<unsigned int, unsigned int> >
00420     hp_line_dof_identities (const FiniteElement<dim,spacedim> &fe_other) const;
00421 
00429     virtual
00430     std::vector<std::pair<unsigned int, unsigned int> >
00431     hp_quad_dof_identities (const FiniteElement<dim,spacedim> &fe_other) const;
00432 
00445     virtual
00446     FiniteElementDomination::Domination
00447     compare_for_face_domination (const FiniteElement<dim,spacedim> &fe_other) const;
00449 
00461     virtual std::size_t memory_consumption () const;
00462 
00463   protected:
00471     virtual FiniteElement<dim,spacedim> * clone() const;
00472 
00473   private:
00474 
00485     static std::vector<unsigned int> get_dpo_vector(const unsigned int degree);
00486 
00494     static
00495     std::vector<unsigned int>
00496     face_lexicographic_to_hierarchic_numbering (const unsigned int degree);
00497 
00503     void initialize_constraints ();
00504 
00510     void initialize_embedding ();
00511 
00517     void initialize_restriction ();
00518 
00526     void initialize_unit_support_points ();
00527 
00535     void initialize_unit_support_points (const Quadrature<1> &points);
00536 
00544     void initialize_unit_face_support_points ();
00545 
00553     void initialize_unit_face_support_points (const Quadrature<1> &points);
00554 
00562     void initialize_quad_dof_index_permutation ();
00563 
00571     const std::vector<unsigned int> face_index_map;
00572 
00573 
00582     struct Implementation;
00583 
00593     template <int, int> friend class FE_Q;
00594 
00595     friend struct FE_Q<dim,spacedim>::Implementation;
00596 };
00597 
00598 
00599 
00602 /* -------------- declaration of explicit specializations ------------- */
00603 
00604 template <>
00605 void FE_Q<1>::initialize_unit_face_support_points ();
00606 
00607 template <>
00608 std::vector<unsigned int>
00609 FE_Q<1>::face_lexicographic_to_hierarchic_numbering (const unsigned int);
00610 
00611 
00612 template <>
00613 void FE_Q<1,2>::initialize_unit_face_support_points ();
00614 
00615 template <>
00616 std::vector<unsigned int>
00617 FE_Q<1,2>::face_lexicographic_to_hierarchic_numbering (const unsigned int);
00618 
00619 template <>
00620 void FE_Q<1,3>::initialize_unit_face_support_points ();
00621 
00622 template <>
00623 std::vector<unsigned int>
00624 FE_Q<1,3>::face_lexicographic_to_hierarchic_numbering (const unsigned int);
00625 
00626 DEAL_II_NAMESPACE_CLOSE
00627 
00628 #endif
00629 
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines

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