tpzintrulep3d.h

Go to the documentation of this file.

#ifndef TPZINTRULEP3D_H
#define TPZINTRULEP3D_H

#include "pzreal.h"
#include "pzvec.h"

class TPZIntRuleP3D {
                
                friend class TPZIntRuleList;
                
    int            fNumInt;
    TPZVec<long double>         fLocationKsi;
    TPZVec<long double>         fLocationEta;
    TPZVec<long double>         fLocationZeta;
    TPZVec<long double>         fWeight;
    int fOrder;

                TPZIntRuleP3D(int &order);
                ~TPZIntRuleP3D();
                

public:
                
    TPZIntRuleP3D(const TPZIntRuleP3D &copy) : fNumInt(copy.fNumInt), fLocationKsi(copy.fLocationKsi), fLocationEta(copy.fLocationEta),
        fLocationZeta(copy.fLocationZeta), fWeight(copy.fWeight), fOrder(copy.fOrder)
    {
        
    }
    
    TPZIntRuleP3D &operator=(const TPZIntRuleP3D &copy)
    {
        fNumInt = copy.fNumInt;
        fLocationKsi = copy.fLocationKsi;
        fLocationEta = copy.fLocationEta;
        fLocationZeta = copy.fLocationZeta;
        fWeight = copy.fWeight;
        fOrder = copy.fOrder;
        return *this;
    }
    
    enum {NRULESPYRAMID_ORDER = 37};
                
    int NInt() const{ return fNumInt; }
                
                void Loc(int i, TPZVec<REAL> &pos) const;
                
                REAL W(int i) const;
    
    int Order()
    {
        return fOrder;
    }

                void Print(std::ostream & out = std::cout);

protected:

                int ComputingCubatureRuleForPyramid(int order);

};

#endif