46template <
int dim,
int spacedim>
67template <
int dim,
int spacedim>
75 this->update_each = update_flags;
77 const unsigned int n_q_points = q.
size();
84 if (this->update_each &
86 compute_manifold_quadrature_weights(q);
89 covariant.resize(n_q_points);
92 contravariant.resize(n_q_points);
95 volume_elements.resize(n_q_points);
100template <
int dim,
int spacedim>
105 const unsigned int n_original_q_points)
107 reinit(update_flags, q);
112 covariant.resize(n_original_q_points);
115 contravariant.resize(n_original_q_points);
118 volume_elements.resize(n_original_q_points);
130 unit_tangentials[i].resize(n_original_q_points);
131 std::fill(unit_tangentials[i].begin(),
132 unit_tangentials[i].end(),
137 .resize(n_original_q_points);
152template <
int dim,
int spacedim>
158 vertices[i] = cell->vertex(i);
164template <
int dim,
int spacedim>
169 cell_manifold_quadrature_weights.resize(quad.
size());
170 for (
unsigned int q = 0; q < quad.
size(); ++q)
174 cell_manifold_quadrature_weights[q][i] =
182template <
int dim,
int spacedim>
189template <
int dim,
int spacedim>
190std::unique_ptr<Mapping<dim, spacedim>>
193 return std::make_unique<MappingManifold<dim, spacedim>>(*this);
198template <
int dim,
int spacedim>
210template <
int dim,
int spacedim>
217 std::array<double, GeometryInfo<dim>::vertices_per_cell> weights;
221 vertices[v] = cell->vertex(v);
224 return cell->get_manifold().get_new_point(
244template <
int dim,
int spacedim>
256 for (
unsigned int i = 0; i < 5; ++i)
306template <
int dim,
int spacedim>
307std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
311 std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase> data_ptr =
312 std::make_unique<InternalData>();
320template <
int dim,
int spacedim>
321std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
328 std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase> data_ptr =
329 std::make_unique<InternalData>();
333 ReferenceCells::get_hypercube<dim>(), quadrature[0]),
334 quadrature[0].
size());
341template <
int dim,
int spacedim>
342std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
347 std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase> data_ptr =
348 std::make_unique<InternalData>();
352 ReferenceCells::get_hypercube<dim>(), quadrature),
362 namespace MappingManifoldImplementation
372 template <
int dim,
int spacedim>
374 maybe_compute_q_points(
376 const typename ::MappingManifold<dim, spacedim>::InternalData
384 for (
unsigned int point = 0; point < quadrature_points.size();
387 quadrature_points[point] =
data.manifold->get_new_point(
390 data.cell_manifold_quadrature_weights[point + data_set]));
403 template <
int dim,
int spacedim>
405 maybe_update_Jacobians(
406 const typename ::QProjector<dim>::DataSetDescriptor data_set,
407 const typename ::MappingManifold<dim, spacedim>::InternalData
414 const unsigned int n_q_points =
data.contravariant.size();
416 std::fill(
data.contravariant.begin(),
417 data.contravariant.end(),
420 for (
unsigned int point = 0; point < n_q_points; ++point)
430 data.cell_manifold_quadrature_weights[point + data_set]));
439 for (
unsigned int i = 0; i < dim; ++i)
442 const double pi = p[i];
443 Assert(pi >= 0 && pi <= 1.0,
445 "Was expecting a quadrature point "
446 "inside the unit reference element."));
450 const double L = pi > .5 ? -pi : 1 - pi;
455 for (
const unsigned int j :
457 data.vertex_weights[j] =
465 data.manifold->get_tangent_vector(P, NP);
467 for (
unsigned int d = 0; d < spacedim; ++d)
468 data.contravariant[point][d][i] = T[d] / L;
474 const unsigned int n_q_points =
data.contravariant.size();
475 for (
unsigned int point = 0; point < n_q_points; ++point)
477 data.covariant[point] =
478 (
data.contravariant[point]).covariant_form();
484 const unsigned int n_q_points =
data.contravariant.size();
485 for (
unsigned int point = 0; point < n_q_points; ++point)
486 data.volume_elements[point] =
487 data.contravariant[point].determinant();
497template <
int dim,
int spacedim>
512 const unsigned int n_q_points = quadrature.
size();
514 data.store_vertices(cell);
515 data.manifold = &(cell->get_manifold());
517 internal::MappingManifoldImplementation::maybe_compute_q_points<dim,
523 internal::MappingManifoldImplementation::maybe_update_Jacobians<dim,
528 const std::vector<double> &weights = quadrature.
get_weights();
543 for (
unsigned int point = 0; point < n_q_points; ++point)
547 const double det =
data.contravariant[point].determinant();
554 Assert(det > 1e-12 * Utilities::fixed_power<dim>(
555 cell->diameter() /
std::sqrt(
double(dim))),
557 cell->center(), det, point)));
559 output_data.
JxW_values[point] = weights[point] * det;
567 for (
unsigned int i = 0; i < spacedim; ++i)
568 for (
unsigned int j = 0; j < dim; ++j)
569 DX_t[j][i] =
data.contravariant[point][i][j];
572 for (
unsigned int i = 0; i < dim; ++i)
573 for (
unsigned int j = 0; j < dim; ++j)
574 G[i][j] = DX_t[i] * DX_t[j];
581 Assert(spacedim == dim + 1,
583 "There is no (unique) cell normal for " +
585 "-dimensional cells in " +
587 "-dimensional space. This only works if the "
588 "space dimension is one greater than the "
589 "dimensionality of the mesh cells."));
593 cross_product_2d(-DX_t[0]);
596 cross_product_3d(DX_t[0], DX_t[1]);
601 if (cell->direction_flag() ==
false)
614 for (
unsigned int point = 0; point < n_q_points; ++point)
622 for (
unsigned int point = 0; point < n_q_points; ++point)
624 data.covariant[point].transpose();
634 namespace MappingManifoldImplementation
647 template <
int dim,
int spacedim>
649 maybe_compute_face_data(
650 const ::MappingManifold<dim, spacedim> &mapping,
651 const typename ::Triangulation<dim, spacedim>::cell_iterator
653 const unsigned int face_no,
654 const unsigned int subface_no,
655 const unsigned int n_q_points,
656 const std::vector<double> &weights,
657 const typename ::MappingManifold<dim, spacedim>::InternalData
675 for (
unsigned int d = 0; d != dim - 1; ++d)
678 data.unit_tangentials.size(),
681 data.aux[d].size() <=
683 .unit_tangentials[face_no +
691 .unit_tangentials[face_no +
702 for (
unsigned int i = 0; i < n_q_points; ++i)
711 (face_no == 0 ? -1 : +1);
715 cross_product_2d(
data.aux[0][i]);
719 cross_product_3d(
data.aux[0][i],
data.aux[1][i]);
735 for (
unsigned int point = 0;
point < n_q_points; ++
point)
743 data.contravariant[
point].transpose()[0];
745 (face_no == 0 ? -1. : +1.) *
757 cross_product_3d(DX_t[0], DX_t[1]);
758 cell_normal /= cell_normal.
norm();
763 cross_product_3d(
data.aux[0][point], cell_normal);
785 const double area_ratio =
787 cell->subface_case(face_no), subface_no);
799 for (
unsigned int point = 0;
point < n_q_points; ++
point)
803 for (
unsigned int point = 0;
point < n_q_points; ++
point)
805 data.covariant[point].transpose();
816 template <
int dim,
int spacedim>
819 const ::MappingManifold<dim, spacedim> &mapping,
820 const typename ::Triangulation<dim, spacedim>::cell_iterator
822 const unsigned int face_no,
823 const unsigned int subface_no,
826 const typename ::MappingManifold<dim, spacedim>::InternalData
831 data.store_vertices(cell);
833 data.manifold = &cell->face(face_no)->get_manifold();
835 maybe_compute_q_points<dim, spacedim>(data_set,
838 maybe_update_Jacobians<dim, spacedim>(data_set,
data);
850 template <
int dim,
int spacedim,
int rank>
859 Assert((
dynamic_cast<const typename ::
860 MappingManifold<dim, spacedim>::InternalData *
>(
861 &mapping_data) !=
nullptr),
863 const typename ::MappingManifold<dim, spacedim>::InternalData
865 static_cast<const typename ::MappingManifold<dim, spacedim>::
866 InternalData &
>(mapping_data);
868 switch (mapping_kind)
875 "update_contravariant_transformation"));
877 for (
unsigned int i = 0; i < output.size(); ++i)
889 "update_contravariant_transformation"));
893 "update_volume_elements"));
898 for (
unsigned int i = 0; i < output.size(); ++i)
902 output[i] /=
data.volume_elements[i];
914 "update_covariant_transformation"));
916 for (
unsigned int i = 0; i < output.size(); ++i)
928 template <
int dim,
int spacedim,
int rank>
937 Assert((
dynamic_cast<const typename ::
938 MappingManifold<dim, spacedim>::InternalData *
>(
939 &mapping_data) !=
nullptr),
941 const typename ::MappingManifold<dim, spacedim>::InternalData
943 static_cast<const typename ::MappingManifold<dim, spacedim>::
944 InternalData &
>(mapping_data);
946 switch (mapping_kind)
953 "update_covariant_transformation"));
957 "update_contravariant_transformation"));
960 for (
unsigned int i = 0; i < output.size(); ++i)
977 "update_covariant_transformation"));
980 for (
unsigned int i = 0; i < output.size(); ++i)
997 "update_covariant_transformation"));
1001 "update_contravariant_transformation"));
1005 "update_volume_elements"));
1008 for (
unsigned int i = 0; i < output.size(); ++i)
1017 output[i] /=
data.volume_elements[i];
1030 template <
int dim,
int spacedim>
1039 Assert((
dynamic_cast<const typename ::
1040 MappingManifold<dim, spacedim>::InternalData *
>(
1041 &mapping_data) !=
nullptr),
1043 const typename ::MappingManifold<dim, spacedim>::InternalData
1045 static_cast<const typename ::MappingManifold<dim, spacedim>::
1046 InternalData &
>(mapping_data);
1048 switch (mapping_kind)
1055 "update_covariant_transformation"));
1059 "update_contravariant_transformation"));
1061 for (
unsigned int q = 0; q < output.size(); ++q)
1062 for (
unsigned int i = 0; i < spacedim; ++i)
1064 double tmp1[dim][dim];
1065 for (
unsigned int J = 0; J < dim; ++J)
1066 for (
unsigned int K = 0;
K < dim; ++
K)
1069 data.contravariant[q][i][0] * input[q][0][J][
K];
1070 for (
unsigned int I = 1; I < dim; ++I)
1072 data.contravariant[q][i][I] * input[q][I][J][K];
1074 for (
unsigned int j = 0; j < spacedim; ++j)
1077 for (
unsigned int K = 0;
K < dim; ++
K)
1079 tmp2[
K] =
data.covariant[q][j][0] * tmp1[0][
K];
1080 for (
unsigned int J = 1; J < dim; ++J)
1081 tmp2[K] +=
data.covariant[q][j][J] * tmp1[J][K];
1083 for (
unsigned int k = 0; k < spacedim; ++k)
1085 output[q][i][j][k] =
1086 data.covariant[q][k][0] * tmp2[0];
1087 for (
unsigned int K = 1;
K < dim; ++
K)
1088 output[q][i][j][k] +=
1089 data.covariant[q][k][K] * tmp2[K];
1101 "update_covariant_transformation"));
1103 for (
unsigned int q = 0; q < output.size(); ++q)
1104 for (
unsigned int i = 0; i < spacedim; ++i)
1106 double tmp1[dim][dim];
1107 for (
unsigned int J = 0; J < dim; ++J)
1108 for (
unsigned int K = 0;
K < dim; ++
K)
1111 data.covariant[q][i][0] * input[q][0][J][
K];
1112 for (
unsigned int I = 1; I < dim; ++I)
1114 data.covariant[q][i][I] * input[q][I][J][K];
1116 for (
unsigned int j = 0; j < spacedim; ++j)
1119 for (
unsigned int K = 0;
K < dim; ++
K)
1121 tmp2[
K] =
data.covariant[q][j][0] * tmp1[0][
K];
1122 for (
unsigned int J = 1; J < dim; ++J)
1123 tmp2[K] +=
data.covariant[q][j][J] * tmp1[J][K];
1125 for (
unsigned int k = 0; k < spacedim; ++k)
1127 output[q][i][j][k] =
1128 data.covariant[q][k][0] * tmp2[0];
1129 for (
unsigned int K = 1;
K < dim; ++
K)
1130 output[q][i][j][k] +=
1131 data.covariant[q][k][K] * tmp2[K];
1144 "update_covariant_transformation"));
1148 "update_contravariant_transformation"));
1152 "update_volume_elements"));
1154 for (
unsigned int q = 0; q < output.size(); ++q)
1155 for (
unsigned int i = 0; i < spacedim; ++i)
1158 for (
unsigned int I = 0; I < dim; ++I)
1160 data.contravariant[q][i][I] /
data.volume_elements[q];
1161 double tmp1[dim][dim];
1162 for (
unsigned int J = 0; J < dim; ++J)
1163 for (
unsigned int K = 0;
K < dim; ++
K)
1165 tmp1[J][
K] = factor[0] * input[q][0][J][
K];
1166 for (
unsigned int I = 1; I < dim; ++I)
1167 tmp1[J][K] += factor[I] * input[q][I][J][K];
1169 for (
unsigned int j = 0; j < spacedim; ++j)
1172 for (
unsigned int K = 0;
K < dim; ++
K)
1174 tmp2[
K] =
data.covariant[q][j][0] * tmp1[0][
K];
1175 for (
unsigned int J = 1; J < dim; ++J)
1176 tmp2[K] +=
data.covariant[q][j][J] * tmp1[J][K];
1178 for (
unsigned int k = 0; k < spacedim; ++k)
1180 output[q][i][j][k] =
1181 data.covariant[q][k][0] * tmp2[0];
1182 for (
unsigned int K = 1;
K < dim; ++
K)
1183 output[q][i][j][k] +=
1184 data.covariant[q][k][K] * tmp2[K];
1199 template <
int dim,
int spacedim,
int rank>
1208 Assert((
dynamic_cast<const typename ::
1209 MappingManifold<dim, spacedim>::InternalData *
>(
1210 &mapping_data) !=
nullptr),
1212 const typename ::MappingManifold<dim, spacedim>::InternalData
1214 static_cast<const typename ::MappingManifold<dim, spacedim>::
1215 InternalData &
>(mapping_data);
1217 switch (mapping_kind)
1224 "update_covariant_transformation"));
1226 for (
unsigned int i = 0; i < output.size(); ++i)
1241template <
int dim,
int spacedim>
1245 const unsigned int face_no,
1258 internal::MappingManifoldImplementation::do_fill_fe_face_values(
1264 ReferenceCells::get_hypercube<dim>(),
1266 cell->combined_face_orientation(face_no),
1267 quadrature[0].
size()),
1275template <
int dim,
int spacedim>
1279 const unsigned int face_no,
1280 const unsigned int subface_no,
1291 internal::MappingManifoldImplementation::do_fill_fe_face_values(
1297 ReferenceCells::get_hypercube<dim>(),
1300 cell->combined_face_orientation(face_no),
1302 cell->subface_case(face_no)),
1310template <
int dim,
int spacedim>
1318 internal::MappingManifoldImplementation::transform_fields(input,
1326template <
int dim,
int spacedim>
1334 internal::MappingManifoldImplementation::transform_differential_forms(
1335 input, mapping_kind, mapping_data, output);
1340template <
int dim,
int spacedim>
1348 switch (mapping_kind)
1351 internal::MappingManifoldImplementation::transform_fields(input,
1360 internal::MappingManifoldImplementation::transform_gradients(
1361 input, mapping_kind, mapping_data, output);
1370template <
int dim,
int spacedim>
1383 switch (mapping_kind)
1389 "update_covariant_transformation"));
1391 for (
unsigned int q = 0; q < output.size(); ++q)
1392 for (
unsigned int i = 0; i < spacedim; ++i)
1393 for (
unsigned int j = 0; j < spacedim; ++j)
1396 for (
unsigned int K = 0; K < dim; ++K)
1398 tmp[K] =
data.covariant[q][j][0] * input[q][i][0][K];
1399 for (
unsigned int J = 1; J < dim; ++J)
1400 tmp[K] +=
data.covariant[q][j][J] * input[q][i][J][K];
1402 for (
unsigned int k = 0; k < spacedim; ++k)
1404 output[q][i][j][k] =
data.covariant[q][k][0] * tmp[0];
1405 for (
unsigned int K = 1; K < dim; ++K)
1406 output[q][i][j][k] +=
data.covariant[q][k][K] * tmp[K];
1419template <
int dim,
int spacedim>
1427 switch (mapping_kind)
1432 internal::MappingManifoldImplementation::transform_hessians(
1433 input, mapping_kind, mapping_data, output);
1441#include "mapping_manifold.inst"
ArrayView< std::remove_reference_t< typename std::iterator_traits< Iterator >::reference >, MemorySpaceType > make_array_view(const Iterator begin, const Iterator end)
void store_vertices(const typename Triangulation< dim, spacedim >::cell_iterator &cell) const
std::vector< double > volume_elements
std::array< std::vector< Tensor< 1, dim > >, GeometryInfo< dim >::faces_per_cell *(dim - 1)> unit_tangentials
std::vector< DerivativeForm< 1, dim, spacedim > > contravariant
std::array< Point< spacedim >, GeometryInfo< dim >::vertices_per_cell > vertices
std::vector< std::array< double, GeometryInfo< dim >::vertices_per_cell > > cell_manifold_quadrature_weights
std::vector< std::vector< Tensor< 1, spacedim > > > aux
virtual std::size_t memory_consumption() const override
void initialize_face(const UpdateFlags update_flags, const Quadrature< dim > &quadrature, const unsigned int n_original_q_points)
virtual void reinit(const UpdateFlags update_flags, const Quadrature< dim > &quadrature) override
ObserverPointer< const Manifold< dim, spacedim > > manifold
std::array< double, GeometryInfo< dim >::vertices_per_cell > vertex_weights
void compute_manifold_quadrature_weights(const Quadrature< dim > &quadrature)
std::vector< DerivativeForm< 1, dim, spacedim > > covariant
Triangulation< dim, spacedim >::cell_iterator cell
virtual Point< spacedim > transform_unit_to_real_cell(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const Point< dim > &p) const override
MappingManifold()=default
virtual Point< dim > transform_real_to_unit_cell(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const Point< spacedim > &p) const override
virtual void transform(const ArrayView< const Tensor< 1, dim > > &input, const MappingKind kind, const typename Mapping< dim, spacedim >::InternalDataBase &internal, const ArrayView< Tensor< 1, spacedim > > &output) const override
virtual void fill_fe_subface_values(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const unsigned int face_no, const unsigned int subface_no, const Quadrature< dim - 1 > &quadrature, const typename Mapping< dim, spacedim >::InternalDataBase &internal_data, internal::FEValuesImplementation::MappingRelatedData< dim, spacedim > &output_data) 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 typename Mapping< dim, spacedim >::InternalDataBase &internal_data, internal::FEValuesImplementation::MappingRelatedData< dim, spacedim > &output_data) const override
virtual UpdateFlags requires_update_flags(const UpdateFlags update_flags) const override
virtual std::unique_ptr< typename Mapping< dim, spacedim >::InternalDataBase > get_data(const UpdateFlags, const Quadrature< dim > &quadrature) const override
virtual std::unique_ptr< Mapping< dim, spacedim > > clone() const override
virtual std::unique_ptr< typename Mapping< dim, spacedim >::InternalDataBase > get_face_data(const UpdateFlags flags, const hp::QCollection< dim - 1 > &quadrature) const override
virtual std::unique_ptr< typename Mapping< dim, spacedim >::InternalDataBase > get_subface_data(const UpdateFlags flags, const Quadrature< dim - 1 > &quadrature) const override
virtual CellSimilarity::Similarity fill_fe_values(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const CellSimilarity::Similarity cell_similarity, const Quadrature< dim > &quadrature, const typename Mapping< dim, spacedim >::InternalDataBase &internal_data, internal::FEValuesImplementation::MappingRelatedData< dim, spacedim > &output_data) const override
Abstract base class for mapping classes.
static constexpr Point< dim, Number > unit_vector(const unsigned int i)
static DataSetDescriptor cell()
const Point< dim > & point(const unsigned int i) const
const std::vector< double > & get_weights() const
const std::vector< Point< dim > > & get_points() const
unsigned int size() const
numbers::NumberTraits< Number >::real_type norm() const
unsigned int size() const
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
static ::ExceptionBase & ExcNotImplemented()
#define Assert(cond, exc)
#define AssertDimension(dim1, dim2)
static ::ExceptionBase & ExcInternalError()
static ::ExceptionBase & ExcDimensionMismatch(std::size_t arg1, std::size_t arg2)
static ::ExceptionBase & ExcMessage(std::string arg1)
@ update_jacobian_pushed_forward_2nd_derivatives
@ update_volume_elements
Determinant of the Jacobian.
@ update_contravariant_transformation
Contravariant transformation.
@ update_jacobian_pushed_forward_grads
@ update_jacobian_grads
Gradient of volume element.
@ update_normal_vectors
Normal vectors.
@ update_JxW_values
Transformed quadrature weights.
@ update_covariant_transformation
Covariant transformation.
@ update_jacobians
Volume element.
@ update_inverse_jacobians
Volume element.
@ update_quadrature_points
Transformed quadrature points.
@ update_jacobian_pushed_forward_3rd_derivatives
@ update_boundary_forms
Outer normal vector, not normalized.
@ mapping_covariant_gradient
@ mapping_contravariant_hessian
@ mapping_covariant_hessian
@ mapping_contravariant_gradient
#define DEAL_II_NOT_IMPLEMENTED()
std::vector< index_type > data
CGAL::Exact_predicates_exact_constructions_kernel_with_sqrt K
std::enable_if_t< std::is_fundamental_v< T >, std::size_t > memory_consumption(const T &t)
Point< spacedim > point(const gp_Pnt &p, const double tolerance=1e-10)
std::string int_to_string(const unsigned int value, const unsigned int digits=numbers::invalid_unsigned_int)
void transform_differential_forms(const ArrayView< const DerivativeForm< rank, dim, spacedim > > &input, const MappingKind mapping_kind, const typename Mapping< dim, spacedim >::InternalDataBase &mapping_data, const ArrayView< Tensor< rank+1, spacedim > > &output)
void do_fill_fe_face_values(const ::MappingQ< dim, spacedim > &mapping, const typename ::Triangulation< dim, spacedim >::cell_iterator &cell, const unsigned int face_no, const unsigned int subface_no, const typename QProjector< dim >::DataSetDescriptor data_set, const Quadrature< dim - 1 > &quadrature, const typename ::MappingQ< dim, spacedim >::InternalData &data, const std::vector< Polynomials::Polynomial< double > > &polynomials_1d, const std::vector< unsigned int > &renumber_lexicographic_to_hierarchic, internal::FEValuesImplementation::MappingRelatedData< dim, spacedim > &output_data)
void transform_fields(const ArrayView< const Tensor< rank, dim > > &input, const MappingKind mapping_kind, const typename Mapping< dim, spacedim >::InternalDataBase &mapping_data, const ArrayView< Tensor< rank, spacedim > > &output)
void transform_gradients(const ArrayView< const Tensor< rank, dim > > &input, const MappingKind mapping_kind, const typename Mapping< dim, spacedim >::InternalDataBase &mapping_data, const ArrayView< Tensor< rank, spacedim > > &output)
void transform_hessians(const ArrayView< const Tensor< 3, dim > > &input, const MappingKind mapping_kind, const typename Mapping< dim, spacedim >::InternalDataBase &mapping_data, const ArrayView< Tensor< 3, spacedim > > &output)
void maybe_compute_face_data(const ::MappingQ< dim, spacedim > &mapping, const typename ::Triangulation< dim, spacedim >::cell_iterator &cell, const unsigned int face_no, const unsigned int subface_no, const unsigned int n_q_points, const std::vector< double > &weights, const typename ::MappingQ< dim, spacedim >::InternalData &data, internal::FEValuesImplementation::MappingRelatedData< dim, spacedim > &output_data)
static const unsigned int invalid_unsigned_int
::VectorizedArray< Number, width > sqrt(const ::VectorizedArray< Number, width > &)
static std_cxx20::ranges::iota_view< unsigned int, unsigned int > face_indices()
static double subface_ratio(const internal::SubfaceCase< dim > &subface_case, const unsigned int subface_no)
static double d_linear_shape_function(const Point< dim > &xi, const unsigned int i)
static std_cxx20::ranges::iota_view< unsigned int, unsigned int > vertex_indices()
DEAL_II_HOST constexpr Number determinant(const SymmetricTensor< 2, dim, Number > &)