16 #ifndef dealii_mapping_manifold_h 17 #define dealii_mapping_manifold_h 20 #include <deal.II/base/config.h> 22 #include <deal.II/base/derivative_form.h> 23 #include <deal.II/base/qprojector.h> 24 #include <deal.II/base/quadrature_lib.h> 25 #include <deal.II/base/table.h> 27 #include <deal.II/dofs/dof_accessor.h> 29 #include <deal.II/fe/mapping.h> 31 #include <deal.II/grid/tria_iterator.h> 35 DEAL_II_NAMESPACE_OPEN
64 template <
int dim,
int spacedim = dim>
79 virtual std::unique_ptr<Mapping<dim, spacedim>>
80 clone()
const override;
190 const unsigned int n_original_q_points);
200 const unsigned int n_original_q_points);
296 std::array<std::vector<Tensor<1, dim>>,
309 mutable std::vector<DerivativeForm<1, dim, spacedim>>
covariant;
323 mutable std::vector<std::vector<Tensor<1, spacedim>>>
aux;
345 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
349 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
354 virtual std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
366 &output_data)
const override;
372 const unsigned int face_no,
376 &output_data)
const override;
382 const unsigned int face_no,
383 const unsigned int subface_no,
387 &output_data)
const override;
402 template <
int dim,
int spacedim>
408 for (
unsigned int i = 0; i < GeometryInfo<dim>::vertices_per_cell; ++i)
414 template <
int dim,
int spacedim>
421 for (
unsigned int q = 0; q < quad.
size(); ++q)
423 for (
unsigned int i = 0; i < GeometryInfo<dim>::vertices_per_cell; ++i)
433 template <
int dim,
int spacedim>
445 DEAL_II_NAMESPACE_CLOSE
std::vector< DerivativeForm< 1, dim, spacedim > > covariant
virtual bool preserves_vertex_locations() const override
virtual void transform(const ArrayView< const Tensor< 1, dim >> &input, const MappingType type, const typename Mapping< dim, spacedim >::InternalDataBase &internal, const ArrayView< Tensor< 1, spacedim >> &output) const override
virtual std::unique_ptr< Mapping< dim, spacedim > > clone() 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
SmartPointer< const Manifold< dim, spacedim > > manifold
virtual void fill_fe_face_values(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const unsigned int face_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 > > contravariant
std::array< std::vector< Tensor< 1, dim > >, GeometryInfo< dim >::faces_per_cell *(dim-1)> unit_tangentials
virtual std::size_t memory_consumption() const override
Triangulation< dim, spacedim >::cell_iterator cell
const Point< dim > & point(const unsigned int i) const
std::vector< Point< spacedim > > vertices
virtual std::unique_ptr< typename Mapping< dim, spacedim >::InternalDataBase > get_subface_data(const UpdateFlags flags, const Quadrature< dim-1 > &quadrature) const override
virtual UpdateFlags requires_update_flags(const UpdateFlags update_flags) const override
void store_vertices(const typename Triangulation< dim, spacedim >::cell_iterator &cell) 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< std::vector< Tensor< 1, spacedim > > > aux
unsigned int size() const
Abstract base class for mapping classes.
void initialize(const UpdateFlags update_flags, const Quadrature< dim > &quadrature, const unsigned int n_original_q_points)
std::vector< double > volume_elements
std::vector< double > vertex_weights
void compute_manifold_quadrature_weights(const Quadrature< dim > &quadrature)
static double d_linear_shape_function(const Point< dim > &xi, const unsigned int i)
virtual Point< spacedim > transform_unit_to_real_cell(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const Point< dim > &p) 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 std::unique_ptr< typename Mapping< dim, spacedim >::InternalDataBase > get_data(const UpdateFlags, const Quadrature< dim > &quadrature) const override
void initialize_face(const UpdateFlags update_flags, const Quadrature< dim > &quadrature, const unsigned int n_original_q_points)
virtual std::unique_ptr< typename Mapping< dim, spacedim >::InternalDataBase > get_face_data(const UpdateFlags flags, const Quadrature< dim-1 > &quadrature) const override
std::vector< std::vector< double > > cell_manifold_quadrature_weights
MappingManifold()=default