tria.h

// ---------------------------------------------------------------------
//
// Copyright (C) 1998 - 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__tria_h
#define dealii__tria_h


#include <deal.II/base/config.h>
#include <deal.II/base/point.h>
#include <deal.II/base/subscriptor.h>
#include <deal.II/base/smartpointer.h>
#include <deal.II/base/geometry_info.h>
#include <deal.II/base/iterator_range.h>
#include <deal.II/base/std_cxx11/function.h>
#include <deal.II/base/std_cxx11/unique_ptr.h>
#include <deal.II/grid/tria_iterator_selector.h>

// Ignore deprecation warnings for auto_ptr.
DEAL_II_DISABLE_EXTRA_DIAGNOSTICS
#include <boost/signals2.hpp>
#include <boost/serialization/vector.hpp>
#include <boost/serialization/map.hpp>
#ifdef DEAL_II_WITH_CXX11
#  include <boost/serialization/unique_ptr.hpp>
#else
#  include <boost/serialization/scoped_ptr.hpp>
#endif
#include <boost/serialization/split_member.hpp>
DEAL_II_ENABLE_EXTRA_DIAGNOSTICS

#include <vector>
#include <list>
#include <map>
#include <numeric>
#include <bitset>


DEAL_II_NAMESPACE_OPEN

template <int dim, int spacedim> class Boundary;
template <int dim, int spacedim> class StraightBoundary;
template <int dim, int spacedim> class Manifold;

namespace GridTools
{
  template<typename CellIterator>  struct PeriodicFacePair;
}

template <int, int, int> class TriaAccessor;
template <int spacedim> class TriaAccessor<0,1,spacedim>;
template <int, int, int> class TriaAccessorBase;

namespace internal
{
  namespace Triangulation
  {
    template <int dim> class TriaLevel;
    template <int dim> class TriaFaces;

    template <typename> class TriaObjects;

    struct Implementation;
  }

  namespace TriaAccessor
  {
    struct Implementation;
  }
}

namespace hp
{
  template <int dim, int spacedim> class DoFHandler;
}


/*------------------------------------------------------------------------*/

template <int structdim>
struct CellData
{
  unsigned int vertices[GeometryInfo<structdim>::vertices_per_cell];

  union
  {
    types::material_id material_id;

    types::boundary_id boundary_id;
  };

  types::manifold_id manifold_id;

  CellData ();
};



struct SubCellData
{
  std::vector<CellData<1> > boundary_lines;

  std::vector<CellData<2> > boundary_quads;

  bool check_consistency (const unsigned int dim) const;
};


/*------------------------------------------------------------------------*/


namespace internal
{
  namespace Triangulation
  {

    template <int dim>
    struct NumberCache
    {
    };

    template <>
    struct NumberCache<1>
    {
      unsigned int n_levels;

      unsigned int n_lines;

      std::vector<unsigned int> n_lines_level;

      unsigned int n_active_lines;

      std::vector<unsigned int> n_active_lines_level;

      NumberCache ();

      std::size_t memory_consumption () const;

      template <class Archive>
      void serialize (Archive &ar,
                      const unsigned int version);
    };


    template <>
    struct NumberCache<2> : public NumberCache<1>
    {
      unsigned int n_quads;

      std::vector<unsigned int> n_quads_level;

      unsigned int n_active_quads;

      std::vector<unsigned int> n_active_quads_level;

      NumberCache ();

      std::size_t memory_consumption () const;

      template <class Archive>
      void serialize (Archive &ar,
                      const unsigned int version);
    };


    template <>
    struct NumberCache<3> : public NumberCache<2>
    {
      unsigned int n_hexes;

      std::vector<unsigned int> n_hexes_level;

      unsigned int n_active_hexes;

      std::vector<unsigned int> n_active_hexes_level;

      NumberCache ();

      std::size_t memory_consumption () const;

      template <class Archive>
      void serialize (Archive &ar,
                      const unsigned int version);
    };
  }
}


/*------------------------------------------------------------------------*/


template <int dim, int spacedim=dim>
class Triangulation : public Subscriptor
{
private:

