GaussianIF.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 _GAUSSIANIF_H_
#define _GAUSSIANIF_H_

#include "MayDay.H"
#include "RealVect.H"

#include "BaseIF.H"

#include "NamespaceHeader.H"


class GaussianIF: public BaseIF
{
public:

  GaussianIF(const RealVect& a_origin,
             const Real&     a_height,
             const RealVect& a_sig2,
             const int&      a_upDir,
             const bool&     a_inside);

  GaussianIF(const GaussianIF& a_inputIF);

  virtual ~GaussianIF();


  virtual void GetParams(RealVect& a_origin,
                         Real&     a_height,
                         RealVect& a_sig2,
                         int&      a_upDir,
                         bool&     a_inside) const;


  virtual void SetParams(const RealVect& a_origin,
                         const Real&     a_height,
                         const RealVect& a_sig2,
                         const int&      a_upDir,
                         const bool&     a_inside);


  virtual Real value(const RealVect& a_point) const;

  virtual BaseIF* newImplicitFunction() const;

protected:
  RealVect m_origin;    // gaussian origin
  Real     m_height;    // height of the gaussian
  RealVect m_sig2;      // standard deviation of gaussian
  int      m_upDir;     // the up direction for the gaussian (the symmetry direction)
  bool     m_inside;    // inside flag

private:
  GaussianIF()
  {
    MayDay::Error("GaussianIF uses strong construction");
  }

  void operator=(const GaussianIF& a_inputIF)
  {
    MayDay::Error("GaussianIF doesn't allow assignment");
  }
};

#include "NamespaceFooter.H"
#endif

---