36#include <boost/container/small_vector.hpp>
45template <
int dim,
typename VectorType,
int spacedim>
47 const unsigned int degree,
49 const VectorType &euler_vector,
50 const unsigned int level)
52 , euler_vector(&euler_vector)
53 , euler_dof_handler(&euler_dof_handler)
55 , support_quadrature(degree)
57 , fe_values(mapping_q,
58 euler_dof_handler.get_fe(),
65template <
int dim,
typename VectorType,
int spacedim>
66std::unique_ptr<Mapping<dim, spacedim>>
69 return std::make_unique<MappingQEulerian<dim, VectorType, spacedim>>(
70 this->get_degree(), *euler_dof_handler, *euler_vector, this->
level);
75template <
int dim,
typename VectorType,
int spacedim>
84 const unsigned int n_q_points = q_iterated.size();
88 const std::vector<unsigned int> renumber =
89 FETools::lexicographic_to_hierarchic_numbering<dim>(map_degree);
92 for (
unsigned int q = 0; q < n_q_points; ++q)
93 this->quadrature_points[renumber[q]] = q_iterated.point(q);
98template <
int dim,
typename VectorType,
int spacedim>
99boost::container::small_vector<Point<spacedim>,
101 ReferenceCells::max_n_vertices<dim>()
112 boost::container::small_vector<Point<spacedim>,
114 ReferenceCells::max_n_vertices<dim>()
119 vertex_locations(a.begin(), a.begin() + cell->n_vertices());
121 return vertex_locations;
126template <
int dim,
typename VectorType,
int spacedim>
127std::vector<Point<spacedim>>
161 Assert(n_components >= spacedim,
164 std::vector<Vector<typename VectorType::value_type>> shift_vector(
167 std::vector<types::global_dof_index> dof_indices(
176 dof_cell->get_mg_dof_indices(dof_indices);
184 std::vector<Point<spacedim>> a(n_support_pts);
185 for (
unsigned int q = 0; q < n_support_pts; ++q)
188 for (
unsigned int d = 0; d < spacedim; ++d)
189 a[q][d] += shift_vector[q][d];
197template <
int dim,
typename VectorType,
int spacedim>
223#include "fe/mapping_q_eulerian.inst"
SupportQuadrature(const unsigned int map_degree)
FEValues< dim, spacedim > fe_values
virtual std::vector< Point< spacedim > > compute_mapping_support_points(const typename Triangulation< dim, spacedim >::cell_iterator &cell) 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
virtual std::unique_ptr< Mapping< dim, spacedim > > clone() const override
MappingQEulerian(const unsigned int degree, const DoFHandler< dim, spacedim > &euler_dof_handler, const VectorType &euler_vector, const unsigned int level=numbers::invalid_unsigned_int)
ObserverPointer< const VectorType, MappingQEulerian< dim, VectorType, spacedim > > euler_vector
ObserverPointer< const DoFHandler< dim, spacedim >, MappingQEulerian< dim, VectorType, spacedim > > euler_dof_handler
Threads::Mutex fe_values_mutex
const SupportQuadrature support_quadrature
virtual boost::container::small_vector< Point< spacedim >, ReferenceCells::max_n_vertices< dim >() > get_vertices(const typename Triangulation< dim, spacedim >::cell_iterator &cell) 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
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
#define Assert(cond, exc)
#define AssertDimension(dim1, dim2)
static ::ExceptionBase & ExcInactiveCell()
static ::ExceptionBase & ExcDimensionMismatch(std::size_t arg1, std::size_t arg2)
typename ActiveSelector::cell_iterator cell_iterator
@ update_values
Shape function values.
@ update_quadrature_points
Transformed quadrature points.
constexpr unsigned int invalid_unsigned_int