16 #include <deal.II/base/memory_consumption.h> 17 #include <deal.II/base/qprojector.h> 18 #include <deal.II/base/quadrature.h> 20 #include <deal.II/dofs/dof_accessor.h> 22 #include <deal.II/fe/fe.h> 23 #include <deal.II/fe/fe_values.h> 24 #include <deal.II/fe/mapping.h> 26 #include <deal.II/grid/tria.h> 27 #include <deal.II/grid/tria_iterator.h> 34 DEAL_II_NAMESPACE_OPEN
40 template <
int dim,
int spacedim>
47 template <
int dim,
int spacedim>
56 template <
int dim,
int spacedim>
59 const std::vector<bool> & r_i_a_f,
60 const std::vector<ComponentMask> &nonzero_c)
71 std::make_pair(
std::make_pair(0U, 0U), 0U))
80 nonzero_c.size() == 1 ?
86 std::bind(
std::not_equal_to<unsigned int>(),
132 ref < RefinementCase<dim>::isotropic_refinement + 1;
148 template <
int dim,
int spacedim>
149 std::pair<std::unique_ptr<FiniteElement<dim, spacedim>>,
unsigned int>
152 return {this->
clone(), multiplicity};
157 template <
int dim,
int spacedim>
168 template <
int dim,
int spacedim>
172 const unsigned int)
const 180 template <
int dim,
int spacedim>
191 template <
int dim,
int spacedim>
195 const unsigned int)
const 203 template <
int dim,
int spacedim>
214 template <
int dim,
int spacedim>
219 const unsigned int)
const 227 template <
int dim,
int spacedim>
238 template <
int dim,
int spacedim>
243 const unsigned int)
const 251 template <
int dim,
int spacedim>
262 template <
int dim,
int spacedim>
267 const unsigned int)
const 274 template <
int dim,
int spacedim>
277 const bool isotropic_restriction_only,
278 const bool isotropic_prolongation_only)
281 ref_case <= RefinementCase<dim>::isotropic_refinement;
284 const unsigned int nc =
287 for (
unsigned int i = 0; i < nc; ++i)
290 (!isotropic_restriction_only ||
295 (!isotropic_prolongation_only ||
304 template <
int dim,
int spacedim>
307 const unsigned int child,
316 "Restriction matrices are only available for refined cells!"));
333 template <
int dim,
int spacedim>
336 const unsigned int child,
345 "Prolongation matrices are only available for refined cells!"));
364 template <
int dim,
int spacedim>
367 const unsigned int index)
const 377 template <
int dim,
int spacedim>
395 template <
int dim,
int spacedim>
401 this->n_components());
417 template <
int dim,
int spacedim>
424 this->n_components());
442 template <
int dim,
int spacedim>
456 component_mask[c] =
true;
463 template <
int dim,
int spacedim>
474 template <
int dim,
int spacedim>
485 template <
int dim,
int spacedim>
497 template <
int dim,
int spacedim>
522 if (component_mask[c] ==
true)
523 block_mask[block] =
true;
532 Assert(component_mask[c] == block_mask[block],
534 "The component mask argument given to this function " 535 "is not a mask where the individual components belonging " 536 "to one block of the finite element are either all " 537 "selected or not selected. You can't call this function " 538 "with a component mask that splits blocks."));
549 template <
int dim,
int spacedim>
552 const unsigned int face,
553 const bool face_orientation,
554 const bool face_flip,
555 const bool face_rotation)
const 575 if ((face_orientation !=
true) || (face_flip !=
false) ||
576 (face_rotation !=
false))
579 "The function in this base class can not handle this case. " 580 "Rather, the derived class you are using must provide " 581 "an overloaded version but apparently hasn't done so. See " 582 "the documentation of this function for more information."));
592 const unsigned int dof_index_on_vertex =
598 face, face_vertex, face_orientation, face_flip, face_rotation) *
600 dof_index_on_vertex);
610 const unsigned int dof_index_on_line = index % this->
dofs_per_line;
614 face, face_line, face_orientation, face_flip, face_rotation) *
632 template <
int dim,
int spacedim>
635 const unsigned int index,
636 const bool face_orientation,
637 const bool face_flip,
638 const bool face_rotation)
const 663 index, 4 * face_orientation + 2 * face_flip + face_rotation);
668 template <
int dim,
int spacedim>
671 const unsigned int index,
672 const bool line_orientation)
const 686 if (line_orientation)
694 template <
int dim,
int spacedim>
699 ref_case < RefinementCase<dim>::isotropic_refinement + 1;
701 for (
unsigned int c = 0;
722 template <
int dim,
int spacedim>
727 ref_case < RefinementCase<dim>::isotropic_refinement + 1;
729 for (
unsigned int c = 0;
750 template <
int dim,
int spacedim>
757 for (
unsigned int c = 0;
778 template <
int dim,
int spacedim>
785 for (
unsigned int c = 0;
806 template <
int dim,
int spacedim>
819 template <
int dim,
int spacedim>
828 template <
int dim,
int spacedim>
835 ExcMessage(
"Constraints for this element are only implemented " 836 "for the case that faces are refined isotropically " 837 "(which is always the case in 2d, and in 3d requires " 838 "that the neighboring cell of a coarse cell presents " 839 "exactly four children on the common face)."));
841 ExcMessage(
"The finite element for which you try to obtain " 842 "hanging node constraints does not appear to " 855 template <
int dim,
int spacedim>
868 12 * this->dofs_per_line +
880 template <
int dim,
int spacedim>
896 template <
int dim,
int spacedim>
912 template <
int dim,
int spacedim>
929 template <
int dim,
int spacedim>
930 std::vector<std::pair<unsigned int, unsigned int>>
935 return std::vector<std::pair<unsigned int, unsigned int>>();
940 template <
int dim,
int spacedim>
941 std::vector<std::pair<unsigned int, unsigned int>>
946 return std::vector<std::pair<unsigned int, unsigned int>>();
951 template <
int dim,
int spacedim>
952 std::vector<std::pair<unsigned int, unsigned int>>
957 return std::vector<std::pair<unsigned int, unsigned int>>();
962 template <
int dim,
int spacedim>
973 template <
int dim,
int spacedim>
980 return ((
typeid(*
this) ==
typeid(f)) && (this->
get_name() == f.
get_name()) &&
988 template <
int dim,
int spacedim>
993 return !(*
this == f);
998 template <
int dim,
int spacedim>
999 const std::vector<Point<dim>> &
1014 template <
int dim,
int spacedim>
1023 template <
int dim,
int spacedim>
1024 const std::vector<Point<dim>> &
1036 template <
int dim,
int spacedim>
1045 template <
int dim,
int spacedim>
1058 template <
int dim,
int spacedim>
1059 const std::vector<
Point<dim - 1>> &
1074 template <
int dim,
int spacedim>
1083 template <
int dim,
int spacedim>
1084 const std::vector<
Point<dim - 1>> &
1098 template <
int dim,
int spacedim>
1107 template <
int dim,
int spacedim>
1110 const unsigned int index)
const 1121 template <
int dim,
int spacedim>
1124 const unsigned int)
const 1131 template <
int dim,
int spacedim>
1137 const unsigned int n_total_components = this->
n_components();
1138 Assert((n_total_components == mask.
size()) || (mask.
size() == 0),
1139 ExcMessage(
"The given ComponentMask has the wrong size."));
1141 const unsigned int n_selected =
1144 ExcMessage(
"You need at least one selected component."));
1146 const unsigned int first_selected =
1151 for (
unsigned int c = 0; c < n_total_components; ++c)
1152 Assert((c < first_selected && (!mask[c])) ||
1153 (c >= first_selected && c < first_selected + n_selected &&
1155 (c >= first_selected + n_selected && !mask[c]),
1156 ExcMessage(
"Error: the given ComponentMask is not contiguous!"));
1159 return get_sub_fe(first_selected, n_selected);
1164 template <
int dim,
int spacedim>
1167 const unsigned int first_component,
1168 const unsigned int n_selected_components)
const 1174 "You can only select a whole FiniteElement, not a part of one."));
1176 (void)first_component;
1177 (void)n_selected_components;
1183 template <
int dim,
int spacedim>
1184 std::pair<Table<2, bool>, std::vector<unsigned int>>
1188 return std::pair<Table<2, bool>, std::vector<unsigned int>>(
1195 template <
int dim,
int spacedim>
1200 std::vector<double> &)
const 1203 ExcMessage(
"The element for which you are calling the current " 1204 "function does not have generalized support points (see " 1205 "the glossary for a definition of generalized support " 1206 "points). Consequently, the current function can not " 1207 "be defined and is not implemented by the element."));
1213 template <
int dim,
int spacedim>
1234 template <
int dim,
int spacedim>
1235 std::vector<unsigned int>
1239 std::vector<unsigned int> retval(nonzero_components.size());
1240 for (
unsigned int i = 0; i < nonzero_components.size(); ++i)
1241 retval[i] = nonzero_components[i].n_selected_components();
1249 template <
int dim,
int spacedim>
1250 std::unique_ptr<typename FiniteElement<dim, spacedim>::InternalDataBase>
1267 template <
int dim,
int spacedim>
1268 std::unique_ptr<typename FiniteElement<dim, spacedim>::InternalDataBase>
1285 template <
int dim,
int spacedim>
1303 DEAL_II_NAMESPACE_CLOSE
static::ExceptionBase & ExcFEHasNoSupportPoints()
virtual double shape_value(const unsigned int i, const Point< dim > &p) const
bool prolongation_is_implemented() const
virtual std::pair< Table< 2, bool >, std::vector< unsigned int > > get_constant_modes() const
static const unsigned int invalid_unsigned_int
const std::vector< Point< dim > > & get_unit_support_points() const
std::vector< std::vector< FullMatrix< double > > > restriction
const unsigned int components
#define AssertDimension(dim1, dim2)
virtual std::vector< std::pair< unsigned int, unsigned int > > hp_line_dof_identities(const FiniteElement< dim, spacedim > &fe_other) const
std::vector< Point< dim > > generalized_support_points
FullMatrix< double > interface_constraints
FiniteElement(const FiniteElementData< dim > &fe_data, const std::vector< bool > &restriction_is_additive_flags, const std::vector< ComponentMask > &nonzero_components)
virtual Tensor< 3, dim > shape_3rd_derivative(const unsigned int i, const Point< dim > &p) const
static::ExceptionBase & ExcUnitShapeValuesDoNotExist()
bool isotropic_restriction_is_implemented() const
unsigned int n_selected_components(const unsigned int overall_number_of_components=numbers::invalid_unsigned_int) const
const std::vector< Point< dim-1 > > & get_generalized_face_support_points() const
virtual void get_interpolation_matrix(const FiniteElement< dim, spacedim > &source, FullMatrix< double > &matrix) const
bool constraints_are_implemented(const ::internal::SubfaceCase< dim > &subface_case=::internal::SubfaceCase< dim >::case_isotropic) const
#define AssertIndexRange(index, range)
const unsigned int dofs_per_quad
virtual std::vector< std::pair< unsigned int, unsigned int > > hp_quad_dof_identities(const FiniteElement< dim, spacedim > &fe_other) const
std::vector< std::pair< std::pair< unsigned int, unsigned int >, unsigned int > > component_to_base_table
std::vector< Point< dim-1 > > unit_face_support_points
virtual bool hp_constraints_are_implemented() const
bool isotropic_prolongation_is_implemented() const
const FiniteElement< dim, spacedim > & get_sub_fe(const ComponentMask &mask) const
unsigned int size() const
TableIndices< 2 > interface_constraints_size() const
virtual std::unique_ptr< FiniteElement< dim, spacedim > > clone() const =0
#define AssertThrow(cond, exc)
virtual Tensor< 2, dim > shape_grad_grad(const unsigned int i, const Point< dim > &p) const
static std::vector< unsigned int > compute_n_nonzero_components(const std::vector< ComponentMask > &nonzero_components)
bool operator!=(const FiniteElement< dim, spacedim > &) const
bool is_primitive() const
virtual bool has_support_on_face(const unsigned int shape_index, const unsigned int face_index) const
std::vector< Point< dim-1 > > generalized_face_support_points
const std::vector< unsigned int > n_nonzero_components_table
static::ExceptionBase & ExcIndexRange(int arg1, int arg2, int arg3)
const unsigned int dofs_per_line
BlockMask block_mask(const FEValuesExtractors::Scalar &scalar) const
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 std::vector< std::pair< unsigned int, unsigned int > > hp_vertex_dof_identities(const FiniteElement< dim, spacedim > &fe_other) const
std::vector< Point< dim > > unit_support_points
const unsigned int first_face_line_index
bool restriction_is_implemented() const
virtual const FullMatrix< double > & get_prolongation_matrix(const unsigned int child, const RefinementCase< dim > &refinement_case=RefinementCase< dim >::isotropic_refinement) const
static::ExceptionBase & ExcMessage(std::string arg1)
const unsigned int first_quad_index
std::vector< std::vector< FullMatrix< double > > > prolongation
#define Assert(cond, exc)
unsigned int component_to_block_index(const unsigned int component) const
virtual Tensor< 1, dim > shape_grad_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const
static::ExceptionBase & ExcDimensionMismatch(std::size_t arg1, std::size_t arg2)
virtual double shape_value_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const
unsigned int adjust_quad_dof_index_for_face_orientation(const unsigned int index, const bool face_orientation, const bool face_flip, const bool face_rotation) const
Abstract base class for mapping classes.
void reinit(const unsigned int n_blocks, const size_type n_elements_per_block)
unsigned int size() const
unsigned int adjust_line_dof_index_for_line_orientation(const unsigned int index, const bool line_orientation) const
static::ExceptionBase & ExcWrongInterfaceMatrixSize(int arg1, int arg2)
bool represents_the_all_selected_mask() const
virtual Point< dim > unit_support_point(const unsigned int index) const
const FullMatrix< double > & constraints(const ::internal::SubfaceCase< dim > &subface_case=::internal::SubfaceCase< dim >::case_isotropic) const
virtual std::string get_name() const =0
bool represents_the_all_selected_mask() const
static::ExceptionBase & ExcProjectionVoid()
bool has_generalized_support_points() const
virtual Point< dim-1 > unit_face_support_point(const unsigned int index) const
virtual bool operator==(const FiniteElement< dim, spacedim > &fe) const
unsigned int n_components() const
virtual std::size_t memory_consumption() const
virtual Tensor< 1, dim > shape_grad(const unsigned int i, const Point< dim > &p) const
virtual void get_face_interpolation_matrix(const FiniteElement< dim, spacedim > &source, FullMatrix< double > &matrix) const
virtual void get_subface_interpolation_matrix(const FiniteElement< dim, spacedim > &source, const unsigned int subface, FullMatrix< double > &matrix) const
unsigned int first_selected_component(const unsigned int overall_number_of_components=numbers::invalid_unsigned_int) const
const unsigned int dofs_per_cell
const std::vector< Point< dim > > & get_generalized_support_points() const
static unsigned int n_children(const RefinementCase< dim > &refinement_case)
size_type n_elements() const
virtual std::unique_ptr< InternalDataBase > get_face_data(const UpdateFlags update_flags, const Mapping< dim, spacedim > &mapping, const Quadrature< dim-1 > &quadrature,::internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim > &output_data) const
virtual Tensor< 3, dim > shape_3rd_derivative_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const
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
Table< 2, int > adjust_quad_dof_index_for_face_orientation_table
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
virtual const FullMatrix< double > & get_restriction_matrix(const unsigned int child, const RefinementCase< dim > &refinement_case=RefinementCase< dim >::isotropic_refinement) const
static::ExceptionBase & ExcEmbeddingVoid()
virtual std::size_t memory_consumption() const
virtual FiniteElementDomination::Domination compare_for_face_domination(const FiniteElement< dim, spacedim > &fe_other) const
const unsigned int first_face_quad_index
std::vector< std::pair< unsigned int, unsigned int > > system_to_component_table
virtual Tensor< 4, dim > shape_4th_derivative(const unsigned int i, const Point< dim > &p) const
unsigned int n_blocks() const
ComponentMask component_mask(const FEValuesExtractors::Scalar &scalar) const
const unsigned int dofs_per_face
const unsigned int first_line_index
const bool cached_primitivity
const std::vector< ComponentMask > nonzero_components
void reinit_restriction_and_prolongation_matrices(const bool isotropic_restriction_only=false, const bool isotropic_prolongation_only=false)
virtual Tensor< 4, dim > shape_4th_derivative_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const
std::vector< std::pair< std::pair< unsigned int, unsigned int >, unsigned int > > system_to_base_table
static::ExceptionBase & ExcNotImplemented()
bool has_generalized_face_support_points() const
unsigned int size() const
const unsigned int dofs_per_vertex
std::vector< std::pair< unsigned int, unsigned int > > face_system_to_component_table
std::pair< std::unique_ptr< FiniteElement< dim, spacedim > >, unsigned int > operator^(const unsigned int multiplicity) const
const std::vector< Point< dim-1 > > & get_unit_face_support_points() const
BlockIndices base_to_block_indices
std::vector< std::pair< std::pair< unsigned int, unsigned int >, unsigned int > > face_system_to_base_table
bool has_support_points() const
types::global_dof_index first_block_of_base(const unsigned int b) const
bool has_face_support_points() const
virtual const FiniteElement< dim, spacedim > & base_element(const unsigned int index) const
std::vector< int > adjust_line_dof_index_for_line_orientation_table
void fill(InputIterator entries, const bool C_style_indexing=true)
const std::vector< bool > restriction_is_additive_flags
virtual void convert_generalized_support_point_values_to_dof_values(const std::vector< Vector< double >> &support_point_values, std::vector< double > &nodal_values) const
std::enable_if< std::is_fundamental< T >::value, std::size_t >::type memory_consumption(const T &t)
virtual Tensor< 2, dim > shape_grad_grad_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const
static::ExceptionBase & ExcInternalError()