16 #ifndef dealii_mapping_q_h
17 #define dealii_mapping_q_h
110 template <
int dim,
int spacedim = dim>
128 const bool use_mapping_q_on_all_cells);
136 virtual std::unique_ptr<Mapping<dim, spacedim>>
137 clone()
const override;
156 &cell)
const override;
294 const unsigned int n_original_q_points);
304 const unsigned int n_original_q_points);
329 shape(
const unsigned int qpoint,
const unsigned int shape_nr)
const;
335 shape(
const unsigned int qpoint,
const unsigned int shape_nr);
341 derivative(
const unsigned int qpoint,
const unsigned int shape_nr)
const;
347 derivative(
const unsigned int qpoint,
const unsigned int shape_nr);
354 const unsigned int shape_nr)
const;
367 const unsigned int shape_nr)
const;
380 const unsigned int shape_nr)
const;
445 std::array<std::vector<Tensor<1, dim>>,
493 (dim <= 2 ? 2 : 4))>;
536 mutable std::vector<AlignedVector<Tensor<1, spacedim>>>
aux;
563 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
569 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
574 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
586 &output_data)
const override;
594 const unsigned int face_no,
598 &output_data)
const override;
604 const unsigned int face_no,
605 const unsigned int subface_no,
609 &output_data)
const override;
618 &output_data)
const override;
685 const std::vector<Table<2, double>>
731 virtual std::vector<Point<spacedim>>
795 template <
int dim,
int spacedim = dim>
805 template <
int dim,
int spacedim>
806 inline const double &
808 const unsigned int shape_nr)
const
816 template <
int dim,
int spacedim>
819 const unsigned int shape_nr)
822 return shape_values[qpoint * n_shape_functions + shape_nr];
826 template <
int dim,
int spacedim>
829 const unsigned int qpoint,
830 const unsigned int shape_nr)
const
833 shape_derivatives.
size());
834 return shape_derivatives[qpoint * n_shape_functions + shape_nr];
839 template <
int dim,
int spacedim>
842 const unsigned int shape_nr)
845 shape_derivatives.
size());
846 return shape_derivatives[qpoint * n_shape_functions + shape_nr];
850 template <
int dim,
int spacedim>
853 const unsigned int qpoint,
854 const unsigned int shape_nr)
const
857 shape_second_derivatives.
size());
858 return shape_second_derivatives[qpoint * n_shape_functions + shape_nr];
862 template <
int dim,
int spacedim>
865 const unsigned int qpoint,
866 const unsigned int shape_nr)
869 shape_second_derivatives.
size());
870 return shape_second_derivatives[qpoint * n_shape_functions + shape_nr];
873 template <
int dim,
int spacedim>
876 const unsigned int qpoint,
877 const unsigned int shape_nr)
const
880 shape_third_derivatives.
size());
881 return shape_third_derivatives[qpoint * n_shape_functions + shape_nr];
885 template <
int dim,
int spacedim>
888 const unsigned int qpoint,
889 const unsigned int shape_nr)
892 shape_third_derivatives.
size());
893 return shape_third_derivatives[qpoint * n_shape_functions + shape_nr];
897 template <
int dim,
int spacedim>
900 const unsigned int qpoint,
901 const unsigned int shape_nr)
const
904 shape_fourth_derivatives.
size());
905 return shape_fourth_derivatives[qpoint * n_shape_functions + shape_nr];
909 template <
int dim,
int spacedim>
912 const unsigned int qpoint,
913 const unsigned int shape_nr)
916 shape_fourth_derivatives.
size());
917 return shape_fourth_derivatives[qpoint * n_shape_functions + shape_nr];
922 template <
int dim,
int spacedim>
Triangulation< dim, spacedim >::cell_iterator cell_of_current_support_points
AlignedVector< DerivativeForm< 1, dim, spacedim > > covariant
QGaussLobatto< 1 > line_support_points
AlignedVector< Tensor< 3, dim > > shape_third_derivatives
Tensor< 2, dim > & second_derivative(const unsigned int qpoint, const unsigned int shape_nr)
virtual std::size_t memory_consumption() const override
std::vector< Point< spacedim > > mapping_support_points
const Tensor< 1, dim > & derivative(const unsigned int qpoint, const unsigned int shape_nr) const
AlignedVector< Tensor< 4, dim > > shape_fourth_derivatives
const double & shape(const unsigned int qpoint, const unsigned int shape_nr) const
const unsigned int polynomial_degree
AlignedVector< double > shape_values
std::array< std::vector< Tensor< 1, dim > >, GeometryInfo< dim >::faces_per_cell *(dim - 1)> unit_tangentials
AlignedVector< DerivativeForm< 1, dim, spacedim > > contravariant
AlignedVector< Tensor< 1, dim > > shape_derivatives
bool tensor_product_quadrature
const Tensor< 3, dim > & third_derivative(const unsigned int qpoint, const unsigned int shape_nr) const
const Tensor< 2, dim > & second_derivative(const unsigned int qpoint, const unsigned int shape_nr) const
AlignedVector< Tensor< 2, dim > > shape_second_derivatives
std::vector< AlignedVector< Tensor< 1, spacedim > > > aux
double & shape(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)
InternalData(const unsigned int polynomial_degree)
void compute_shape_function_values(const std::vector< Point< dim >> &unit_points)
AlignedVector< double > volume_elements
void initialize(const UpdateFlags update_flags, const Quadrature< dim > &quadrature, const unsigned int n_original_q_points)
Tensor< 4, dim > & fourth_derivative(const unsigned int qpoint, const unsigned int shape_nr)
Tensor< 3, dim > & third_derivative(const unsigned int qpoint, const unsigned int shape_nr)
Tensor< 1, dim > & derivative(const unsigned int qpoint, const unsigned int shape_nr)
const Tensor< 4, dim > & fourth_derivative(const unsigned int qpoint, const unsigned int shape_nr) const
const unsigned int n_shape_functions
internal::MatrixFreeFunctions::ShapeInfo< VectorizedArrayType > shape_info
AlignedVector< VectorizedArrayType > scratch
virtual void add_line_support_points(const typename Triangulation< dim, spacedim >::cell_iterator &cell, std::vector< Point< spacedim >> &a) const
const std::vector< unsigned int > renumber_lexicographic_to_hierarchic
const Table< 2, double > support_point_weights_cell
virtual std::vector< Point< spacedim > > compute_mapping_support_points(const typename Triangulation< dim, spacedim >::cell_iterator &cell) const
virtual std::unique_ptr< typename Mapping< dim, spacedim >::InternalDataBase > get_data(const UpdateFlags, const Quadrature< dim > &quadrature) const override
void fill_mapping_data_for_generic_points(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const ArrayView< const Point< dim >> &unit_points, const UpdateFlags update_flags, internal::FEValuesImplementation::MappingRelatedData< dim, spacedim > &output_data) 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_subface_data(const UpdateFlags flags, const Quadrature< dim - 1 > &quadrature) 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 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
virtual void fill_fe_immersed_surface_values(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const NonMatching::ImmersedSurfaceQuadrature< 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 bool preserves_vertex_locations() 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 std::unique_ptr< typename Mapping< dim, spacedim >::InternalDataBase > get_face_data(const UpdateFlags flags, const hp::QCollection< dim - 1 > &quadrature) const override
Point< dim > transform_real_to_unit_cell_internal(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const Point< spacedim > &p, const Point< dim > &initial_p_unit) const
const std::vector< Table< 2, double > > support_point_weights_perimeter_to_interior
const std::vector< Point< 1 > > line_support_points
virtual Point< spacedim > transform_unit_to_real_cell(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const Point< dim > &p) const override
const std::vector< Polynomials::Polynomial< double > > polynomials_1d
virtual std::unique_ptr< Mapping< dim, spacedim > > clone() const override
const std::vector< Point< dim > > unit_cell_support_points
virtual void transform_points_real_to_unit_cell(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const ArrayView< const Point< spacedim >> &real_points, const ArrayView< Point< dim >> &unit_points) const override
virtual Point< dim > transform_real_to_unit_cell(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const Point< spacedim > &p) const override
MappingQ(const unsigned int polynomial_degree)
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
virtual bool is_compatible_with(const ReferenceCell &reference_cell) const override
virtual void add_quad_support_points(const typename Triangulation< dim, spacedim >::cell_iterator &cell, std::vector< Point< spacedim >> &a) const
unsigned int get_degree() const
Abstract base class for mapping classes.
#define DEAL_II_DEPRECATED
#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)