EBCoarsen.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 _EBCOARSEN_H_
#define _EBCOARSEN_H_

#include "REAL.H"
#include "FArrayBox.H"
#include "LevelData.H"
#include "EBLevelGrid.H"
#include "EBCellFAB.H"
#include "Interval.H"
#include "VoFIterator.H"
#include "NamespaceHeader.H"
class EBIndexSpace;
class VoFStencil;


class EBCoarsen
{
public:

  EBCoarsen();

  ~EBCoarsen();


  EBCoarsen(const EBLevelGrid& eblgFine,
            const EBLevelGrid& eblgCoar,
            const int& nref,
            const int& nvar);


  void define(const EBLevelGrid& eblgFine,
              const EBLevelGrid& eblgCoar,
              const int& nref,
              const int& nvar);


  bool isDefined() const;


  void
  coarsenFine(LevelData<EBCellFAB>&       a_coarData,
              const LevelData<EBCellFAB>& a_fineData,
              const Interval&             a_variables);
protected:
  void defineStencil(const LayoutData<IntVectSet>&  a_cfivs);

  void coarsenIrreg(EBCellFAB&       a_coar,
                    const EBCellFAB& a_fine,
                    const DataIndex& a_dit,
                    const Interval&  a_variables);

  void checkStencil(VoFStencil& a_stencil,
                    const Real& a_sum) const;

  bool getVofsSideDir(Vector<VolIndex>&     a_vofList,
                      const VolIndex&       a_vof,
                      const IntVectSet&     a_cfivs,
                      const EBISBox&        a_ebisBox,
                      const int&            a_dir,
                      const Side::LoHiSide& a_side,
                      const int&            a_stepSize) const;

  void addWeight(VoFStencil&             a_stencil,
                 const Vector<VolIndex>& a_vofList,
                 const Real&             a_weight,
                 const bool&             a_hasVof) const;

  void getCoarsenVoFStencil(VoFStencil&       a_stencil,
                            const EBISBox&    a_ebisBoxCoar,
                            const EBISBox&    a_ebisBoxFine,
                            const VolIndex&   a_vofCoar,
                            const DataIndex&  a_datInd,
                            const IntVectSet& a_cfivs);

  void coarsenFAB(EBCellFAB&       a_coar,
                  const EBCellFAB& a_fine,
                  const DataIndex& a_datInd,
                  const Interval&  a_variables);

  void
  setDefaultValues();

  bool m_isDefined;

  DisjointBoxLayout m_gridsCoar;
  DisjointBoxLayout m_gridsFine;
  DisjointBoxLayout m_coarsenedFineGrids;
  ProblemDomain     m_domainCoar;
  ProblemDomain     m_domainFine;

  EBISLayout m_ebislCoar;
  EBISLayout m_ebislFine;

  EBISLayout m_coarsenedFineEBISL;
  RefCountedPtr<LayoutData<IntVectSet> > m_cfivsPtr;

  int m_refRat;
  int m_nComp;
  LevelData<EBCellFAB> m_coarsenedFineData;

  LayoutData<BaseIVFAB<VoFStencil> > m_coarsenStencil;
  LayoutData<VoFIterator>           m_vofIt;

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

};

#include "NamespaceFooter.H"
#endif

---