  typedef ::internal::Triangulation::Iterators<dim, spacedim> IteratorSelector;

public:
  static const StraightBoundary<dim,spacedim> straight_boundary;

  enum MeshSmoothing
  {
    none                               = 0x0,
    limit_level_difference_at_vertices = 0x1,
    eliminate_unrefined_islands        = 0x2,
    patch_level_1                      = 0x4,
    coarsest_level_1                   = 0x8,
    allow_anisotropic_smoothing        = 0x10,
    eliminate_refined_inner_islands    = 0x100,
    eliminate_refined_boundary_islands = 0x200,
    do_not_produce_unrefined_islands   = 0x400,

    smoothing_on_refinement            = (limit_level_difference_at_vertices |
                                          eliminate_unrefined_islands),
    smoothing_on_coarsening            = (eliminate_refined_inner_islands |
                                          eliminate_refined_boundary_islands |
                                          do_not_produce_unrefined_islands),

    maximum_smoothing                  = 0xffff ^ allow_anisotropic_smoothing
  };

  typedef TriaIterator      <CellAccessor<dim,spacedim>         > cell_iterator;

  typedef TriaActiveIterator<CellAccessor<dim,spacedim>         > active_cell_iterator;

  typedef TriaIterator      <TriaAccessor<dim-1, dim, spacedim> > face_iterator;
  typedef TriaActiveIterator<TriaAccessor<dim-1, dim, spacedim> > active_face_iterator;

  typedef typename IteratorSelector::vertex_iterator        vertex_iterator;
  typedef typename IteratorSelector::active_vertex_iterator active_vertex_iterator;

  typedef typename IteratorSelector::line_iterator        line_iterator;
  typedef typename IteratorSelector::active_line_iterator active_line_iterator;

  typedef typename IteratorSelector::quad_iterator        quad_iterator;
  typedef typename IteratorSelector::active_quad_iterator active_quad_iterator;

  typedef typename IteratorSelector::hex_iterator         hex_iterator;
  typedef typename IteratorSelector::active_hex_iterator  active_hex_iterator;

  struct DistortedCellList : public ::ExceptionBase
  {
    virtual ~DistortedCellList () DEAL_II_NOEXCEPT;

    std::list<typename Triangulation<dim,spacedim>::cell_iterator>
    distorted_cells;
  };

  static const unsigned int dimension = dim;

  static const unsigned int space_dimension = spacedim;

  Triangulation (const MeshSmoothing smooth_grid = none,
                 const bool check_for_distorted_cells = false);

  Triangulation (const Triangulation<dim, spacedim> &t);

#ifdef DEAL_II_WITH_CXX11

  Triangulation (Triangulation<dim, spacedim> &&tria);

  Triangulation &operator = (Triangulation<dim, spacedim> &&tria);
#endif

  virtual ~Triangulation ();

  virtual void clear ();

  virtual void set_mesh_smoothing (const MeshSmoothing mesh_smoothing);

  virtual const MeshSmoothing &get_mesh_smoothing() const;

  void set_boundary (const types::manifold_id   number,
                     const Boundary<dim,spacedim> &boundary_object);


  void set_boundary (const types::manifold_id number);

  void set_manifold (const types::manifold_id   number,
                     const Manifold<dim,spacedim> &manifold_object);


  void set_manifold (const types::manifold_id number);

  void set_all_manifold_ids (const types::manifold_id number);

  void set_all_manifold_ids_on_boundary (const types::manifold_id number);

  void set_all_manifold_ids_on_boundary (const types::boundary_id b_id,
                                         const types::manifold_id number);


  const Boundary<dim,spacedim> &get_boundary (const types::manifold_id number) const;

  const Manifold<dim,spacedim> &get_manifold (const types::manifold_id number) const;

  std::vector<types::boundary_id> get_boundary_ids() const;

  std::vector<types::manifold_id> get_manifold_ids() const;

  virtual void copy_triangulation (const Triangulation<dim, spacedim> &other_tria);

  virtual void create_triangulation (const std::vector<Point<spacedim> >    &vertices,
                                     const std::vector<CellData<dim> > &cells,
                                     const SubCellData                 &subcelldata);

