persistent_tria.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_persistent_tria_h
#define dealii_persistent_tria_h


#include <deal.II/base/config.h>
#include <deal.II/base/smartpointer.h>
#include <deal.II/grid/tria.h>
#include <vector>

DEAL_II_NAMESPACE_OPEN

template <int dim, int spacedim=dim>
class PersistentTriangulation : public Triangulation<dim, spacedim>
{
public:
  static const unsigned int dimension = dim;
  static const unsigned int spacedimension = spacedim;

  PersistentTriangulation (const Triangulation<dim, spacedim> &coarse_grid);

  PersistentTriangulation (const PersistentTriangulation<dim, spacedim> &old_tria);

  virtual ~PersistentTriangulation () = default;

  virtual void execute_coarsening_and_refinement ();

  void restore ();

  void restore (const unsigned int step_no);

  unsigned int n_refinement_steps () const;

  virtual void copy_triangulation (const Triangulation<dim, spacedim> &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);

  virtual void write_flags(std::ostream &out) const;

  virtual void read_flags(std::istream &in);

  virtual void clear_flags();

  virtual std::size_t memory_consumption () const;

  DeclException0 (ExcTriaNotEmpty);
  DeclException0 (ExcFlagsNotCleared);

private:
  SmartPointer<const Triangulation<dim,spacedim>,PersistentTriangulation<dim,spacedim> > coarse_grid;

  std::vector<std::vector<bool> >   refine_flags;

  std::vector<std::vector<bool> >   coarsen_flags;
};


DEAL_II_NAMESPACE_CLOSE

#endif