22#ifdef DEAL_II_BOOST_HAS_BROKEN_HEADER_DEPRECATIONS
23# define BOOST_ALLOW_DEPRECATED_HEADERS
25#include <boost/geometry.hpp>
26#ifdef DEAL_II_BOOST_HAS_BROKEN_HEADER_DEPRECATIONS
27# undef BOOST_ALLOW_DEPRECATED_HEADERS
35template <
int dim,
int spacedim>
36boost::container::small_vector<Point<spacedim>,
41 boost::container::small_vector<Point<spacedim>,
45 vertices.push_back(cell->
vertex(i));
52template <
int dim,
int spacedim>
53boost::container::small_vector<Point<spacedim>,
57 const unsigned int face_no)
const
59 boost::container::small_vector<Point<spacedim>,
63 const auto &cell_vertices = get_vertices(cell);
67 for (
const unsigned int v :
70 face_vertices.push_back(
72 face_no, v, face_orientation)]);
80template <
int dim,
int spacedim>
84 const bool map_barycenter_of_reference_cell)
const
86 if (map_barycenter_of_reference_cell)
88 return transform_unit_to_real_cell(
93 const auto vertices = get_vertices(cell);
95 for (
const auto &v : vertices)
103template <
int dim,
int spacedim>
108 if (preserves_vertex_locations())
116template <
int dim,
int spacedim>
129template <
int dim,
int spacedim>
137 for (
unsigned int i = 0; i < real_points.size(); ++i)
146 unit_points[i][0] = std::numeric_limits<double>::infinity();
153template <
int dim,
int spacedim>
157 const unsigned int face_no,
173 return Point<dim - 1>{unit_cell_pt[1]};
175 return Point<dim - 1>{unit_cell_pt[0]};
180 return Point<dim - 1>{unit_cell_pt[1], unit_cell_pt[2]};
182 return Point<dim - 1>{unit_cell_pt[0], unit_cell_pt[2]};
184 return Point<dim - 1>{unit_cell_pt[0], unit_cell_pt[1]};
195template <
int dim,
int spacedim>
199 const unsigned int face_no,
207 fill_fe_face_values(cell, face_no, quadrature[0], internal_data, output_data);
212template <
int dim,
int spacedim>
216 const unsigned int face_no,
223 ExcMessage(
"Use of a deprecated interface, please implement "
224 "fill_fe_face_values taking a hp::QCollection argument"));
234template <
int dim,
int spacedim>
235std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
241 return get_face_data(update_flags, quadrature[0]);
246template <
int dim,
int spacedim>
247std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
253 ExcMessage(
"Use of a deprecated interface, please implement "
254 "fill_fe_face_values taking a hp::QCollection argument"));
258 return std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>();
265template <
int dim,
int spacedim>
272template <
int dim,
int spacedim>
282template <
int dim,
int spacedim>
286 return sizeof(*this);
292template <
int dim,
int spacedim>
293const Mapping<dim, spacedim> &
297 Assert(reference_cells.size() == 1,
299 "This function can only work for triangulations that "
300 "use only a single cell type -- for example, only triangles "
301 "or only quadrilaterals. For mixed meshes, there is no "
302 "single linear mapping object that can be used for all "
303 "cells of the triangulation. The triangulation you are "
304 "passing to this function uses multiple cell types."));
306 return reference_cells.front()
313#include "mapping.inst"
virtual std::size_t memory_consumption() const
virtual void reinit(const UpdateFlags update_flags, const Quadrature< dim > &quadrature)
Point< dim - 1 > project_real_point_to_unit_point_on_face(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const unsigned int face_no, const Point< spacedim > &p) const
virtual BoundingBox< spacedim > get_bounding_box(const typename Triangulation< dim, spacedim >::cell_iterator &cell) const
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
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
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
virtual boost::container::small_vector< Point< spacedim >, GeometryInfo< dim >::vertices_per_cell > get_vertices(const typename Triangulation< dim, spacedim >::cell_iterator &cell) const
virtual Point< spacedim > get_center(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const bool map_barycenter_of_reference_cell=true) const
virtual std::unique_ptr< InternalDataBase > get_face_data(const UpdateFlags update_flags, const hp::QCollection< dim - 1 > &quadrature) const
std_cxx20::ranges::iota_view< unsigned int, unsigned int > vertex_indices() const
unsigned int n_vertices() const
BoundingBox< spacedim > bounding_box() const
Point< spacedim > & vertex(const unsigned int i) const
ReferenceCell reference_cell() const
const std::vector< ReferenceCell > & get_reference_cells() const
unsigned int size() const
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
constexpr std::array< unsigned int, GeometryInfo< dim >::faces_per_cell > GeometryInfo< dim >::unit_normal_direction
static ::ExceptionBase & ExcTransformationFailed()
static ::ExceptionBase & ExcNotImplemented()
#define Assert(cond, exc)
#define AssertDimension(dim1, dim2)
static ::ExceptionBase & ExcMessage(std::string arg1)
#define AssertThrow(cond, exc)
TriaIterator< CellAccessor< dim, spacedim > > cell_iterator
@ update_default
No update.
unsigned char combined_face_orientation(const unsigned int face) const
#define DEAL_II_ASSERT_UNREACHABLE()
const Mapping< dim, spacedim > & get_default_linear_mapping(const Triangulation< dim, spacedim > &triangulation)
void reference_cell(Triangulation< dim, spacedim > &tria, const ReferenceCell &reference_cell)
Point< 1 > transform_real_to_unit_cell(const std::array< Point< spacedim >, GeometryInfo< 1 >::vertices_per_cell > &vertices, const Point< spacedim > &p)
static constexpr unsigned int vertices_per_cell
static constexpr std::array< unsigned int, faces_per_cell > unit_normal_direction
static constexpr unsigned int vertices_per_face