Reference documentation for deal.II version Git 1206453fa0 2020-07-09 18:20:27 -0600
\(\newcommand{\dealvcentcolon}{\mathrel{\mathop{:}}}\) \(\newcommand{\dealcoloneq}{\dealvcentcolon\mathrel{\mkern-1.2mu}=}\) \(\newcommand{\jump}[1]{\left[\!\left[ #1 \right]\!\right]}\) \(\newcommand{\average}[1]{\left\{\!\left\{ #1 \right\}\!\right\}}\)
mapping_q_eulerian.h
Go to the documentation of this file.
1 // ---------------------------------------------------------------------
2 //
3 // Copyright (C) 2001 - 2018 by the deal.II authors
4 //
5 // This file is part of the deal.II library.
6 //
7 // The deal.II library is free software; you can use it, redistribute
8 // it, and/or modify it under the terms of the GNU Lesser General
9 // Public License as published by the Free Software Foundation; either
10 // version 2.1 of the License, or (at your option) any later version.
11 // The full text of the license can be found in the file LICENSE.md at
12 // the top level directory of deal.II.
13 //
14 // ---------------------------------------------------------------------
15 
16 
17 #ifndef dealii_mapping_q_eulerian_h
18 #define dealii_mapping_q_eulerian_h
19 
20 #include <deal.II/base/config.h>
21 
24 
27 
28 #include <deal.II/fe/fe_values.h>
29 #include <deal.II/fe/mapping_q.h>
30 
32 
33 
35 
36 template <typename>
37 class Vector;
38 
39 
42 
92 template <int dim, typename VectorType = Vector<double>, int spacedim = dim>
93 class MappingQEulerian : public MappingQ<dim, spacedim>
94 {
95 public:
112  MappingQEulerian(const unsigned int degree,
114  const VectorType & euler_vector,
115  const unsigned int level = numbers::invalid_unsigned_int);
116 
123  virtual std::array<Point<spacedim>, GeometryInfo<dim>::vertices_per_cell>
125  const override;
126 
131  virtual std::unique_ptr<Mapping<dim, spacedim>>
132  clone() const override;
133 
139  virtual bool
140  preserves_vertex_locations() const override;
141 
147 
148 protected:
159  const typename Triangulation<dim, spacedim>::cell_iterator &cell,
160  const CellSimilarity::Similarity cell_similarity,
161  const Quadrature<dim> & quadrature,
162  const typename Mapping<dim, spacedim>::InternalDataBase & internal_data,
164  &output_data) const override;
165 
171 
178 
179 
180 private:
184  const unsigned int level;
185 
191  class MappingQEulerianGeneric : public MappingQGeneric<dim, spacedim>
192  {
193  public:
198  const unsigned int degree,
200 
207  virtual std::array<Point<spacedim>, GeometryInfo<dim>::vertices_per_cell>
209  &cell) const override;
210 
216  virtual std::vector<Point<spacedim>>
218  const typename Triangulation<dim, spacedim>::cell_iterator &cell)
219  const override;
220 
226  virtual bool
227  preserves_vertex_locations() const override;
228 
229  private:
234 
235 
240  class SupportQuadrature : public Quadrature<dim>
241  {
242  public:
246  SupportQuadrature(const unsigned int map_degree);
247  };
248 
253 
263 
268  };
269 };
270 
274 /*----------------------------------------------------------------------*/
275 
276 #ifndef DOXYGEN
277 
278 template <int dim, typename VectorType, int spacedim>
279 inline bool
281 {
282  return false;
283 }
284 
285 #endif // DOXYGEN
286 
287 
289 
290 
291 #endif // dealii_mapping_q_eulerian_h
static const unsigned int invalid_unsigned_int
Definition: types.h:191
SmartPointer< const DoFHandler< dim, spacedim >, MappingQEulerian< dim, VectorType, spacedim > > euler_dof_handler
const unsigned int level
static ::ExceptionBase & ExcInactiveCell()
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
MappingQEulerianGeneric(const unsigned int degree, const MappingQEulerian< dim, VectorType, spacedim > &mapping_q_eulerian)
virtual std::array< Point< spacedim >, GeometryInfo< dim >::vertices_per_cell > get_vertices(const typename Triangulation< dim, spacedim >::cell_iterator &cell) const override
virtual std::vector< Point< spacedim > > compute_mapping_support_points(const typename Triangulation< dim, spacedim >::cell_iterator &cell) const override
#define DeclException0(Exception0)
Definition: exceptions.h:470
#define DEAL_II_NAMESPACE_CLOSE
Definition: config.h:363
MappingQEulerian(const unsigned int degree, const DoFHandler< dim, spacedim > &euler_dof_handler, const VectorType &euler_vector, const unsigned int level=numbers::invalid_unsigned_int)
virtual std::array< Point< spacedim >, GeometryInfo< dim >::vertices_per_cell > get_vertices(const typename Triangulation< dim, spacedim >::cell_iterator &cell) const override
virtual bool preserves_vertex_locations() const override
virtual std::unique_ptr< Mapping< dim, spacedim > > clone() const override
#define DEAL_II_NAMESPACE_OPEN
Definition: config.h:362
const MappingQEulerian< dim, VectorType, spacedim > & mapping_q_eulerian
virtual bool preserves_vertex_locations() const override
SmartPointer< const VectorType, MappingQEulerian< dim, VectorType, spacedim > > euler_vector