---

QuadCFStencil.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.
//

#ifndef _QuadCFStencil_H_
#define _QuadCFStencil_H_

#include <iostream>
#include <math.h>
#include <assert.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"




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



  bool isDefined() const;



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


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



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


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

};




#endif

---