CD_EBCoarseToFineInterp.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_EBCOARSETOFINEINTERP_H
#define CD_EBCOARSETOFINEINTERP_H
 
// Chombo includes
#include <EBLevelGrid.H>
#include <LevelData.H>
#include <EBCellFAB.H>
#include <BaseIVFAB.H>
 
// Our includes
#include <CD_NamespaceHeader.H>
 
class EBCoarseToFineInterp
{
public:
  enum Type
  {
    PWC,
    ConservativePWC,
    ConservativeMinMod,
    ConservativeMonotonizedCentral,
    ConservativeSuperbee,
  };
 
  EBCoarseToFineInterp() noexcept;
 
  EBCoarseToFineInterp(const EBLevelGrid& a_eblgFine,
                       const EBLevelGrid& a_eblgCoFi,
                       const EBLevelGrid& a_eblgCoar,
                       const int&         a_refRat) noexcept;
 
  EBCoarseToFineInterp(const EBCoarseToFineInterp&) = delete;
 
  EBCoarseToFineInterp(const EBCoarseToFineInterp&&) = delete;
 
  virtual ~EBCoarseToFineInterp() noexcept;
 
  EBCoarseToFineInterp&
  operator=(const EBCoarseToFineInterp&) = delete;
 
  EBCoarseToFineInterp&
  operator=(const EBCoarseToFineInterp&&) = delete;
 
  virtual void
  define(const EBLevelGrid& a_eblgFine,
         const EBLevelGrid& a_eblgCoFi,
         const EBLevelGrid& a_eblgCoar,
         const int&         a_refRat) noexcept;
 
  virtual void
  interpolate(LevelData<EBCellFAB>&             a_fineData,
              const LevelData<EBCellFAB>&       a_coarData,
              const Interval&                   a_variables,
              const EBCoarseToFineInterp::Type& a_interpType) const noexcept;
 
  virtual void
  interpolate(LevelData<BaseIVFAB<Real>>&       a_fineData,
              const LevelData<BaseIVFAB<Real>>& a_coarData,
              const Interval&                   a_variables,
              const EBCoarseToFineInterp::Type& a_interpType) const noexcept;
 
protected:
  enum class SlopeLimiter
  {
    MinMod,
    MonotonizedCentral,
    Superbee
  };
 
  bool m_isDefined;
 
  int m_refRat;
 
  mutable Copier m_ebCopier;
 
  mutable std::map<IntVect, Copier> m_cellCopiers;
 
  EBLevelGrid m_eblgFine;
 
  EBLevelGrid m_eblgCoar;
 
  EBLevelGrid m_eblgCoFi;
 
  LevelData<BaseIVFAB<Real>> m_volumeWeights;
 
  LevelData<BaseIVFAB<Real>> m_areaWeights;
 
  mutable LevelData<BaseIVFAB<Real>> m_irregCoFi;
 
  mutable LayoutData<VoFIterator> m_fineVoFs;
 
  mutable LayoutData<VoFIterator> m_coarVoFs;
 
  virtual void
  defineWeights() noexcept;
 
  virtual void
  interpolatePWC(EBCellFAB&       a_fineData,
                 const EBCellFAB& a_coarData,
                 const DataIndex& a_dit,
                 const int&       a_fineVar,
                 const int&       a_coarVar) const noexcept;
 
  virtual void
  interpolateConservativePWC(EBCellFAB&       a_fineData,
                             const EBCellFAB& a_coarData,
                             const DataIndex& a_dit,
                             const int&       a_fineVar,
                             const int&       a_coarVar) const noexcept;
 
  virtual void
  interpolateConservativeSlope(EBCellFAB&          a_fineData,
                               const EBCellFAB&    a_coarData,
                               const DataIndex&    a_dit,
                               const int&          a_fineVar,
                               const int&          a_coarVar,
                               const SlopeLimiter& a_limiter) const noexcept;
 
  virtual void
  interpolatePWC(BaseIVFAB<Real>&       a_fineData,
                 const BaseIVFAB<Real>& a_coarData,
                 const DataIndex&       a_dit,
                 const int&             a_fineVar,
                 const int&             a_coarVar) const noexcept;
 
  virtual void
  interpolateConservativePWC(BaseIVFAB<Real>&       a_fineData,
                             const BaseIVFAB<Real>& a_coarData,
                             const DataIndex&       a_dit,
                             const int&             a_fineVar,
                             const int&             a_coarVar) const noexcept;
 
  virtual void
  checkConservation(const EBCellFAB& a_fineData,
                    const EBCellFAB& a_coarData,
                    const DataIndex& a_dit,
                    const int        a_fineVar,
                    const int        a_coarVar) const noexcept;
 
  virtual void
  checkConservation(const BaseIVFAB<Real>& a_fineData,
                    const BaseIVFAB<Real>& a_coarData,
                    const DataIndex&       a_dit,
                    const int              a_fineVar,
                    const int              a_coarVar) const noexcept;
};
 
#include <CD_NamespaceFooter.H>
 
#endif