  virtual void create_triangulation_compatibility (
    const std::vector<Point<spacedim> >    &vertices,
    const std::vector<CellData<dim> > &cells,
    const SubCellData                 &subcelldata);

  void flip_all_direction_flags();

  void set_all_refine_flags ();

  void refine_global (const unsigned int times = 1);

  virtual void execute_coarsening_and_refinement ();

  virtual bool prepare_coarsening_and_refinement ();

  /*
   * @}
   */

  enum CellStatus
  {
    CELL_PERSIST,
    CELL_REFINE,
    CELL_COARSEN,
    CELL_INVALID
  };

  template<typename T>
  struct CellWeightSum
  {
    typedef T result_type;

    template<typename InputIterator>
    T operator()(InputIterator first, InputIterator last) const
    {
      return std::accumulate (first, last, T());
    }
  };

  struct Signals
  {
    boost::signals2::signal<void ()> create;

    boost::signals2::signal<void ()> pre_refinement;

    boost::signals2::signal<void ()> post_refinement;

    boost::signals2::signal<void ()> mesh_movement;

    boost::signals2::signal<void (const typename Triangulation<dim, spacedim>::cell_iterator &cell)> pre_coarsening_on_cell;

    boost::signals2::signal<void (const typename Triangulation<dim, spacedim>::cell_iterator &cell)> post_refinement_on_cell;

    boost::signals2::signal<void (const Triangulation<dim, spacedim> &destination_tria)> copy;

    boost::signals2::signal<void ()> clear;

    boost::signals2::signal<void ()> any_change;

    boost::signals2::signal<unsigned int (const cell_iterator &,
                                          const CellStatus),
                                                CellWeightSum<unsigned int> > cell_weight;
  };

  mutable Signals signals;

  /*
   * @}
   */

  void save_refine_flags (std::ostream &out) const;

  void save_refine_flags (std::vector<bool> &v) const;

  void load_refine_flags (std::istream &in);

  void load_refine_flags (const std::vector<bool> &v);

  void save_coarsen_flags (std::ostream &out) const;

  void save_coarsen_flags (std::vector<bool> &v) const;

  void load_coarsen_flags (std::istream &out);

  void load_coarsen_flags (const std::vector<bool> &v);

  bool get_anisotropic_refinement_flag() const;

  /*
   * @}
   */

  void clear_user_flags ();

  void save_user_flags (std::ostream &out) const;

  void save_user_flags (std::vector<bool> &v) const;

  void load_user_flags (std::istream &in);

  void load_user_flags (const std::vector<bool> &v);

  void clear_user_flags_line ();

  void save_user_flags_line (std::ostream &out) const;

  void save_user_flags_line (std::vector<bool> &v) const;

  void load_user_flags_line (std::istream &in);

  void load_user_flags_line (const std::vector<bool> &v);

  void clear_user_flags_quad ();

  void save_user_flags_quad (std::ostream &out) const;

  void save_user_flags_quad (std::vector<bool> &v) const;

  void load_user_flags_quad (std::istream &in);

  void load_user_flags_quad (const std::vector<bool> &v);


  void clear_user_flags_hex ();

  void save_user_flags_hex (std::ostream &out) const;

  void save_user_flags_hex (std::vector<bool> &v) const;

  void load_user_flags_hex (std::istream &in);

  void load_user_flags_hex (const std::vector<bool> &v);

  void clear_user_data ();

  void save_user_indices (std::vector<unsigned int> &v) const;

  void load_user_indices (const std::vector<unsigned int> &v);

  void save_user_pointers (std::vector<void *> &v) const;

  void load_user_pointers (const std::vector<void *> &v);

  void save_user_indices_line (std::vector<unsigned int> &v) const;

  void load_user_indices_line (const std::vector<unsigned int> &v);

  void save_user_indices_quad (std::vector<unsigned int> &v) const;

  void load_user_indices_quad (const std::vector<unsigned int> &v);

  void save_user_indices_hex (std::vector<unsigned int> &v) const;

  void load_user_indices_hex (const std::vector<unsigned int> &v);
  void save_user_pointers_line (std::vector<void *> &v) const;

  void load_user_pointers_line (const std::vector<void *> &v);

