deal.II version GIT relicensing-2167-g9622207b8f 2024-11-21 12:40:00+00:00
|
#include <deal.II/fe/fe_update_flags.h>
Public Types | |
using | ShapeVector = ::Table< 2, double > |
using | GradientVector = ::Table< 2, Tensor< 1, spacedim > > |
using | HessianVector = ::Table< 2, Tensor< 2, spacedim > > |
using | ThirdDerivativeVector = ::Table< 2, Tensor< 3, spacedim > > |
Public Member Functions | |
void | initialize (const unsigned int n_quadrature_points, const FiniteElement< dim, spacedim > &fe, const UpdateFlags flags) |
std::size_t | memory_consumption () const |
Public Attributes | |
ShapeVector | shape_values |
GradientVector | shape_gradients |
HessianVector | shape_hessians |
ThirdDerivativeVector | shape_3rd_derivatives |
std::vector< unsigned int > | shape_function_to_row_table |
A class that stores all of the shape function related data used in FEValues, FEFaceValues, and FESubfaceValues objects. Objects of this kind will be given as output argument when FEValues::reinit() calls FiniteElement::fill_fe_values().
Definition at line 410 of file fe_update_flags.h.
using internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim >::ShapeVector = ::Table<2, double> |
Storage type for shape values. Each row in the matrix denotes the values of a single shape function at the different points, columns are for a single point with the different shape functions.
If a shape function has more than one non-zero component (in deal.II diction: it is non-primitive), then we allocate one row per non-zero component, and shift subsequent rows backward. Lookup of the correct row for a shape function is thus simple in case the entire finite element is primitive (i.e. all shape functions are primitive), since then the shape function number equals the row number. Otherwise, use the shape_function_to_row_table array to get at the first row that belongs to this particular shape function, and navigate among all the rows for this shape function using the FiniteElement::get_nonzero_components() function which tells us which components are non-zero and thus have a row in the array presently under discussion.
Definition at line 446 of file fe_update_flags.h.
using internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim >::GradientVector = ::Table<2, Tensor<1, spacedim> > |
Storage type for gradients. The layout of data is the same as for the ShapeVector data type.
Definition at line 452 of file fe_update_flags.h.
using internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim >::HessianVector = ::Table<2, Tensor<2, spacedim> > |
Likewise for second order derivatives.
Definition at line 457 of file fe_update_flags.h.
using internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim >::ThirdDerivativeVector = ::Table<2, Tensor<3, spacedim> > |
And the same also applies to the third order derivatives.
Definition at line 462 of file fe_update_flags.h.
void internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim >::initialize | ( | const unsigned int | n_quadrature_points, |
const FiniteElement< dim, spacedim > & | fe, | ||
const UpdateFlags | flags | ||
) |
Initialize all vectors to correct size.
Definition at line 83 of file fe_values.cc.
std::size_t internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim >::memory_consumption | ( | ) | const |
Compute and return an estimate for the memory consumption (in bytes) of this object.
Definition at line 140 of file fe_values.cc.
ShapeVector internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim >::shape_values |
Store the values of the shape functions at the quadrature points. See the description of the data type for the layout of the data in this field.
Definition at line 469 of file fe_update_flags.h.
GradientVector internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim >::shape_gradients |
Store the gradients of the shape functions at the quadrature points. See the description of the data type for the layout of the data in this field.
Definition at line 476 of file fe_update_flags.h.
HessianVector internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim >::shape_hessians |
Store the 2nd derivatives of the shape functions at the quadrature points. See the description of the data type for the layout of the data in this field.
Definition at line 483 of file fe_update_flags.h.
ThirdDerivativeVector internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim >::shape_3rd_derivatives |
Store the 3rd derivatives of the shape functions at the quadrature points. See the description of the data type for the layout of the data in this field.
Definition at line 490 of file fe_update_flags.h.
std::vector<unsigned int> internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim >::shape_function_to_row_table |
When asked for the value (or gradient, or Hessian) of shape function i's c-th vector component, we need to look it up in the shape_values, shape_gradients and shape_hessians arrays. The question is where in this array does the data for shape function i, component c reside. This is what this table answers.
The format of the table is as follows: - It has dofs_per_cell times n_components entries. - The entry that corresponds to shape function i, component c is i * n_components + c
. - The value stored at this position indicates the row in shape_values and the other tables where the corresponding datum is stored for all the quadrature points.
In the general, vector-valued context, the number of components is larger than one, but for a given shape function, not all vector components may be nonzero (e.g., if a shape function is primitive, then exactly one vector component is non-zero, while the others are all zero). For such zero components, shape_values and friends do not have a row. Consequently, for vector components for which shape function i is zero, the entry in the current table is numbers::invalid_unsigned_int.
On the other hand, the table is guaranteed to have at least one valid index for each shape function. In particular, for a primitive finite element, each shape function has exactly one nonzero component and so for each i, there is exactly one valid index within the range [i*n_components, (i+1)*n_components)
.
Definition at line 521 of file fe_update_flags.h.