40template <
int dim,
int spacedim>
47template <
int dim,
int spacedim>
56template <
int dim,
int spacedim>
59 const std::vector<bool> & r_i_a_f,
60 const std::vector<ComponentMask> &nonzero_c)
62 , adjust_line_dof_index_for_line_orientation_table(
63 dim == 3 ? this->n_dofs_per_line() : 0)
64 , system_to_base_table(this->n_dofs_per_cell())
65 , component_to_base_table(this->components,
66 std::make_pair(
std::make_pair(0
U, 0
U), 0
U))
71 restriction_is_additive_flags(
73 std::vector<
bool>(fe_data.n_dofs_per_cell(), r_i_a_f[0]) :
76 nonzero_c.size() == 1 ?
79 , n_nonzero_components_table(compute_n_nonzero_components(nonzero_components))
80 , cached_primitivity(
std::find_if(n_nonzero_components_table.
begin(),
81 n_nonzero_components_table.
end(),
82 [](const unsigned
int n_components) {
83 return n_components != 1U;
84 }) == n_nonzero_components_table.end())
86 Assert(restriction_is_additive_flags.size() == this->n_dofs_per_cell(),
88 this->n_dofs_per_cell()));
90 for (
unsigned int i = 0; i < nonzero_components.size(); ++i)
92 Assert(nonzero_components[i].size() == this->n_components(),
94 Assert(nonzero_components[i].n_selected_components() >= 1,
97 Assert(n_nonzero_components_table[i] <= this->n_components(),
104 if (this->is_primitive())
106 system_to_component_table.resize(this->n_dofs_per_cell());
107 for (
unsigned int j = 0; j < this->n_dofs_per_cell(); ++j)
108 system_to_component_table[j] = std::pair<unsigned, unsigned>(0, j);
110 face_system_to_component_table.resize(this->n_unique_faces());
111 for (
unsigned int f = 0; f < this->n_unique_faces(); ++f)
113 face_system_to_component_table[f].resize(this->n_dofs_per_face(f));
114 for (
unsigned int j = 0; j < this->n_dofs_per_face(f); ++j)
115 face_system_to_component_table[f][j] =
116 std::pair<unsigned, unsigned>(0, j);
120 for (
unsigned int j = 0; j < this->n_dofs_per_cell(); ++j)
121 system_to_base_table[j] = std::make_pair(std::make_pair(0
U, 0
U), j);
123 face_system_to_base_table.resize(this->n_unique_faces());
124 for (
unsigned int f = 0; f < this->n_unique_faces(); ++f)
126 face_system_to_base_table[f].resize(this->n_dofs_per_face(f));
127 for (
unsigned int j = 0; j < this->n_dofs_per_face(f); ++j)
128 face_system_to_base_table[f][j] =
129 std::make_pair(std::make_pair(0
U, 0
U), j);
133 base_to_block_indices.reinit(1, 1);
140 ref < RefinementCase<dim>::isotropic_refinement + 1;
154 adjust_quad_dof_index_for_face_orientation_table.resize(
155 this->n_unique_quads());
157 for (
unsigned int f = 0; f < this->n_unique_quads(); ++f)
159 adjust_quad_dof_index_for_face_orientation_table[f] =
165 adjust_quad_dof_index_for_face_orientation_table[f].fill(0);
169 unit_face_support_points.resize(this->n_unique_faces());
170 generalized_face_support_points.resize(this->n_unique_faces());
175template <
int dim,
int spacedim>
176std::pair<std::unique_ptr<FiniteElement<dim, spacedim>>,
unsigned int>
179 return {this->
clone(), multiplicity};
184template <
int dim,
int spacedim>
195template <
int dim,
int spacedim>
199 const unsigned int)
const
207template <
int dim,
int spacedim>
218template <
int dim,
int spacedim>
222 const unsigned int)
const
230template <
int dim,
int spacedim>
241template <
int dim,
int spacedim>
246 const unsigned int)
const
254template <
int dim,
int spacedim>
265template <
int dim,
int spacedim>
270 const unsigned int)
const
278template <
int dim,
int spacedim>
289template <
int dim,
int spacedim>
294 const unsigned int)
const
301template <
int dim,
int spacedim>
304 const bool isotropic_restriction_only,
305 const bool isotropic_prolongation_only)
308 ref_case <= RefinementCase<dim>::isotropic_refinement;
311 const unsigned int nc =
314 for (
unsigned int i = 0; i < nc; ++i)
318 (!isotropic_restriction_only ||
324 (!isotropic_prolongation_only ||
333template <
int dim,
int spacedim>
336 const unsigned int child,
343 "Restriction matrices are only available for refined cells!"));
356template <
int dim,
int spacedim>
359 const unsigned int child,
366 "Prolongation matrices are only available for refined cells!"));
382template <
int dim,
int spacedim>
385 const unsigned int index)
const
394template <
int dim,
int spacedim>
412template <
int dim,
int spacedim>
418 this->n_components());
434template <
int dim,
int spacedim>
441 this->n_components());
459template <
int dim,
int spacedim>
480template <
int dim,
int spacedim>
491template <
int dim,
int spacedim>
502template <
int dim,
int spacedim>
514template <
int dim,
int spacedim>
551 "The component mask argument given to this function "
552 "is not a mask where the individual components belonging "
553 "to one block of the finite element are either all "
554 "selected or not selected. You can't call this function "
555 "with a component mask that splits blocks."));
566template <
int dim,
int spacedim>
569 const unsigned int face,
570 const bool face_orientation,
571 const bool face_flip,
572 const bool face_rotation)
const
590 if ((face_orientation !=
true) || (face_flip !=
false) ||
591 (face_rotation !=
false))
594 "The function in this base class can not handle this case. "
595 "Rather, the derived class you are using must provide "
596 "an overloaded version but apparently hasn't done so. See "
597 "the documentation of this function for more information."));
607 const unsigned int dof_index_on_vertex =
618 dof_index_on_vertex);
625 const unsigned int index =
629 const unsigned int dof_index_on_line = index % this->
n_dofs_per_line();
646 const unsigned int index =
655template <
int dim,
int spacedim>
658 const unsigned int index,
659 const unsigned int face,
660 const bool face_orientation,
661 const bool face_flip,
662 const bool face_rotation)
const
693 index, 4 * face_orientation + 2 * face_flip + face_rotation);
698template <
int dim,
int spacedim>
701 const unsigned int index,
702 const bool line_orientation)
const
713 this->n_dofs_per_line(),
715 if (line_orientation)
723template <
int dim,
int spacedim>
728 ref_case < RefinementCase<dim>::isotropic_refinement + 1;
730 for (
unsigned int c = 0;
751template <
int dim,
int spacedim>
756 ref_case < RefinementCase<dim>::isotropic_refinement + 1;
758 for (
unsigned int c = 0;
779template <
int dim,
int spacedim>
786 for (
unsigned int c = 0;
807template <
int dim,
int spacedim>
814 for (
unsigned int c = 0;
835template <
int dim,
int spacedim>
843 const unsigned int face_no = 0;
854template <
int dim,
int spacedim>
863template <
int dim,
int spacedim>
871 const unsigned int face_no = 0;
876 ExcMessage(
"Constraints for this element are only implemented "
877 "for the case that faces are refined isotropically "
878 "(which is always the case in 2d, and in 3d requires "
879 "that the neighboring cell of a coarse cell presents "
880 "exactly four children on the common face)."));
883 ExcMessage(
"The finite element for which you try to obtain "
884 "hanging node constraints does not appear to "
897template <
int dim,
int spacedim>
904 const unsigned int face_no = 0;
925template <
int dim,
int spacedim>
941template <
int dim,
int spacedim>
946 const unsigned int)
const
958template <
int dim,
int spacedim>
964 const unsigned int)
const
976template <
int dim,
int spacedim>
977std::vector<std::pair<unsigned int, unsigned int>>
982 return std::vector<std::pair<unsigned int, unsigned int>>();
987template <
int dim,
int spacedim>
988std::vector<std::pair<unsigned int, unsigned int>>
993 return std::vector<std::pair<unsigned int, unsigned int>>();
998template <
int dim,
int spacedim>
999std::vector<std::pair<unsigned int, unsigned int>>
1002 const unsigned int)
const
1005 return std::vector<std::pair<unsigned int, unsigned int>>();
1010template <
int dim,
int spacedim>
1014 const unsigned int)
const
1022template <
int dim,
int spacedim>
1029 return ((
typeid(*
this) ==
typeid(f)) && (this->
get_name() == f.
get_name()) &&
1037template <
int dim,
int spacedim>
1042 return !(*
this == f);
1047template <
int dim,
int spacedim>
1048const std::vector<Point<dim>> &
1063template <
int dim,
int spacedim>
1072template <
int dim,
int spacedim>
1073const std::vector<Point<dim>> &
1085template <
int dim,
int spacedim>
1094template <
int dim,
int spacedim>
1106template <
int dim,
int spacedim>
1107const std::vector<
Point<dim - 1>> &
1109 const unsigned int face_no)
const
1125template <
int dim,
int spacedim>
1128 const unsigned int face_no)
const
1136template <
int dim,
int spacedim>
1139 const unsigned int index,
1140 const unsigned int face_no)
const
1152template <
int dim,
int spacedim>
1155 const unsigned int)
const
1162template <
int dim,
int spacedim>
1168 const unsigned int n_total_components = this->
n_components();
1169 Assert((n_total_components == mask.
size()) || (mask.
size() == 0),
1170 ExcMessage(
"The given ComponentMask has the wrong size."));
1172 const unsigned int n_selected =
1175 ExcMessage(
"You need at least one selected component."));
1177 const unsigned int first_selected =
1182 for (
unsigned int c = 0; c < n_total_components; ++c)
1183 Assert((c < first_selected && (!mask[c])) ||
1184 (c >= first_selected && c < first_selected + n_selected &&
1186 (c >= first_selected + n_selected && !mask[c]),
1187 ExcMessage(
"Error: the given ComponentMask is not contiguous!"));
1190 return get_sub_fe(first_selected, n_selected);
1195template <
int dim,
int spacedim>
1198 const unsigned int first_component,
1199 const unsigned int n_selected_components)
const
1205 "You can only select a whole FiniteElement, not a part of one."));
1207 (void)first_component;
1208 (void)n_selected_components;
1214template <
int dim,
int spacedim>
1215std::pair<Table<2, bool>, std::vector<unsigned int>>
1219 return std::pair<Table<2, bool>, std::vector<unsigned int>>(
1226template <
int dim,
int spacedim>
1231 std::vector<double> &)
const
1234 ExcMessage(
"The element for which you are calling the current "
1235 "function does not have generalized support points (see "
1236 "the glossary for a definition of generalized support "
1237 "points). Consequently, the current function can not "
1238 "be defined and is not implemented by the element."));
1244template <
int dim,
int spacedim>
1265template <
int dim,
int spacedim>
1266std::vector<unsigned int>
1268 const std::vector<ComponentMask> &nonzero_components)
1281template <
int dim,
int spacedim>
1282std::unique_ptr<typename FiniteElement<dim, spacedim>::InternalDataBase>
1300template <
int dim,
int spacedim>
1301std::unique_ptr<typename FiniteElement<dim, spacedim>::InternalDataBase>
1319template <
int dim,
int spacedim>
1323 const unsigned int face_no,
1327 const ::internal::FEValuesImplementation::MappingRelatedData<dim,
1349template <
int dim,
int spacedim>
1353 const unsigned int face_no,
1357 const ::internal::FEValuesImplementation::MappingRelatedData<dim,
1366 ExcMessage(
"Use of a deprecated interface, please implement "
1367 "fill_fe_face_values taking a hp::QCollection argument"));
1372 (void)mapping_internal;
1381template <
int dim,
int spacedim>
1382std::unique_ptr<typename FiniteElement<dim, spacedim>::InternalDataBase>
1400template <
int dim,
int spacedim>
unsigned int size() const
bool represents_the_all_selected_mask() const
unsigned int size() const
bool represents_the_all_selected_mask() const
unsigned int size() const
unsigned int n_selected_components(const unsigned int overall_number_of_components=numbers::invalid_unsigned_int) const
unsigned int first_selected_component(const unsigned int overall_number_of_components=numbers::invalid_unsigned_int) const
unsigned int get_first_line_index() const
unsigned int n_dofs_per_vertex() const
unsigned int n_dofs_per_cell() const
unsigned int n_dofs_per_line() const
unsigned int get_first_quad_index(const unsigned int quad_no=0) const
unsigned int n_dofs_per_face(unsigned int face_no=0, unsigned int child=0) const
unsigned int n_blocks() const
unsigned int n_components() const
unsigned int n_unique_quads() const
ReferenceCell reference_cell() const
unsigned int get_first_face_line_index(const unsigned int face_no=0) const
unsigned int get_first_face_quad_index(const unsigned int face_no=0) const
unsigned int n_unique_faces() const
unsigned int n_dofs_per_quad(unsigned int face_no=0) const
virtual std::size_t memory_consumption() const
bool constraints_are_implemented(const ::internal::SubfaceCase< dim > &subface_case=::internal::SubfaceCase< dim >::case_isotropic) const
bool isotropic_prolongation_is_implemented() const
virtual std::string get_name() const =0
virtual Point< dim > unit_support_point(const unsigned int index) const
virtual FiniteElementDomination::Domination compare_for_domination(const FiniteElement< dim, spacedim > &fe_other, const unsigned int codim=0) const
const std::vector< unsigned int > n_nonzero_components_table
virtual std::unique_ptr< InternalDataBase > get_data(const UpdateFlags update_flags, const Mapping< dim, spacedim > &mapping, const Quadrature< dim > &quadrature, ::internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim > &output_data) const =0
virtual const FullMatrix< double > & get_restriction_matrix(const unsigned int child, const RefinementCase< dim > &refinement_case=RefinementCase< dim >::isotropic_refinement) const
std::vector< std::vector< Point< dim - 1 > > > unit_face_support_points
bool prolongation_is_implemented() const
virtual std::pair< Table< 2, bool >, std::vector< unsigned int > > get_constant_modes() const
virtual const FiniteElement< dim, spacedim > & base_element(const unsigned int index) const
virtual Tensor< 1, dim > shape_grad(const unsigned int i, const Point< dim > &p) const
virtual std::unique_ptr< InternalDataBase > get_face_data(const UpdateFlags update_flags, const Mapping< dim, spacedim > &mapping, const hp::QCollection< dim - 1 > &quadrature, ::internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim > &output_data) const
virtual Point< dim - 1 > unit_face_support_point(const unsigned int index, const unsigned int face_no=0) const
static std::vector< unsigned int > compute_n_nonzero_components(const std::vector< ComponentMask > &nonzero_components)
virtual std::unique_ptr< InternalDataBase > get_subface_data(const UpdateFlags update_flags, const Mapping< dim, spacedim > &mapping, const Quadrature< dim - 1 > &quadrature, ::internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim > &output_data) const
bool has_face_support_points(const unsigned int face_no=0) const
virtual std::vector< std::pair< unsigned int, unsigned int > > hp_quad_dof_identities(const FiniteElement< dim, spacedim > &fe_other, const unsigned int face_no=0) const
ComponentMask component_mask(const FEValuesExtractors::Scalar &scalar) const
bool operator!=(const FiniteElement< dim, spacedim > &) const
bool has_support_points() const
std::vector< std::vector< FullMatrix< double > > > restriction
bool has_generalized_support_points() const
virtual std::unique_ptr< FiniteElement< dim, spacedim > > clone() const =0
virtual Tensor< 2, dim > shape_grad_grad(const unsigned int i, const Point< dim > &p) const
std::vector< Table< 2, int > > adjust_quad_dof_index_for_face_orientation_table
virtual bool operator==(const FiniteElement< dim, spacedim > &fe) const
unsigned int adjust_line_dof_index_for_line_orientation(const unsigned int index, const bool line_orientation) const
const std::vector< Point< dim > > & get_unit_support_points() const
BlockIndices base_to_block_indices
void reinit_restriction_and_prolongation_matrices(const bool isotropic_restriction_only=false, const bool isotropic_prolongation_only=false)
const std::vector< Point< dim - 1 > > & get_unit_face_support_points(const unsigned int face_no=0) const
bool isotropic_restriction_is_implemented() const
std::vector< int > adjust_line_dof_index_for_line_orientation_table
virtual bool has_support_on_face(const unsigned int shape_index, const unsigned int face_index) 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 FiniteElement< dim, spacedim > & get_sub_fe(const ComponentMask &mask) const
virtual Tensor< 4, dim > shape_4th_derivative_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const
virtual double shape_value_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const
virtual void get_subface_interpolation_matrix(const FiniteElement< dim, spacedim > &source, const unsigned int subface, FullMatrix< double > &matrix, const unsigned int face_no=0) const
unsigned int component_to_block_index(const unsigned int component) const
virtual void fill_fe_face_values(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const unsigned int face_no, const hp::QCollection< dim - 1 > &quadrature, const Mapping< dim, spacedim > &mapping, const typename Mapping< dim, spacedim >::InternalDataBase &mapping_internal, const ::internal::FEValuesImplementation::MappingRelatedData< dim, spacedim > &mapping_data, const InternalDataBase &fe_internal, ::internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim > &output_data) const
BlockMask block_mask(const FEValuesExtractors::Scalar &scalar) const
virtual void get_interpolation_matrix(const FiniteElement< dim, spacedim > &source, FullMatrix< double > &matrix) const
virtual Tensor< 3, dim > shape_3rd_derivative(const unsigned int i, const Point< dim > &p) const
std::vector< std::pair< std::pair< unsigned int, unsigned int >, unsigned int > > system_to_base_table
virtual std::vector< std::pair< unsigned int, unsigned int > > hp_line_dof_identities(const FiniteElement< dim, spacedim > &fe_other) const
virtual void convert_generalized_support_point_values_to_dof_values(const std::vector< Vector< double > > &support_point_values, std::vector< double > &nodal_values) const
FiniteElement(const FiniteElementData< dim > &fe_data, const std::vector< bool > &restriction_is_additive_flags, const std::vector< ComponentMask > &nonzero_components)
std::vector< std::pair< unsigned int, unsigned int > > system_to_component_table
virtual Tensor< 3, dim > shape_3rd_derivative_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const
const FullMatrix< double > & constraints(const ::internal::SubfaceCase< dim > &subface_case=::internal::SubfaceCase< dim >::case_isotropic) const
std::pair< std::unique_ptr< FiniteElement< dim, spacedim > >, unsigned int > operator^(const unsigned int multiplicity) const
const std::vector< bool > restriction_is_additive_flags
std::vector< std::pair< std::pair< unsigned int, unsigned int >, unsigned int > > component_to_base_table
TableIndices< 2 > interface_constraints_size() const
std::vector< std::vector< std::pair< std::pair< unsigned int, unsigned int >, unsigned int > > > face_system_to_base_table
std::vector< Point< dim > > unit_support_points
virtual Tensor< 1, dim > shape_grad_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const
bool restriction_is_implemented() const
std::vector< std::vector< std::pair< unsigned int, unsigned int > > > face_system_to_component_table
FullMatrix< double > interface_constraints
const std::vector< Point< dim > > & get_generalized_support_points() const
virtual Tensor< 2, dim > shape_grad_grad_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const
virtual const FullMatrix< double > & get_prolongation_matrix(const unsigned int child, const RefinementCase< dim > &refinement_case=RefinementCase< dim >::isotropic_refinement) const
virtual double shape_value(const unsigned int i, const Point< dim > &p) const
virtual void get_face_interpolation_matrix(const FiniteElement< dim, spacedim > &source, FullMatrix< double > &matrix, const unsigned int face_no=0) const
virtual std::size_t memory_consumption() const
virtual Tensor< 4, dim > shape_4th_derivative(const unsigned int i, const Point< dim > &p) const
std::vector< Point< dim > > generalized_support_points
unsigned int adjust_quad_dof_index_for_face_orientation(const unsigned int index, const unsigned int face_no, const bool face_orientation, const bool face_flip, const bool face_rotation) const
const std::vector< ComponentMask > nonzero_components
virtual std::vector< std::pair< unsigned int, unsigned int > > hp_vertex_dof_identities(const FiniteElement< dim, spacedim > &fe_other) const
virtual bool hp_constraints_are_implemented() const
std::vector< std::vector< FullMatrix< double > > > prolongation
types::global_dof_index first_block_of_base(const unsigned int b) const
unsigned int size() const
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
@ update_default
No update.
static ::ExceptionBase & ExcNotImplemented()
#define Assert(cond, exc)
static ::ExceptionBase & ExcWrongInterfaceMatrixSize(int arg1, int arg2)
#define AssertDimension(dim1, dim2)
static ::ExceptionBase & ExcUnitShapeValuesDoNotExist()
static ::ExceptionBase & ExcProjectionVoid()
#define AssertIndexRange(index, range)
static ::ExceptionBase & ExcInternalError()
static ::ExceptionBase & ExcEmbeddingVoid()
static ::ExceptionBase & ExcDimensionMismatch(std::size_t arg1, std::size_t arg2)
static ::ExceptionBase & ExcMessage(std::string arg1)
static ::ExceptionBase & ExcFEHasNoSupportPoints()
#define AssertThrow(cond, exc)
@ neither_element_dominates
void reference_cell(Triangulation< dim, spacedim > &tria, const ReferenceCell &reference_cell)
std::enable_if< std::is_fundamental< T >::value, std::size_t >::type memory_consumption(const T &t)
constexpr const ReferenceCell Quadrilateral
VectorType::value_type * end(VectorType &V)
VectorType::value_type * begin(VectorType &V)
static const unsigned int invalid_unsigned_int
static unsigned int n_children(const RefinementCase< dim > &refinement_case)