ReferenceHeightIF.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 _REFERENCEHEIGHTIF_H_
#define _REFERENCEHEIGHTIF_H_

#include "Notation.H"
#include "RealVect.H"
#include "IndexTM.H"
#include "BaseIF.H"
#include "NamespaceHeader.H"


class ReferenceHeightIF : public BaseIF
{
public:

  ReferenceHeightIF(const  BaseIF & a_implicitFunction,
                    const  Real   & a_referenceHeight):
    m_implicitFunction(a_implicitFunction.newImplicitFunction()),
    m_referenceHeight(a_referenceHeight){;}

  virtual ~ReferenceHeightIF()
  {
    delete m_implicitFunction;
  }


  virtual Real value(const RealVect& a_point) const
  {
    IndexTM<Real,GLOBALDIM> pt;

    for (int idir = 0; idir < SpaceDim; ++idir)
      {
        pt[idir] = a_point[idir];
      }

    if (GLOBALDIM == 3 && SpaceDim == 2)
    {
      pt[SpaceDim] = m_referenceHeight;
    }

    return value(pt);
  }



  virtual Real value(const IndexTM<Real,GLOBALDIM>& a_point) const
  {
    return m_implicitFunction->value(a_point);
  };


  virtual IndexTM<Real,GLOBALDIM> grad(const IndexTM<Real,GLOBALDIM>& a_point) const
  {
    return m_implicitFunction->grad(a_point);
  };


  virtual IndexTM<Real,GLOBALDIM> normal(const IndexTM<Real,GLOBALDIM>& a_point) const
  {
    return m_implicitFunction->normal(a_point);
  };


  virtual Vector<IndexTM <Real,GLOBALDIM> > gradGrad(const IndexTM<Real,GLOBALDIM>& a_point) const
  {
    return m_implicitFunction->gradGrad(a_point);
  };


  virtual Vector<IndexTM <Real,GLOBALDIM> > gradNormal(const IndexTM<Real,GLOBALDIM>& a_point) const
  {
    return m_implicitFunction->gradNormal(a_point);
  };


  virtual ReferenceHeightIF* newImplicitFunction() const
  {
    return new ReferenceHeightIF(*m_implicitFunction, m_referenceHeight);
  }


  Real getReferenceHeight() const
  {
    return m_referenceHeight;
  }

private:

  // Default constructor is disabled in order to prevent initialization
  // without reference height. Base classes should all use the constructor
  // version that takes a reference height
  ReferenceHeightIF()
  {
  }

  BaseIF* m_implicitFunction;
  Real    m_referenceHeight;
};

#include "NamespaceFooter.H"
#endif

---