Reference documentation for deal.II version 8.1.0
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Modules Pages
fe_q_base.h
1 // ---------------------------------------------------------------------
2 // @f$Id: fe_q_base.h 30276 2013-08-10 10:52:57Z bangerth @f$
3 //
4 // Copyright (C) 2000 - 2013 by the deal.II authors
5 //
6 // This file is part of the deal.II library.
7 //
8 // The deal.II library is free software; you can use it, redistribute
9 // it, and/or modify it under the terms of the GNU Lesser General
10 // Public License as published by the Free Software Foundation; either
11 // version 2.1 of the License, or (at your option) any later version.
12 // The full text of the license can be found in the file LICENSE at
13 // the top level of the deal.II distribution.
14 //
15 // ---------------------------------------------------------------------
16 
17 #ifndef __deal2__fe_q_base_h
18 #define __deal2__fe_q_base_h
19 
20 #include <deal.II/base/config.h>
21 #include <deal.II/fe/fe_poly.h>
22 #include <deal.II/base/thread_management.h>
23 
25 
26 
29 
37 template <class POLY, int dim=POLY::dimension, int spacedim=dim>
38 class FE_Q_Base : public FE_Poly<POLY,dim,spacedim>
39 {
40 public:
44  FE_Q_Base (const POLY &poly_space,
45  const FiniteElementData<dim> &fe_data,
46  const std::vector<bool> &restriction_is_additive_flags);
47 
57  virtual void
59  FullMatrix<double> &matrix) const;
60 
61 
71  virtual void
73  FullMatrix<double> &matrix) const;
74 
84  virtual void
86  const unsigned int subface,
87  FullMatrix<double> &matrix) const;
88 
97  virtual bool has_support_on_face (const unsigned int shape_index,
98  const unsigned int face_index) const;
99 
122  virtual const FullMatrix<double> &
123  get_restriction_matrix (const unsigned int child,
125 
152  virtual const FullMatrix<double> &
153  get_prolongation_matrix (const unsigned int child,
155 
191  virtual
192  unsigned int face_to_cell_index (const unsigned int face_dof_index,
193  const unsigned int face,
194  const bool face_orientation = true,
195  const bool face_flip = false,
196  const bool face_rotation = false) const;
197 
211  virtual bool hp_constraints_are_implemented () const;
212 
230  virtual
231  std::vector<std::pair<unsigned int, unsigned int> >
233 
238  virtual
239  std::vector<std::pair<unsigned int, unsigned int> >
240  hp_line_dof_identities (const FiniteElement<dim,spacedim> &fe_other) const;
241 
246  virtual
247  std::vector<std::pair<unsigned int, unsigned int> >
248  hp_quad_dof_identities (const FiniteElement<dim,spacedim> &fe_other) const;
249 
258  virtual
262 
263 protected:
270  static std::vector<unsigned int> get_dpo_vector(const unsigned int degree);
271 
277  void initialize (const std::vector<Point<1> > &support_points_1d);
278 
283  void initialize_constraints (const std::vector<Point<1> > &points);
284 
289  void initialize_unit_support_points (const std::vector<Point<1> > &points);
290 
295  void initialize_unit_face_support_points (const std::vector<Point<1> > &points);
296 
302 
309  struct Implementation;
310 
311  /*
312  * Declare implementation friend.
313  */
314  friend struct FE_Q_Base<POLY,dim,spacedim>::Implementation;
315 
316 private:
317  /*
318  * Mutex for protecting initialization of restriction and embedding matrix.
319  */
320  mutable Threads::Mutex mutex;
321 };
322 
323 
326 DEAL_II_NAMESPACE_CLOSE
327 
328 #endif
virtual bool has_support_on_face(const unsigned int shape_index, const unsigned int face_index) const
virtual std::vector< std::pair< unsigned int, unsigned int > > hp_line_dof_identities(const FiniteElement< dim, spacedim > &fe_other) const
virtual unsigned int face_to_cell_index(const unsigned int face_dof_index, const unsigned int face, const bool face_orientation=true, const bool face_flip=false, const bool face_rotation=false) const
const std::vector< bool > restriction_is_additive_flags
Definition: fe.h:1909
virtual const FullMatrix< double > & get_prolongation_matrix(const unsigned int child, const RefinementCase< dim > &refinement_case=RefinementCase< dim >::isotropic_refinement) const
virtual std::vector< std::pair< unsigned int, unsigned int > > hp_quad_dof_identities(const FiniteElement< dim, spacedim > &fe_other) const
const unsigned int degree
Definition: fe_base.h:287
virtual const FullMatrix< double > & get_restriction_matrix(const unsigned int child, const RefinementCase< dim > &refinement_case=RefinementCase< dim >::isotropic_refinement) const
virtual bool hp_constraints_are_implemented() const
virtual void get_subface_interpolation_matrix(const FiniteElement< dim, spacedim > &source, const unsigned int subface, FullMatrix< double > &matrix) const
void initialize_constraints(const std::vector< Point< 1 > > &points)
virtual FiniteElementDomination::Domination compare_for_face_domination(const FiniteElement< dim, spacedim > &fe_other) const
FE_Q_Base(const POLY &poly_space, const FiniteElementData< dim > &fe_data, const std::vector< bool > &restriction_is_additive_flags)
static std::vector< unsigned int > get_dpo_vector(const unsigned int degree)
void initialize_quad_dof_index_permutation()
virtual void get_interpolation_matrix(const FiniteElement< dim, spacedim > &source, FullMatrix< double > &matrix) const
void initialize(const std::vector< Point< 1 > > &support_points_1d)
void initialize_unit_face_support_points(const std::vector< Point< 1 > > &points)
void initialize_unit_support_points(const std::vector< Point< 1 > > &points)
virtual std::vector< std::pair< unsigned int, unsigned int > > hp_vertex_dof_identities(const FiniteElement< dim, spacedim > &fe_other) const
virtual void get_face_interpolation_matrix(const FiniteElement< dim, spacedim > &source, FullMatrix< double > &matrix) const
POLY poly_space
Definition: fe_poly.h:284