tpzintrulet.h

Go to the documentation of this file.

#ifndef TPZINTRULET_H
#define TPZINTRULET_H

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

class TPZIntRuleT {
                
                friend class TPZIntRuleList;
                
    int            fNumInt;
    TPZManVector<long double>   fLocationKsi;
    TPZManVector<long double>   fLocationEta;
    TPZManVector<long double>   fWeight;
    int fOrder;

                TPZIntRuleT(int order);
    
    TPZIntRuleT(const TPZIntRuleT &copy) : fNumInt(copy.fNumInt), fLocationKsi(copy.fLocationKsi), fLocationEta(copy.fLocationEta),
        fWeight(copy.fWeight), fOrder(copy.fOrder)
    {
        
    }
    
    TPZIntRuleT &operator=(const TPZIntRuleT &copy)
    {
        fNumInt = copy.fNumInt;
        fLocationKsi = copy.fLocationKsi;
        fLocationEta = copy.fLocationEta;
        fWeight = copy.fWeight;
        fOrder = copy.fOrder;
        return *this;
    }
                ~TPZIntRuleT();

                int ComputingSymmetricCubatureRule(int order);

                void TransformBarycentricCoordInCartesianCoord(long double baryvec[],long double weightvec[]);

public:
    enum {NRULESTRIANGLE_ORDER = 21};

    int NInt() const { return fNumInt;}

    void Loc(int i, TPZVec<REAL> &pos) const;

    REAL W(int i) const;
    
    int Order() const
    {
        return fOrder;
    }

};

#endif