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>,
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);
65 const auto face_orientation = cell->combined_face_orientation(face_no);
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(
89 cell, cell->reference_cell().template barycenter<dim>());
93 const auto vertices = get_vertices(cell);
97 return center / cell->n_vertices();
103template <
int dim,
int spacedim>
108 if (preserves_vertex_locations())
109 return cell->bounding_box();
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,
167 const unsigned int unit_normal_direction =
172 if (unit_normal_direction == 0)
173 return Point<dim - 1>{unit_cell_pt[1]};
174 else if (unit_normal_direction == 1)
175 return Point<dim - 1>{unit_cell_pt[0]};
179 if (unit_normal_direction == 0)
180 return Point<dim - 1>{unit_cell_pt[1], unit_cell_pt[2]};
181 else if (unit_normal_direction == 1)
182 return Point<dim - 1>{unit_cell_pt[0], unit_cell_pt[2]};
183 else if (unit_normal_direction == 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>
296 const auto &reference_cells =
triangulation.get_reference_cells();
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()
307 .template get_default_linear_mapping<dim, spacedim>();
313#include "mapping.inst"
virtual std::size_t memory_consumption() const
virtual void reinit(const UpdateFlags update_flags, const Quadrature< dim > &quadrature)
Abstract base class for mapping classes.
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
unsigned int size() const
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
static ::ExceptionBase & ExcNotImplemented()
#define Assert(cond, exc)
#define AssertDimension(dim1, dim2)
static ::ExceptionBase & ExcMessage(std::string arg1)
#define AssertThrow(cond, exc)
@ update_default
No update.
#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)
const ::parallel::distributed::Triangulation< dim, spacedim > * triangulation