mapping_q_cache.h

Go to the documentation of this file.

// ---------------------------------------------------------------------
//
// Copyright (C) 2019 - 2023 by the deal.II authors
//
// This file is part of the deal.II library.
//
// The deal.II library is free software; you can use it, redistribute
// it, and/or modify it under the terms of the GNU Lesser General
// Public License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.
// The full text of the license can be found in the file LICENSE.md at
// the top level directory of deal.II.
//
// ---------------------------------------------------------------------
 
#ifndef dealii_mapping_q_cache_h
#define dealii_mapping_q_cache_h
 
 
#include <deal.II/base/config.h>
 
#include <deal.II/base/function.h>
#include <deal.II/base/mg_level_object.h>
 
#include <deal.II/fe/mapping_q.h>
 
#include <deal.II/grid/tria.h>
 
 
DEAL_II_NAMESPACE_OPEN
 
// Forward declarations
#ifndef DOXYGEN
template <int dim, int spacedim>
DEAL_II_CXX20_REQUIRES((concepts::is_valid_dim_spacedim<dim, spacedim>))
class DoFHandler;
#endif
 
 
template <int dim, int spacedim = dim>
class MappingQCache : public MappingQ<dim, spacedim>
{
public:
  explicit MappingQCache(const unsigned int polynomial_degree);
 
  explicit MappingQCache(const MappingQCache<dim, spacedim> &mapping);
 
  ~MappingQCache();
 
  virtual std::unique_ptr<Mapping<dim, spacedim>>
  clone() const override;
 
  virtual bool
  preserves_vertex_locations() const override;
 
  void
  initialize(const Mapping<dim, spacedim> &      mapping,
             const Triangulation<dim, spacedim> &triangulation);
 
  DEAL_II_DEPRECATED void
  initialize(const Triangulation<dim, spacedim> &triangulation,
             const MappingQ<dim, spacedim> &     mapping);
 
  void
  initialize(const Triangulation<dim, spacedim> &triangulation,
             const std::function<std::vector<Point<spacedim>>(
               const typename Triangulation<dim, spacedim>::cell_iterator &)>
               &compute_points_on_cell);
 
  void
  initialize(const Mapping<dim, spacedim> &      mapping,
             const Triangulation<dim, spacedim> &tria,
             const std::function<Point<spacedim>(
               const typename Triangulation<dim, spacedim>::cell_iterator &,
               const Point<spacedim> &)> &       transformation_function,
             const bool function_describes_relative_displacement);
 
  void
  initialize(const Mapping<dim, spacedim> &      mapping,
             const Triangulation<dim, spacedim> &tria,
             const Function<spacedim> &          transformation_function,
             const bool function_describes_relative_displacement);
 
  template <typename VectorType>
  void
  initialize(const Mapping<dim, spacedim> &   mapping,
             const DoFHandler<dim, spacedim> &dof_handler,
             const VectorType &               vector,
             const bool vector_describes_relative_displacement);
 
  template <typename VectorType>
  void
  initialize(const Mapping<dim, spacedim> &   mapping,
             const DoFHandler<dim, spacedim> &dof_handler,
             const MGLevelObject<VectorType> &vectors,
             const bool vector_describes_relative_displacement);
 
  virtual boost::container::small_vector<Point<spacedim>,
                                         GeometryInfo<dim>::vertices_per_cell>
  get_vertices(const typename Triangulation<dim, spacedim>::cell_iterator &cell)
    const override;
 
  std::size_t
  memory_consumption() const;
 
protected:
  virtual std::vector<Point<spacedim>>
  compute_mapping_support_points(
    const typename Triangulation<dim, spacedim>::cell_iterator &cell)
    const override;
 
private:
  std::shared_ptr<std::vector<std::vector<std::vector<Point<spacedim>>>>>
    support_point_cache;
 
  boost::signals2::connection clear_signal;
 
  bool uses_level_info;
};
 
DEAL_II_NAMESPACE_CLOSE
 
#endif