17 #ifndef dealii_matrix_free_fe_evaluation_data_h
18 #define dealii_matrix_free_fe_evaluation_data_h
50 <<
"You are requesting information from an FEEvaluation/FEFaceEvaluation "
51 <<
"object for which this kind of information has not been computed. What "
52 <<
"information these objects compute is determined by the update_* flags "
53 <<
"you pass to MatrixFree::reinit() via MatrixFree::AdditionalData. "
54 <<
"Here, the operation you are attempting requires the <" << arg1
55 <<
"> flag to be set, but it was apparently not specified "
56 <<
"upon initialization.");
62 int n_q_points_1d = fe_degree + 1,
63 int n_components_ = 1,
77 namespace MatrixFreeFunctions
79 template <
int,
typename>
80 class MappingDataOnTheFly;
113 template <
int dim,
typename Number,
bool is_face>
118 MappingInfoStorage<(is_face ? dim - 1 : dim), dim, Number>;
165 const unsigned int n_components);
194 JxW(
const unsigned int q_point)
const;
396 const std::vector<unsigned int> &
494 const std::array<unsigned int, n_lanes> &
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>
1023 template <
int dim,
typename Number,
bool is_face>
1028 InitializationData{&shape_info, nullptr, nullptr, 0, 0, nullptr},
1035 template <
int dim,
typename Number,
bool is_face>
1037 const InitializationData &initialization_data,
1041 :
data(initialization_data.shape_info)
1086 template <
int dim,
typename Number,
bool is_face>
1088 const std::shared_ptr<
1126 .jacobian_gradients_non_inverse[0]
1134 template <
int dim,
typename Number,
bool is_face>
1187 template <
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>()));
1241 template <
int dim,
typename Number,
bool is_face>
1247 ExcMessage(
"Faces can only be set if the is_face template parameter "
1272 template <
int dim,
typename Number,
bool is_face>
1275 const unsigned int q_point)
const
1280 "update_normal_vectors"));
1290 template <
int dim,
typename Number,
bool is_face>
1293 const unsigned int q_point)
const
1300 template <
int dim,
typename Number,
bool is_face>
1307 "update_values|update_gradients"));
1319 template <
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
>(
1356 template <
int dim,
typename Number,
bool is_face>
1359 const unsigned int q_point)
const
1364 "update_gradients"));
1373 template <
int dim,
typename Number,
bool is_face>
1374 inline const Number *
1382 template <
int dim,
typename Number,
bool is_face>
1394 template <
int dim,
typename Number,
bool is_face>
1395 inline const Number *
1406 template <
int dim,
typename Number,
bool is_face>
1419 template <
int dim,
typename Number,
bool is_face>
1420 inline const Number *
1432 template <
int dim,
typename Number,
bool is_face>
1445 template <
int dim,
typename Number,
bool is_face>
1446 inline const Number *
1457 template <
int dim,
typename Number,
bool is_face>
1469 template <
int dim,
typename Number,
bool is_face>
1486 template <
int dim,
typename Number,
bool is_face>
1498 template <
int dim,
typename Number,
bool is_face>
1508 template <
int dim,
typename Number,
bool is_face>
1514 "FEEvaluation was not initialized with a MatrixFree object!"));
1520 template <
int dim,
typename Number,
bool is_face>
1521 inline const std::vector<unsigned int> &
1529 template <
int dim,
typename Number,
bool is_face>
1538 template <
int dim,
typename Number,
bool is_face>
1551 template <
int dim,
typename Number,
bool is_face>
1560 template <
int dim,
typename Number,
bool is_face>
1569 template <
int dim,
typename Number,
bool is_face>
1578 template <
int dim,
typename Number,
bool is_face>
1587 template <
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."));
1604 template <
int dim,
typename Number,
bool is_face>
1613 template <
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."));
1631 template <
int dim,
typename Number,
bool is_face>
1640 template <
int dim,
typename Number,
bool is_face>
1649 template <
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];
1698 template <
int dim,
typename Number,
bool is_face>
1699 template <
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) {
1717 template <
int dim,
typename Number,
bool is_face>
1718 template <
typename T>
1723 internal::process_cell_or_face_data(this->get_cell_ids(),
1726 [](
const auto &local,
auto &global) {
1733 template <
int dim,
typename Number,
bool is_face>
1734 template <
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) {
1752 template <
int dim,
typename Number,
bool is_face>
1753 template <
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)
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
Point< dim, Number > quadrature_point(const unsigned int q) const
std_cxx20::ranges::iota_view< unsigned int, unsigned int > quadrature_point_indices() const
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
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
FEEvaluationData & operator=(const FEEvaluationData &other)
const std::array< unsigned int, n_lanes > & get_cell_or_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
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)
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
Tensor< 1, dim, Number > get_normal_vector(const unsigned int q_point) const
const Number * begin_values() const
internal::MatrixFreeFunctions::GeometryType cell_type
const Number * begin_hessians() const
const std::array< unsigned int, n_lanes > & get_cell_ids() const
Number * begin_dof_values()
std::array< unsigned int, n_lanes > cell_ids
unsigned int get_current_cell_index() const
unsigned int get_subface_index() const
Number * begin_gradients()
bool divergence_is_requested
Number shape_info_number_type
const unsigned int active_quad_index
unsigned int get_active_fe_index() const
Number * values_from_gradients_quad
const Number * begin_gradients() const
void set_cell_data(AlignedVector< T > &array, const T &value) const
bool gradients_quad_submitted
const std::vector< unsigned int > & get_internal_dof_numbering() const
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
const internal::MatrixFreeFunctions::DoFInfo & get_dof_info() const
bool is_interior_face() const
const std::array< unsigned int, n_lanes > & get_face_ids() const
virtual ~FEEvaluationData()=default
T read_cell_data(const AlignedVector< T > &array) const
ArrayView< Number > get_scratch_data() 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)
Tensor< 2, dim, Number > inverse_jacobian(const unsigned int q_point) const
bool values_quad_initialized
FEEvaluationData(const ShapeInfoType &shape_info, const bool is_interior_face=true)
bool dof_values_initialized
typename internal::VectorizedArrayTrait< Number >::value_type ScalarNumber
std::array< unsigned int, n_lanes > face_ids
const ShapeInfoType & get_shape_info() const
FEEvaluationData(const InitializationData &initialization_data, const bool is_interior_face, const unsigned int quad_no, const unsigned int first_selected_component)
Number * begin_hessians()
std::uint8_t get_face_orientation(const unsigned int v=0) const
bool hessians_quad_initialized
unsigned int get_cell_or_face_batch_id() const
Tensor< 1, dim, Number > normal_vector(const unsigned int q_point) const
const unsigned int first_selected_component
const Number * begin_dof_values() const
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)
static ::ExceptionBase & ExcAccessToUninitializedField()
static ::ExceptionBase & ExcInternalError()
static ::ExceptionBase & ExcNotInitialized()
#define Assert(cond, exc)
static ::ExceptionBase & ExcMatrixFreeAccessToUninitializedMappingField(std::string arg1)
#define AssertDimension(dim1, dim2)
#define AssertIndexRange(index, range)
#define DeclException1(Exception1, type1, outsequence)
static ::ExceptionBase & ExcMessage(std::string arg1)
@ general
No special properties.
Point< spacedim > point(const gp_Pnt &p, const double tolerance=1e-10)
void quadrature_points(const Triangulation< dim, spacedim > &triangulation, const Quadrature< dim > &quadrature, const std::vector< std::vector< BoundingBox< spacedim >>> &global_bounding_boxes, ParticleHandler< dim, spacedim > &particle_handler, const Mapping< dim, spacedim > &mapping=(ReferenceCells::get_hypercube< dim >() .template get_default_linear_mapping< dim, spacedim >()), const std::vector< std::vector< double >> &properties={})
SymmetricTensor< 2, dim, Number > d(const Tensor< 2, dim, Number > &F, const Tensor< 2, dim, Number > &dF_dt)
SymmetricTensor< 2, dim, Number > e(const Tensor< 2, dim, Number > &F)
VectorType::value_type * begin(VectorType &V)
static const unsigned int invalid_unsigned_int
boost::integer_range< IncrementableType > iota_view
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()