tpzintrulet3d.h

Go to the documentation of this file.

#ifndef TPZINTRULET3D_H
#define TPZINTRULET3D_H

#include "pzreal.h"
#include "pzmanvector.h"

class TPZIntRuleT3D {
                
    friend class TPZIntRuleList;
    
    int            fNumInt;
    TPZManVector<long double>   fLocationKsi;
    TPZManVector<long double>   fLocationEta;
    TPZManVector<long double>   fLocationZeta;
    TPZManVector<long double>   fWeight;
    int fOrder;
                
    TPZIntRuleT3D(int i = 2);
    ~TPZIntRuleT3D();

    TPZIntRuleT3D(const TPZIntRuleT3D &copy) : fNumInt(copy.fNumInt), fLocationKsi(copy.fLocationKsi), fLocationEta(copy.fLocationEta),
        fLocationZeta(copy.fLocationZeta), fWeight(copy.fWeight), fOrder(copy.fOrder)
    {
        
    }
    TPZIntRuleT3D &operator=(const TPZIntRuleT3D &copy)
    {
        fNumInt = copy.fNumInt;
        fLocationKsi = copy.fLocationKsi;
        fLocationEta = copy.fLocationEta;
        fLocationZeta = copy.fLocationZeta;
        fWeight = copy.fWeight;
        fOrder = copy.fOrder;
        return *this;
    }
                int ComputingSymmetricCubatureRule(int order);
                
                void TransformBarycentricCoordInCartesianCoord(long double baryvec[],long double weightvec[]);

public:
                enum {NRULESTETRAHEDRA_ORDER = 14};
                
                int NInt() const { return fNumInt;}
                
                void Loc(int i, TPZVec<REAL> &pos) const;
                
                REAL W(int i) const;
    
    int Order()
    {
        return fOrder;
    }
};

#endif