|
Reference documentation for deal.II version 9.2.0
|
\(\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\}}\)
Go to the documentation of this file.
16 #ifndef dealii_fe_poly_h
17 #define dealii_fe_poly_h
74 int dim = PolynomialType::dimension,
103 std::vector<unsigned int>
110 std::vector<unsigned int>
134 const unsigned int component)
const override;
157 const unsigned int component)
const override;
180 const unsigned int component)
const override;
204 const unsigned int component)
const override;
228 const unsigned int component)
const override;
244 virtual std::unique_ptr<
252 &output_data)
const override
256 std::unique_ptr<typename FiniteElement<dim, spacedim>::InternalDataBase>
257 data_ptr = std_cxx14::make_unique<InternalData>();
258 auto &data =
dynamic_cast<InternalData &
>(*data_ptr);
261 const unsigned int n_q_points = quadrature.
size();
266 std::vector<double> values(
268 std::vector<Tensor<1, dim>> grads(
270 std::vector<Tensor<2, dim>> grad_grads(
272 std::vector<Tensor<3, dim>> third_derivatives(
274 std::vector<Tensor<4, dim>>
294 !((output_data.shape_values.n_rows() > 0) &&
295 (output_data.shape_values.n_cols() == n_q_points)))
296 data.shape_values.reinit(this->dofs_per_cell, n_q_points);
299 data.shape_gradients.reinit(this->dofs_per_cell, n_q_points);
302 data.shape_hessians.reinit(this->dofs_per_cell, n_q_points);
305 data.shape_3rd_derivatives.reinit(this->dofs_per_cell, n_q_points);
312 for (
unsigned int i = 0; i < n_q_points; ++i)
329 if (output_data.shape_values.n_rows() > 0)
331 if (output_data.shape_values.n_cols() == n_q_points)
333 output_data.shape_values[k][i] = values[k];
336 data.shape_values[k][i] = values[k];
344 data.shape_gradients[k][i] = grads[k];
348 data.shape_hessians[k][i] = grad_grads[k];
352 data.shape_3rd_derivatives[k][i] = third_derivatives[k];
364 const ::internal::FEValuesImplementation::MappingRelatedData<dim,
370 &output_data)
const override;
375 const unsigned int face_no,
379 const ::internal::FEValuesImplementation::MappingRelatedData<dim,
385 &output_data)
const override;
390 const unsigned int face_no,
391 const unsigned int sub_no,
395 const ::internal::FEValuesImplementation::MappingRelatedData<dim,
401 &output_data)
const override;
479 const unsigned int n_q_points)
const;
509 const unsigned int n_q_points)
const;
void correct_hessians(internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim > &output_data, const internal::FEValuesImplementation::MappingRelatedData< dim, spacedim > &mapping_data, const unsigned int n_q_points) const
virtual Tensor< 1, dim > shape_grad(const unsigned int i, const Point< dim > &p) const override
virtual double shape_value(const unsigned int i, const Point< dim > &p) const override
PolynomialType poly_space
unsigned int get_degree() const
Table< 2, Tensor< 3, dim > > shape_3rd_derivatives
@ update_3rd_derivatives
Third derivatives of shape functions.
virtual std::unique_ptr< typename FiniteElement< dim, spacedim >::InternalDataBase > get_data(const UpdateFlags update_flags, const Mapping< dim, spacedim > &, const Quadrature< dim > &quadrature, ::internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim > &output_data) const override
void correct_third_derivatives(internal::FEValuesImplementation::FiniteElementRelatedData< dim, spacedim > &output_data, const internal::FEValuesImplementation::MappingRelatedData< dim, spacedim > &mapping_data, const unsigned int n_q_points) const
virtual Tensor< 1, dim > shape_grad_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const override
friend class InternalDataBase
virtual UpdateFlags requires_update_flags(const UpdateFlags update_flags) const override
Table< 2, double > shape_values
@ update_values
Shape function values.
virtual Tensor< 2, dim > shape_grad_grad_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const override
virtual Tensor< 4, dim > shape_4th_derivative_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const override
virtual Tensor< 3, dim > shape_3rd_derivative(const unsigned int i, const Point< dim > &p) const override
const Point< dim > & point(const unsigned int i) const
Abstract base class for mapping classes.
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
virtual Tensor< 4, dim > shape_4th_derivative(const unsigned int i, const Point< dim > &p) const override
virtual double shape_value_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const override
const unsigned int dofs_per_cell
const std::vector< ComponentMask > nonzero_components
std::vector< unsigned int > get_poly_space_numbering() const
#define DEAL_II_NAMESPACE_OPEN
@ update_gradients
Shape function gradients.
virtual std::size_t memory_consumption() const override
Table< 2, Tensor< 2, dim > > shape_hessians
Table< 2, Tensor< 1, dim > > shape_gradients
virtual Tensor< 3, dim > shape_3rd_derivative_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const override
FE_Poly(const PolynomialType &poly_space, const FiniteElementData< dim > &fe_data, const std::vector< bool > &restriction_is_additive_flags, const std::vector< ComponentMask > &nonzero_components)
@ update_hessians
Second derivatives of shape functions.
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
const std::vector< bool > restriction_is_additive_flags
virtual void fill_fe_face_values(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const unsigned int face_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
unsigned int size() const
#define DEAL_II_NAMESPACE_CLOSE
virtual Tensor< 2, dim > shape_grad_grad(const unsigned int i, const Point< dim > &p) const override
std::vector< unsigned int > get_poly_space_numbering_inverse() const