00001
00002
00003
00004
00005
00006
00007
00008
00009
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
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