16 #ifndef dealii_mapping_q_generic_h 17 #define dealii_mapping_q_generic_h 134 template <
int dim,
int spacedim = dim>
151 virtual std::unique_ptr<Mapping<dim, spacedim>>
152 clone()
const override;
171 &cell)
const override;
252 template <
typename Number,
typename Number2>
301 const unsigned int n_original_q_points);
311 const unsigned int n_original_q_points);
336 shape(
const unsigned int qpoint,
const unsigned int shape_nr)
const;
342 shape(
const unsigned int qpoint,
const unsigned int shape_nr);
348 derivative(
const unsigned int qpoint,
const unsigned int shape_nr)
const;
354 derivative(
const unsigned int qpoint,
const unsigned int shape_nr);
361 const unsigned int shape_nr)
const;
374 const unsigned int shape_nr)
const;
387 const unsigned int shape_nr)
const;
452 std::array<std::vector<Tensor<1, dim>>,
535 mutable std::vector<DerivativeForm<1, dim, spacedim>>
covariant;
549 mutable std::vector<std::vector<Tensor<1, spacedim>>>
aux;
576 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
582 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
587 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
599 &output_data)
const override;
607 const unsigned int face_no,
611 &output_data)
const override;
617 const unsigned int face_no,
618 const unsigned int subface_no,
622 &output_data)
const override;
690 const std::vector<Table<2, double>>
736 virtual std::vector<Point<spacedim>>
806 template <
int dim,
int spacedim>
807 inline const double &
809 const unsigned int qpoint,
810 const unsigned int shape_nr)
const 818 template <
int dim,
int spacedim>
821 const unsigned int shape_nr)
828 template <
int dim,
int spacedim>
831 const unsigned int qpoint,
832 const unsigned int shape_nr)
const 841 template <
int dim,
int spacedim>
844 const unsigned int qpoint,
845 const unsigned int shape_nr)
853 template <
int dim,
int spacedim>
856 const unsigned int qpoint,
857 const unsigned int shape_nr)
const 865 template <
int dim,
int spacedim>
868 const unsigned int qpoint,
869 const unsigned int shape_nr)
876 template <
int dim,
int spacedim>
879 const unsigned int qpoint,
880 const unsigned int shape_nr)
const 888 template <
int dim,
int spacedim>
891 const unsigned int qpoint,
892 const unsigned int shape_nr)
900 template <
int dim,
int spacedim>
903 const unsigned int qpoint,
904 const unsigned int shape_nr)
const 912 template <
int dim,
int spacedim>
915 const unsigned int qpoint,
916 const unsigned int shape_nr)
925 template <
int dim,
int spacedim>
934 template <
int dim,
int spacedim>
935 template <
typename Number,
typename Number2>
947 const std::vector<Point<spacedim>> support_points =
950 const unsigned int n_points = unit_points.size();
958 for (
unsigned int i = 0; i < n_points; i += n_lanes)
959 if (n_points - i > 1)
962 for (
unsigned int j = 0; j < n_lanes; ++j)
963 if (i + j < n_points)
964 for (
unsigned int d = 0;
d < spacedim; ++
d)
965 p_vec[
d][j] = unit_points[i + j][
d];
967 for (
unsigned int d = 0;
d < spacedim; ++
d)
968 p_vec[
d][j] = unit_points[i][
d];
981 for (
unsigned int j = 0; j < n_lanes && i + j < n_points; ++j)
984 for (
unsigned int d = 0;
d < spacedim; ++
d)
985 for (
unsigned int e = 0;
e < dim; ++
e)
986 jac_j[
d][
e] = jac[
d][
e][j];
std::vector< Tensor< 2, dim > > shape_second_derivatives
#define AssertDimension(dim1, dim2)
const unsigned int polynomial_degree
virtual void add_quad_support_points(const typename Triangulation< dim, spacedim >::cell_iterator &cell, std::vector< Point< spacedim >> &a) const
AlignedVector< VectorizedArray< double > > values_quad
SymmetricTensor< 2, dim, Number > e(const Tensor< 2, dim, Number > &F)
const Tensor< 1, dim > & derivative(const unsigned int qpoint, const unsigned int shape_nr) const
virtual UpdateFlags requires_update_flags(const UpdateFlags update_flags) const override
virtual void add_line_support_points(const typename Triangulation< dim, spacedim >::cell_iterator &cell, std::vector< Point< spacedim >> &a) const
const Tensor< 3, dim > & third_derivative(const unsigned int qpoint, const unsigned int shape_nr) const
virtual Point< spacedim > transform_unit_to_real_cell(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const Point< dim > &p) const override
#define AssertIndexRange(index, range)
virtual std::unique_ptr< Mapping< dim, spacedim > > clone() const override
const Table< 2, double > support_point_weights_cell
AlignedVector< VectorizedArray< double > > gradients_quad
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
MappingQGeneric(const unsigned int polynomial_degree)
bool tensor_product_quadrature
AlignedVector< VectorizedArray< double > > scratch
InternalData(const unsigned int polynomial_degree)
const unsigned int polynomial_degree
const std::vector< Polynomials::Polynomial< double > > polynomials_1d
QGaussLobatto< 1 > line_support_points
virtual std::unique_ptr< typename Mapping< dim, spacedim >::InternalDataBase > get_face_data(const UpdateFlags flags, const hp::QCollection< dim - 1 > &quadrature) const override
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 bool is_compatible_with(const ReferenceCell &reference_cell) const override
void compute_shape_function_values(const std::vector< Point< dim >> &unit_points)
std::vector< std::vector< Tensor< 1, spacedim > > > aux
#define Assert(cond, exc)
const std::vector< Point< dim > > unit_cell_support_points
void reference_cell(const ReferenceCell &reference_cell, Triangulation< dim, spacedim > &tria)
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
Abstract base class for mapping classes.
void transform_variable(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const MappingKind kind, const ArrayView< const Point< dim >> &unit_points, const ArrayView< const Number > &input, const ArrayView< Number2 > &output) const
#define DEAL_II_NAMESPACE_CLOSE
std::vector< Point< spacedim > > mapping_support_points
std::vector< Tensor< 3, dim > > shape_third_derivatives
virtual std::unique_ptr< typename Mapping< dim, spacedim >::InternalDataBase > get_data(const UpdateFlags, const Quadrature< dim > &quadrature) const override
static constexpr std::size_t size()
void initialize(const UpdateFlags update_flags, const Quadrature< dim > &quadrature, const unsigned int n_original_q_points)
std::vector< double > volume_elements
SymmetricTensor< 2, dim, Number > d(const Tensor< 2, dim, Number > &F, const Tensor< 2, dim, Number > &dF_dt)
const unsigned int n_shape_functions
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< 1, dim > > shape_derivatives
const std::vector< Table< 2, double > > support_point_weights_perimeter_to_interior
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
AlignedVector< VectorizedArray< double > > hessians_quad
AlignedVector< VectorizedArray< double > > values_dofs
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
unsigned int get_degree() const
const double & shape(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 std::vector< Point< 1 > > line_support_points
std::pair< typename ProductTypeNoPoint< Number, Number2 >::type, Tensor< 1, dim, typename ProductTypeNoPoint< Number, Number2 >::type > > evaluate_tensor_product_value_and_gradient(const std::vector< Polynomials::Polynomial< double >> &poly, const std::vector< Number > &values, const Point< dim, Number2 > &p, const bool d_linear=false, const std::vector< unsigned int > &renumber={})
#define DEAL_II_NAMESPACE_OPEN
virtual bool preserves_vertex_locations() const override
const Tensor< 2, dim > & second_derivative(const unsigned int qpoint, const unsigned int shape_nr) const
Triangulation< dim, spacedim >::cell_iterator cell_of_current_support_points
virtual std::size_t memory_consumption() const override
virtual std::vector< Point< spacedim > > compute_mapping_support_points(const typename Triangulation< dim, spacedim >::cell_iterator &cell) const
const std::vector< unsigned int > renumber_lexicographic_to_hierarchic
static ::ExceptionBase & ExcNotImplemented()
std::vector< DerivativeForm< 1, dim, spacedim > > contravariant
virtual BoundingBox< spacedim > get_bounding_box(const typename Triangulation< dim, spacedim >::cell_iterator &cell) const override
std::vector< Tensor< 4, dim > > shape_fourth_derivatives
void initialize_face(const UpdateFlags update_flags, const Quadrature< dim > &quadrature, const unsigned int n_original_q_points)
std::vector< double > shape_values
virtual Point< dim > transform_real_to_unit_cell(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const Point< spacedim > &p) const override
const Tensor< 4, dim > & fourth_derivative(const unsigned int qpoint, const unsigned int shape_nr) const
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
std::vector< DerivativeForm< 1, dim, spacedim > > covariant
internal::MatrixFreeFunctions::ShapeInfo< VectorizedArray< double > > shape_info