  void save_user_pointers_quad (std::vector<void *> &v) const;

  void load_user_pointers_quad (const std::vector<void *> &v);

  void save_user_pointers_hex (std::vector<void *> &v) const;

  void load_user_pointers_hex (const std::vector<void *> &v);

  /*
   * @}
   */

  cell_iterator        begin       (const unsigned int level = 0) const;

  active_cell_iterator begin_active(const unsigned int level = 0) const;

  cell_iterator        end () const;

  cell_iterator        end (const unsigned int level) const;

  active_cell_iterator end_active (const unsigned int level) const;


  cell_iterator        last () const;

  active_cell_iterator last_active () const;

  IteratorRange<cell_iterator>        cell_iterators () const;

  IteratorRange<active_cell_iterator> active_cell_iterators () const;

  IteratorRange<cell_iterator>        cell_iterators_on_level (const unsigned int level) const;

  IteratorRange<active_cell_iterator> active_cell_iterators_on_level (const unsigned int level) const;

  /*
   * @}
   */

  /*---------------------------------------*/
  /*---------------------------------------*/

  face_iterator        begin_face       () const;

  active_face_iterator begin_active_face() const;

  face_iterator        end_face () const;

  /*
   * @}
   */

  /*---------------------------------------*/
  /*---------------------------------------*/

  vertex_iterator        begin_vertex() const;

  active_vertex_iterator begin_active_vertex() const;

  vertex_iterator        end_vertex() const;

  /*
   * @}
   */

  unsigned int n_lines () const;

  unsigned int n_lines (const unsigned int level) const;

  unsigned int n_active_lines () const;

  unsigned int n_active_lines (const unsigned int level) const;

  unsigned int n_quads () const;

  unsigned int n_quads (const unsigned int level) const;

  unsigned int n_active_quads () const;

  unsigned int n_active_quads (const unsigned int level) const;

  unsigned int n_hexs() const;

  unsigned int n_hexs(const unsigned int level) const;

  unsigned int n_active_hexs() const;

  unsigned int n_active_hexs(const unsigned int level) const;

  unsigned int n_cells () const;

  unsigned int n_cells (const unsigned int level) const;

  unsigned int n_active_cells () const;

  virtual types::global_dof_index n_global_active_cells () const;


  unsigned int n_active_cells (const unsigned int level) const;

  unsigned int n_faces () const;

  unsigned int n_active_faces () const;

  unsigned int n_levels () const;

  virtual
  unsigned int n_global_levels () const;

  virtual
  bool has_hanging_nodes() const;

  unsigned int n_vertices () const;

  const std::vector<Point<spacedim> > &
  get_vertices () const;

  unsigned int n_used_vertices () const;

  bool vertex_used (const unsigned int index) const;

  const std::vector<bool> &
  get_used_vertices () const;

  unsigned int max_adjacent_cells () const;

  virtual types::subdomain_id locally_owned_subdomain () const;

  Triangulation<dim,spacedim> &
  get_triangulation ();

  const Triangulation<dim,spacedim> &
  get_triangulation () const;


  /*
   * @}
   */

  unsigned int n_raw_lines () const;

  unsigned int n_raw_lines (const unsigned int level) const;

  unsigned int n_raw_quads () const;

  unsigned int n_raw_quads (const unsigned int level) const;

  unsigned int n_raw_hexs () const;

  unsigned int n_raw_hexs (const unsigned int level) const;

  unsigned int n_raw_cells (const unsigned int level) const;

  unsigned int n_raw_faces () const;

  /*
   * @}
   */

  virtual std::size_t memory_consumption () const;

  template <class Archive>
  void save (Archive &ar,
             const unsigned int version) const;

  template <class Archive>
  void load (Archive &ar,
             const unsigned int version);


  virtual void
  add_periodicity
  (const std::vector<GridTools::PeriodicFacePair<cell_iterator> > &);

  const std::map<std::pair<cell_iterator, unsigned int>,std::pair<std::pair<cell_iterator, unsigned int>, std::bitset<3> > > &
  get_periodic_face_map() const;


