16#ifndef dealii_matrix_free_fe_evaluation_data_h
17#define dealii_matrix_free_fe_evaluation_data_h
49 <<
"You are requesting information from an FEEvaluation/FEFaceEvaluation "
50 <<
"object for which this kind of information has not been computed. What "
51 <<
"information these objects compute is determined by the update_* flags "
52 <<
"you pass to MatrixFree::reinit() via MatrixFree::AdditionalData. "
53 <<
"Here, the operation you are attempting requires the <" << arg1
54 <<
"> flag to be set, but it was apparently not specified "
55 <<
"upon initialization.");
61 int n_q_points_1d = fe_degree + 1,
62 int n_components_ = 1,
63 typename Number = double,
76 namespace MatrixFreeFunctions
78 template <
int,
typename>
79 class MappingDataOnTheFly;
112template <
int dim,
typename Number,
bool is_face>
117 MappingInfoStorage<(is_face ? dim - 1 : dim), dim, Number>;
164 const unsigned int n_components);
193 JxW(
const unsigned int q_point)
const;
236 DEAL_II_DEPRECATED_EARLY_WITH_COMMENT(
"Use normal_vector() instead.")
376 internal::MatrixFreeFunctions::GeometryType
396 const
std::vector<
unsigned int> &
508 const std::array<unsigned int, n_lanes> &
537 const std::array<unsigned int, n_lanes> &
577 template <
typename T>
587 template <
typename T>
599 template <
typename T>
611 template <
typename T>
648 const std::shared_ptr<
999 template <
int,
int,
typename,
bool,
typename>
1002 template <
int,
int,
int,
int,
typename,
typename>
1023template <
int dim,
typename Number,
bool is_face>
1028 InitializationData{&shape_info, nullptr, nullptr, 0, 0, nullptr},
1035template <
int dim,
typename Number,
bool is_face>
1037 const InitializationData &initialization_data,
1041 :
data(initialization_data.shape_info)
1086template <
int dim,
typename Number,
bool is_face>
1088 const std::shared_ptr<
1126 .jacobian_gradients_non_inverse[0]
1134template <
int dim,
typename Number,
bool is_face>
1187template <
int dim,
typename Number,
bool is_face>
1195 const unsigned int tensor_dofs_per_component =
1196 Utilities::fixed_power<dim>(
data->
data.front().fe_degree + 1);
1199 const unsigned int size_scratch_data =
1203 const unsigned int size_data_arrays =
1205 (
n_components * ((dim * (dim + 1)) / 2 + 2 * dim + 2) *
1210 const unsigned int allocated_size = size_scratch_data + size_data_arrays + 12;
1214 Number(numbers::signaling_nan<ScalarNumber>()));
1241template <
int dim,
typename Number,
bool is_face>
1247 ExcMessage(
"Faces can only be set if the is_face template parameter "
1272template <
int dim,
typename Number,
bool is_face>
1275 const unsigned int q_point)
const
1280 "update_normal_vectors"));
1290template <
int dim,
typename Number,
bool is_face>
1293 const unsigned int q_point)
const
1300template <
int dim,
typename Number,
bool is_face>
1307 "update_values|update_gradients"));
1319template <
int dim,
typename Number,
bool is_face>
1322 const unsigned int q)
const
1327 "update_quadrature_points"));
1332 if (is_face ==
false &&
1340 for (
unsigned int d = 0;
d < dim; ++
d)
1341 point[d] += jac[d][d] *
static_cast<typename Number::value_type
>(
1344 for (
unsigned int d = 0;
d < dim; ++
d)
1345 for (
unsigned int e = 0;
e < dim; ++
e)
1346 point[d] += jac[d][e] *
static_cast<typename Number::value_type
>(
1356template <
int dim,
typename Number,
bool is_face>
1359 const unsigned int q_point)
const
1364 "update_gradients"));
1373template <
int dim,
typename Number,
bool is_face>
1374inline const Number *
1382template <
int dim,
typename Number,
bool is_face>
1394template <
int dim,
typename Number,
bool is_face>
1395inline const Number *
1406template <
int dim,
typename Number,
bool is_face>
1419template <
int dim,
typename Number,
bool is_face>
1420inline const Number *
1432template <
int dim,
typename Number,
bool is_face>
1445template <
int dim,
typename Number,
bool is_face>
1446inline const Number *
1457template <
int dim,
typename Number,
bool is_face>
1469template <
int dim,
typename Number,
bool is_face>
1486template <
int dim,
typename Number,
bool is_face>
1498template <
int dim,
typename Number,
bool is_face>
1508template <
int dim,
typename Number,
bool is_face>
1514 "FEEvaluation was not initialized with a MatrixFree object!"));
1520template <
int dim,
typename Number,
bool is_face>
1521inline const std::vector<unsigned int> &
1529template <
int dim,
typename Number,
bool is_face>
1538template <
int dim,
typename Number,
bool is_face>
1551template <
int dim,
typename Number,
bool is_face>
1560template <
int dim,
typename Number,
bool is_face>
1569template <
int dim,
typename Number,
bool is_face>
1578template <
int dim,
typename Number,
bool is_face>
1587template <
int dim,
typename Number,
bool is_face>
1596 ExcMessage(
"All face numbers can only be queried for ECL at exterior "
1597 "faces. Use get_face_no() in other cases."));
1604template <
int dim,
typename Number,
bool is_face>
1613template <
int dim,
typename Number,
bool is_face>
1616 const unsigned int v)
const
1623 ExcMessage(
"All face numbers can only be queried for ECL at exterior "
1624 "faces. Use get_face_no() in other cases."));
1631template <
int dim,
typename Number,
bool is_face>
1640template <
int dim,
typename Number,
bool is_face>
1649template <
int dim,
typename Number,
bool is_face>
1660 template <std::size_t N,
1661 typename VectorOfArrayType,
1665 process_cell_or_face_data(
const std::array<unsigned int, N> indices,
1666 VectorOfArrayType &array,
1670 for (
unsigned int i = 0; i < N; ++i)
1674 fu(out[i], array[indices[i] / N][indices[i] % N]);
1678 template <std::
size_t N,
typename VectorOfArrayType,
typename ArrayType>
1680 set_valid_element_to_array(
const std::array<unsigned int, N> indices,
1681 const VectorOfArrayType &array,
1691 for (
unsigned int i = 0; i < N; ++i)
1692 out[i] = array[indices[
index] / N][indices[
index] % N];
1698template <
int dim,
typename Number,
bool is_face>
1699template <
typename T>
1705 internal::set_valid_element_to_array(this->get_cell_ids(), array, out);
1706 internal::process_cell_or_face_data(this->get_cell_ids(),
1709 [](
auto &local,
const auto &global) {
1717template <
int dim,
typename Number,
bool is_face>
1718template <
typename T>
1723 internal::process_cell_or_face_data(this->get_cell_ids(),
1726 [](
const auto &local,
auto &global) {
1733template <
int dim,
typename Number,
bool is_face>
1734template <
typename T>
1740 internal::set_valid_element_to_array(this->get_cell_ids(), array, out);
1741 internal::process_cell_or_face_data(this->get_face_ids(),
1744 [](
auto &local,
const auto &global) {
1752template <
int dim,
typename Number,
bool is_face>
1753template <
typename T>
1758 internal::process_cell_or_face_data(this->get_face_ids(),
1761 [](
const auto &local,
auto &global) {
void resize_fast(const size_type new_size)
void resize(const size_type new_size)
ArrayView(const std::initializer_list< std::remove_cv_t< value_type > > &initializer_list) DEAL_II_CXX20_REQUIRES(std void reinit(value_type *starting_element, const std::size_t n_elements)
AlignedVector< VectorizedArrayType > * scratch_data_array
const unsigned int quad_no
bool hessians_quad_submitted
internal::MatrixFreeFunctions::GeometryType get_cell_type() const
const Tensor< 2, dim, Number > * jacobian
const MappingInfoStorageType::QuadratureDescriptor * descriptor
const unsigned int n_quadrature_points
const MappingInfoStorageType * mapping_data
internal::MatrixFreeFunctions::DoFInfo::DoFAccessIndex dof_access_index
std::uint8_t get_face_no(const unsigned int v=0) const
ArrayView< Number > scratch_data
internal::MatrixFreeFunctions::DoFInfo::DoFAccessIndex get_dof_access_index() const
const Point< dim, Number > * quadrature_points
unsigned int subface_index
const Tensor< 1, dim *(dim+1)/2, Tensor< 1, dim, Number > > * jacobian_gradients_non_inverse
bool values_quad_submitted
FEEvaluationData(const FEEvaluationData &other)=default
const ShapeInfoType & get_shape_info() const
void set_face_data(AlignedVector< T > &array, const T &value) const
void reinit_face(const internal::MatrixFreeFunctions::FaceToCellTopology< n_lanes > &face)
Number JxW(const unsigned int q_point) const
const std::array< unsigned int, n_lanes > & get_face_ids() const
unsigned int get_first_selected_component() const
const unsigned int n_fe_components
const ShapeInfoType * data
T read_face_data(const AlignedVector< T > &array) const
Number * gradients_from_hessians_quad
std::shared_ptr< internal::MatrixFreeFunctions::MappingDataOnTheFly< dim, Number > > mapped_geometry
internal::MatrixFreeFunctions::DoFInfo DoFInfo
unsigned int get_active_quadrature_index() const
const internal::MatrixFreeFunctions::DoFInfo & get_dof_info() const
const std::array< unsigned int, n_lanes > & get_cell_ids() const
unsigned int get_mapping_data_index_offset() const
bool gradients_quad_initialized
void set_data_pointers(AlignedVector< Number > *scratch_data, const unsigned int n_components)
Tensor< 2, dim, Number > inverse_jacobian(const unsigned int q_point) const
std::array< std::uint8_t, n_lanes > face_numbers
const unsigned int active_fe_index
static constexpr unsigned int n_lanes
const Tensor< 1, dim, Number > * normal_vectors
internal::MatrixFreeFunctions::GeometryType cell_type
const Number * begin_gradients() const
const std::array< unsigned int, n_lanes > & get_cell_or_face_ids() const
Tensor< 1, dim, Number > normal_vector(const unsigned int q_point) const
std::array< unsigned int, n_lanes > cell_ids
unsigned int get_current_cell_index() const
unsigned int get_subface_index() const
bool divergence_is_requested
Number shape_info_number_type
Tensor< 1, dim, Number > get_normal_vector(const unsigned int q_point) const
const unsigned int active_quad_index
unsigned int get_active_fe_index() const
Number * values_from_gradients_quad
void set_cell_data(AlignedVector< T > &array, const T &value) const
bool gradients_quad_submitted
unsigned int get_quadrature_index() const
std::array< std::uint8_t, n_lanes > face_orientations
const ScalarNumber * quadrature_weights
const Tensor< 1, dim *(dim+1)/2, Tensor< 1, dim, Number > > * jacobian_gradients
const Tensor< 1, dim, Number > * normal_x_jacobian
bool is_interior_face() const
std_cxx20::ranges::iota_view< unsigned int, unsigned int > quadrature_point_indices() const
const std::vector< unsigned int > & get_internal_dof_numbering() const
FEEvaluationData & operator=(const FEEvaluationData &other)
virtual ~FEEvaluationData()=default
ArrayView< Number > get_scratch_data() const
T read_cell_data(const AlignedVector< T > &array) const
static constexpr unsigned int dimension
internal::MatrixFreeFunctions::ShapeInfo< VectorizedArrayType > ShapeInfoType
FEEvaluationData(const std::shared_ptr< internal::MatrixFreeFunctions::MappingDataOnTheFly< dim, Number > > &mapping_data, const unsigned int n_fe_components, const unsigned int first_selected_component)
bool values_quad_initialized
FEEvaluationData(const ShapeInfoType &shape_info, const bool is_interior_face=true)
Point< dim, Number > quadrature_point(const unsigned int q) const
const Number * begin_values() const
bool dof_values_initialized
typename internal::VectorizedArrayTrait< Number >::value_type ScalarNumber
std::array< unsigned int, n_lanes > face_ids
FEEvaluationData(const InitializationData &initialization_data, const bool is_interior_face, const unsigned int quad_no, const unsigned int first_selected_component)
const Number * begin_dof_values() const
std::uint8_t get_face_orientation(const unsigned int v=0) const
const Number * begin_hessians() const
bool hessians_quad_initialized
unsigned int get_cell_or_face_batch_id() const
const unsigned int first_selected_component
const unsigned int dofs_per_component
const unsigned int n_q_points
static constexpr unsigned int n_components
const Point< dim > & point(const unsigned int i) const
#define DEAL_II_ALWAYS_INLINE
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
#define DeclException0(Exception0)
#define Assert(cond, exc)
#define AssertDimension(dim1, dim2)
#define AssertIndexRange(index, range)
static ::ExceptionBase & ExcInternalError()
static ::ExceptionBase & ExcAccessToUninitializedField()
static ::ExceptionBase & ExcNotInitialized()
#define DeclException1(Exception1, type1, outsequence)
static ::ExceptionBase & ExcMatrixFreeAccessToUninitializedMappingField(std::string arg1)
static ::ExceptionBase & ExcMessage(std::string arg1)
@ general
No special properties.
Point< spacedim > point(const gp_Pnt &p, const double tolerance=1e-10)
SymmetricTensor< 2, dim, Number > e(const Tensor< 2, dim, Number > &F)
SymmetricTensor< 2, dim, Number > d(const Tensor< 2, dim, Number > &F, const Tensor< 2, dim, Number > &dF_dt)
static const unsigned int invalid_unsigned_int
const InputIterator OutputIterator out
boost::integer_range< IncrementableType > iota_view
::VectorizedArray< Number, width > max(const ::VectorizedArray< Number, width > &, const ::VectorizedArray< Number, width > &)
const MappingInfoStorageType * mapping_data
unsigned int active_quad_index
const MappingInfoStorageType::QuadratureDescriptor * descriptor
const ShapeInfoType * shape_info
unsigned int active_fe_index
@ dof_access_face_exterior
@ dof_access_face_interior
unsigned char subface_index
unsigned char interior_face_no
std::array< unsigned int, vectorization_width > cells_interior
types::boundary_id exterior_face_no
std::array< unsigned int, vectorization_width > cells_exterior
unsigned char face_orientation
Quadrature< structdim > quadrature
AlignedVector< unsigned int > data_index_offsets
unsigned int dofs_per_component_on_cell
std::vector< UnivariateShapeData< Number > > data
std::vector< unsigned int > lexicographic_numbering
static constexpr std::size_t width()