CD_EBGhostCellInterpolator.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_EBGHOSTCELLINTERPOLATOR_H
#define CD_EBGHOSTCELLINTERPOLATOR_H
 
// Std includes
#include <map>
 
// Chombo includes
#include <EBLevelGrid.H>
 
// Our includes
#include <CD_Location.H>
#include <CD_NamespaceHeader.H>
 
class EBGhostCellInterpolator
{
public:
  enum Type
  {
    PWC,
    MinMod,
    MonotonizedCentral,
    Superbee
  };
 
  EBGhostCellInterpolator() noexcept;
 
  EBGhostCellInterpolator(const EBGhostCellInterpolator& a_other) = delete;
 
  EBGhostCellInterpolator(const EBLevelGrid& a_eblgFine,
                          const EBLevelGrid& a_eblgCoFi,
                          const EBLevelGrid& a_eblgCoar,
                          const IntVect&     a_ghostVector,
                          const int          a_refRat,
                          const int          a_ghostCF) noexcept;
 
  virtual ~EBGhostCellInterpolator() noexcept;
 
  virtual void
  define(const EBLevelGrid& a_eblgFine,
         const EBLevelGrid& a_eblgCoFi,
         const EBLevelGrid& a_eblgCoar,
         const IntVect&     a_ghostVector,
         const int          a_refRat,
         const int          a_ghostCF) noexcept;
 
  virtual void
  interpolate(LevelData<EBCellFAB>&       a_phiFine,
              const LevelData<EBCellFAB>& a_phiCoar,
              const Interval              a_variables,
              const Type                  a_interpType) const noexcept;
 
protected:
  bool m_isDefined;
 
  EBLevelGrid m_eblgFine;
 
  EBLevelGrid m_eblgCoar;
 
  EBLevelGrid m_eblgCoFi;
 
  IntVect m_ghostVector;
 
  int m_refRat;
 
  int m_ghostCF;
 
  LayoutData<std::map<std::pair<int, Side::LoHiSide>, Box>> m_regularGhostRegions;
 
  Copier m_copier;
 
  LayoutData<BaseIVFAB<VolIndex>> m_coarsenedFineGhosts;
 
  mutable LayoutData<VoFIterator> m_fineIrregCells;
 
  mutable LayoutData<VoFIterator> m_coarIrregCells;
 
  mutable LayoutData<BaseIVFAB<Real>> m_coarIrregSlopes;
 
  virtual void
  defineGhostRegions() noexcept;
 
  virtual void
  defineBuffers() noexcept;
 
  virtual void
  interpolateRegular(FArrayBox&       a_phiFine,
                     const FArrayBox& a_phiCoar,
                     const DataIndex& a_dit,
                     const int        a_fineVar,
                     const int        a_coarVar,
                     const Type       a_interpType) const noexcept;
 
  virtual void
  interpolateIrregular(EBCellFAB&       a_phiFine,
                       const EBCellFAB& a_phiCoar,
                       const DataIndex& a_dit,
                       const int        a_fineVar,
                       const int        a_coarVar,
                       const Type       a_interpType) const noexcept;
 
  static inline Real
  minmod(const Real& dwl, const Real& dwr) noexcept;
 
  static Real
  superbee(const Real& dwl, const Real& dwr) noexcept;
 
  static inline Real
  monotonizedCentral(const Real& dwl, const Real& dwr) noexcept;
};
 
#include <CD_NamespaceFooter.H>
 
#endif