mapping_q_eulerian.h

Go to the documentation of this file.

// ------------------------------------------------------------------------
//
// SPDX-License-Identifier: LGPL-2.1-or-later
// Copyright (C) 2008 - 2025 by the deal.II authors
//
// This file is part of the deal.II library.
//
// Part of the source code is dual licensed under Apache-2.0 WITH
// LLVM-exception OR LGPL-2.1-or-later. Detailed license information
// governing the source code and code contributions can be found in
// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II.
//
// ------------------------------------------------------------------------
 
 
#ifndef dealii_mapping_q_eulerian_h
#define dealii_mapping_q_eulerian_h
 
#include <deal.II/base/config.h>
 
#include <deal.II/base/mutex.h>
#include <deal.II/base/observer_pointer.h>
 
#include <deal.II/dofs/dof_accessor.h>
#include <deal.II/dofs/dof_handler.h>
 
#include <deal.II/fe/fe_values.h>
#include <deal.II/fe/mapping_q.h>
 
#include <deal.II/grid/tria_iterator.h>
 
#include <boost/container/small_vector.hpp>
 
 
DEAL_II_NAMESPACE_OPEN
 
template <typename>
class Vector;
 


template <int dim, typename VectorType = Vector<double>, int spacedim = dim>
class MappingQEulerian : public MappingQ<dim, spacedim>
{
public:
  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 boost::container::small_vector<Point<spacedim>,
#ifndef _MSC_VER
                                         ReferenceCells::max_n_vertices<dim>()
#else
                                         GeometryInfo<dim>::vertices_per_cell
#endif
                                         >
  get_vertices(const typename Triangulation<dim, spacedim>::cell_iterator &cell)
    const override;

  virtual std::unique_ptr<Mapping<dim, spacedim>>
  clone() const override;

  virtual bool
  preserves_vertex_locations() const override;
 
  // for documentation, see the Mapping base class
  virtual std::vector<Point<spacedim>>
  compute_mapping_support_points(
    const typename Triangulation<dim, spacedim>::cell_iterator &cell)
    const override;

  DeclException0(ExcInactiveCell);
 
protected:
  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;

  ObserverPointer<const VectorType, MappingQEulerian<dim, VectorType, spacedim>>
    euler_vector;

  ObserverPointer<const DoFHandler<dim, spacedim>,
                  MappingQEulerian<dim, VectorType, spacedim>>
    euler_dof_handler;
 
private:
  const unsigned int level;

  class SupportQuadrature : public Quadrature<dim>
  {
  public:
    SupportQuadrature(const unsigned int map_degree);
  };

  const SupportQuadrature support_quadrature;

  const MappingQ<dim, spacedim> mapping_q;

  mutable FEValues<dim, spacedim> fe_values;

  mutable Threads::Mutex fe_values_mutex;
};

 
 
/*----------------------------------------------------------------------*/
 
#ifndef DOXYGEN
 
template <int dim, typename VectorType, int spacedim>
inline bool
MappingQEulerian<dim, VectorType, spacedim>::preserves_vertex_locations() const
{
  return false;
}
 
#endif // DOXYGEN
 
 
DEAL_II_NAMESPACE_CLOSE
 
 
#endif // dealii_mapping_q_eulerian_h