TiltedCylinderIF.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 _TILTEDCYLINDERIF_H_
#define _TILTEDCYLINDERIF_H_

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

#include "BaseIF.H"

#include "NamespaceHeader.H"


class TiltedCylinderIF: public BaseIF
{
public:

  TiltedCylinderIF(const Real&     a_radius,
                   const RealVect& a_direction,
                   const RealVect& a_point,
                   const bool&     a_inside);

  TiltedCylinderIF(const TiltedCylinderIF& a_inputIF);

  virtual ~TiltedCylinderIF();


  virtual void GetParams(Real&     a_radius,
                         RealVect& a_direction,
                         RealVect& a_point,
                         bool&     a_inside) const;


  virtual void SetParams(const Real&     a_radius,
                         const RealVect& a_direction,
                         const RealVect& a_point,
                         const bool&     a_inside);


  virtual Real value(const RealVect& a_point) const;

  virtual BaseIF* newImplicitFunction() const;

  virtual bool fastIntersection(const RealVect& a_low, const RealVect& a_high) const
  {return !m_tiltedAxis;}


  virtual GeometryService::InOut InsideOutside(const RealVect& a_low, const RealVect& a_high) const ;

protected:
  Real     m_radius;      // radius
  RealVect m_direction;   // axis direction
  bool     m_tiltedAxis;  // false if straight cylinder
  int      m_coordDir;    // coordinate direction of axis for a straight cylinder
  RealVect m_point;       // point on the axis
  bool     m_inside;      // inside flag

  Real     m_radius2;     // precomputed radius squared
  Real     m_length2;     // precomputed length squared of direction vector

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

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

#include "NamespaceFooter.H"
#endif

---