GaussianSphereIF.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 _GAUSSIANSPHEREIF_H_
#define _GAUSSIANSPHEREIF_H_

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

#include "BaseIF.H"

#include "NamespaceHeader.H"


class GaussianSphereIF: public BaseIF
{
public:

  GaussianSphereIF(const Real&     a_sigma,
                   const RealVect& a_center,
                   const bool&     a_inside);

  GaussianSphereIF(const GaussianSphereIF& a_inputIF);

  virtual ~GaussianSphereIF();


  virtual void GetParams(Real&     a_sigma,
                         RealVect& a_center,
                         bool&     a_inside) const;


  virtual void SetParams(const Real&     a_sigma,
                         const RealVect& a_center,
                         const bool&     a_inside);


  virtual Real value(const RealVect& a_point) const;


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


  virtual IndexTM<Real,GLOBALDIM> grad(const IndexTM<Real,GLOBALDIM>& a_point) const;


  virtual Vector<IndexTM <Real,GLOBALDIM> > gradGrad(const IndexTM<Real,GLOBALDIM>& a_point) const;

  virtual BaseIF* newImplicitFunction() const;

protected:
  Real     m_sigma;     // sigma
  RealVect m_center;    // center
  bool     m_inside;    // inside flag

  Real     m_sigma2;    // precomputed sigma squared
  Real     m_normal;    // precomputed normalizing factor

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

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

#include "NamespaceFooter.H"
#endif

---