persistent_tria.h

// ---------------------------------------------------------------------
//
// Copyright (C) 1999 - 2018 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_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() override = default;

  virtual void
  execute_coarsening_and_refinement() override;

  void
  restore();

  void
  restore(const unsigned int step_no);

  unsigned int
  n_refinement_steps() const;

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

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

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

  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 override;

  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