16 #ifndef dealii_fe_poly_h 17 #define dealii_fe_poly_h 20 #include <deal.II/base/quadrature.h> 21 #include <deal.II/base/std_cxx14/memory.h> 23 #include <deal.II/fe/fe.h> 25 DEAL_II_NAMESPACE_OPEN
69 template <
class PolynomialType,
70 int dim = PolynomialType::dimension,
99 std::vector<unsigned int>
106 std::vector<unsigned int>
130 const unsigned int component)
const override;
153 const unsigned int component)
const override;
176 const unsigned int component)
const override;
200 const unsigned int component)
const override;
224 const unsigned int component)
const override;
234 virtual std::unique_ptr<
242 &output_data)
const override 246 auto data = std_cxx14::make_unique<InternalData>();
249 const unsigned int n_q_points = quadrature.
size();
254 std::vector<double> values(
256 std::vector<Tensor<1, dim>> grads(
258 std::vector<Tensor<2, dim>> grad_grads(
260 std::vector<Tensor<3, dim>> third_derivatives(
262 std::vector<Tensor<4, dim>>
282 !((output_data.shape_values.n_rows() > 0) &&
283 (output_data.shape_values.n_cols() == n_q_points)))
284 data->shape_values.reinit(this->dofs_per_cell, n_q_points);
287 data->shape_gradients.reinit(this->dofs_per_cell, n_q_points);
290 data->shape_hessians.reinit(this->dofs_per_cell, n_q_points);
293 data->shape_3rd_derivatives.reinit(this->dofs_per_cell, n_q_points);
298 if (update_flags & (update_values | update_gradients | update_hessians |
299 update_3rd_derivatives))
300 for (
unsigned int i = 0; i < n_q_points; ++i)
302 poly_space.compute(quadrature.
point(i),
316 if (update_flags & update_values)
317 if (output_data.shape_values.n_rows() > 0)
319 if (output_data.shape_values.n_cols() == n_q_points)
321 output_data.shape_values[k][i] = values[k];
324 data->shape_values[k][i] = values[k];
330 if (update_flags & update_gradients)
332 data->shape_gradients[k][i] = grads[k];
334 if (update_flags & update_hessians)
336 data->shape_hessians[k][i] = grad_grads[k];
338 if (update_flags & update_3rd_derivatives)
340 data->shape_3rd_derivatives[k][i] = third_derivatives[k];
342 return std::move(data);
352 const ::internal::FEValuesImplementation::MappingRelatedData<dim,
358 &output_data)
const override;
363 const unsigned int face_no,
367 const ::internal::FEValuesImplementation::MappingRelatedData<dim,
373 &output_data)
const override;
376 fill_fe_subface_values(
378 const unsigned int face_no,
379 const unsigned int sub_no,
383 const ::internal::FEValuesImplementation::MappingRelatedData<dim,
389 &output_data)
const override;
473 const unsigned int n_q_points,
474 const unsigned int dof)
const;
485 DEAL_II_NAMESPACE_CLOSE
virtual double shape_value(const unsigned int i, const Point< dim > &p) const override
Table< 2, Tensor< 2, dim > > shape_hessians
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 unsigned int dof) const
virtual Tensor< 4, dim > shape_4th_derivative(const unsigned int i, const Point< dim > &p) const override
virtual Tensor< 1, dim > shape_grad(const unsigned int i, const Point< dim > &p) const override
virtual Tensor< 2, dim > shape_grad_grad(const unsigned int i, const Point< dim > &p) const override
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
const Point< dim > & point(const unsigned int i) const
virtual Tensor< 2, dim > shape_grad_grad_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const override
Table< 2, Tensor< 1, dim > > shape_gradients
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)
unsigned int size() const
Third derivatives of shape functions.
Abstract base class for mapping classes.
Table< 2, Tensor< 3, dim > > shape_3rd_derivatives
unsigned int get_degree() const
Second derivatives of shape functions.
const unsigned int dofs_per_cell
std::vector< unsigned int > get_poly_space_numbering_inverse() const
virtual Tensor< 3, dim > shape_3rd_derivative_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const override
Table< 2, double > shape_values
std::vector< unsigned int > get_poly_space_numbering() const
virtual UpdateFlags requires_update_flags(const UpdateFlags update_flags) const override
virtual Tensor< 4, dim > shape_4th_derivative_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const override
Shape function gradients.
virtual Tensor< 3, dim > shape_3rd_derivative(const unsigned int i, const Point< dim > &p) const override
const std::vector< ComponentMask > nonzero_components
virtual double shape_value_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const override
PolynomialType poly_space
const std::vector< bool > restriction_is_additive_flags
virtual Tensor< 1, dim > shape_grad_component(const unsigned int i, const Point< dim > &p, const unsigned int component) const override