40template <
int dim,
int spacedim>
42 : use_mapping_q1_on_current_cell(false)
47template <
int dim,
int spacedim>
60template <
int dim,
int spacedim>
62 const bool use_mapping_q_on_all_cells)
63 : polynomial_degree(degree)
71 use_mapping_q_on_all_cells(degree == 1 || use_mapping_q_on_all_cells ||
81 qp_mapping(this->polynomial_degree > 1 ?
88template <
int dim,
int spacedim>
90 : polynomial_degree(mapping.polynomial_degree)
91 , use_mapping_q_on_all_cells(mapping.use_mapping_q_on_all_cells)
95 std::shared_ptr<const Mapping<dim, spacedim>> other_q1_map =
97 q1_mapping = std::dynamic_pointer_cast<const MappingQGeneric<dim, spacedim>>(
103 if (this->polynomial_degree == 1)
105 qp_mapping = q1_mapping;
109 std::shared_ptr<const Mapping<dim, spacedim>> other_qp_map =
112 std::dynamic_pointer_cast<const MappingQGeneric<dim, spacedim>>(
120template <
int dim,
int spacedim>
129template <
int dim,
int spacedim>
138template <
int dim,
int spacedim>
142 return (
q1_mapping->requires_update_flags(in) |
148template <
int dim,
int spacedim>
149std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
153 std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase> data_ptr =
154 std::make_unique<InternalData>();
155 auto &data =
dynamic_cast<InternalData &
>(*data_ptr);
159 std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>>
168 std::move(
q1_mapping->get_data(update_flags, quadrature)));
173 std::move(do_get_data.return_value()));
179template <
int dim,
int spacedim>
180std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
185 std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase> data_ptr =
186 std::make_unique<InternalData>();
187 auto &data =
dynamic_cast<InternalData &
>(*data_ptr);
189 std::unique_ptr<typename MappingQGeneric<dim, spacedim>::InternalDataBase> (
196 std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>>
205 std::move(
q1_mapping->get_face_data(update_flags, quadrature)));
210 std::move(do_get_data.return_value()));
216template <
int dim,
int spacedim>
217std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>
222 std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase> data_ptr =
223 std::make_unique<InternalData>();
224 auto &data =
dynamic_cast<InternalData &
>(*data_ptr);
228 std::unique_ptr<typename Mapping<dim, spacedim>::InternalDataBase>>
238 std::move(
q1_mapping->get_subface_data(update_flags, quadrature)));
243 std::move(do_get_data.return_value()));
250template <
int dim,
int spacedim>
262 Assert(
dynamic_cast<const InternalData *
>(&internal_data) !=
nullptr,
264 const InternalData &data =
static_cast<const InternalData &
>(internal_data);
268 data.use_mapping_q1_on_current_cell =
279 ((data.use_mapping_q1_on_current_cell ==
false) &&
286 if (data.use_mapping_q1_on_current_cell)
288 updated_cell_similarity,
290 *data.mapping_q1_data,
294 updated_cell_similarity,
296 *data.mapping_qp_data,
299 return updated_cell_similarity;
304template <
int dim,
int spacedim>
308 const unsigned int face_no,
316 Assert(
dynamic_cast<const InternalData *
>(&internal_data) !=
nullptr,
318 const InternalData &data =
static_cast<const InternalData &
>(internal_data);
326 data.use_mapping_q1_on_current_cell =
331 if (data.use_mapping_q1_on_current_cell)
333 cell, face_no, quadrature, *data.mapping_q1_data, output_data);
336 cell, face_no, quadrature, *data.mapping_qp_data, output_data);
340template <
int dim,
int spacedim>
344 const unsigned int face_no,
345 const unsigned int subface_no,
353 Assert(
dynamic_cast<const InternalData *
>(&internal_data) !=
nullptr,
355 const InternalData &data =
static_cast<const InternalData &
>(internal_data);
363 data.use_mapping_q1_on_current_cell =
368 if (data.use_mapping_q1_on_current_cell)
373 *data.mapping_q1_data,
380 *data.mapping_qp_data,
386template <
int dim,
int spacedim>
396 const InternalData *data =
dynamic_cast<const InternalData *
>(&mapping_data);
400 if (data->use_mapping_q1_on_current_cell)
401 q1_mapping->transform(input, mapping_kind, *data->mapping_q1_data, output);
404 qp_mapping->transform(input, mapping_kind, *data->mapping_qp_data, output);
409template <
int dim,
int spacedim>
419 &mapping_data) !=
nullptr),
421 const InternalData *data =
dynamic_cast<const InternalData *
>(&mapping_data);
424 if (data->use_mapping_q1_on_current_cell)
425 q1_mapping->transform(input, mapping_kind, *data->mapping_q1_data, output);
428 qp_mapping->transform(input, mapping_kind, *data->mapping_qp_data, output);
433template <
int dim,
int spacedim>
443 &mapping_data) !=
nullptr),
445 const InternalData *data =
dynamic_cast<const InternalData *
>(&mapping_data);
448 if (data->use_mapping_q1_on_current_cell)
449 q1_mapping->transform(input, mapping_kind, *data->mapping_q1_data, output);
452 qp_mapping->transform(input, mapping_kind, *data->mapping_qp_data, output);
457template <
int dim,
int spacedim>
467 &mapping_data) !=
nullptr),
469 const InternalData *data =
dynamic_cast<const InternalData *
>(&mapping_data);
472 if (data->use_mapping_q1_on_current_cell)
473 q1_mapping->transform(input, mapping_kind, *data->mapping_q1_data, output);
476 qp_mapping->transform(input, mapping_kind, *data->mapping_qp_data, output);
481template <
int dim,
int spacedim>
491 &mapping_data) !=
nullptr),
493 const InternalData *data =
dynamic_cast<const InternalData *
>(&mapping_data);
496 if (data->use_mapping_q1_on_current_cell)
497 q1_mapping->transform(input, mapping_kind, *data->mapping_q1_data, output);
500 qp_mapping->transform(input, mapping_kind, *data->mapping_qp_data, output);
505template <
int dim,
int spacedim>
515 return qp_mapping->transform_unit_to_real_cell(cell, p);
517 return q1_mapping->transform_unit_to_real_cell(cell, p);
522template <
int dim,
int spacedim>
530 return qp_mapping->transform_real_to_unit_cell(cell, p);
532 return q1_mapping->transform_real_to_unit_cell(cell, p);
537template <
int dim,
int spacedim>
538std::unique_ptr<Mapping<dim, spacedim>>
541 return std::make_unique<MappingQ<dim, spacedim>>(
547template <
int dim,
int spacedim>
555 qp_mapping->compute_mapping_support_points(cell));
562template <
int dim,
int spacedim>
568 ExcMessage(
"The dimension of your mapping (" +
570 ") and the reference cell cell_type (" +
572 " ) do not agree."));
580#include "mapping_q.inst"
virtual std::size_t memory_consumption() const override
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
virtual std::unique_ptr< typename Mapping< dim, spacedim >::InternalDataBase > get_data(const UpdateFlags, const Quadrature< dim > &quadrature) const override
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
std::shared_ptr< const MappingQGeneric< dim, spacedim > > q1_mapping
MappingQ(const unsigned int polynomial_degree, const bool use_mapping_q_on_all_cells=false)
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
virtual BoundingBox< spacedim > get_bounding_box(const typename Triangulation< dim, spacedim >::cell_iterator &cell) const override
virtual bool is_compatible_with(const ReferenceCell &reference_cell) const override
virtual Point< dim > transform_real_to_unit_cell(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const Point< spacedim > &p) const override
const unsigned int polynomial_degree
virtual bool preserves_vertex_locations() const override
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
const bool use_mapping_q_on_all_cells
virtual Point< spacedim > transform_unit_to_real_cell(const typename Triangulation< dim, spacedim >::cell_iterator &cell, const Point< dim > &p) 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
std::shared_ptr< const MappingQGeneric< dim, spacedim > > qp_mapping
virtual std::unique_ptr< typename Mapping< dim, spacedim >::InternalDataBase > get_face_data(const UpdateFlags flags, const hp::QCollection< dim - 1 > &quadrature) const override
virtual std::unique_ptr< Mapping< dim, spacedim > > clone() const override
unsigned int get_degree() const
Abstract base class for mapping classes.
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
#define Assert(cond, exc)
#define AssertDimension(dim1, dim2)
static ::ExceptionBase & ExcInternalError()
static ::ExceptionBase & ExcMessage(std::string arg1)
Task< RT > new_task(const std::function< RT()> &function)
void reference_cell(Triangulation< dim, spacedim > &tria, const ReferenceCell &reference_cell)
std::enable_if< std::is_fundamental< T >::value, std::size_t >::type memory_consumption(const T &t)
std::unique_ptr< To > dynamic_unique_cast(std::unique_ptr< From > &&p)
std::string to_string(const number value, const unsigned int digits=numbers::invalid_unsigned_int)