33#include <boost/container/small_vector.hpp>
42template <
int dim,
typename VectorType,
int spacedim>
45 const VectorType &euler_transform_vectors)
47 , euler_transform_vectors(&euler_transform_vectors)
48 , shiftmap_dof_handler(&shiftmap_dof_handler)
53template <
int dim,
typename VectorType,
int spacedim>
54boost::container::small_vector<Point<spacedim>,
56 ReferenceCells::max_n_vertices<dim>()
64 boost::container::small_vector<Point<spacedim>,
66 ReferenceCells::max_n_vertices<dim>()
71 vertices(cell->n_vertices());
77 AssertDimension(spacedim, shiftmap_dof_handler->get_fe().n_dofs_per_vertex());
78 AssertDimension(shiftmap_dof_handler->get_fe(0).n_components(), spacedim);
81 euler_transform_vectors->size());
87 *cell, shiftmap_dof_handler);
91 Assert(dof_cell->is_active() ==
true, ExcInactiveCell());
95 shiftmap_dof_handler->get_fe().n_dofs_per_cell());
96 dof_cell->get_dof_values(*euler_transform_vectors, mapping_values);
105 for (
unsigned int j = 0; j < spacedim; ++j)
106 shift_vector[j] = mapping_values(i * spacedim + j);
110 vertices[i] = cell->vertex(i) + shift_vector;
117template <
int dim,
typename VectorType,
int spacedim>
118std::vector<Point<spacedim>>
122 const auto vertices = this->get_vertices(cell);
133template <
int dim,
typename VectorType,
int spacedim>
134std::unique_ptr<Mapping<dim, spacedim>>
137 return std::make_unique<MappingQ1Eulerian<dim, VectorType, spacedim>>(*this);
142template <
int dim,
typename VectorType,
int spacedim>
169#include "fe/mapping_q1_eulerian.inst"
virtual std::vector< Point< spacedim > > compute_mapping_support_points(const typename Triangulation< dim, spacedim >::cell_iterator &cell) const override
virtual std::unique_ptr< Mapping< dim, spacedim > > clone() 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
MappingQ1Eulerian(const DoFHandler< dim, spacedim > &euler_dof_handler, const VectorType &euler_vector)
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)
typename ActiveSelector::cell_iterator cell_iterator
static std_cxx20::ranges::iota_view< unsigned int, unsigned int > vertex_indices()