00001 #ifdef CH_LANG_CC 00002 /* 00003 * _______ __ 00004 * / ___/ / ___ __ _ / / ___ 00005 * / /__/ _ \/ _ \/ V \/ _ \/ _ \ 00006 * \___/_//_/\___/_/_/_/_.__/\___/ 00007 * Please refer to Copyright.txt, in Chombo's root directory. 00008 */ 00009 #endif 00010 00011 #ifndef _TORUSIF_H_ 00012 #define _TORUSIF_H_ 00013 00014 #include "MayDay.H" 00015 #include "RealVect.H" 00016 00017 #include "BaseIF.H" 00018 00019 #include "NamespaceHeader.H" 00020 00022 00026 class TorusIF: public BaseIF 00027 { 00028 public: 00030 00035 TorusIF(const Real& a_majorRadius, 00036 const Real& a_minorRadius, 00037 const RealVect& a_center, 00038 const bool& a_inside); 00039 00041 TorusIF(const TorusIF& a_inputIF); 00042 00044 virtual ~TorusIF(); 00045 00047 00050 virtual void GetParams(Real& a_majorRadius, 00051 Real& a_minorRadius, 00052 RealVect& a_center, 00053 bool& a_inside) const; 00054 00056 00059 virtual void SetParams(const Real& a_majorRadius, 00060 const Real& a_minorRadius, 00061 const RealVect& a_center, 00062 const bool& a_inside); 00063 00065 00068 virtual Real value(const RealVect& a_point) const; 00069 00070 virtual BaseIF* newImplicitFunction() const; 00071 00072 protected: 00073 Real m_majorRadius; // major radius 00074 Real m_minorRadius; // minor radius 00075 00076 RealVect m_center; // center 00077 00078 bool m_inside; // inside flag 00079 00080 Real m_minorRadius2; // precomputed minor radius squared 00081 00082 private: 00083 TorusIF() 00084 { 00085 MayDay::Error("TorusIF uses strong construction"); 00086 } 00087 00088 void operator=(const TorusIF& a_inputIF) 00089 { 00090 MayDay::Error("TorusIF doesn't allow assignment"); 00091 } 00092 }; 00093 00094 #include "NamespaceFooter.H" 00095 #endif