CD_EBReflux.H

Go to the documentation of this file.

/* chombo-discharge
 * Copyright © 2023 SINTEF Energy Research.
 * Please refer to Copyright.txt and LICENSE in the chombo-discharge root directory.
 */
 
#ifndef CD_EBReflux_H
#define CD_EBReflux_H
 
// Std includes
#include <map>
 
// Chombo includes
#include <Copier.H>
#include <EBLevelGrid.H>
 
// Our includes
#include <CD_NamespaceHeader.H>
 
class EBReflux
{
public:
  EBReflux() noexcept;
 
  EBReflux(const EBReflux& a_other) = delete;
 
  EBReflux(const EBLevelGrid& a_eblg,
           const EBLevelGrid& a_eblgFine,
           const EBLevelGrid& a_eblgCoFi,
           const int          a_refRat) noexcept;
 
  virtual ~EBReflux() noexcept;
 
  virtual void
  define(const EBLevelGrid& a_eblg,
         const EBLevelGrid& a_eblgFine,
         const EBLevelGrid& a_eblgCoFi,
         const int          a_refRat) noexcept;
 
  virtual void
  reflux(LevelData<EBCellFAB>&       a_Lphi,
         const LevelData<EBFluxFAB>& a_flux,
         const LevelData<EBFluxFAB>& a_fineFlux,
         const Interval              a_variables,
         const Real                  a_scaleCoarFlux,
         const Real                  a_scaleFineFlux) const noexcept;
 
protected:
  bool m_isDefined;
 
  int m_refRat;
 
  EBLevelGrid m_eblgFine;
 
  EBLevelGrid m_eblg;
 
  EBLevelGrid m_eblgCoFi;
 
  Copier m_copier;
 
  LayoutData<std::map<std::pair<int, Side::LoHiSide>, BaseIFFAB<FaceStencil>>> m_fluxCoarseningStencils;
 
  LayoutData<std::map<std::pair<int, Side::LoHiSide>, DenseIntVectSet>> m_regularCoarseFineRegions;
 
  mutable LayoutData<std::map<std::pair<int, Side::LoHiSide>, FaceIterator>> m_fluxCoarseningRegions;
 
  mutable LayoutData<std::map<std::pair<int, Side::LoHiSide>, VoFIterator>> m_irregularCoarseFineRegions;
 
  virtual void
  defineRegionsCF() noexcept;
 
  virtual void
  defineStencils() noexcept;
 
  virtual void
  defineBuffers() noexcept;
 
  virtual void
  coarsenFluxesCF(LevelData<EBFluxFAB>&       a_coarFluxes,
                  const LevelData<EBFluxFAB>& a_fineFluxes,
                  const int                   a_coarVar,
                  const int                   a_fineVar) const noexcept;
 
  virtual void
  refluxIntoCoarse(LevelData<EBCellFAB>&       a_Lphi,
                   const LevelData<EBFluxFAB>& a_oldFluxes,
                   const LevelData<EBFluxFAB>& a_newFluxes,
                   const int                   a_phiVar,
                   const int                   a_oldFluxVar,
                   const int                   a_newFluxVar,
                   const Real                  a_scaleCoarFlux,
                   const Real                  a_scaleFineFlux) const noexcept;
};
 
#include <CD_NamespaceFooter.H>
 
#endif