16 #ifndef dealii_fe_system_h 17 # define dealii_fe_system_h 31 # include <type_traits> 40 template <
int dim,
int spacedim>
215 template <
int dim,
int spacedim = dim>
250 const unsigned int n_elements);
259 const unsigned int n1,
261 const unsigned int n2);
270 const unsigned int n1,
272 const unsigned int n2,
274 const unsigned int n3);
283 const unsigned int n1,
285 const unsigned int n2,
287 const unsigned int n3,
289 const unsigned int n4);
298 const unsigned int n1,
300 const unsigned int n2,
302 const unsigned int n3,
304 const unsigned int n4,
306 const unsigned int n5);
494 const std::vector<unsigned int> &multiplicities);
496 # if !defined(__INTEL_COMPILER) || __INTEL_COMPILER >= 1900 542 (std::is_same<typename std::decay<FEPairs>::type,
543 std::pair<std::unique_ptr<FiniteElement<dim, spacedim>>,
544 unsigned int>>::value ||
545 std::is_base_of<FiniteElement<dim, spacedim>,
546 typename std::decay<FEPairs>::type>::value)...>::type>
562 const std::initializer_list<
581 std::move(other_fe_system.generalized_support_points_index_table);
600 virtual std::unique_ptr<FiniteElement<dim, spacedim>>
601 clone()
const override;
613 get_sub_fe(
const unsigned int first_component,
614 const unsigned int n_selected_components)
const override;
642 const unsigned int component)
const override;
672 const unsigned int component)
const override;
702 const unsigned int component)
const override;
733 const unsigned int component)
const override;
764 const unsigned int component)
const override;
796 const unsigned int face_index)
const override;
823 const unsigned int child,
856 const unsigned int child,
900 const unsigned int face,
901 const bool face_orientation =
true,
902 const bool face_flip =
false,
903 const bool face_rotation =
false)
const override;
914 virtual Point<dim - 1>
916 const unsigned int face_no = 0)
const override;
925 virtual std::pair<Table<2, bool>, std::vector<unsigned int>>
958 const unsigned int face_no = 0)
const override;
976 const unsigned int subface,
978 const unsigned int face_no = 0)
const override;
995 virtual std::vector<std::pair<unsigned int, unsigned int>>
1003 virtual std::vector<std::pair<unsigned int, unsigned int>>
1011 virtual std::vector<std::pair<unsigned int, unsigned int>>
1013 const unsigned int face_no = 0)
const override;
1020 const unsigned int codim = 0)
const override final;
1040 const std::vector<Vector<double>> &support_point_values,
1041 std::vector<double> & dof_values)
const override;
1055 virtual std::unique_ptr<
1063 &output_data)
const override;
1067 virtual std::unique_ptr<
1075 &output_data)
const override;
1077 virtual std::unique_ptr<
1085 &output_data)
const override;
1094 const ::internal::FEValuesImplementation::MappingRelatedData<dim,
1100 &output_data)
const override;
1107 const unsigned int face_no,
1111 const ::internal::FEValuesImplementation::MappingRelatedData<dim,
1117 &output_data)
const override;
1122 const unsigned int face_no,
1123 const unsigned int sub_no,
1127 const ::internal::FEValuesImplementation::MappingRelatedData<dim,
1133 &output_data)
const override;
1145 template <
int dim_1>
1150 const unsigned int face_no,
1151 const unsigned int sub_no,
1159 &output_data)
const;
1174 std::vector<std::pair<std::unique_ptr<const FiniteElement<dim, spacedim>>,
1199 const std::vector<unsigned int> &multiplicities);
1212 template <
int structdim>
1213 std::vector<std::pair<unsigned int, unsigned int>>
1215 const unsigned int face_no = 0)
const;
1244 const unsigned int base_no,
1273 typename std::vector<
1274 std::unique_ptr<typename FiniteElement<dim, spacedim>::InternalDataBase>>
1285 mutable std::vector<
1303 namespace FESystemImplementation
1305 template <
int dim,
int spacedim>
1308 const std::initializer_list<
1312 return std::count_if(
1315 [](
const std::pair<std::unique_ptr<FiniteElement<dim, spacedim>>,
1316 unsigned int> &fe_system) {
1317 return fe_system.second > 0;
1323 template <
int dim,
int spacedim>
1324 std::pair<std::unique_ptr<FiniteElement<dim, spacedim>>,
unsigned int>
1327 return std::make_pair(std::move(fe.
clone()), 1u);
1332 template <
int dim,
int spacedim>
1340 return std::forward<
1341 std::pair<std::unique_ptr<FiniteElement<dim, spacedim>>,
unsigned int>>(
1349 # if !defined(__INTEL_COMPILER) || __INTEL_COMPILER >= 1900 1354 template <
int dim,
int spacedim>
1355 template <
class... FEPairs,
typename>
1358 {internal::FESystemImplementation::promote_to_fe_pair<dim, spacedim>(
1359 std::forward<FEPairs>(fe_pairs))...})
1364 template <
int dim,
int spacedim>
1366 const std::initializer_list<
1376 , base_elements(
internal::FESystemImplementation::count_nonzeros(fe_systems))
1378 std::vector<const FiniteElement<dim, spacedim> *> fes;
1379 std::vector<unsigned int> multiplicities;
1382 [&fes, &multiplicities](
1383 const std::pair<std::unique_ptr<FiniteElement<dim, spacedim>>,
1384 unsigned int> &fe_system) {
1385 fes.push_back(fe_system.first.get());
1386 multiplicities.push_back(fe_system.second);
1389 for (
const auto &p : fe_systems)
virtual std::size_t memory_consumption() const override
virtual Point< dim > unit_support_point(const unsigned int index) const override
std::vector< internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim > > base_fe_output_objects
static const unsigned int invalid_unsigned_int
virtual double shape_value_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const override
void build_interface_constraints()
FESystem(FESystem< dim, spacedim > &&other_fe_system) noexcept
void initialize(const std::vector< const FiniteElement< dim, spacedim > *> &fes, const std::vector< unsigned int > &multiplicities)
virtual std::string get_name() const override
virtual std::vector< std::pair< unsigned int, unsigned int > > hp_line_dof_identities(const FiniteElement< dim, spacedim > &fe_other) const override
Contents is actually a matrix.
virtual Tensor< 4, dim > shape_4th_derivative(const unsigned int i, const Point< dim > &p) const override
virtual std::pair< Table< 2, bool >, std::vector< unsigned int > > get_constant_modes() const override
virtual ~FESystem() override=default
std::vector< std::pair< unsigned int, unsigned int > > hp_object_dof_identities(const FiniteElement< dim, spacedim > &fe_other, const unsigned int face_no=0) const
virtual void convert_generalized_support_point_values_to_dof_values(const std::vector< Vector< double >> &support_point_values, std::vector< double > &dof_values) const override
virtual bool has_support_on_face(const unsigned int shape_index, const unsigned int face_index) const override
virtual std::unique_ptr< typename FiniteElement< dim, spacedim >::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 override
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 override
virtual Point< dim - 1 > unit_face_support_point(const unsigned int index, const unsigned int face_no=0) const override
virtual Tensor< 1, dim > shape_grad_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const override
virtual const FiniteElement< dim, spacedim > & base_element(const unsigned int index) const override
virtual std::unique_ptr< FiniteElement< dim, spacedim > > clone() const =0
virtual double shape_value(const unsigned int i, const Point< dim > &p) const override
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 override
void compute_fill(const Mapping< dim, spacedim > &mapping, const typename Triangulation< dim, spacedim >::cell_iterator &cell, const unsigned int face_no, const unsigned int sub_no, const hp::QCollection< dim_1 > &quadrature, const CellSimilarity::Similarity cell_similarity, const typename Mapping< dim, spacedim >::InternalDataBase &mapping_internal, const typename FiniteElement< dim, spacedim >::InternalDataBase &fe_data, const internal::FEValuesImplementation::MappingRelatedData< dim, spacedim > &mapping_data, internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim > &output_data) const
virtual const FiniteElement< dim, spacedim > & get_sub_fe(const unsigned int first_component, const unsigned int n_selected_components) const override
virtual Tensor< 4, dim > shape_4th_derivative_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const override
InternalData(const unsigned int n_base_elements)
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 override
virtual Tensor< 3, dim > shape_3rd_derivative(const unsigned int i, const Point< dim > &p) const override
std::vector< std::vector< std::size_t > > generalized_support_points_index_table
virtual std::unique_ptr< FiniteElement< dim, spacedim > > clone() const override
virtual UpdateFlags requires_update_flags(const UpdateFlags update_flags) const override
virtual void fill_fe_values(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const CellSimilarity::Similarity cell_similarity, const Quadrature< dim > &quadrature, const Mapping< dim, spacedim > &mapping, const typename Mapping< dim, spacedim >::InternalDataBase &mapping_internal, const ::internal::FEValuesImplementation::MappingRelatedData< dim, spacedim > &mapping_data, const typename FiniteElement< dim, spacedim >::InternalDataBase &fe_internal, ::internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim > &output_data) const override
constexpr ReturnType< rank, T >::value_type & extract(T &t, const ArrayType &indices)
static const unsigned int invalid_face_number
virtual const FullMatrix< double > & get_prolongation_matrix(const unsigned int child, const RefinementCase< dim > &refinement_case=RefinementCase< dim >::isotropic_refinement) const override
Abstract base class for mapping classes.
#define DEAL_II_NAMESPACE_CLOSE
void set_fe_data(const unsigned int base_no, std::unique_ptr< typename FiniteElement< dim, spacedim >::InternalDataBase >)
virtual std::unique_ptr< typename FiniteElement< dim, spacedim >::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 override
virtual Tensor< 2, dim > shape_grad_grad(const unsigned int i, const Point< dim > &p) const override
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 typename FiniteElement< dim, spacedim >::InternalDataBase &fe_internal, ::internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim > &output_data) const override
virtual bool hp_constraints_are_implemented() const override
std::vector< std::unique_ptr< typename FiniteElement< dim, spacedim >::InternalDataBase > > base_fe_datas
virtual std::vector< std::pair< unsigned int, unsigned int > > hp_vertex_dof_identities(const FiniteElement< dim, spacedim > &fe_other) const override
virtual FiniteElementDomination::Domination compare_for_domination(const FiniteElement< dim, spacedim > &fe_other, const unsigned int codim=0) const override final
virtual Tensor< 3, dim > shape_3rd_derivative_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const override
virtual void get_interpolation_matrix(const FiniteElement< dim, spacedim > &source, FullMatrix< double > &matrix) const override
FiniteElement< dim, spacedim >::InternalDataBase & get_fe_data(const unsigned int base_no) const
virtual const FullMatrix< double > & get_restriction_matrix(const unsigned int child, const RefinementCase< dim > &refinement_case=RefinementCase< dim >::isotropic_refinement) const override
#define DEAL_II_NAMESPACE_OPEN
virtual void get_face_interpolation_matrix(const FiniteElement< dim, spacedim > &source, FullMatrix< double > &matrix, const unsigned int face_no=0) const override
FESystem(const FiniteElement< dim, spacedim > &fe, const unsigned int n_elements)
internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim > & get_fe_output_object(const unsigned int base_no) const
virtual Tensor< 1, dim > shape_grad(const unsigned int i, const Point< dim > &p) const override
virtual Tensor< 2, dim > shape_grad_grad_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const override
unsigned int n_base_elements() const
virtual std::unique_ptr< typename FiniteElement< dim, spacedim >::InternalDataBase > get_data(const UpdateFlags update_flags, const Mapping< dim, spacedim > &mapping, const Quadrature< dim > &quadrature, ::internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim > &output_data) const override
std::vector< std::pair< std::unique_ptr< const FiniteElement< dim, spacedim > >, unsigned int > > base_elements
virtual void fill_fe_subface_values(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const unsigned int face_no, const unsigned int sub_no, const Quadrature< 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 typename FiniteElement< dim, spacedim >::InternalDataBase &fe_internal, ::internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim > &output_data) const override