tria_levels.h

Go to the documentation of this file.

// ------------------------------------------------------------------------
//
// SPDX-License-Identifier: LGPL-2.1-or-later
// Copyright (C) 1998 - 2024 by the deal.II authors
//
// This file is part of the deal.II library.
//
// Part of the source code is dual licensed under Apache-2.0 WITH
// LLVM-exception OR LGPL-2.1-or-later. Detailed license information
// governing the source code and code contributions can be found in
// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II.
//
// ------------------------------------------------------------------------
 
#ifndef dealii_tria_levels_h
#define dealii_tria_levels_h
 
 
#include <deal.II/base/config.h>
 
#include <deal.II/base/point.h>
 
#include <deal.II/grid/reference_cell.h>
#include <deal.II/grid/tria_objects.h>
#include <deal.II/grid/tria_objects_orientations.h>
 
#include <boost/serialization/utility.hpp>
 
#include <cstdint>
#include <vector>
 
DEAL_II_NAMESPACE_OPEN
 
namespace internal
{
  namespace TriangulationImplementation
  {
    class TriaLevel
    {
    public:
      TriaLevel(const unsigned int dim)
        : dim(dim)
        , cells(dim)
      {}
 
      TriaLevel()
        : dim(numbers::invalid_unsigned_int)
        , cells(numbers::invalid_unsigned_int)
      {}
 
      unsigned int dim;
 
      std::vector<std::uint8_t> refine_flags;
 
      std::vector<std::uint8_t> refine_choice;
 
      std::vector<bool> coarsen_flags;
 
 
      std::vector<unsigned int> active_cell_indices;
 
      std::vector<types::global_cell_index> global_active_cell_indices;
 
      std::vector<types::global_cell_index> global_level_cell_indices;
 
      std::vector<std::pair<int, int>> neighbors;
 
      std::vector<types::subdomain_id> subdomain_ids;
 
      std::vector<types::subdomain_id> level_subdomain_ids;
 
      std::vector<int> parents;
 
      std::vector<bool> direction_flags;
 
      TriaObjects cells;
 
      TriaObjectsOrientations face_orientations;
 
      std::vector<ReferenceCell> reference_cell;
 
      std::vector<unsigned int> cell_vertex_indices_cache;
 
      std::size_t
      memory_consumption() const;
 
      template <class Archive>
      void
      serialize(Archive &ar, const unsigned int version);
    };
 
 
    template <class Archive>
    void
    TriaLevel::serialize(Archive &ar, const unsigned int)
    {
      ar &dim;
 
      ar &refine_flags &coarsen_flags;
 
      // do not serialize `active_cell_indices` and `vertex_indices_cache`
      // here. instead of storing them to the stream and re-reading them again
      // later, we just rebuild them in Triangulation::load()
 
      ar &neighbors;
      ar &subdomain_ids;
      ar &level_subdomain_ids;
      ar &parents;
      ar &direction_flags;
      ar &cells;
      ar &face_orientations;
      ar &reference_cell;
    }
 
  } // namespace TriangulationImplementation
} // namespace internal
 
DEAL_II_NAMESPACE_CLOSE
 
#endif