Stencils.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 _STENCILS_H_
#define _STENCILS_H_

#include "Vector.H"
#include "VolIndex.H"
#include "FaceIndex.H"
#include "REAL.H"
#include "NamespaceHeader.H"


class VoFStencil
{
public:

  VoFStencil();

  ~VoFStencil();

  void clear();


  void
  add(const VolIndex& vof,const Real& weight, int ivar = 0);


  int
  size() const;


  const VolIndex&
  vof(int i) const;


  const Real&
  weight(int i) const;


  Real& weight(int i);


  const int& variable(int i) const;

  int& variable(int i);


  VoFStencil&
  operator+=(const VoFStencil& a_vofsten);


  void operator*=(const Real& scaling);

  VoFStencil&
  operator=(const VoFStencil& a_vofsten);

  VoFStencil(const VoFStencil& a_vofstenin);

  void setAllVariables(int a_var)
  {
    if(vofs.size() > 0)
      {
        variables = Vector<int>(vofs.size(), a_var);
      }
  }
protected:

  Vector<VolIndex> vofs;
  Vector<Real> weights;

  Vector<int> variables;
};


class FaceStencil
{
public:
  FaceStencil();

  ~FaceStencil();

  void clear();


  void
  add(const FaceIndex& face,const Real& weight, int variable=0);


  int
  size() const;


  const FaceIndex&
  face(int i) const;


  const Real&
  weight(int i) const;

  const int& variable(int i) const;

  int& variable(int i);


  FaceStencil&
  operator+=(const FaceStencil& a_facesten);

  FaceStencil&
  operator=(const FaceStencil& a_facesten);

  FaceStencil(const FaceStencil& a_facesten);


  void operator*=(const Real& scaling);

private:
  Vector<FaceIndex> faces;
  Vector<Real> weights;
  Vector<int> variables;

};

class EBCellFAB;
class EBFaceFAB;
extern Real applyVoFStencil(const VoFStencil& a_vofSten, const EBCellFAB& a_fab, const int& a_comp);

extern Real applyFaceStencil(const FaceStencil& a_faceSten, const EBFaceFAB& a_fab, const int& a_comp);

/**************/
inline int
VoFStencil::size() const
{
  return weights.size();
}
/**************/
/**************/
inline const  VolIndex&
VoFStencil::vof(int i) const
{
  return vofs[i];
}
/**************/
/**************/
inline const Real&
VoFStencil::weight(int i) const
{
  return weights[i];
}

inline Real&
VoFStencil::weight(int i)
{
  return weights[i];
}

#ifdef CH_EXPLICIT_TEMPLATES
//
// Extra no-op stuff required for successful explicit template instantiation.
//
inline int linearSize( const FaceStencil& fs )
{
    CH_assert(0);
    return -1;
}
template<> inline
void linearIn<FaceStencil>(FaceStencil& a_outputT, const void* const inBuf)
{
    CH_assert(0);
}
template<> inline
void linearOut<FaceStencil>(void* const a_outBuf, const FaceStencil& a_inputT)
{
    CH_assert(0);
}


inline int linearSize( const VoFStencil& fs )
{
    CH_assert(0);
    return -1;
}
template<> inline
void linearIn<VoFStencil>(VoFStencil& a_outputT, const void* const inBuf)
{
    CH_assert(0);
}
template<> inline
void linearOut<VoFStencil>(void* const a_outBuf, const VoFStencil& a_inputT)
{
    CH_assert(0);
}


inline int linearSize( const bool& fs )
{
    CH_assert(0);
    return -1;
}
template<> inline
void linearIn<bool>(bool& a_outputT, const void* const inBuf)
{
    CH_assert(0);
}
template<> inline
void linearOut<bool>(void* const a_outBuf, const bool& a_inputT)
{
    CH_assert(0);
}


#endif // CH_EXPLICIT_TEMPLATES


#include "NamespaceFooter.H"
#endif

---