CD_EBLeastSquaresMultigridInterpolator.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_EBLEASTSQUARESMULTIGRIDINTERPOLATOR_H
#define CD_EBLEASTSQUARESMULTIGRIDINTERPOLATOR_H
 
// Std includes
#include <map>
 
// Chombo includes
#include <EBLevelGrid.H>
#include <QuadCFInterp.H>
#include <AggStencil.H>
 
// Our includes
#include <CD_Location.H>
#include <CD_EBMultigridInterpolator.H>
#include <CD_CoarseInterpQuadCF.H>
#include <CD_NamespaceHeader.H>
 
class EBLeastSquaresMultigridInterpolator : public EBMultigridInterpolator
{
public:
  using CellLocation = Location::Cell;
 
  EBLeastSquaresMultigridInterpolator() = delete;
 
  EBLeastSquaresMultigridInterpolator(const EBLeastSquaresMultigridInterpolator& a_other) = delete;
 
  EBLeastSquaresMultigridInterpolator(const EBLevelGrid& a_eblgFine,
                                      const EBLevelGrid& a_eblgCoFi,
                                      const EBLevelGrid& a_eblgCoar,
                                      const CellLocation a_dataLocation,
                                      const IntVect&     a_ghostVector,
                                      const int          a_refRat,
                                      const int          a_ghostCF,
                                      const int          a_order,
                                      const int          a_weighting = 0) noexcept;
 
  virtual ~EBLeastSquaresMultigridInterpolator() noexcept;
 
  virtual int
  getGhostCF() const noexcept override;
 
  virtual std::pair<DerivStencil, DerivStencil>
  getInterpolationStencilRegular(const IntVect&       a_fineGhost,
                                 const DataIndex&     a_din,
                                 const int            a_dir,
                                 const Side::LoHiSide a_side) const noexcept override;
 
  virtual std::pair<VoFStencil, VoFStencil>
  getInterpolationStencilEB(const VolIndex& a_fineGhost, const DataIndex& a_din) const noexcept override;
 
  virtual void
  coarseFineInterp(LevelData<EBCellFAB>&       a_phiFine,
                   const LevelData<EBCellFAB>& a_phiCoar,
                   const Interval              a_variables) const noexcept override;
 
  virtual void
  coarseFineInterpH(LevelData<EBCellFAB>& a_phiFine, const Interval a_variables) const noexcept override;
 
  virtual void
  coarseFineInterpH(EBCellFAB& a_phiFine, const Interval a_variables, const DataIndex& a_din) const noexcept override;
 
protected:
  static constexpr int m_stenComp = 0;
 
  static constexpr int m_numStenComp = 1;
 
  static constexpr int m_comp = 0;
 
  CellLocation m_dataLocation;
 
  EBLevelGrid m_eblgFine;
 
  EBLevelGrid m_eblgCoar;
 
  EBLevelGrid m_eblgCoFi;
 
  Copier m_copier;
 
  mutable LevelData<EBCellFAB> m_coarsenedFineData;
 
  mutable LayoutData<VoFIterator> m_ghostIterFine;
 
  bool m_isDefined;
 
  IntVect m_ghostVectorFine;
 
  IntVect m_ghostVectorCoFi;
 
  int m_refRat;
 
  int m_ghostCF;
 
  int m_order;
 
  int m_weight;
 
  LayoutData<IntVectSet> m_ghostCells;
 
  LayoutData<BaseIVFAB<VoFStencil>> m_fineStencils;
 
  LayoutData<BaseIVFAB<VoFStencil>> m_coarStencils;
 
  LayoutData<RefCountedPtr<AggStencil<EBCellFAB, EBCellFAB>>> m_aggFineStencils;
 
  LayoutData<RefCountedPtr<AggStencil<EBCellFAB, EBCellFAB>>> m_aggCoarStencils;
 
  LayoutData<std::map<std::pair<int, Side::LoHiSide>, Box>> m_cfivs;
 
  LayoutData<CoarseInterpQuadCF> m_loCoarseInterpCF[SpaceDim];
 
  LayoutData<CoarseInterpQuadCF> m_hiCoarseInterpCF[SpaceDim];
 
  virtual void
  defineGhostRegions() noexcept;
 
  virtual void
  defineBuffers() noexcept;
 
  virtual void
  defineCoarseInterp() noexcept;
 
  virtual void
  defineStencilsEBCF() noexcept;
 
  virtual void
  makeAggStencils() noexcept;
 
  virtual bool
  getStencil(VoFStencil&            a_stencilFine,
             VoFStencil&            a_stencilCoar,
             const CellLocation&    a_dataLocation,
             const VolIndex&        a_ghostVofFine,
             const VolIndex&        a_ghostVofCoar,
             const EBISBox&         a_ebisboxFine,
             const EBISBox&         a_ebisboxCoar,
             const DenseIntVectSet& a_validFineCells,
             const DenseIntVectSet& a_validCoarCells,
             const Real&            a_dxFine,
             const Real&            a_dxCoar,
             const int&             a_order,
             const int&             a_weight) const noexcept;
 
  virtual void
  regularCoarseFineInterp(LevelData<EBCellFAB>&       a_finePhi,
                          const LevelData<EBCellFAB>& a_coarPhi,
                          const int                   a_fineVar,
                          const int                   a_coarVar) const noexcept;
};
 
#include <CD_NamespaceFooter.H>
 
#endif