CD_EBMGLeastSquaresInterpolator.H

/* chombo-discharge
 * Copyright © 2021 SINTEF Energy Research.
 * Please refer to Copyright.txt and LICENSE in the chombo-discharge root directory.
 */
 
#ifndef CD_EBMultigridInterpolator_H
#define CD_EBMultigridInterpolator_H
 
// Std includes
#include <map>
 
// Chombo includes
#include <EBLevelGrid.H>
#include <QuadCFInterp.H>
#include <AggStencil.H>
 
// Our includes
#include <CD_Location.H>
#include <CD_NamespaceHeader.H>
 
class EBMultigridInterpolator
{
public:
  using CellLocation = Location::Cell;
 
  EBMultigridInterpolator() = delete;
 
  EBMultigridInterpolator(const EBMultigridInterpolator& a_other) = delete;
 
  EBMultigridInterpolator(const EBLevelGrid& a_eblgFine,
                          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);
 
  ~EBMultigridInterpolator();
 
  int
  getGhostCF() const;
 
  void
  coarseFineInterp(LevelData<EBCellFAB>& a_phiFine, const LevelData<EBCellFAB>& a_phiCoar, const Interval a_variables);
 
  void
  slowCoarseFineInterp(LevelData<EBCellFAB>&       a_phiFine,
                       const LevelData<EBCellFAB>& a_phiCoar,
                       const Interval              a_variables);
 
  void
  coarseFineInterpH(LevelData<EBCellFAB>& a_phiFine, const Interval a_variables) const;
 
  void
  slowCoarseFineInterpH(LevelData<EBCellFAB>& a_phiFine, const Interval a_variables) const;
 
  void
  coarseFineInterpH(EBCellFAB& a_phiFine, const Interval a_variables, const DataIndex& a_dit) const;
 
  void
  slowCoarseFineInterpH(EBCellFAB& a_phiFine, const Interval a_variables, const DataIndex& a_dit) const;
 
protected:
  static constexpr int m_stenComp = 0;
 
  static constexpr int m_numStenComp = 1;
 
  static constexpr int m_comp = 0;
 
  RefCountedPtr<QuadCFInterp> m_scalarInterpolator;
 
  RefCountedPtr<QuadCFInterp> m_vectorInterpolator;
 
  CellLocation m_dataLocation;
 
  EBLevelGrid m_eblgFine;
 
  EBLevelGrid m_eblgCoar;
 
  EBLevelGrid m_eblgCoFi;
 
  BoxLayout m_grownCoarBoxesLayout;
 
  mutable BoxLayoutData<EBCellFAB> m_grownCoarData;
 
  mutable LayoutData<VoFIterator> m_vofIterFine;
 
  IntVect m_ghostVector;
 
  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;
 
  void
  defineGhostRegions();
 
  void
  defineGrids(const EBLevelGrid& a_eblgFine, const EBLevelGrid& a_eblgCoar);
 
  void
  defineBuffers();
 
  void
  defineStencilsEBCF();
 
  void
  makeAggStencils();
 
  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);
};
 
#include <CD_NamespaceFooter.H>
 
#endif