00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029 #ifndef _ChannelBL_H_
00030 #define _ChannelBL_H_
00031
00032 #include "MayDay.H"
00033 #include "LoHiSide.H"
00034 #include "IntVect.H"
00035 #include "RealVect.H"
00036 #include "BaseLevelSet.H"
00037 #include "REAL.H"
00038
00039 using std::pair;
00040
00042
00048 class ChannelBL: public BaseLevelSet
00049 {
00050 public:
00051
00053
00056 ChannelBL(const RealVect& a_centerPlaneNormal,
00057 const Real& a_radius);
00058
00060 virtual ~ChannelBL();
00061
00063
00067 virtual bool isRegular(const Box& region,
00068 const Box& domain,
00069 const RealVect& origin,
00070 const Real& dx) const;
00071
00073
00077 virtual bool isCovered(const Box& region,
00078 const Box& domain,
00079 const RealVect& origin,
00080 const Real& dx) const;
00081
00083
00093 pair<int, Side::LoHiSide> upDirection(const RealVect& a_midpt,const IntVect& a_iv) const;
00094
00096
00101 virtual Real localFuncValue(const RealVect& a_independentCoords,
00102 const int& a_upDirection,
00103 const IntVect& a_iv,
00104 const Box& a_domain,
00105 const RealVect& a_origin,
00106 const Real& a_dx) const;
00107
00108
00110
00114 virtual
00115 BaseLevelSet* new_baseLevelSet() const;
00116 virtual Real LevelSurface(const RealVect& a_vec) const;
00117 Real getCenterPlanePt(const RealVect& a_indCoords) const;
00118
00119 protected:
00120
00121 int m_upDir;
00122
00123 RealVect m_centerPlaneNormal;
00124 Real m_radius;
00125 Real m_upDirDistance;
00126
00127 private:
00129 ChannelBL()
00130 {
00131 MayDay::Error("ChannelBL uses strong construction");
00132 }
00133 ChannelBL(const ChannelBL& a_ramplg)
00134 {
00135 MayDay::Error("ChannelBL disallows copy construction");
00136 }
00137 void operator=(const ChannelBL& a_ramplg)
00138 {
00139 MayDay::Error("ChannelBL disallows assignment");
00140 }
00141
00142 };
00143 #endif