  BOOST_SERIALIZATION_SPLIT_MEMBER()

  
  DeclException1 (ExcInvalidLevel,
                  int,
                  << "The given level " << arg1
                  << " is not in the valid range!");
  DeclException2 (ExcTriangulationNotEmpty,
                  int, int,
                  << "You are trying to perform an operation on a triangulation "
                  << "that is only allowed if the triangulation is currently empty. "
                  << "However, it currently stores " << arg1 << " vertices and has "
                  << "cells on " << arg2 << " levels.");
  DeclException0 (ExcGridReadError);
  DeclException0 (ExcFacesHaveNoLevel);
  DeclException1 (ExcEmptyLevel,
                  int,
                  << "You tried to do something on level " << arg1
                  << ", but this level is empty.");
  DeclException0 (ExcNonOrientableTriangulation);

  DeclException1 (ExcBoundaryIdNotFound,
                  types::boundary_id,
                  << "The given boundary_id " << arg1
                  << " is not defined in this Triangulation!");

  /*
   * @}
   */

protected:
  MeshSmoothing                    smooth_grid;

  static void write_bool_vector (const unsigned int       magic_number1,
                                 const std::vector<bool> &v,
                                 const unsigned int       magic_number2,
                                 std::ostream            &out);

  static void read_bool_vector (const unsigned int       magic_number1,
                                std::vector<bool>       &v,
                                const unsigned int       magic_number2,
                                std::istream            &in);

  void update_periodic_face_map ();


private:
  std::vector<GridTools::PeriodicFacePair<cell_iterator> > periodic_face_pairs_level_0;

  std::map<std::pair<cell_iterator, unsigned int>, std::pair<std::pair<cell_iterator, unsigned int>, std::bitset<3> > > periodic_face_map;

  typedef TriaRawIterator   <CellAccessor<dim,spacedim>         > raw_cell_iterator;
  typedef TriaRawIterator   <TriaAccessor<dim-1, dim, spacedim> > raw_face_iterator;
  typedef typename IteratorSelector::raw_vertex_iterator          raw_vertex_iterator;
  typedef typename IteratorSelector::raw_line_iterator            raw_line_iterator;
  typedef typename IteratorSelector::raw_quad_iterator            raw_quad_iterator;
  typedef typename IteratorSelector::raw_hex_iterator             raw_hex_iterator;

  raw_cell_iterator    begin_raw   (const unsigned int level = 0) const;

  raw_cell_iterator    end_raw (const unsigned int level) const;

  /*
   * @}
   */

  raw_line_iterator
  begin_raw_line   (const unsigned int level = 0) const;

  line_iterator
  begin_line       (const unsigned int level = 0) const;

  active_line_iterator
  begin_active_line(const unsigned int level = 0) const;

  line_iterator        end_line () const;

  /*
   * @}
   */

  raw_quad_iterator
  begin_raw_quad   (const unsigned int level = 0) const;

  quad_iterator
  begin_quad       (const unsigned int level = 0) const;

  active_quad_iterator
  begin_active_quad (const unsigned int level = 0) const;

  quad_iterator
  end_quad () const;

  /*
   * @}
   */

  raw_hex_iterator
  begin_raw_hex   (const unsigned int level = 0) const;

  hex_iterator
  begin_hex       (const unsigned int level = 0) const;

  active_hex_iterator
  begin_active_hex (const unsigned int level = 0) const;

  hex_iterator
  end_hex () const;

  /*
   * @}
   */


  void clear_despite_subscriptions ();

  void reset_active_cell_indices ();

  DistortedCellList execute_refinement ();

  void execute_coarsening ();

  void fix_coarsen_flags ();

  std::vector<::internal::Triangulation::TriaLevel<dim>*> levels;

  std_cxx11::unique_ptr<::internal::Triangulation::TriaFaces<dim> > faces;


  std::vector<Point<spacedim> >              vertices;

  std::vector<bool>                     vertices_used;

  std::map<types::manifold_id, SmartPointer<const Manifold<dim,spacedim> , Triangulation<dim, spacedim> > >  manifold;


  bool                             anisotropic_refinement;


  const bool check_for_distorted_cells;

  ::internal::Triangulation::NumberCache<dim> number_cache;

