EBStencil.H

Go to the documentation of this file.

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

#ifndef _EBSTENCIL_H_
#define _EBSTENCIL_H_

#include "Stencils.H"
#include "Vector.H"
#include "VolIndex.H"
#include "REAL.H"
#include "Box.H"
#include "FArrayBox.H"
#include "EBIndexSpace.H"
#include "EBCellFAB.H"
#include "NamespaceHeader.H"


class EBStencil
{
public:

  ~EBStencil();


  EBStencil(const Vector<VolIndex>&      a_srcVofs,
            const BaseIVFAB<VoFStencil>& a_vofstencil,
            const Box&                   a_box,
            const EBISBox&               a_ebisbox,
            const IntVect&               a_ghostVectPhi,
            const IntVect&               a_ghostVectLph,
            int a_varDest = 0,
            bool a_doRelaxOpt = false);



  EBStencil(const Vector<VolIndex>& a_srcVofs,
            const Vector<VoFStencil>& a_vofStencil,
            const Box& a_boxLph,
            const Box& a_boxPhi,
            const EBISBox& a_ebisBoxLph,
            const EBISBox& a_ebisBoxPhi,
            const IntVect& a_ghostVectLph,
            const IntVect& a_ghostVectPhi,
            int a_varDest = 0);



  void
  apply(EBCellFAB& a_lofphi, const EBCellFAB& a_phi, bool incrementOnly = false) const;



  void
  apply(EBCellFAB& a_lofphi, const EBCellFAB& a_phi,
        const BaseIVFAB<Real>& a_alphaWeight,
        Real a_alpha, Real a_beta, bool incrementOnly = false) const;

  void
  relax(EBCellFAB& a_phi,
        const EBCellFAB& a_rhs,
        const BaseIVFAB<Real>& a_alphaWeight,
        const BaseIVFAB<Real>&  a_betaWeight,
        Real a_alpha, Real a_beta, Real a_safety) const;


  void
  cache(const EBCellFAB& a_lphi) const;


  void
  cachePhi(const EBCellFAB& a_lphi) const;



  void
  uncache(EBCellFAB& a_lphi) const;


  void
  uncachePhi(EBCellFAB& a_phi) const;



  void
  computeOffsets(const Vector<VolIndex>& a_srcVoFs, const BaseIVFAB<VoFStencil>& a_vofstencil);


  struct {
    int offset;
    bool multiValued;
  } typedef destTerm_t;

  struct {
    int offset;
    Real weight;
    //only used when debugging
    //VolIndex vof;
  } typedef stencilTerm;

  struct {
    Vector<stencilTerm> single;
    Vector<stencilTerm> multi;
  } typedef ebstencil_t;

protected:

  Box m_box;
  EBISBox m_ebisBox;
  Box m_lphBox;
  Box m_phiBox;
  IntVect m_ghostVectPhi;
  IntVect m_ghostVectLph;
  int m_destVar;
  Vector<ebstencil_t> m_ebstencil;
  Vector<destTerm_t>  m_destTerms;
  Vector<destTerm_t>  m_sourTerms;
  Vector<int>         m_alphaBeta;
  mutable Vector<Real> m_cache;
  mutable Vector<Real> m_cachePhi;
  //only used when debugging
  //Vector<VolIndex> m_srcVoFs;
  bool m_doRelaxOpt;
private:
  /*
     Default constructor.
  */
  EBStencil()
  {
    MayDay::Error("invalid operator");
  }

  /*
     Operator =.
  */
  void
  operator=(const EBStencil& stenin)
  {
    MayDay::Error("invalid operator");
  }

  /*
     Copy constructor.
  */
  EBStencil(const EBStencil& stenin)
  {
    MayDay::Error("invalid operator");
  }

};


#include "NamespaceFooter.H"
#endif

---