rectangular_mesh.hh

Go to the documentation of this file.

#pragma once
 
// external libraries
#include <deal.II/base/config.h>
#include <deal.II/grid/tria.h>
#ifdef DEAL_II_WITH_MPI
#include <deal.II/distributed/shared_tria.h>
#endif
 
// standard library
#include <type_traits>
 
// DiFfRG
#include <DiFfRG/common/linear_algebra.hh>
#include <DiFfRG/common/utils.hh>
#include <DiFfRG/discretization/mesh/configuration_mesh.hh>
 
namespace DiFfRG
{
  using namespace dealii;
 
  struct RectangularMeshOptions {
    bool origin_cell_centered = false;
  };
 
  namespace internal
  {
    template <typename Tria> struct TriangulationHolder {
      Tria triangulation;
    };
 
#ifdef DEAL_II_WITH_MPI
    template <int dim> struct TriangulationHolder<dealii::parallel::shared::Triangulation<dim>> {
      dealii::parallel::shared::Triangulation<dim> triangulation{MPI_COMM_WORLD, dealii::Triangulation<dim>::none,
                                                                 /*allow_artificial_cells=*/false};
    };
#endif
  } // namespace internal
 
  template <uint dim_, typename TriangulationType_ = DefaultTriangulation<dim_>>
  class RectangularMesh : protected internal::TriangulationHolder<TriangulationType_>
  {
  public:
    static constexpr uint dim = dim_;
    using TriangulationType = TriangulationType_;
    static constexpr bool is_rectangular = true;
 
    static constexpr bool is_parallel = !std::is_same_v<TriangulationType, dealii::Triangulation<dim_>>;
 
    [[deprecated("Please use RectangularMesh(const Config::ConfigurationMesh<dim> &mesh_config) instead")]]
    RectangularMesh(const ConfigTree &config);
 
    [[deprecated("Please use RectangularMesh(const Config::ConfigurationMesh<dim> &mesh_config, "
                 "RectangularMeshOptions options) instead")]]
    RectangularMesh(const ConfigTree &config, RectangularMeshOptions options);
 
    RectangularMesh(const Config::ConfigurationMesh<dim> &mesh_config);
 
    RectangularMesh(const Config::ConfigurationMesh<dim> &mesh_config, RectangularMeshOptions options);
 
    TriangulationType &get_triangulation() { return this->triangulation; }
    const TriangulationType &get_triangulation() const { return this->triangulation; }
 
  protected:
    virtual void make_grid();
 
    const Config::ConfigurationMesh<dim> mesh_config;
    const RectangularMeshOptions options;
  };
 
  template <uint dim> using RectangularMeshSerial = RectangularMesh<dim, dealii::Triangulation<dim>>;
 
#ifdef DEAL_II_WITH_MPI
  template <uint dim>
  using RectangularMeshParallel = RectangularMesh<dim, dealii::parallel::shared::Triangulation<dim>>;
#endif
} // namespace DiFfRG