CD_EBCoarseFineParticleMesh.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_EBCOARSEFINEPARTICLEMESH_H
#define CD_EBCOARSEFINEPARTICLEMESH_H
 
// Chombo includes
#include <DisjointBoxLayout.H>
#include <EBLevelGrid.H>
#include <EBCellFAB.H>
#include <ProblemDomain.H>
 
// Our includes
#include <CD_EBCoarAve.H>
#include <CD_NamespaceHeader.H>
 
class EBCoarseFineParticleMesh
{
public:
  enum Average
  {
    Arithmetic,
    Conservative
  };
 
  EBCoarseFineParticleMesh() noexcept;
 
  EBCoarseFineParticleMesh(const EBLevelGrid& a_eblgCoar,
                           const EBLevelGrid& a_eblgFine,
                           const int          a_refRat,
                           const IntVect      a_ghost) noexcept;
 
  virtual ~EBCoarseFineParticleMesh() noexcept;
 
  virtual void
  define(const EBLevelGrid& a_eblgCoar,
         const EBLevelGrid& a_eblgFine,
         const int          a_refRat,
         const IntVect      a_ghost) noexcept;
 
  virtual void
  addFineGhostsToCoarse(LevelData<EBCellFAB>& a_coarData, const LevelData<EBCellFAB>& a_fineData) const noexcept;
 
  virtual void
  addFiCoDataToFine(LevelData<EBCellFAB>& a_fineData, const LevelData<EBCellFAB>& a_fiCoData) const noexcept;
 
  virtual void
  addInvalidCoarseToFine(LevelData<EBCellFAB>& a_fineData, const LevelData<EBCellFAB>& a_coarData) const noexcept;
 
  virtual void
  exchangeAndAddFiCoData(LevelData<EBCellFAB>& a_fiCoData) const noexcept;
 
  virtual void
  restrictAndAddFiCoDataToCoar(LevelData<EBCellFAB>&                   a_coarData,
                               const LevelData<EBCellFAB>&             a_fiCoData,
                               const EBCoarseFineParticleMesh::Average a_average) const noexcept;
 
  const EBLevelGrid&
  getEblgFiCo() const;
 
  template <int NCOMP>
  LevelData<EBCellFAB>&
  getBufferFiCo() const noexcept;
 
protected:
  static constexpr int m_comp = 0;
 
  static constexpr int m_nComp = 1;
 
  bool m_isDefined;
 
  bool m_verbose;
 
  EBLevelGrid m_eblgCoar;
 
  EBLevelGrid m_eblgFine;
 
  EBLevelGrid m_eblgFiCo;
 
  EBLevelGrid m_eblgCoFi;
 
  Copier m_copierFiCoToFineIncludeGhosts;
 
  Copier m_copierCoFiToCoarIncludeGhosts;
 
  Copier m_copierFiCoToFineNoGhosts;
 
  Copier m_copierFiCoToFiCo;
 
  mutable LevelData<EBCellFAB> m_bufferFiCoReal;
 
  mutable LevelData<EBCellFAB> m_bufferFiCoRealVect;
 
  mutable LayoutData<VoFIterator> m_vofIterFineGhosts;
 
  mutable LayoutData<VoFIterator> m_vofIterCoFiGhosts;
 
  mutable LayoutData<VoFIterator> m_vofIterCoar;
 
  LayoutData<BaseIVFAB<VoFStencil>> m_cellArithmeticStencils;
 
  LayoutData<BaseIVFAB<VoFStencil>> m_cellConservativeStencils;
 
  int m_refRat;
 
  IntVect m_ghost;
 
  void
  defineVoFIterators() noexcept;
 
  void
  defineStencils() noexcept;
 
  void
  defineBufferFiCo() noexcept;
 
  void
  conservativeAverageAndAdd(EBCellFAB& a_coarData, const EBCellFAB& a_fineData, const DataIndex& a_din) const noexcept;
 
  void
  arithmeticAverageAndAdd(EBCellFAB& a_coarData, const EBCellFAB& a_fineData, const DataIndex& a_din) const noexcept;
};
 
#include <CD_NamespaceFooter.H>
 
#endif