ComplementIF.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 _COMPLEMENTIF_H_
#define _COMPLEMENTIF_H_

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

#include "BaseIF.H"

#include "NamespaceHeader.H"


class ComplementIF: public BaseIF
{
public:

  ComplementIF(const BaseIF& a_impFunc,
               const bool&   a_complement = true);

  ComplementIF(const ComplementIF& a_inputIF,
               const bool&         a_complement = true);

  virtual ~ComplementIF();


  virtual void GetParams(bool& a_complement) const;


  virtual void SetParams(const bool& a_complement);


  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 IndexTM<Real,GLOBALDIM> normal(const IndexTM<Real,GLOBALDIM>& a_point) const;


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


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

  virtual BaseIF* newImplicitFunction() const;

  virtual bool fastIntersection(const RealVect& a_low, const RealVect& a_high) const
  {return m_impFunc->fastIntersection(a_low, a_high);}

  virtual GeometryService::InOut InsideOutside(const RealVect& a_low, const RealVect& a_high) const;
protected:
  bool m_complement;               // complement the function (default = true)

  BaseIF* m_impFunc; // implicit function to complement

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

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

#include "NamespaceFooter.H"
#endif

---