  std_cxx11::unique_ptr<std::map<unsigned int, types::boundary_id> > vertex_to_boundary_id_map_1d;


  std_cxx11::unique_ptr<std::map<unsigned int, types::manifold_id> > vertex_to_manifold_id_map_1d;

  // make a couple of classes friends
  template <int,int,int> friend class TriaAccessorBase;
  template <int,int,int> friend class TriaAccessor;
  friend class TriaAccessor<0, 1, spacedim>;

  friend class CellAccessor<dim, spacedim>;

  friend struct ::internal::TriaAccessor::Implementation;

  friend class hp::DoFHandler<dim,spacedim>;

  friend struct ::internal::Triangulation::Implementation;

  template <typename>
  friend class ::internal::Triangulation::TriaObjects;

  // explicitly check for sensible template arguments, but not on windows
  // because MSVC creates bogus warnings during normal compilation
#ifdef DEAL_II_WITH_CXX11
#ifndef DEAL_II_MSVC
  static_assert (dim<=spacedim,
                 "The dimension <dim> of a Triangulation must be less than or "
                 "equal to the space dimension <spacedim> in which it lives.");
#endif
#endif

};


#ifndef DOXYGEN


template <int structdim>
inline
CellData<structdim>::CellData ()
{
  for (unsigned int i=0; i<GeometryInfo<structdim>::vertices_per_cell; ++i)
    vertices[i] = numbers::invalid_unsigned_int;

  material_id = 0;

  // And the manifold to be invalid
  manifold_id = numbers::invalid_manifold_id;
}



namespace internal
{
  namespace Triangulation
  {
    template <class Archive>
    void NumberCache<1>::serialize (Archive &ar,
                                    const unsigned int)
    {
      ar &n_levels;
      ar &n_lines &n_lines_level;
      ar &n_active_lines &n_active_lines_level;
    }


    template <class Archive>
    void NumberCache<2>::serialize (Archive &ar,
                                    const unsigned int version)
    {
      this->NumberCache<1>::serialize (ar, version);

      ar &n_quads &n_quads_level;
      ar &n_active_quads &n_active_quads_level;
    }


    template <class Archive>
    void NumberCache<3>::serialize (Archive &ar,
                                    const unsigned int version)
    {
      this->NumberCache<2>::serialize (ar, version);

      ar &n_hexes &n_hexes_level;
      ar &n_active_hexes &n_active_hexes_level;
    }

  }
}


template <int dim, int spacedim>
inline
bool
Triangulation<dim,spacedim>::vertex_used(const unsigned int index) const
{
  Assert (index < vertices_used.size(),
          ExcIndexRange(index, 0, vertices_used.size()));
  return vertices_used[index];
}



template <int dim, int spacedim>
inline
unsigned int Triangulation<dim, spacedim>::n_levels () const
{
  return number_cache.n_levels;
}

template <int dim, int spacedim>
inline
unsigned int Triangulation<dim, spacedim>::n_global_levels () const
{
  return number_cache.n_levels;
}


template <int dim, int spacedim>
inline
unsigned int
Triangulation<dim, spacedim>::n_vertices () const
{
  return vertices.size();
}



template <int dim, int spacedim>
inline
const std::vector<Point<spacedim> > &
Triangulation<dim, spacedim>::get_vertices () const
{
  return vertices;
}


template <int dim, int spacedim>
template <class Archive>
void
Triangulation<dim,spacedim>::save (Archive &ar,
                                   const unsigned int) const
{
  // as discussed in the documentation, do not store the signals as
  // well as boundary and manifold description but everything else
  ar &smooth_grid;
  ar &levels;
  ar &faces;
  ar &vertices;
  ar &vertices_used;

  ar &anisotropic_refinement;
  ar &number_cache;

  ar &check_for_distorted_cells;

  if (dim == 1)
    {
      ar &vertex_to_boundary_id_map_1d;
      ar &vertex_to_manifold_id_map_1d;
    }
}



