Reference documentation for deal.II version GIT relicensing-878-g25c11e0287 2024-06-21 14:50:02+00:00
\(\newcommand{\dealvcentcolon}{\mathrel{\mathop{:}}}\) \(\newcommand{\dealcoloneq}{\dealvcentcolon\mathrel{\mkern-1.2mu}=}\) \(\newcommand{\jump}[1]{\left[\!\left[ #1 \right]\!\right]}\) \(\newcommand{\average}[1]{\left\{\!\left\{ #1 \right\}\!\right\}}\)
Loading...
Searching...
No Matches
fe_system.h
Go to the documentation of this file.
1// ------------------------------------------------------------------------
2//
3// SPDX-License-Identifier: LGPL-2.1-or-later
4// Copyright (C) 1999 - 2023 by the deal.II authors
5//
6// This file is part of the deal.II library.
7//
8// Part of the source code is dual licensed under Apache-2.0 WITH
9// LLVM-exception OR LGPL-2.1-or-later. Detailed license information
10// governing the source code and code contributions can be found in
11// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II.
12//
13// ------------------------------------------------------------------------
14
15#ifndef dealii_fe_system_h
16#define dealii_fe_system_h
17
18#include <deal.II/base/config.h>
19
20#include <deal.II/fe/fe.h>
21#include <deal.II/fe/fe_tools.h>
22
23#include <memory>
24#include <type_traits>
25#include <utility>
26#include <vector>
27
29
30// Forward declaration
31#ifndef DOXYGEN
32template <int dim, int spacedim>
33class FE_Enriched;
34#endif
35
206template <int dim, int spacedim = dim>
207class FESystem : public FiniteElement<dim, spacedim>
208{
209public:
214 FESystem() = delete;
215
247 const unsigned int n_elements);
248
256 const unsigned int n1,
258 const unsigned int n2);
259
267 const unsigned int n1,
269 const unsigned int n2,
271 const unsigned int n3);
272
280 const unsigned int n1,
282 const unsigned int n2,
284 const unsigned int n3,
286 const unsigned int n4);
287
295 const unsigned int n1,
297 const unsigned int n2,
299 const unsigned int n3,
301 const unsigned int n4,
303 const unsigned int n5);
304
490 FESystem(const std::vector<const FiniteElement<dim, spacedim> *> &fes,
491 const std::vector<unsigned int> &multiplicities);
492
493#if !defined(__INTEL_COMPILER) || __INTEL_COMPILER >= 1900
536 template <
537 class... FEPairs,
538 typename = enable_if_all_t<
539 (std::is_same_v<std::decay_t<FEPairs>,
540 std::pair<std::unique_ptr<FiniteElement<dim, spacedim>>,
541 unsigned int>> ||
542 std::is_base_of_v<FiniteElement<dim, spacedim>,
543 std::decay_t<FEPairs>>)...>>
544 FESystem(FEPairs &&...fe_pairs);
545
559 const std::initializer_list<
560 std::pair<std::unique_ptr<FiniteElement<dim, spacedim>>, unsigned int>>
561 &fe_systems);
562#endif
563
569
573 FESystem(FESystem<dim, spacedim> &&other_fe_system) noexcept
574 : FiniteElement<dim, spacedim>(std::move(other_fe_system))
575 {
576 base_elements = std::move(other_fe_system.base_elements);
578 std::move(other_fe_system.generalized_support_points_index_table);
579 }
580
584 virtual ~FESystem() override = default;
585
594 virtual std::string
595 get_name() const override;
596
597 virtual std::unique_ptr<FiniteElement<dim, spacedim>>
598 clone() const override;
599
600 virtual UpdateFlags
601 requires_update_flags(const UpdateFlags update_flags) const override;
602
603 // make variant with ComponentMask also available:
604 using FiniteElement<dim, spacedim>::get_sub_fe;
605
609 virtual const FiniteElement<dim, spacedim> &
610 get_sub_fe(const unsigned int first_component,
611 const unsigned int n_selected_components) const override;
612
625 virtual double
626 shape_value(const unsigned int i, const Point<dim> &p) const override;
627
636 virtual double
637 shape_value_component(const unsigned int i,
638 const Point<dim> &p,
639 const unsigned int component) const override;
640
655 virtual Tensor<1, dim>
656 shape_grad(const unsigned int i, const Point<dim> &p) const override;
657
666 virtual Tensor<1, dim>
667 shape_grad_component(const unsigned int i,
668 const Point<dim> &p,
669 const unsigned int component) const override;
670
686 virtual Tensor<2, dim>
687 shape_grad_grad(const unsigned int i, const Point<dim> &p) const override;
688
697 virtual Tensor<2, dim>
698 shape_grad_grad_component(const unsigned int i,
699 const Point<dim> &p,
700 const unsigned int component) const override;
701
716 virtual Tensor<3, dim>
717 shape_3rd_derivative(const unsigned int i,
718 const Point<dim> &p) const override;
719
728 virtual Tensor<3, dim>
729 shape_3rd_derivative_component(const unsigned int i,
730 const Point<dim> &p,
731 const unsigned int component) const override;
732
747 virtual Tensor<4, dim>
748 shape_4th_derivative(const unsigned int i,
749 const Point<dim> &p) const override;
750
759 virtual Tensor<4, dim>
760 shape_4th_derivative_component(const unsigned int i,
761 const Point<dim> &p,
762 const unsigned int component) const override;
763
775 virtual void
777 FullMatrix<double> &matrix) const override;
778
785 virtual const FiniteElement<dim, spacedim> &
786 base_element(const unsigned int index) const override;
787
792 virtual bool
793 has_support_on_face(const unsigned int shape_index,
794 const unsigned int face_index) const override;
795
819 virtual const FullMatrix<double> &
821 const unsigned int child,
822 const RefinementCase<dim> &refinement_case =
824
852 virtual const FullMatrix<double> &
854 const unsigned int child,
855 const RefinementCase<dim> &refinement_case =
857
889 virtual unsigned int
891 const unsigned int face_dof_index,
892 const unsigned int face,
893 const unsigned char combined_orientation =
895
899 virtual Point<dim>
900 unit_support_point(const unsigned int index) const override;
901
905 virtual Point<dim - 1>
906 unit_face_support_point(const unsigned int index,
907 const unsigned int face_no = 0) const override;
908
916 virtual std::pair<Table<2, bool>, std::vector<unsigned int>>
917 get_constant_modes() const override;
918
931 virtual bool
933
946 virtual void
948 FullMatrix<double> &matrix,
949 const unsigned int face_no = 0) const override;
950
951
964 virtual void
966 const FiniteElement<dim, spacedim> &source,
967 const unsigned int subface,
968 FullMatrix<double> &matrix,
969 const unsigned int face_no = 0) const override;
970
986 virtual std::vector<std::pair<unsigned int, unsigned int>>
988 const FiniteElement<dim, spacedim> &fe_other) const override;
989
994 virtual std::vector<std::pair<unsigned int, unsigned int>>
996 const FiniteElement<dim, spacedim> &fe_other) const override;
997
1002 virtual std::vector<std::pair<unsigned int, unsigned int>>
1004 const unsigned int face_no = 0) const override;
1005
1011 const unsigned int codim = 0) const override final;
1012
1029 virtual void
1031 const std::vector<Vector<double>> &support_point_values,
1032 std::vector<double> &dof_values) const override;
1033
1042 virtual std::size_t
1043 memory_consumption() const override;
1044
1057 {
1062
1068 unsigned int in_index;
1069
1075 unsigned int out_index;
1076 };
1077
1078protected:
1079 virtual std::unique_ptr<
1082 const UpdateFlags update_flags,
1083 const Mapping<dim, spacedim> &mapping,
1084 const Quadrature<dim> &quadrature,
1086 spacedim>
1087 &output_data) const override;
1088
1089 using FiniteElement<dim, spacedim>::get_face_data;
1090
1091 virtual std::unique_ptr<
1094 const UpdateFlags update_flags,
1095 const Mapping<dim, spacedim> &mapping,
1096 const hp::QCollection<dim - 1> &quadrature,
1098 spacedim>
1099 &output_data) const override;
1100
1101 virtual std::unique_ptr<
1104 const UpdateFlags update_flags,
1105 const Mapping<dim, spacedim> &mapping,
1106 const Quadrature<dim - 1> &quadrature,
1108 spacedim>
1109 &output_data) const override;
1110
1111 virtual void
1114 const CellSimilarity::Similarity cell_similarity,
1115 const Quadrature<dim> &quadrature,
1116 const Mapping<dim, spacedim> &mapping,
1117 const typename Mapping<dim, spacedim>::InternalDataBase &mapping_internal,
1119 &mapping_data,
1120 const typename FiniteElement<dim, spacedim>::InternalDataBase &fe_internal,
1122 spacedim>
1123 &output_data) const override;
1124
1125 using FiniteElement<dim, spacedim>::fill_fe_face_values;
1126
1127 virtual void
1130 const unsigned int face_no,
1131 const hp::QCollection<dim - 1> &quadrature,
1132 const Mapping<dim, spacedim> &mapping,
1133 const typename Mapping<dim, spacedim>::InternalDataBase &mapping_internal,
1135 &mapping_data,
1136 const typename FiniteElement<dim, spacedim>::InternalDataBase &fe_internal,
1138 spacedim>
1139 &output_data) const override;
1140
1141 virtual void
1144 const unsigned int face_no,
1145 const unsigned int sub_no,
1146 const Quadrature<dim - 1> &quadrature,
1147 const Mapping<dim, spacedim> &mapping,
1148 const typename Mapping<dim, spacedim>::InternalDataBase &mapping_internal,
1150 &mapping_data,
1151 const typename FiniteElement<dim, spacedim>::InternalDataBase &fe_internal,
1153 spacedim>
1154 &output_data) const override;
1155
1166 template <class Q_or_QC>
1167 void
1169 const Mapping<dim, spacedim> &mapping,
1171 const unsigned int face_no,
1172 const unsigned int sub_no,
1173 const Q_or_QC &quadrature,
1174 const CellSimilarity::Similarity cell_similarity,
1175 const typename Mapping<dim, spacedim>::InternalDataBase &mapping_internal,
1176 const typename FiniteElement<dim, spacedim>::InternalDataBase &fe_data,
1178 &mapping_data,
1180 &output_data) const;
1181
1182private:
1187
1197 std::vector<Table<2, unsigned int>> primitive_offset_tables;
1198
1202 std::vector<std::vector<BaseOffsets>> nonprimitive_offset_tables;
1203
1211 std::vector<std::pair<std::unique_ptr<const FiniteElement<dim, spacedim>>,
1212 unsigned int>>
1214
1227 std::vector<std::vector<std::size_t>> generalized_support_points_index_table;
1228
1234 void
1235 initialize(const std::vector<const FiniteElement<dim, spacedim> *> &fes,
1236 const std::vector<unsigned int> &multiplicities);
1237
1241 void
1243
1249 template <int structdim>
1250 std::vector<std::pair<unsigned int, unsigned int>>
1252 const unsigned int face_no = 0) const;
1253
1260 class InternalData : public FiniteElement<dim, spacedim>::InternalDataBase
1261 {
1262 public:
1267 InternalData(const unsigned int n_base_elements);
1268
1273 ~InternalData() override;
1274
1279 void
1281 const unsigned int base_no,
1282 std::unique_ptr<typename FiniteElement<dim, spacedim>::InternalDataBase>);
1283
1289 get_fe_data(const unsigned int base_no) const;
1290
1297 get_fe_output_object(const unsigned int base_no) const;
1298
1299 private:
1310 typename std::vector<
1311 std::unique_ptr<typename FiniteElement<dim, spacedim>::InternalDataBase>>
1313
1322 mutable std::vector<
1325 };
1326
1331 mutable std::mutex restriction_matrix_mutex;
1332 mutable std::mutex prolongation_matrix_mutex;
1333
1334 friend class FE_Enriched<dim, spacedim>;
1335};
1336
1337//------------------------variadic template constructor------------------------
1338
1339#ifndef DOXYGEN
1340namespace internal
1341{
1342 namespace FESystemImplementation
1343 {
1344 template <int dim, int spacedim>
1345 unsigned int
1346 count_nonzeros(
1347 const std::initializer_list<
1348 std::pair<std::unique_ptr<FiniteElement<dim, spacedim>>, unsigned int>>
1349 &fe_systems)
1350 {
1351 return std::count_if(
1352 fe_systems.begin(),
1353 fe_systems.end(),
1354 [](const std::pair<std::unique_ptr<FiniteElement<dim, spacedim>>,
1355 unsigned int> &fe_system) {
1356 return fe_system.second > 0;
1357 });
1358 }
1359
1360
1361
1362 template <int dim, int spacedim>
1363 std::pair<std::unique_ptr<FiniteElement<dim, spacedim>>, unsigned int>
1364 promote_to_fe_pair(const FiniteElement<dim, spacedim> &fe)
1365 {
1366 return std::make_pair(std::move(fe.clone()), 1u);
1367 }
1368
1369
1370
1371 template <int dim, int spacedim>
1372 auto
1373 promote_to_fe_pair(std::pair<std::unique_ptr<FiniteElement<dim, spacedim>>,
1374 unsigned int> &&p)
1375 -> decltype(std::forward<
1376 std::pair<std::unique_ptr<FiniteElement<dim, spacedim>>,
1377 unsigned int>>(p))
1378 {
1379 return std::forward<
1380 std::pair<std::unique_ptr<FiniteElement<dim, spacedim>>, unsigned int>>(
1381 p);
1382 }
1383 } // namespace FESystemImplementation
1384} // namespace internal
1385
1386
1387
1388# if !defined(__INTEL_COMPILER) || __INTEL_COMPILER >= 1900
1389// We are just forwarding/delegating to the constructor taking a
1390// std::initializer_list. If we decide to remove the deprecated constructors, we
1391// might just use the variadic constructor with a suitable static_assert instead
1392// of the std::enable_if.
1393template <int dim, int spacedim>
1394template <class... FEPairs, typename>
1395FESystem<dim, spacedim>::FESystem(FEPairs &&...fe_pairs)
1396 : FESystem<dim, spacedim>(
1397 {internal::FESystemImplementation::promote_to_fe_pair<dim, spacedim>(
1398 std::forward<FEPairs>(fe_pairs))...})
1399{}
1400
1401
1402
1403template <int dim, int spacedim>
1405 const std::initializer_list<
1406 std::pair<std::unique_ptr<FiniteElement<dim, spacedim>>, unsigned int>>
1407 &fe_systems)
1408 : FiniteElement<dim, spacedim>(
1409 FETools::Compositing::multiply_dof_numbers<dim, spacedim>(fe_systems),
1411 spacedim>(
1412 fe_systems),
1413 FETools::Compositing::compute_nonzero_components<dim, spacedim>(
1414 fe_systems))
1415 , base_elements(internal::FESystemImplementation::count_nonzeros(fe_systems))
1416{
1417 std::vector<const FiniteElement<dim, spacedim> *> fes;
1418 std::vector<unsigned int> multiplicities;
1419
1420 const auto extract =
1421 [&fes, &multiplicities](
1422 const std::pair<std::unique_ptr<FiniteElement<dim, spacedim>>,
1423 unsigned int> &fe_system) {
1424 fes.push_back(fe_system.first.get());
1425 multiplicities.push_back(fe_system.second);
1426 };
1427
1428 for (const auto &p : fe_systems)
1429 extract(p);
1430
1431 initialize(fes, multiplicities);
1432}
1433# endif
1434#endif // DOXYGEN
1435
1437
1438#endif
void set_fe_data(const unsigned int base_no, std::unique_ptr< typename FiniteElement< dim, spacedim >::InternalDataBase >)
InternalData(const unsigned int n_base_elements)
std::vector< std::unique_ptr< typename FiniteElement< dim, spacedim >::InternalDataBase > > base_fe_datas
Definition fe_system.h:1312
std::vector< internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim > > base_fe_output_objects
Definition fe_system.h:1324
FiniteElement< dim, spacedim >::InternalDataBase & get_fe_data(const unsigned int base_no) const
internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim > & get_fe_output_object(const unsigned int base_no) const
std::mutex restriction_matrix_mutex
Definition fe_system.h:1331
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
FESystem(const std::initializer_list< std::pair< std::unique_ptr< FiniteElement< dim, spacedim > >, unsigned int > > &fe_systems)
virtual Tensor< 2, dim > shape_grad_grad(const unsigned int i, const Point< dim > &p) const override
FESystem(const FESystem< dim, spacedim > &)=delete
virtual const FiniteElement< dim, spacedim > & base_element(const unsigned int index) const override
FESystem(FEPairs &&...fe_pairs)
FESystem()=delete
virtual Point< dim > unit_support_point(const unsigned int index) const override
std::vector< std::vector< BaseOffsets > > nonprimitive_offset_tables
Definition fe_system.h:1202
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
std::vector< std::pair< std::unique_ptr< const FiniteElement< dim, spacedim > >, unsigned int > > base_elements
Definition fe_system.h:1213
virtual Tensor< 3, dim > shape_3rd_derivative(const unsigned int i, const Point< dim > &p) const override
std::vector< Table< 2, unsigned int > > primitive_offset_tables
Definition fe_system.h:1197
virtual std::string get_name() const override
virtual const FullMatrix< double > & get_prolongation_matrix(const unsigned int child, const RefinementCase< dim > &refinement_case=RefinementCase< dim >::isotropic_refinement) const override
virtual const FiniteElement< dim, spacedim > & get_sub_fe(const unsigned int first_component, const unsigned int n_selected_components) const override
FESystem(const FiniteElement< dim, spacedim > &fe1, const unsigned int n1, const FiniteElement< dim, spacedim > &fe2, const unsigned int n2, const FiniteElement< dim, spacedim > &fe3, const unsigned int n3)
virtual Tensor< 1, dim > shape_grad_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const override
virtual UpdateFlags requires_update_flags(const UpdateFlags update_flags) const override
virtual FiniteElementDomination::Domination compare_for_domination(const FiniteElement< dim, spacedim > &fe_other, const unsigned int codim=0) const override final
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 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
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 Q_or_QC &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 Tensor< 3, dim > shape_3rd_derivative_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const override
virtual const FullMatrix< double > & get_restriction_matrix(const unsigned int child, const RefinementCase< dim > &refinement_case=RefinementCase< dim >::isotropic_refinement) const override
virtual bool has_support_on_face(const unsigned int shape_index, const unsigned int face_index) const override
virtual void get_interpolation_matrix(const FiniteElement< dim, spacedim > &source, FullMatrix< double > &matrix) const override
virtual Tensor< 1, dim > shape_grad(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
FESystem(const std::vector< const FiniteElement< dim, spacedim > * > &fes, const std::vector< unsigned int > &multiplicities)
virtual std::vector< std::pair< unsigned int, unsigned int > > hp_line_dof_identities(const FiniteElement< dim, spacedim > &fe_other) 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
static const unsigned int invalid_face_number
Definition fe_system.h:1186
virtual bool hp_constraints_are_implemented() const override
virtual std::vector< std::pair< unsigned int, unsigned int > > hp_vertex_dof_identities(const FiniteElement< dim, spacedim > &fe_other) const override
FESystem(const FiniteElement< dim, spacedim > &fe1, const unsigned int n1, const FiniteElement< dim, spacedim > &fe2, const unsigned int n2, const FiniteElement< dim, spacedim > &fe3, const unsigned int n3, const FiniteElement< dim, spacedim > &fe4, const unsigned int n4, const FiniteElement< dim, spacedim > &fe5, const unsigned int n5)
void build_interface_constraints()
virtual Tensor< 4, dim > shape_4th_derivative_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const override
virtual unsigned int face_to_cell_index(const unsigned int face_dof_index, const unsigned int face, const unsigned char combined_orientation=ReferenceCell::default_combined_face_orientation()) const override
virtual Point< dim - 1 > unit_face_support_point(const unsigned int index, const unsigned int face_no=0) const override
void initialize(const std::vector< const FiniteElement< dim, spacedim > * > &fes, const std::vector< unsigned int > &multiplicities)
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
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 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
virtual double shape_value_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const override
virtual double shape_value(const unsigned int i, const Point< dim > &p) const override
FESystem(const FiniteElement< dim, spacedim > &fe1, const unsigned int n1, const FiniteElement< dim, spacedim > &fe2, const unsigned int n2, const FiniteElement< dim, spacedim > &fe3, const unsigned int n3, const FiniteElement< dim, spacedim > &fe4, const unsigned int n4)
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< 2, dim > shape_grad_grad_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const override
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)
std::vector< std::vector< std::size_t > > generalized_support_points_index_table
Definition fe_system.h:1227
virtual std::unique_ptr< FiniteElement< dim, spacedim > > clone() const override
FESystem(FESystem< dim, spacedim > &&other_fe_system) noexcept
Definition fe_system.h:573
std::mutex prolongation_matrix_mutex
Definition fe_system.h:1332
virtual std::size_t memory_consumption() const override
virtual ~FESystem() override=default
FESystem(const FiniteElement< dim, spacedim > &fe1, const unsigned int n1, const FiniteElement< dim, spacedim > &fe2, const unsigned int n2)
virtual Tensor< 4, dim > shape_4th_derivative(const unsigned int i, const Point< dim > &p) const override
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 std::unique_ptr< FiniteElement< dim, spacedim > > clone() const =0
unsigned int n_base_elements() const
Abstract base class for mapping classes.
Definition mapping.h:318
Definition point.h:111
static constexpr unsigned char default_combined_face_orientation()
#define DEAL_II_NAMESPACE_OPEN
Definition config.h:502
#define DEAL_II_NAMESPACE_CLOSE
Definition config.h:503
UpdateFlags
std::vector< ComponentMask > compute_nonzero_components(const std::vector< const FiniteElement< dim, spacedim > * > &fes, const std::vector< unsigned int > &multiplicities, const bool do_tensor_product=true)
std::vector< bool > compute_restriction_is_additive_flags(const std::vector< const FiniteElement< dim, spacedim > * > &fes, const std::vector< unsigned int > &multiplicities)
FiniteElementData< dim > multiply_dof_numbers(const std::vector< const FiniteElement< dim, spacedim > * > &fes, const std::vector< unsigned int > &multiplicities, const bool do_tensor_product=true)
constexpr ReturnType< rank, T >::value_type & extract(T &t, const ArrayType &indices)
static const unsigned int invalid_unsigned_int
Definition types.h:220
STL namespace.
unsigned int n_nonzero_components
Definition fe_system.h:1061
unsigned int out_index
Definition fe_system.h:1075
unsigned int in_index
Definition fe_system.h:1068
typename enable_if_all< Values... >::type enable_if_all_t