Chombo + EB: MonomialIF.H Source File

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 _MONOMIALIF_H_
00012 #define _MONOMIALIF_H_
00013 
00014 #include "MayDay.H"
00015 #include "RealVect.H"
00016 
00017 #include "BaseIF.H"
00018 #include "OneDIntegrator.H"
00019 
00020 #include "NamespaceHeader.H"
00021 
00022 ///
00023 /**
00024    class to implement 1-D function that gets integrated 
00025    to compute EBMoments for monomials
00026 */
00027 class MonomialEBMomentFunc : public BaseOneDFunc
00028 {
00029 public:
00030   ///
00031   /**
00032        Q = power of monomial
00033        p = moment being calculated
00034    */
00035   MonomialEBMomentFunc(const IvSDMinOne& a_Q,
00036                        const IvSpaceDim& a_p)
00037   {
00038     m_Q = a_Q;
00039     m_p = a_p;
00040   }
00041 
00042 
00043   virtual ~MonomialEBMomentFunc()
00044   {; }
00045 
00046   virtual Real value(const Real& location) const;
00047 
00048 private:
00049   //power of monomial
00050   IvSDMinOne m_Q;
00051   //moment being calculated
00052   IvSpaceDim m_p;
00053   //strong construction is best when it can be used
00054   MonomialEBMomentFunc()
00055   { ;  }
00056   
00057 };
00058 ///
00059 /**
00060    This implicit function specifies a Monomial.
00061 */
00062 class MonomialIF: public BaseIF
00063 {
00064 public:
00065   ///
00066   /**
00067      in 2d
00068      y = (x)^p[0] 
00069      in 3d
00070      z = [(x)^p[0]*(y)^p[1]]
00071 
00072      p == a_power
00073   */
00074   MonomialIF(const IndexTM<int, SpaceDim-1> & a_power,
00075              const bool                     & a_inside);
00076 
00077   /// Destructor
00078   virtual ~MonomialIF()
00079   {
00080   }
00081 
00082   virtual Real value(const IndexTM<int,SpaceDim> & a_partialDerivative,
00083                      const IndexTM<Real,SpaceDim>& a_point) const;
00084 
00085   virtual Real value(const RealVect& a_point) const;
00086 
00087   virtual Real value(const IndexTM<Real,SpaceDim>& a_point) const;
00088 
00089   virtual BaseIF* newImplicitFunction() const;
00090 
00091   ///return int x^p dV for the vof
00092   virtual IndMomSpaceDim getExactVolumeMoments(const VolIndex    & a_vof,
00093                                                const Real        & a_dx) const;
00094                                           
00095   
00096   ///return int_eb x^p dA for the eb
00097   virtual IndMomSpaceDim getExactEBMoments(const VolIndex   & a_vof,
00098                                            const Real       & a_dx) const;
00099                                           
00100                                           
00101   ///return int_eb x^p n_i dA for the eb
00102   virtual IndMomSpaceDim getExactEBNormalMoments(const VolIndex   & a_vof,
00103                                                  const Real       & a_dx,
00104                                                  const int        & a_ni) const;
00105                                           
00106   ///
00107   virtual IndMomSpaceDim getExactEBNormalPartialDerivs(const VolIndex   & a_vof,
00108                                                        const Real       & a_dx,
00109                                                        const int        & a_ni) const;
00110 
00111   ///
00112   Real
00113   getZerothDerivativeOfXNormal(const Real& a_xloc) const;
00114   ///
00115   Real
00116   getFirstDerivativeOfXNormal(const Real& a_xloc) const;
00117   ///
00118   Real
00119   getSecondDerivativeOfXNormal(const Real& a_xloc) const;
00120   ///
00121   Real
00122   getThirdDerivativeOfXNormal(const Real& a_xloc) const;
00123   ///
00124   Real
00125   getFourthDerivativeOfXNormal(const Real& a_xloc) const;
00126   ///
00127   Real
00128   getZerothDerivativeOfYNormal(const Real& a_xloc) const; 
00129   ///
00130   Real
00131   getFirstDerivativeOfYNormal(const Real& a_xloc) const;
00132   ///
00133   Real
00134   getSecondDerivativeOfYNormal(const Real& a_xloc) const;
00135   ///
00136   Real
00137   getThirdDerivativeOfYNormal(const Real& a_xloc) const;
00138   ///
00139   Real
00140   getFourthDerivativeOfYNormal(const Real& a_xloc) const;
00141 
00142 
00143   ///return int x^p  dA for the face
00144   virtual IndMomSDMinOne getExactFaceMoments(const FaceIndex  & a_face,
00145                                              const Real       & a_dx) const;
00146                                           
00147   void getTloThi(bool& a_noEB,
00148                  Real& a_tlo, 
00149                  Real& a_thi,
00150                  const VolIndex   & a_vof,
00151                  const Real       & a_dx) const;
00152 
00153   Real getEBNormMomX(const Real        & a_tlo,
00154                      const Real        & a_thi,
00155                      const IvSpaceDim  & a_p) const;
00156 
00157   Real getEBNormMomY(const Real        & a_tlo,
00158                      const Real        & a_thi,
00159                      const IvSpaceDim  & a_p) const;
00160 
00161   Real xToTheQ(const IndexTM<Real,SpaceDim>& a_point) const;
00162 
00163   IndMomSpaceDim 
00164   getUncutVolumeSection(const RealVect& a_xlo, 
00165                         const RealVect& a_xhi, 
00166                         const RealVect& a_xbar,
00167                           const Real    & a_dx) const;
00168 
00169   IndMomSpaceDim 
00170   getCutVolumeSection(const RealVect & a_xlo, 
00171                       const RealVect & a_xhi, 
00172                       const RealVect & a_xbar,
00173                       const Real     & a_dx) const;
00174 
00175   void shiftToXBar(IndMomSpaceDim& a_moment,
00176                    const RealVect& a_xbar) const;
00177 protected:
00178 
00179   IndexTM<int, SpaceDim-1>  m_power;
00180   bool                      m_inside;
00181 
00182 private:
00183   MonomialIF();
00184 
00185 };
00186 
00187 #include "NamespaceFooter.H"
00188 #endif