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;
244 template <
typename Number,
typename Number2>
293 const unsigned int n_original_q_points);
303 const unsigned int n_original_q_points);
328 shape(
const unsigned int qpoint,
const unsigned int shape_nr)
const;
334 shape(
const unsigned int qpoint,
const unsigned int shape_nr);
340 derivative(
const unsigned int qpoint,
const unsigned int shape_nr)
const;
346 derivative(
const unsigned int qpoint,
const unsigned int shape_nr);
353 const unsigned int shape_nr)
const;
366 const unsigned int shape_nr)
const;
379 const unsigned int shape_nr)
const;
444 std::array<std::vector<Tensor<1, dim>>,
527 mutable std::vector<DerivativeForm<1, dim, spacedim>>
covariant;
541 mutable std::vector<std::vector<Tensor<1, spacedim>>>
aux;
568 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
574 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
579 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
591 &output_data)
const override;
599 const unsigned int face_no,
603 &output_data)
const override;
609 const unsigned int face_no,
610 const unsigned int subface_no,
614 &output_data)
const override;
682 const std::vector<Table<2, double>>
728 virtual std::vector<Point<spacedim>>
798 template <
int dim,
int spacedim>
799 inline const double &
801 const unsigned int qpoint,
802 const unsigned int shape_nr)
const 810 template <
int dim,
int spacedim>
813 const unsigned int shape_nr)
820 template <
int dim,
int spacedim>
823 const unsigned int qpoint,
824 const unsigned int shape_nr)
const 833 template <
int dim,
int spacedim>
836 const unsigned int qpoint,
837 const unsigned int shape_nr)
845 template <
int dim,
int spacedim>
848 const unsigned int qpoint,
849 const unsigned int shape_nr)
const 857 template <
int dim,
int spacedim>
860 const unsigned int qpoint,
861 const unsigned int shape_nr)
868 template <
int dim,
int spacedim>
871 const unsigned int qpoint,
872 const unsigned int shape_nr)
const 880 template <
int dim,
int spacedim>
883 const unsigned int qpoint,
884 const unsigned int shape_nr)
892 template <
int dim,
int spacedim>
895 const unsigned int qpoint,
896 const unsigned int shape_nr)
const 904 template <
int dim,
int spacedim>
907 const unsigned int qpoint,
908 const unsigned int shape_nr)
917 template <
int dim,
int spacedim>
926 template <
int dim,
int spacedim>
927 template <
typename Number,
typename Number2>
939 const std::vector<Point<spacedim>> support_points =
942 const unsigned int n_points = unit_points.size();
950 for (
unsigned int i = 0; i < n_points; i += n_lanes)
951 if (n_points - i > 1)
954 for (
unsigned int j = 0; j < n_lanes; ++j)
955 if (i + j < n_points)
956 for (
unsigned int d = 0;
d < spacedim; ++
d)
957 p_vec[
d][j] = unit_points[i + j][
d];
959 for (
unsigned int d = 0;
d < spacedim; ++
d)
960 p_vec[
d][j] = unit_points[i][
d];
973 for (
unsigned int j = 0; j < n_lanes && i + j < n_points; ++j)
976 for (
unsigned int d = 0;
d < spacedim; ++
d)
977 for (
unsigned int e = 0;
e < dim; ++
e)
978 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
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
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
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