tria_base.h

Go to the documentation of this file.

// ---------------------------------------------------------------------
//
// Copyright (C) 2008 - 2020 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_distributed_tria_base_h
#define dealii_distributed_tria_base_h


#include <deal.II/base/config.h>

#include <deal.II/base/mpi.h>
#include <deal.II/base/smartpointer.h>
#include <deal.II/base/subscriptor.h>
#include <deal.II/base/template_constraints.h>

#include <deal.II/grid/tria.h>

#include <functional>
#include <list>
#include <set>
#include <tuple>
#include <utility>
#include <vector>


DEAL_II_NAMESPACE_OPEN

namespace parallel
{
  template <int dim, int spacedim = dim>
  class TriangulationBase : public ::Triangulation<dim, spacedim>
  {
  public:
    TriangulationBase(
      MPI_Comm mpi_communicator,
      const typename ::Triangulation<dim, spacedim>::MeshSmoothing
                 smooth_grid = (::Triangulation<dim, spacedim>::none),
      const bool check_for_distorted_cells = false);

    virtual ~TriangulationBase() override;

    virtual const MPI_Comm &
    get_communicator() const;

    virtual bool
    is_multilevel_hierarchy_constructed() const = 0;

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

    unsigned int
    n_locally_owned_active_cells() const;

    virtual types::global_cell_index
    n_global_active_cells() const override;

    virtual std::size_t
    memory_consumption() const override;


    virtual unsigned int
    n_global_levels() const override;

    types::subdomain_id
    locally_owned_subdomain() const override;

    const std::set<types::subdomain_id> &
    ghost_owners() const;

    const std::set<types::subdomain_id> &
    level_ghost_owners() const;

    DEAL_II_DEPRECATED virtual std::map<unsigned int,
                                        std::set<::types::subdomain_id>>
    compute_vertices_with_ghost_neighbors() const;

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

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

  protected:
    MPI_Comm mpi_communicator;

    types::subdomain_id my_subdomain;

    types::subdomain_id n_subdomains;

    struct NumberCache
    {
      unsigned int n_locally_owned_active_cells;
      types::global_cell_index n_global_active_cells;
      unsigned int n_global_levels;
      std::set<types::subdomain_id> ghost_owners;
      std::set<types::subdomain_id> level_ghost_owners;

      NumberCache();
    };

    NumberCache number_cache;

    virtual void
    update_number_cache();
  };

  template <int dim, int spacedim = dim>
  using Triangulation DEAL_II_DEPRECATED = TriangulationBase<dim, spacedim>;



  template <int dim, int spacedim = dim>
  class DistributedTriangulationBase
    : public ::parallel::TriangulationBase<dim, spacedim>
  {
  public:
    DistributedTriangulationBase(
      MPI_Comm mpi_communicator,
      const typename ::Triangulation<dim, spacedim>::MeshSmoothing
                 smooth_grid = (::Triangulation<dim, spacedim>::none),
      const bool check_for_distorted_cells = false);
  };

} // namespace parallel

DEAL_II_NAMESPACE_CLOSE

#endif