16#ifndef dealii_mapping_fe_h
17#define dealii_mapping_fe_h
59template <
int dim,
int spacedim = dim>
74 virtual std::unique_ptr<Mapping<dim, spacedim>>
75 clone()
const override;
87 &cell)
const override;
200 const unsigned int n_original_q_points);
210 const unsigned int n_original_q_points);
225 shape(
const unsigned int qpoint,
const unsigned int shape_nr)
const;
231 shape(
const unsigned int qpoint,
const unsigned int shape_nr);
237 derivative(
const unsigned int qpoint,
const unsigned int shape_nr)
const;
243 derivative(
const unsigned int qpoint,
const unsigned int shape_nr);
250 const unsigned int shape_nr)
const;
263 const unsigned int shape_nr)
const;
276 const unsigned int shape_nr)
const;
334 std::array<std::vector<Tensor<1, dim>>,
362 mutable std::vector<DerivativeForm<1, dim, spacedim>>
covariant;
376 mutable std::vector<std::vector<Tensor<1, spacedim>>>
aux;
408 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
414 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
419 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
431 &output_data)
const override;
439 const unsigned int face_no,
443 &output_data)
const override;
449 const unsigned int face_no,
450 const unsigned int subface_no,
454 &output_data)
const override;
461 const std::unique_ptr<FiniteElement<dim, spacedim>>
fe;
472 virtual std::vector<Point<spacedim>>
488template <
int dim,
int spacedim>
491 const unsigned int shape_nr)
const
499template <
int dim,
int spacedim>
502 const unsigned int shape_nr)
504 AssertIndexRange(qpoint * n_shape_functions + shape_nr, shape_values.size());
505 return shape_values[qpoint * n_shape_functions + shape_nr];
509template <
int dim,
int spacedim>
512 const unsigned int qpoint,
513 const unsigned int shape_nr)
const
516 shape_derivatives.size());
517 return shape_derivatives[qpoint * n_shape_functions + shape_nr];
522template <
int dim,
int spacedim>
525 const unsigned int shape_nr)
528 shape_derivatives.size());
529 return shape_derivatives[qpoint * n_shape_functions + shape_nr];
533template <
int dim,
int spacedim>
536 const unsigned int qpoint,
537 const unsigned int shape_nr)
const
540 shape_second_derivatives.size());
541 return shape_second_derivatives[qpoint * n_shape_functions + shape_nr];
545template <
int dim,
int spacedim>
548 const unsigned int qpoint,
549 const unsigned int shape_nr)
552 shape_second_derivatives.size());
553 return shape_second_derivatives[qpoint * n_shape_functions + shape_nr];
556template <
int dim,
int spacedim>
559 const unsigned int qpoint,
560 const unsigned int shape_nr)
const
563 shape_third_derivatives.size());
564 return shape_third_derivatives[qpoint * n_shape_functions + shape_nr];
568template <
int dim,
int spacedim>
571 const unsigned int qpoint,
572 const unsigned int shape_nr)
575 shape_third_derivatives.size());
576 return shape_third_derivatives[qpoint * n_shape_functions + shape_nr];
580template <
int dim,
int spacedim>
583 const unsigned int qpoint,
584 const unsigned int shape_nr)
const
587 shape_fourth_derivatives.size());
588 return shape_fourth_derivatives[qpoint * n_shape_functions + shape_nr];
592template <
int dim,
int spacedim>
595 const unsigned int qpoint,
596 const unsigned int shape_nr)
599 shape_fourth_derivatives.size());
600 return shape_fourth_derivatives[qpoint * n_shape_functions + shape_nr];
605template <
int dim,
int spacedim>
const Tensor< 3, dim > & third_derivative(const unsigned int qpoint, const unsigned int shape_nr) const
Tensor< 2, dim > & second_derivative(const unsigned int qpoint, const unsigned int shape_nr)
void initialize_face(const UpdateFlags update_flags, const Quadrature< dim > &quadrature, const unsigned int n_original_q_points)
std::vector< Point< spacedim > > mapping_support_points
const Tensor< 4, dim > & fourth_derivative(const unsigned int qpoint, const unsigned int shape_nr) const
std::vector< std::vector< Tensor< 1, spacedim > > > aux
std::vector< Tensor< 3, dim > > shape_third_derivatives
std::vector< Tensor< 2, dim > > shape_second_derivatives
Tensor< 1, dim > & derivative(const unsigned int qpoint, const unsigned int shape_nr)
std::vector< Tensor< 4, dim > > shape_fourth_derivatives
Tensor< 3, dim > & third_derivative(const unsigned int qpoint, const unsigned int shape_nr)
void initialize(const UpdateFlags update_flags, const Quadrature< dim > &quadrature, const unsigned int n_original_q_points)
Triangulation< dim, spacedim >::cell_iterator cell_of_current_support_points
Tensor< 4, dim > & fourth_derivative(const unsigned int qpoint, const unsigned int shape_nr)
void compute_shape_function_values(const std::vector< Point< dim > > &unit_points)
std::array< std::vector< Tensor< 1, dim > >, GeometryInfo< dim >::faces_per_cell *(dim - 1)> unit_tangentials
std::vector< double > quadrature_weights
const FiniteElement< dim, spacedim > & fe
std::vector< DerivativeForm< 1, dim, spacedim > > contravariant
std::vector< DerivativeForm< 1, dim, spacedim > > covariant
const Tensor< 1, dim > & derivative(const unsigned int qpoint, const unsigned int shape_nr) const
double & shape(const unsigned int qpoint, const unsigned int shape_nr)
std::vector< double > shape_values
std::vector< Tensor< 1, dim > > shape_derivatives
virtual std::size_t memory_consumption() const override
const Tensor< 2, dim > & second_derivative(const unsigned int qpoint, const unsigned int shape_nr) const
const unsigned int n_shape_functions
std::vector< double > volume_elements
const double & shape(const unsigned int qpoint, const unsigned int shape_nr) const
const unsigned int polynomial_degree
virtual CellSimilarity::Similarity fill_fe_values(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const CellSimilarity::Similarity cell_similarity, const Quadrature< dim > &quadrature, const typename Mapping< dim, spacedim >::InternalDataBase &internal_data, internal::FEValuesImplementation::MappingRelatedData< dim, spacedim > &output_data) const override
const unsigned int polynomial_degree
virtual std::unique_ptr< typename Mapping< dim, spacedim >::InternalDataBase > get_face_data(const UpdateFlags flags, const hp::QCollection< dim - 1 > &quadrature) const override
virtual Point< dim > transform_real_to_unit_cell(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const Point< spacedim > &p) const override
virtual UpdateFlags requires_update_flags(const UpdateFlags update_flags) const override
virtual void fill_fe_face_values(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const unsigned int face_no, const hp::QCollection< dim - 1 > &quadrature, const typename Mapping< dim, spacedim >::InternalDataBase &internal_data, internal::FEValuesImplementation::MappingRelatedData< dim, spacedim > &output_data) const override
Table< 2, double > mapping_support_point_weights
virtual std::vector< Point< spacedim > > compute_mapping_support_points(const typename Triangulation< dim, spacedim >::cell_iterator &cell) const
virtual void transform(const ArrayView< const Tensor< 1, dim > > &input, const MappingKind kind, const typename Mapping< dim, spacedim >::InternalDataBase &internal, const ArrayView< Tensor< 1, spacedim > > &output) const override
virtual BoundingBox< spacedim > get_bounding_box(const typename Triangulation< dim, spacedim >::cell_iterator &cell) const override
virtual std::unique_ptr< typename Mapping< dim, spacedim >::InternalDataBase > get_data(const UpdateFlags, const Quadrature< dim > &quadrature) const override
virtual std::unique_ptr< typename Mapping< dim, spacedim >::InternalDataBase > get_subface_data(const UpdateFlags flags, const Quadrature< dim - 1 > &quadrature) const override
unsigned int get_degree() const
virtual bool is_compatible_with(const ReferenceCell &reference_cell) const override
virtual void fill_fe_subface_values(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const unsigned int face_no, const unsigned int subface_no, const Quadrature< dim - 1 > &quadrature, const typename Mapping< dim, spacedim >::InternalDataBase &internal_data, internal::FEValuesImplementation::MappingRelatedData< dim, spacedim > &output_data) const override
virtual bool preserves_vertex_locations() const override
virtual std::unique_ptr< Mapping< dim, spacedim > > clone() const override
virtual Point< spacedim > transform_unit_to_real_cell(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const Point< dim > &p) const override
const std::unique_ptr< FiniteElement< dim, spacedim > > fe
Abstract base class for mapping classes.
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
#define AssertIndexRange(index, range)