16 #ifndef dealii_mapping_fe_h 17 #define dealii_mapping_fe_h 57 template <
int dim,
int spacedim = dim>
72 virtual std::unique_ptr<Mapping<dim, spacedim>>
73 clone()
const override;
85 &cell)
const override;
195 const unsigned int n_original_q_points);
205 const unsigned int n_original_q_points);
220 shape(
const unsigned int qpoint,
const unsigned int shape_nr)
const;
226 shape(
const unsigned int qpoint,
const unsigned int shape_nr);
232 derivative(
const unsigned int qpoint,
const unsigned int shape_nr)
const;
238 derivative(
const unsigned int qpoint,
const unsigned int shape_nr);
245 const unsigned int shape_nr)
const;
258 const unsigned int shape_nr)
const;
271 const unsigned int shape_nr)
const;
329 std::array<std::vector<Tensor<1, dim>>,
357 mutable std::vector<DerivativeForm<1, dim, spacedim>>
covariant;
371 mutable std::vector<std::vector<Tensor<1, spacedim>>>
aux;
403 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
409 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
414 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
426 &output_data)
const override;
434 const unsigned int face_no,
438 &output_data)
const override;
444 const unsigned int face_no,
445 const unsigned int subface_no,
449 &output_data)
const override;
456 const std::unique_ptr<FiniteElement<dim, spacedim>>
fe;
467 virtual std::vector<Point<spacedim>>
483 template <
int dim,
int spacedim>
484 inline const double &
486 const unsigned int shape_nr)
const 494 template <
int dim,
int spacedim>
497 const unsigned int shape_nr)
504 template <
int dim,
int spacedim>
507 const unsigned int qpoint,
508 const unsigned int shape_nr)
const 517 template <
int dim,
int spacedim>
520 const unsigned int shape_nr)
528 template <
int dim,
int spacedim>
531 const unsigned int qpoint,
532 const unsigned int shape_nr)
const 540 template <
int dim,
int spacedim>
543 const unsigned int qpoint,
544 const unsigned int shape_nr)
551 template <
int dim,
int spacedim>
554 const unsigned int qpoint,
555 const unsigned int shape_nr)
const 563 template <
int dim,
int spacedim>
566 const unsigned int qpoint,
567 const unsigned int shape_nr)
575 template <
int dim,
int spacedim>
578 const unsigned int qpoint,
579 const unsigned int shape_nr)
const 587 template <
int dim,
int spacedim>
590 const unsigned int qpoint,
591 const unsigned int shape_nr)
600 template <
int dim,
int spacedim>
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 std::unique_ptr< typename Mapping< dim, spacedim >::InternalDataBase > get_data(const UpdateFlags, const Quadrature< dim > &quadrature) const override
std::vector< Tensor< 4, dim > > shape_fourth_derivatives
virtual std::unique_ptr< Mapping< dim, spacedim > > clone() const override
MappingFE(const FiniteElement< dim, spacedim > &fe)
std::vector< double > volume_elements
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
std::vector< std::vector< Tensor< 1, spacedim > > > aux
std::vector< double > shape_values
InternalData(const FiniteElement< dim, spacedim > &fe)
#define AssertIndexRange(index, range)
unsigned int get_degree() const
virtual BoundingBox< spacedim > get_bounding_box(const typename Triangulation< dim, spacedim >::cell_iterator &cell) const override
virtual bool preserves_vertex_locations() 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 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 Point< spacedim > transform_unit_to_real_cell(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const Point< dim > &p) const override
virtual std::unique_ptr< typename Mapping< dim, spacedim >::InternalDataBase > get_subface_data(const UpdateFlags flags, const Quadrature< dim - 1 > &quadrature) const override
std::vector< Tensor< 3, dim > > shape_third_derivatives
void initialize_face(const UpdateFlags update_flags, const Quadrature< dim > &quadrature, const unsigned int n_original_q_points)
std::vector< Tensor< 1, dim > > shape_derivatives
void initialize(const UpdateFlags update_flags, const Quadrature< dim > &quadrature, const unsigned int n_original_q_points)
Abstract base class for mapping classes.
virtual std::size_t memory_consumption() const override
#define DEAL_II_NAMESPACE_CLOSE
std::vector< Tensor< 2, dim > > shape_second_derivatives
Triangulation< dim, spacedim >::cell_iterator cell_of_current_support_points
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
std::vector< Point< spacedim > > mapping_support_points
const double & shape(const unsigned int qpoint, const unsigned int shape_nr) const
std::vector< DerivativeForm< 1, dim, spacedim > > covariant
const std::unique_ptr< FiniteElement< dim, spacedim > > fe
const Tensor< 3, dim > & third_derivative(const unsigned int qpoint, const unsigned int shape_nr) const
virtual std::vector< Point< spacedim > > compute_mapping_support_points(const typename Triangulation< dim, spacedim >::cell_iterator &cell) const
const Tensor< 2, dim > & second_derivative(const unsigned int qpoint, const unsigned int shape_nr) const
std::array< std::vector< Tensor< 1, dim > >, GeometryInfo< dim >::faces_per_cell *(dim - 1)> unit_tangentials
const Tensor< 1, dim > & derivative(const unsigned int qpoint, const unsigned int shape_nr) const
#define DEAL_II_NAMESPACE_OPEN
void compute_shape_function_values(const std::vector< Point< dim >> &unit_points)
virtual std::unique_ptr< typename Mapping< dim, spacedim >::InternalDataBase > get_face_data(const UpdateFlags flags, const hp::QCollection< dim - 1 > &quadrature) const override
const unsigned int polynomial_degree
std::vector< double > quadrature_weights
Table< 2, double > mapping_support_point_weights
const unsigned int n_shape_functions
const unsigned int polynomial_degree
const Tensor< 4, dim > & fourth_derivative(const unsigned int qpoint, const unsigned int shape_nr) const
std::vector< DerivativeForm< 1, dim, spacedim > > contravariant
virtual UpdateFlags requires_update_flags(const UpdateFlags update_flags) const override
const FiniteElement< dim, spacedim > & fe