QuadCFStencil.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 _QUADCFSTENCIL_H_
#define _QUADCFSTENCIL_H_

#include <iostream>
#include <math.h>
#include "SPACE.H"
#include <stdlib.h>
#include "REAL.H"
#include "IntVect.H"
#include "Box.H"
#include "DisjointBoxLayout.H"
#include "FArrayBox.H"
#include "IntVectSet.H"
#include "CFStencil.H"
#include "DerivStencil.H"
#include "NamespaceHeader.H"


class QuadCFStencil
{

public:


  Real computeSecondDerivative(
                               const BaseFab<Real> & a_phic,
                               int a_derivdir,
                               int a_ivar,
                               const IntVect& a_ivin,
                               Real a_dx) const;


  Real computeFirstDerivative(
                              const BaseFab<Real> & a_phic,
                              int a_derivdir,
                              int a_ivar,
                              const IntVect& a_ivin,
                              Real a_dx) const;


  Real computeMixedDerivative(
                              const BaseFab<Real> & a_phic,
                              int a_ivar,
                              const IntVect& a_ivin,
                              Real a_dx) const;

  QuadCFStencil();

  ~QuadCFStencil();

  QuadCFStencil(
                const Box& a_fine_domain,
                const Box& a_grid,
                const DisjointBoxLayout& a_fineBoxes,
                const DisjointBoxLayout& a_coarBoxes,
                int a_refRatio,
                int a_direction,
                Side::LoHiSide a_hiorlo);

  QuadCFStencil(
                const ProblemDomain& a_fine_domain,
                const Box& a_grid,
                const DisjointBoxLayout& a_fineBoxes,
                const DisjointBoxLayout& a_coarBoxes,
                int a_refRatio,
                int a_direction,
                Side::LoHiSide a_hiorlo);

  void   define(
                const Box& a_fine_domain,
                const Box& a_grid,
                const DisjointBoxLayout& a_fineBoxes,
                const DisjointBoxLayout& a_coarBoxes,
                int a_refRatio,
                int a_direction,
                Side::LoHiSide a_hiorlo);

  void   define(
                const ProblemDomain& a_fine_domain,
                const Box& a_grid,
                const DisjointBoxLayout& a_fineBoxes,
                const DisjointBoxLayout& a_coarBoxes,
                int a_refRatio,
                int a_direction,
                Side::LoHiSide a_hiorlo);


  void   define(
                const ProblemDomain& a_fine_domain,
                const Box& a_grid,
                const Vector<Box>& a_periodicBoxes,
                const Vector<Box>& a_coarsenedPeriodicBoxes,
                const DisjointBoxLayout& a_coarBoxes,
                int a_refRatio,
                int a_direction,
                Side::LoHiSide a_hiorlo);



  bool isDefined() const;


  bool isEmpty() const
  {
    CH_assert(m_baseCFS.isDefined());
    return m_baseCFS.isEmpty();
  }


  const IntVectSet& getFineIVS() const
  {
    CH_assert(m_baseCFS.isDefined());
    return m_baseCFS.getFineIVS();
  }


  const IntVectSet& getCoarIVS() const
  {
    CH_assert(m_baseCFS.isDefined());
    return m_baseCFS.getCoarIVS();
  }

  bool finePacked() const
  {
    return m_baseCFS.isPacked();
  }

  const Box& packedBox() const
  {
    return m_baseCFS.packedBox();
  }


protected:

  //has define fcn been called?
  bool m_isDefined;

  //direction of face
  int m_direction;


  IntVectSet m_ivsStandard;


  IntVectSet m_ivsQuadd;



  BaseFab<DerivStencil> m_mixedSten;


  BaseFab<DerivStencil> m_second[SpaceDim];


  BaseFab<DerivStencil> m_firstD[SpaceDim];


  BaseFab<bool> m_dropOrd;

  CFStencil m_baseCFS;


private:
  // internally useful but not for public consumption

  //
  void  setDefaultValues();

  /*
    add fab to the stencil.  For internal use only.
  */
  void addFabToSten(const BaseFab<Real> & fabin_a,
                    DerivStencil&  sten_a);

  // there is no assignment operator for this class
  void operator= (const QuadCFStencil& cfs_in)
  {
  }

  // there is no copy constructor for this class
  QuadCFStencil(const QuadCFStencil& cfs_in)
  {
  }

  void buildStencils(const IntVectSet& ivsAllGood);

};

#include "NamespaceFooter.H"
#endif

---