intergrid_map.h

// ---------------------------------------------------------------------
//
// Copyright (C) 1999 - 2017 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 at
// the top level of the deal.II distribution.
//
// ---------------------------------------------------------------------

#ifndef dealii_intergrid_map_h
#define dealii_intergrid_map_h

#include <deal.II/base/config.h>
#include <deal.II/base/smartpointer.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 <class MeshType>
class InterGridMap : public Subscriptor
{
public:

  typedef typename MeshType::cell_iterator 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;

  SmartPointer<const MeshType,InterGridMap<MeshType> > source_grid;

  SmartPointer<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