intergrid_map.h

Go to the documentation of this file.

// ------------------------------------------------------------------------
//
// SPDX-License-Identifier: LGPL-2.1-or-later
// Copyright (C) 1999 - 2023 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_intergrid_map_h
#define dealii_intergrid_map_h
 
#include <deal.II/base/config.h>
 
#include <deal.II/base/observer_pointer.h>
 
#include <deal.II/dofs/dof_accessor.h>
 
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
 
DEAL_II_NAMESPACE_OPEN
 
 
template <typename MeshType>
DEAL_II_CXX20_REQUIRES(concepts::is_triangulation_or_dof_handler<MeshType>)
class InterGridMap : public EnableObserverPointer
{
public:
  using cell_iterator = typename MeshType::cell_iterator;
 
  InterGridMap();
 
  void
  make_mapping(const MeshType &source_grid, const MeshType &destination_grid);
 
  cell_iterator
  operator[](const cell_iterator &source_cell) const;
 
  void
  clear();
 
  const MeshType &
  get_source_grid() const;
 
  const MeshType &
  get_destination_grid() const;
 
  std::size_t
  memory_consumption() const;
 
  DeclException1(ExcInvalidKey,
                 cell_iterator,
                 << "The iterator " << arg1 << " is not valid as key for "
                 << "this map.");
  DeclException0(ExcIncompatibleGrids);
 
private:
  std::vector<std::vector<cell_iterator>> mapping;
 
  ObserverPointer<const MeshType, InterGridMap<MeshType>> source_grid;
 
  ObserverPointer<const MeshType, InterGridMap<MeshType>> destination_grid;
 
  void
  set_mapping(const cell_iterator &src_cell, const cell_iterator &dst_cell);
 
  void
  set_entries_to_cell(const cell_iterator &src_cell,
                      const cell_iterator &dst_cell);
};
 
 
DEAL_II_NAMESPACE_CLOSE
 
#endif