16 #include <deal.II/base/array_view.h> 17 #include <deal.II/base/memory_consumption.h> 18 #include <deal.II/base/polynomial.h> 19 #include <deal.II/base/quadrature.h> 20 #include <deal.II/base/quadrature_lib.h> 21 #include <deal.II/base/std_cxx14/memory.h> 22 #include <deal.II/base/tensor_product_polynomials.h> 23 #include <deal.II/base/utilities.h> 25 #include <deal.II/dofs/dof_accessor.h> 27 #include <deal.II/fe/fe_q.h> 28 #include <deal.II/fe/fe_values.h> 29 #include <deal.II/fe/mapping_q.h> 31 #include <deal.II/grid/tria_iterator.h> 33 #include <deal.II/lac/full_matrix.h> 38 DEAL_II_NAMESPACE_OPEN
41 template <
int dim,
int spacedim>
43 : use_mapping_q1_on_current_cell(false)
48 template <
int dim,
int spacedim>
61 template <
int dim,
int spacedim>
72 use_mapping_q_on_all_cells(degree == 1 || use_mapping_q_on_all_cells ||
89 template <
int dim,
int spacedim>
96 std::shared_ptr<const Mapping<dim, spacedim>> other_q1_map =
110 std::shared_ptr<const Mapping<dim, spacedim>> other_qp_map =
121 template <
int dim,
int spacedim>
130 template <
int dim,
int spacedim>
139 template <
int dim,
int spacedim>
143 return (
q1_mapping->requires_update_flags(in) |
149 template <
int dim,
int spacedim>
150 std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
154 auto data = std_cxx14::make_unique<InternalData>();
158 std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>>
167 std::move(
q1_mapping->get_data(update_flags, quadrature)));
171 typename MappingQGeneric<dim, spacedim>::InternalData>(
172 std::move(do_get_data.return_value()));
173 return std::move(data);
178 template <
int dim,
int spacedim>
179 std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
184 auto data = std_cxx14::make_unique<InternalData>();
188 std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>>
198 std::move(
q1_mapping->get_face_data(update_flags, quadrature)));
202 typename MappingQGeneric<dim, spacedim>::InternalData>(
203 std::move(do_get_data.return_value()));
204 return std::move(data);
209 template <
int dim,
int spacedim>
210 std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
215 auto data = std_cxx14::make_unique<InternalData>();
219 std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>>
229 std::move(
q1_mapping->get_subface_data(update_flags, quadrature)));
233 typename MappingQGeneric<dim, spacedim>::InternalData>(
234 std::move(do_get_data.return_value()));
235 return std::move(data);
241 template <
int dim,
int spacedim>
253 Assert(dynamic_cast<const InternalData *>(&internal_data) !=
nullptr,
279 updated_cell_similarity,
285 updated_cell_similarity,
290 return updated_cell_similarity;
295 template <
int dim,
int spacedim>
299 const unsigned int face_no,
307 Assert(dynamic_cast<const InternalData *>(&internal_data) !=
nullptr,
331 template <
int dim,
int spacedim>
335 const unsigned int face_no,
336 const unsigned int subface_no,
344 Assert(dynamic_cast<const InternalData *>(&internal_data) !=
nullptr,
377 template <
int dim,
int spacedim>
391 if (data->use_mapping_q1_on_current_cell)
392 q1_mapping->transform(input, mapping_type, *data->mapping_q1_data, output);
395 qp_mapping->transform(input, mapping_type, *data->mapping_qp_data, output);
400 template <
int dim,
int spacedim>
410 &mapping_data) !=
nullptr),
424 template <
int dim,
int spacedim>
434 &mapping_data) !=
nullptr),
448 template <
int dim,
int spacedim>
458 &mapping_data) !=
nullptr),
472 template <
int dim,
int spacedim>
482 &mapping_data) !=
nullptr),
496 template <
int dim,
int spacedim>
506 return qp_mapping->transform_unit_to_real_cell(cell, p);
508 return q1_mapping->transform_unit_to_real_cell(cell, p);
513 template <
int dim,
int spacedim>
521 return qp_mapping->transform_real_to_unit_cell(cell, p);
523 return q1_mapping->transform_real_to_unit_cell(cell, p);
528 template <
int dim,
int spacedim>
529 std::unique_ptr<Mapping<dim, spacedim>>
532 return std_cxx14::make_unique<MappingQ<dim, spacedim>>(
539 #include "mapping_q.inst" 542 DEAL_II_NAMESPACE_CLOSE
virtual std::size_t memory_consumption() const override
#define AssertDimension(dim1, dim2)
std::shared_ptr< const MappingQGeneric< dim, spacedim > > qp_mapping
unsigned int get_degree() const
virtual std::unique_ptr< typename Mapping< dim, spacedim >::InternalDataBase > get_data(const UpdateFlags, const Quadrature< dim > &quadrature) const override
Task< RT > new_task(const std::function< RT()> &function)
bool use_mapping_q1_on_current_cell
const bool use_mapping_q_on_all_cells
std::unique_ptr< To > dynamic_unique_cast(std::unique_ptr< From > &&p)
#define Assert(cond, exc)
Abstract base class for mapping classes.
const unsigned int polynomial_degree
std::unique_ptr< typename MappingQGeneric< dim, spacedim >::InternalData > mapping_qp_data
virtual Point< spacedim > transform_unit_to_real_cell(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const Point< dim > &p) const override
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 std::unique_ptr< typename Mapping< dim, spacedim >::InternalDataBase > get_face_data(const UpdateFlags flags, const Quadrature< dim-1 > &quadrature) const override
std::unique_ptr< typename MappingQGeneric< dim, spacedim >::InternalData > mapping_q1_data
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
MappingQ(const unsigned int polynomial_degree, const bool use_mapping_q_on_all_cells=false)
virtual Point< dim > transform_real_to_unit_cell(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const Point< spacedim > &p) const override
std::enable_if< std::is_fundamental< T >::value, std::size_t >::type memory_consumption(const T &t)
static::ExceptionBase & ExcInternalError()
std::shared_ptr< const MappingQGeneric< dim, spacedim > > q1_mapping