EBQuadCFInterp.H

Go to the documentation of this file.

#ifdef CH_LANG_CC
/*
*      _______              __
*     / ___/ /  ___  __ _  / /  ___
*    / /__/ _ \/ _ \/  V \/ _ \/ _ \
*    \___/_//_/\___/_/_/_/_.__/\___/
*    Please refer to Copyright.txt, in Chombo's root directory.
*/
#endif

// dtgraves fri, sept 14, 2001
// mfbarad made it eb aware in June, 2005

#ifndef _EBQUADCFINTERP_H_
#define _EBQUADCFINTERP_H_

#include "REAL.H"
#include "FArrayBox.H"
#include "LevelData.H"
#include "DisjointBoxLayout.H"
#include "EBIndexSpace.H"
#include "EBISLayout.H"
#include "EBCellFAB.H"
#include "Interval.H"
#include "Stencils.H"
#include "QuadCFInterp.H"
#include "BaseIVFAB.H"
#include "CornerCopier.H"
#include "VoFIterator.H"
#include "NamespaceHeader.H"


class EBQuadCFInterp: public QuadCFInterp
{
public:

  EBQuadCFInterp();

  virtual ~EBQuadCFInterp();


  EBQuadCFInterp(const DisjointBoxLayout& a_dblFine,
                 const DisjointBoxLayout& a_dblCoar,
                 const EBISLayout&        a_ebislFine,
                 const EBISLayout&        a_ebislCoar,
                 const ProblemDomain&     a_domainCoar,
                 const int&               a_nref,
                 const int&               a_nvar,
                 const LayoutData<IntVectSet>&  a_cfivs,
                 const EBIndexSpace* const a_ebisPtr = Chombo_EBIS::instance(),
                 bool a_doEBCFCrossing = true);


  void define(const DisjointBoxLayout& a_dblFine,
              const DisjointBoxLayout& a_dblCoar,
              const EBISLayout&        a_ebislFine,
              const EBISLayout&        a_ebislCoar,
              const ProblemDomain&     a_domainCoar,
              const int&               a_nref,
              const int&               a_nvar,
              const LayoutData<IntVectSet>&  a_cfivs,
              const EBIndexSpace* const a_ebisPtr = Chombo_EBIS::instance(),
              bool a_doEBCFCrossing = true);



  void
  interpolate(LevelData<EBCellFAB>&       a_fineData,
              const LevelData<EBCellFAB>& a_coarData,
              const Interval&             a_variables,
              bool a_doOnlyRegularInterp = false);

  bool isDefined() const;

protected:


  bool getEBCFIVS(const LayoutData<IntVectSet>& a_cfivs);


  void buildEBCFStencils();


  void buildEBCFCornerStencils(const LayoutData<IntVectSet>& a_cfivs);


  void interpEBCFCrossing(LevelData<EBCellFAB>&       a_fineData,
                          const LevelData<EBCellFAB>& a_coarData,
                          const Interval&             a_variables);


  void interpEBCFCorners(LevelData<EBCellFAB>&       a_fineData,
                         const LevelData<EBCellFAB>& a_coarData,
                         const Interval&             a_variables);



  EBISLayout        m_ebislCoar;
  EBISLayout        m_ebislCoarsenedFine;
  EBISLayout        m_ebislFine;
  DisjointBoxLayout m_gridsCoar;
  DisjointBoxLayout m_gridsCoarsenedFine;
  DisjointBoxLayout m_gridsFine;
  ProblemDomain     m_domainFine;
  ProblemDomain     m_domainCoar;
  bool              m_isDefined;
  int               m_nComp;
  int               m_refRat;

  LevelData<EBCellFAB> m_ebBufferCoarsenedFine;

  //cache the EB-CF vofiterator, these are on the fine index space
  LayoutData<VoFIterator> m_vofItEBCFLo[SpaceDim];
  LayoutData<VoFIterator> m_vofItEBCFHi[SpaceDim];
  LayoutData<VoFIterator> m_vofItCorners;
  LayoutData<VoFIterator> m_vofItEdges;

  //cache the EB-CF vofiterator,
  //coarse index space vofs, fine layout
  LayoutData<BaseIVFAB<VoFStencil> > m_coarStencilLo[SpaceDim];
  LayoutData<BaseIVFAB<VoFStencil> > m_coarStencilHi[SpaceDim];

  //These are the stencils that interpolate with phistar to get ghost val
  //phiStar is indexed at the same point as the ghost phi
  //(UNLIKE in regular, non eb quadcfinterp) but lives between
  //coarse cell centers (just like in non-eb quadcfinterp)
  LayoutData<BaseIVFAB<VoFStencil> > m_fineStencilLo[SpaceDim];
  LayoutData<BaseIVFAB<VoFStencil> > m_fineStencilHi[SpaceDim];

  //fine index space fine layout
  LayoutData<IntVectSet> m_ebcfivsLo[SpaceDim];
  LayoutData<IntVectSet> m_ebcfivsHi[SpaceDim];
  LayoutData<BaseIVFAB<VoFStencil> > m_stencilCorners;
  LayoutData<BaseIVFAB<VoFStencil> > m_stencilEdges;

  // Copier object to handle corner cells
  CornerCopier m_cornerCopier;

  bool m_doEBCFCrossing;

private:

  //disallowed for all the usual reasons
  EBQuadCFInterp(const EBQuadCFInterp& ebcin)
  {
    MayDay::Error("ebqcfi 2 invalid operator");
  }
  void operator=(const EBQuadCFInterp& fabin)
  {
    MayDay::Error("ebqcfi 3 invalid operator");
  }
};

#include "NamespaceFooter.H"
#endif

---