---

Moments.H

Go to the documentation of this file.

/* _______              __
  / ___/ /  ___  __ _  / /  ___
 / /__/ _ \/ _ \/  ' \/ _ \/ _ \
 \___/_//_/\___/_/_/_/_.__/\___/ 
*/
//
// This software is copyright (C) by the Lawrence Berkeley
// National Laboratory.  Permission is granted to reproduce
// this software for non-commercial purposes provided that
// this notice is left intact.
// 
// It is acknowledged that the U.S. Government has rights to
// this software under Contract DE-AC03-765F00098 between
// the U.S.  Department of Energy and the University of
// California.
//
// This software is provided as a professional and academic
// contribution for joint exchange. Thus it is experimental,
// is provided ``as is'', with no warranties of any kind
// whatsoever, no support, no promise of updates, or printed
// documentation. By using this software, you acknowledge
// that the Lawrence Berkeley National Laboratory and
// Regents of the University of California shall have no
// liability with respect to the infringement of other
// copyrights by any part of this software.
#ifndef _Moments_H_
#define _Moments_H_

class edgeMo
{
private: 
  RealVect m_Lo;
  RealVect m_Hi;
  Real m_edgeLength;
  RealVect m_edgeCentroid;
  bool m_covered;
  int m_varOfInt;

public:
  void define (const RealVect& a_loPt,
               const RealVect& a_hiPt, 
               const int& a_direction, 
               const bool& a_covered);
  
  edgeMo(){};
  
  //access member data
  RealVect getLo();
  RealVect getHi();
  Real getEdgeLength();
  RealVect getEdgeCentroid();
  bool isCovered();
  int direction();
  
  //integrate monomial
  Real moment(const IntVect& monomial);
  
};
class faceMo
 {
 private: 
   edgeMo m_edges[4]; //lexigraphic order:0Lo,0Hi,1Lo,1Hi,2Ho,2Hi
   Real m_normalVec[2];
   bool m_regular;
   bool m_covered;
   bool m_dontKnow;
   RealVect m_centroid;
   Real m_areaFrac;
   int m_faceNormal;
   Real m_bdLength;
   
   public:
   void setFaceArea(const Real& area);
   Real getFaceArea();
   void setFaceCentroid(const RealVect& centroid);
   RealVect getFaceCentroid();
   bool m_normalConstructed;
   bool isRegular()const;
   bool isCovered()const;
   bool dontKnow() const;
   void define(edgeMo a_edges[4],
               int& faceNormal, 
               bool& covered,
               bool& regular, 
               bool& dontKnow);
   faceMo(){};
   void getNormal(Real normal[2])const;
   Real getBdLength()const;
   void makeNormal();
   void getEdges(edgeMo Edges[4])const;
   bool normalConstructed()const;
   int getFaceNormal()const;
   void normalize(Real normalVec[2]);
 
 private:
    // disallowed for all the usual reasons
  void operator=(const faceMo& a_input)
  {
    MayDay::Error("invalid operator");
  };

  // disallowed for all the usual reasons
  faceMo(const faceMo& a_input)
  {
    MayDay::Error("invalid operator");
  };
};
 
class vofMo
{
private:
  faceMo m_faces[6];
  Real m_normalVec[3];
  Real m_bdArea;

public:
  void define(faceMo a_faces[6]);
  vofMo(){};
  void makeNormal();
  void getNormal(Real normalVec[3])const;
  void getFaces(faceMo faces[6])const;
  
  void setNormal(Real normalVec[3]);
  void normalize(Real normalVec[3]);
  Real getBdArea();
 };

class Moments
{
public: 
  
  Moments();
  
  Vector<Real> momentCalc3D(const int& a_order,
                             vofMo& a_vof);
  
  Vector<Real> momentCalc2D(const int& a_order,
                            const faceMo& a_face);
  
  void  listOfMoments(const int& a_order,
                      Vector<IntVect>& a_exponents);
  
  void makeMatrix(const Vector<IntVect>& a_list,
                  const Vector<IntVect>& a_listPlus,
                  const Vector<Real>& a_normalVec,
                  Real** a_A,
                  const int& a_faceNormal=999
                  );//use default faceNormal for 3D matrix or when SpaceDim==2
};



#endif

---