CD_EBCoarAve.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_EBCoarAve_H
#define CD_EBCoarAve_H
 
// Chombo includes
#include <AggStencil.H>
#include <EBLevelGrid.H>
#include <VoFIterator.H>
#include <FaceIterator.H>
 
// Our includes
#include <CD_Average.H>
#include <CD_NamespaceHeader.H>
 
class EBCoarAve
{
public:
  enum class Type
  {
    Arithmetic,
    Conservative,
    Harmonic
  };
 
  EBCoarAve() noexcept;
 
  EBCoarAve(const EBCoarAve& a_other) = delete;
 
  EBCoarAve(const DisjointBoxLayout& a_dblFine,
            const DisjointBoxLayout& a_dblCoar,
            const EBISLayout&        a_ebislFine,
            const EBISLayout&        a_ebislCoar,
            const ProblemDomain&     a_domainCoar,
            const int&               a_refRat,
            const EBIndexSpace*      a_ebisPtr) noexcept;
 
  EBCoarAve(const EBLevelGrid& a_eblgFine,
            const EBLevelGrid& a_eblgCoar,
            const EBLevelGrid& a_eblgCoFi,
            const int&         a_refRat) noexcept;
 
  virtual ~EBCoarAve() noexcept;
 
  virtual void
  define(const EBLevelGrid& a_eblgFine,
         const EBLevelGrid& a_eblgCoar,
         const EBLevelGrid& a_eblgCoFi,
         const int&         a_refRat) noexcept;
 
  virtual void
  operator=(const EBCoarAve& fabin) = delete;
 
  virtual void
  averageData(LevelData<EBCellFAB>&       a_coarData,
              const LevelData<EBCellFAB>& a_fineData,
              const Interval&             a_variables,
              const Average&              a_average) const noexcept;
 
  virtual void
  averageData(LevelData<EBFluxFAB>&       a_coarData,
              const LevelData<EBFluxFAB>& a_fineData,
              const Interval&             a_variables,
              const Average&              a_average) const noexcept;
 
  virtual void
  averageData(LevelData<BaseIVFAB<Real>>&       a_coarData,
              const LevelData<BaseIVFAB<Real>>& a_fineData,
              const Interval&                   a_variables,
              const Average&                    a_average) const noexcept;
 
protected:
  bool m_isDefined;
 
  int m_refRat;
 
  EBLevelGrid m_eblgFine;
 
  EBLevelGrid m_eblgCoar;
 
  EBLevelGrid m_eblgCoFi;
 
  mutable std::map<IntVect, Copier> m_cellCopiers;
 
  mutable std::map<IntVect, Copier> m_faceCopiers;
 
  mutable std::map<IntVect, Copier> m_ebCopiers;
 
  LayoutData<IntVectSet> m_irregSetsCoFi;
 
  mutable LayoutData<VoFIterator> m_irregCellsCoFi;
 
  mutable LayoutData<std::map<int, FaceIterator>> m_irregFacesCoFi;
 
  LayoutData<BaseIVFAB<VoFStencil>> m_cellArithmeticStencils;
 
  LayoutData<BaseIVFAB<VoFStencil>> m_cellHarmonicStencils;
 
  LayoutData<BaseIVFAB<VoFStencil>> m_cellConservativeStencils;
 
  LayoutData<std::map<int, BaseIFFAB<FaceStencil>>> m_faceArithmeticStencils;
 
  LayoutData<std::map<int, BaseIFFAB<FaceStencil>>> m_faceHarmonicStencils;
 
  LayoutData<std::map<int, BaseIFFAB<FaceStencil>>> m_faceConservativeStencils;
 
  LayoutData<BaseIVFAB<VoFStencil>> m_ebArithmeticStencils;
 
  LayoutData<BaseIVFAB<VoFStencil>> m_ebHarmonicStencils;
 
  LayoutData<BaseIVFAB<VoFStencil>> m_ebConservativeStencils;
 
  virtual void
  defineCellStencils() noexcept;
 
  virtual void
  defineFaceStencils() noexcept;
 
  virtual void
  defineEBStencils() noexcept;
 
  virtual void
  defineBuffers() noexcept;
 
  virtual void
  arithmeticAverage(EBCellFAB&       a_coarData,
                    const EBCellFAB& a_fineData,
                    const DataIndex& a_datInd,
                    const int&       a_coarVar,
                    const int&       a_fineVar) const noexcept;
  virtual void
  harmonicAverage(EBCellFAB&       a_coarData,
                  const EBCellFAB& a_fineData,
                  const DataIndex& a_datInd,
                  const int&       a_coarVar,
                  const int&       a_fineVar) const noexcept;
 
  virtual void
  conservativeAverage(EBCellFAB&       a_coarData,
                      const EBCellFAB& a_fineData,
                      const DataIndex& a_datInd,
                      const int&       a_coarVar,
                      const int&       a_fineVar) const noexcept;
 
  virtual void
  arithmeticAverage(EBFaceFAB&       a_coarData,
                    const EBFaceFAB& a_fineData,
                    const DataIndex& a_datInd,
                    const int&       a_fineVar,
                    const int&       a_coarVar,
                    const int&       a_dir) const noexcept;
 
  virtual void
  harmonicAverage(EBFaceFAB&       a_coarData,
                  const EBFaceFAB& a_fineData,
                  const DataIndex& a_datInd,
                  const int&       a_coarVar,
                  const int&       a_fineVar,
                  const int&       a_dir) const noexcept;
 
  virtual void
  conservativeAverage(EBFaceFAB&       a_coarData,
                      const EBFaceFAB& a_fineData,
                      const DataIndex& a_datInd,
                      const int&       a_coarVar,
                      const int&       a_fineVar,
                      const int&       a_dir) const noexcept;
 
  virtual void
  arithmeticAverage(BaseIVFAB<Real>&       a_coarData,
                    const BaseIVFAB<Real>& a_fineData,
                    const DataIndex&       a_datInd,
                    const int&             a_coarVar,
                    const int&             a_fineVar) const noexcept;
 
  virtual void
  harmonicAverage(BaseIVFAB<Real>&       a_coarData,
                  const BaseIVFAB<Real>& a_fineData,
                  const DataIndex&       a_datInd,
                  const int&             a_coarVar,
                  const int&             a_fineVar) const noexcept;
 
  virtual void
  conservativeAverage(BaseIVFAB<Real>&       a_coarData,
                      const BaseIVFAB<Real>& a_fineData,
                      const DataIndex&       a_datInd,
                      const int&             a_coarVar,
                      const int&             a_fineVar) const noexcept;
};
 
#include <CD_NamespaceFooter.H>
 
#endif