---

EBQuadCFInterp.H

Go to the documentation of this file.

/* _______              __
  / ___/ /  ___  __ _  / /  ___
 / /__/ _ \/ _ \/  ' \/ _ \/ _ \
 \___/_//_/\___/_/_/_/_.__/\___/ 
*/
//
// This software is copyright (C) by the Lawrence Berkeley
// National Laboratory.  Permission is granted to reproduce
// this software for non-commercial purposes provided that
// this notice is left intact.
// 
// It is acknowledged that the U.S. Government has rights to
// this software under Contract DE-AC03-765F00098 between
// the U.S.  Department of Energy and the University of
// California.
//
// This software is provided as a professional and academic
// contribution for joint exchange. Thus it is experimental,
// is provided ``as is'', with no warranties of any kind
// whatsoever, no support, no promise of updates, or printed
// documentation. By using this software, you acknowledge
// that the Lawrence Berkeley National Laboratory and
// Regents of the University of California shall have no
// liability with respect to the infringement of other
// copyrights by any part of this software.
//
// dtgraves fri, sept 14, 2001

#ifndef _EB_QUAD_CF_INTERP_H_
#define _EB_QUAD_CF_INTERP_H_

#include "REAL.H"
#include "FArrayBox.H"
#include "LevelData.H"
#include "DisjointBoxLayout.H"
#include "EBISLayout.H"
#include "EBCellFAB.H"
#include "Interval.H"
#include "Stencils.H"
#include "QuadCFStencil.H"
#include "BaseIVFAB.H"


class EBQuadCFInterp
{
public:

  EBQuadCFInterp();

  ~EBQuadCFInterp();


  EBQuadCFInterp(const DisjointBoxLayout& a_dblFine,
                 const DisjointBoxLayout& a_dblCoar,
                 const EBISLayout& a_ebislFine,
                 const EBISLayout& a_ebislCoar,
                 const Box& a_domainCoar,
                 const int& a_nref,
                 const int& a_nvar);


  void define(const DisjointBoxLayout& a_dblFine,
              const DisjointBoxLayout& a_dblCoar,
              const EBISLayout& a_ebislFine,
              const EBISLayout& a_ebislCoar,
              const Box& a_domainCoar,
              const int& a_nref,
              const int& a_nvar);


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

  bool isDefined() const;
protected:

  /* Coarse / Fine interpolation operator.
     using one sided differencing
     only at a particular face
  */
  void regCFInterp(BaseFab<Real> & a_phif,
                   const BaseFab<Real> & a_phic,
                   const QuadCFStencil& a_qcfs,
                   const Side::LoHiSide a_hiorlo,
                   const int& a_idir,
                   const Interval& a_variables) const;

  void irrCFInterp(EBCellFAB& a_phif,
                   const EBCellFAB& a_phic,
                   const DataIndex& a_datInd,
                   const Interval& a_variables) const;

  //get extended phi (lives next to interpivs)
  void getPhiStar(BaseFab<Real> & a_phistar, 
                  const BaseFab<Real> & a_phic,
                  const QuadCFStencil& a_qcfs,
                  const Side::LoHiSide a_hiorlo,
                  const int& a_idir,
                  const Interval& a_variables) const;

  //interpolate over correct intvectset
  void interpOnIVS(BaseFab<Real> & a_phif,
                   const BaseFab<Real> & a_phiStar, 
                   const QuadCFStencil& a_qcfs,
                   const Side::LoHiSide a_hiorlo,
                   const int& a_idir,
                   const Interval& a_variables) const;

  
  bool extractFromSten(Real& a_value,
                       const VoFStencil& a_vofsten,
                       const EBCellFAB& a_coarData,
                       const int& ivar) const;

  void  makeStencils();

  //number of components
  int m_nComp;

  //refinement ratio between levels
  int m_refRat;

  //cfstencils on low-side faces
  //lives on fine input grid layout
  LayoutData<QuadCFStencil> m_loQCFS[SpaceDim];

  //cfstencils on high-side faces
  //lives on fine input grid layout
  LayoutData<QuadCFStencil> m_hiQCFS[SpaceDim];

  //buffer for coarse data 
  //lives on fine proc layout (coarsened fine grids with ghost)
  LevelData<EBCellFAB> m_coarBuffer;

  // big ebisl commentary 
  //(not a big commentary. a commentary about a big ebisl):
  //this is the ebisl for the buffers.  This could be done
  //with a 2*SpaceDim smaller layouts but that would be
  //more complicated and (i suspect) more expensive
  //in terms of communication ---not obvious since
  //this makes a bigger ebisl 
  EBISLayout        m_coarsenedFineEBISL;
  DisjointBoxLayout m_coarsenedFineGrids;

  //stored inputs
  DisjointBoxLayout m_inputFineGrids;
  DisjointBoxLayout m_inputCoarGrids;
  EBISLayout        m_inputFineEBISL;
  EBISLayout        m_inputCoarEBISL;
  Box m_domainCoar;
  Box m_domainFine;

  bool m_isDefined;

  //these are the stencils for the derivatives
  //only on irregular cells inside the coarse-fine interface.
  //lives in the coarse index space
  LayoutData<BaseIVFAB<VoFStencil> >  m_stencils[SpaceDim];

  //These are the cells on the coarse/fine interface
  //that are irregular.
  //These are in the fine index space.  
  LayoutData<IntVectSet>  m_irregRegions;

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");
  }
};

#endif

---