CD_PolyhedralEBGraph.H

Go to the documentation of this file.

/*
 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
 *
 * SPDX-License-Identifier: GPL-3.0-or-later
 */
 
#ifndef CD_POLYHEDRALEBGRAPH_H
#define CD_POLYHEDRALEBGRAPH_H
 
// Chombo includes
#include <BaseFab.H>
#include <BaseIF.H>
#include <DisjointBoxLayout.H>
#include <IVSFAB.H>
#include <IntVectSet.H>
#include <LayoutData.H>
#include <LevelData.H>
#include <ProblemDomain.H>
#include <RealVect.H>
 
// Our includes
#include <CD_CutCellSurface.H>
#include <CD_NamespaceHeader.H>
 
class PolyhedralEBGraph
{
public:
  static constexpr int s_covered = -1;
 
  static constexpr int s_cut = 0;
 
  static constexpr int s_regular = 1;
 
  static constexpr int s_faceClosed = 0;
 
  static constexpr int s_faceSameLevel = 1;
 
  static constexpr int s_faceFiner = 2;
 
  static constexpr int s_faceCoarser = 3;
 
  static constexpr int s_faceBoundary = 4;
 
  PolyhedralEBGraph();
 
  virtual ~PolyhedralEBGraph();
 
  virtual void
  define(const BaseIF&        a_function,
         const Vector<Box>&   a_cutTiles,
         const ProblemDomain& a_domain,
         const RealVect&      a_probLo,
         const Real           a_dx,
         const int            a_numGhost,
         const Vector<Box>&   a_covered);
 
  virtual void
  define(const PolyhedralEBGraph& a_source, const DisjointBoxLayout& a_grids, const BaseIF& a_function);
 
  bool
  equals(const PolyhedralEBGraph& a_other) const;
 
  static void
  link(PolyhedralEBGraph& a_coarse, const PolyhedralEBGraph& a_fine);
 
  static void
  findUnresolvedCells(const BaseFab<Real>& a_nodeValues,
                      const Box&           a_region,
                      const Vector<Box>&   a_covered,
                      BaseFab<bool>&       a_unresolved);
 
  bool
  isDefined() const noexcept;
 
  const ProblemDomain&
  getDomain() const noexcept;
 
  Real
  getDx() const noexcept;
 
  const DisjointBoxLayout&
  getGrids() const noexcept;
 
  const LayoutData<IntVectSet>&
  getCutCells() const noexcept;
 
  const LevelData<BaseFab<signed char>>&
  getCellStates() const noexcept;
 
  const LevelData<BaseFab<signed char>>&
  getFaceStates() const noexcept;
 
  const LevelData<BaseFab<signed char>>&
  getRefinedMask() const noexcept;
 
  const LevelData<IVSFAB<PolyhedralEB::CutCellSurface>>&
  getSurfaces() const noexcept;
 
protected:
  bool m_isDefined;
 
  int m_numGhost;
 
  Real m_dx;
 
  RealVect m_probLo;
 
  ProblemDomain m_domain;
 
  DisjointBoxLayout m_grids;
 
  Vector<Box> m_covered;
 
  LayoutData<IntVectSet> m_cutCells;
 
  LevelData<BaseFab<signed char>> m_cellStates;
 
  LevelData<BaseFab<signed char>> m_faceStates;
 
  LevelData<BaseFab<signed char>> m_refined;
 
  LevelData<IVSFAB<PolyhedralEB::CutCellSurface>> m_surfaces;
 
  void
  defineGrids(const Vector<Box>& a_cutTiles);
 
  void
  defineData();
 
  void
  defineGhostCells(const BaseIF& a_function, const LevelData<BaseFab<signed char>>& a_carried);
 
  void
  markCarried(LevelData<BaseFab<signed char>>& a_carried) const;
 
  void
  defineCells(const BaseIF& a_function, const LevelData<BaseFab<signed char>>& a_carried);
 
  void
  defineOuterFaces(const LevelData<BaseFab<signed char>>& a_carried);
};
 
#include <CD_NamespaceFooter.H>
 
#endif