pzshapeprism.h

Go to the documentation of this file.

#ifndef SHAPEPRISMHPP
#define SHAPEPRISMHPP

#include "pzfmatrix.h"
#include "pzstack.h"
#include "pztrnsform.h"
#include "tpzprism.h"
#include "pzshtmat.h"

namespace pzshape {
                
                class TPZShapePrism : public pztopology::TPZPrism {
                                
                public:
                                
                                static void Shape(TPZVec<REAL> &pt, TPZVec<int64_t> &id, TPZVec<int> &order,
                                                                                                  TPZFMatrix<REAL> &phi,TPZFMatrix<REAL> &dphi);
                                
                                static void SideShape(int side, TPZVec<REAL> &pt, TPZVec<int64_t> &id, TPZVec<int> &order,
                                                                                                                  TPZFMatrix<REAL> &phi,TPZFMatrix<REAL> &dphi);
        
        static void ShapeOrder(TPZVec<int64_t> &id, TPZVec<int> &order, TPZGenMatrix<int> &shapeorders);//, TPZVec<int64_t> &sides
        
        static void SideShapeOrder(int side,  TPZVec<int64_t> &id, int order, TPZGenMatrix<int> &shapeorders);
        
                                
                                static void CornerShape(TPZVec<REAL> &pt, TPZFMatrix<REAL> &phi, TPZFMatrix<REAL> &dphi);

                                static void ShapeGenerating(TPZVec<REAL> &pt, TPZFMatrix<REAL> &phi, TPZFMatrix<REAL> &dphi);
                                
                                static void ShapeInternal(TPZVec<REAL> &x, int order,TPZFMatrix<REAL> &phi,
                                                                                                                                  TPZFMatrix<REAL> &dphi);
                                
                                static void ProjectPoint3dPrismaToRib(int rib, TPZVec<REAL> &in, REAL &outval);
                                
                                static void ProjectPoint3dPrismaToFace(int face, TPZVec<REAL> &in, TPZVec<REAL> &outval);
                                
                                static void TransformPoint3dPrismaFace(int transid, int face, TPZVec<REAL> &in, TPZVec<REAL> &out);
                                
                                
                                static void TransformDerivativeFromRibToPrisma(int rib,int num,TPZFMatrix<REAL> &dphi);
                                
                                static void TransformDerivativeFromFaceToPrisma(int face,int num,TPZFMatrix<REAL> &dphi);
                                
                                static void TransformDerivativeFace3dPrisma(int transid, int face, int num, TPZFMatrix<REAL> &in);
                                
                                static REAL gFaceTrans3dPrisma2d[5][2][3];
                                
                                static REAL gFaceSum3dPrisma2d[5][2];
                                
                                static REAL gRibTrans3dPrisma1d[9][3];
                                
                                static REAL gRibSum3dPrisma1d[9];
                                
                                static int NConnectShapeF(int side, int order);
                                
                                static int NShapeF(TPZVec<int> &order);
        
       static void ShapeInternal(int side, TPZVec<REAL> &x, int order, TPZFMatrix<REAL> &phi, TPZFMatrix<REAL> &dphi);
       static void ShapeCorner(TPZVec<REAL> &pt, TPZFMatrix<REAL> &phi, TPZFMatrix<REAL> &dphi);
                                
                };
                
};

#endif