persistent_tria.h

Go to the documentation of this file.

// -----------------------------------------------------------------------------
//
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception OR LGPL-2.1-or-later
// Copyright (C) 1999 - 2024 by the deal.II authors
//
// This file is part of the deal.II library.
//
// Detailed license information governing the source code and contributions
// can be found in LICENSE.md and CONTRIBUTING.md at the top level directory.
//
// -----------------------------------------------------------------------------
 
#ifndef dealii_persistent_tria_h
#define dealii_persistent_tria_h
 
 
#include <deal.II/base/config.h>
 
#include <deal.II/base/observer_pointer.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 constexpr 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(
    const TriangulationDescription::Description<dim, spacedim>
      &construction_data) 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:
  ObserverPointer<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