16 #ifndef dealii_mapping_fe_h
17 #define dealii_mapping_fe_h
58 template <
int dim,
int spacedim = dim>
73 virtual std::unique_ptr<Mapping<dim, spacedim>>
74 clone()
const override;
86 &cell)
const override;
199 const unsigned int n_original_q_points);
209 const unsigned int n_original_q_points);
224 shape(
const unsigned int qpoint,
const unsigned int shape_nr)
const;
230 shape(
const unsigned int qpoint,
const unsigned int shape_nr);
236 derivative(
const unsigned int qpoint,
const unsigned int shape_nr)
const;
242 derivative(
const unsigned int qpoint,
const unsigned int shape_nr);
249 const unsigned int shape_nr)
const;
262 const unsigned int shape_nr)
const;
275 const unsigned int shape_nr)
const;
333 std::array<std::vector<Tensor<1, dim>>,
361 mutable std::vector<DerivativeForm<1, dim, spacedim>>
covariant;
375 mutable std::vector<std::vector<Tensor<1, spacedim>>>
aux;
407 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
413 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
418 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
430 &output_data)
const override;
438 const unsigned int face_no,
442 &output_data)
const override;
448 const unsigned int face_no,
449 const unsigned int subface_no,
453 &output_data)
const override;
460 const std::unique_ptr<FiniteElement<dim, spacedim>>
fe;
471 virtual std::vector<Point<spacedim>>
487 template <
int dim,
int spacedim>
488 inline const double &
490 const unsigned int shape_nr)
const
498 template <
int dim,
int spacedim>
501 const unsigned int shape_nr)
503 AssertIndexRange(qpoint * n_shape_functions + shape_nr, shape_values.size());
504 return shape_values[qpoint * n_shape_functions + shape_nr];
508 template <
int dim,
int spacedim>
511 const unsigned int qpoint,
512 const unsigned int shape_nr)
const
515 shape_derivatives.size());
516 return shape_derivatives[qpoint * n_shape_functions + shape_nr];
521 template <
int dim,
int spacedim>
524 const unsigned int shape_nr)
527 shape_derivatives.size());
528 return shape_derivatives[qpoint * n_shape_functions + shape_nr];
532 template <
int dim,
int spacedim>
535 const unsigned int qpoint,
536 const unsigned int shape_nr)
const
539 shape_second_derivatives.size());
540 return shape_second_derivatives[qpoint * n_shape_functions + shape_nr];
544 template <
int dim,
int spacedim>
547 const unsigned int qpoint,
548 const unsigned int shape_nr)
551 shape_second_derivatives.size());
552 return shape_second_derivatives[qpoint * n_shape_functions + shape_nr];
555 template <
int dim,
int spacedim>
558 const unsigned int qpoint,
559 const unsigned int shape_nr)
const
562 shape_third_derivatives.size());
563 return shape_third_derivatives[qpoint * n_shape_functions + shape_nr];
567 template <
int dim,
int spacedim>
570 const unsigned int qpoint,
571 const unsigned int shape_nr)
574 shape_third_derivatives.size());
575 return shape_third_derivatives[qpoint * n_shape_functions + shape_nr];
579 template <
int dim,
int spacedim>
582 const unsigned int qpoint,
583 const unsigned int shape_nr)
const
586 shape_fourth_derivatives.size());
587 return shape_fourth_derivatives[qpoint * n_shape_functions + shape_nr];
591 template <
int dim,
int spacedim>
594 const unsigned int qpoint,
595 const unsigned int shape_nr)
598 shape_fourth_derivatives.size());
599 return shape_fourth_derivatives[qpoint * n_shape_functions + shape_nr];
604 template <
int dim,
int spacedim>
Triangulation< dim, spacedim >::cell_iterator cell_of_current_support_points
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< 1, dim > & 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)
std::vector< std::vector< Tensor< 1, spacedim > > > aux
std::vector< Tensor< 3, dim > > shape_third_derivatives
std::vector< Tensor< 2, dim > > shape_second_derivatives
const Tensor< 3, dim > & third_derivative(const unsigned int qpoint, const unsigned int shape_nr) const
std::vector< Tensor< 4, dim > > shape_fourth_derivatives
void initialize(const UpdateFlags update_flags, const Quadrature< dim > &quadrature, const unsigned int n_original_q_points)
Tensor< 1, dim > & derivative(const unsigned int qpoint, const unsigned int shape_nr)
Tensor< 3, dim > & third_derivative(const unsigned int qpoint, const unsigned int shape_nr)
std::array< std::vector< Tensor< 1, dim > >, GeometryInfo< dim >::faces_per_cell *(dim - 1)> unit_tangentials
std::vector< double > quadrature_weights
const double & shape(const unsigned int qpoint, const unsigned int shape_nr) const
const FiniteElement< dim, spacedim > & fe
std::vector< DerivativeForm< 1, dim, spacedim > > contravariant
std::vector< DerivativeForm< 1, dim, spacedim > > covariant
const Tensor< 4, dim > & fourth_derivative(const unsigned int qpoint, const unsigned int shape_nr) const
void compute_shape_function_values(const std::vector< Point< dim >> &unit_points)
std::vector< double > shape_values
std::vector< Tensor< 1, dim > > shape_derivatives
Tensor< 4, dim > & fourth_derivative(const unsigned int qpoint, const unsigned int shape_nr)
virtual std::size_t memory_consumption() const override
double & shape(const unsigned int qpoint, const unsigned int shape_nr)
InternalData(const FiniteElement< dim, spacedim > &fe)
const unsigned int n_shape_functions
std::vector< double > volume_elements
const Tensor< 2, dim > & second_derivative(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
MappingFE(const FiniteElement< dim, spacedim > &fe)
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 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 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 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)
void reference_cell(Triangulation< dim, spacedim > &tria, const ReferenceCell &reference_cell)