template <int dim, int spacedim>
template <class Archive>
void
Triangulation<dim,spacedim>::load (Archive &ar,
                                   const unsigned int)
{
  // clear previous content. this also calls the respective signal
  clear ();

  // as discussed in the documentation, do not store the signals as
  // well as boundary and manifold description but everything else
  ar &smooth_grid;
  ar &levels;
  ar &faces;
  ar &vertices;
  ar &vertices_used;

  ar &anisotropic_refinement;
  ar &number_cache;

  // the levels do not serialize the active_cell_indices because
  // they are easy enough to rebuild upon re-loading data. do
  // this here. don't forget to first resize the fields appropriately
  {
    for (unsigned int l=0; l<levels.size(); ++l)
      levels[l]->active_cell_indices.resize (levels[l]->refine_flags.size());
    reset_active_cell_indices ();
  }


  bool my_check_for_distorted_cells;
  ar &my_check_for_distorted_cells;

  Assert (my_check_for_distorted_cells == check_for_distorted_cells,
          ExcMessage ("The triangulation loaded into here must have the "
                      "same setting with regard to reporting distorted "
                      "cell as the one previously stored."));

  if (dim == 1)
    {
      ar &vertex_to_boundary_id_map_1d;
      ar &vertex_to_manifold_id_map_1d;
    }

  // trigger the create signal to indicate
  // that new content has been imported into
  // the triangulation
  signals.create();
}


/* -------------- declaration of explicit specializations ------------- */

template <> unsigned int Triangulation<1,1>::n_raw_lines (const unsigned int level) const;
template <> unsigned int Triangulation<1,1>::n_quads () const;
template <> unsigned int Triangulation<1,1>::n_quads (const unsigned int level) const;
template <> unsigned int Triangulation<1,1>::n_raw_quads (const unsigned int level) const;
template <> unsigned int Triangulation<2,2>::n_raw_quads (const unsigned int level) const;
template <> unsigned int Triangulation<1,1>::n_raw_hexs (const unsigned int level) const;
template <> unsigned int Triangulation<1,1>::n_active_quads (const unsigned int level) const;
template <> unsigned int Triangulation<1,1>::n_active_quads () const;
template <> unsigned int Triangulation<1,1>::max_adjacent_cells () const;


// -------------------------------------------------------------------
// -- Explicit specializations for codimension one grids


template <> unsigned int Triangulation<1,2>::n_raw_lines (const unsigned int level) const;
template <> unsigned int Triangulation<1,2>::n_quads () const;
template <> unsigned int Triangulation<1,2>::n_quads (const unsigned int level) const;
template <> unsigned int Triangulation<1,2>::n_raw_quads (const unsigned int level) const;
template <> unsigned int Triangulation<2,3>::n_raw_quads (const unsigned int level) const;
template <> unsigned int Triangulation<1,2>::n_raw_hexs (const unsigned int level) const;
template <> unsigned int Triangulation<1,2>::n_active_quads (const unsigned int level) const;
template <> unsigned int Triangulation<1,2>::n_active_quads () const;
template <> unsigned int Triangulation<1,2>::max_adjacent_cells () const;

// -------------------------------------------------------------------
// -- Explicit specializations for codimension two grids


template <> unsigned int Triangulation<1,3>::n_raw_lines (const unsigned int level) const;
template <> unsigned int Triangulation<1,3>::n_quads () const;
template <> unsigned int Triangulation<1,3>::n_quads (const unsigned int level) const;
template <> unsigned int Triangulation<1,3>::n_raw_quads (const unsigned int level) const;
template <> unsigned int Triangulation<2,3>::n_raw_quads (const unsigned int level) const;
template <> unsigned int Triangulation<1,3>::n_raw_hexs (const unsigned int level) const;
template <> unsigned int Triangulation<1,3>::n_active_quads (const unsigned int level) const;
template <> unsigned int Triangulation<1,3>::n_active_quads () const;
template <> unsigned int Triangulation<1,3>::max_adjacent_cells () const;


#endif // DOXYGEN

DEAL_II_NAMESPACE_CLOSE

// Include tria_accessor.h here, so that it is possible for an end
// user to use the iterators of Triangulation<dim> directly without
// the need to include tria_accessor.h separately. (Otherwise the
// iterators are an 'opaque' or 'incomplete' type.)
#include <deal.II/grid/tria_accessor.h>

#endif