pzshapequad.h

Go to the documentation of this file.

#ifndef SHAPEQUADHPP
#define SHAPEQUADHPP

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

#ifdef _AUTODIFF
#include "fadType.h"
#endif

namespace pzshape{

                class TPZShapeQuad  : public pztopology::TPZQuadrilateral{

                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 ShapeCorner(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,int quad_transformation_index);
        
        static void ShapeInternal(TPZVec<REAL> &x, int order,
                                  TPZFMatrix<REAL> &phi,TPZFMatrix<REAL> &dphi);
                                
#ifdef _AUTODIFF

                                static void Shape2dQuadInternal(TPZVec<FADREAL> &x, int order,
                                                                                                                                                                TPZVec<FADREAL> &phi,int quad_transformation_index);
#endif

                                static void TransformDerivative2dQ(int transid, int num, TPZFMatrix<REAL> &in);
                                
                                static void TransformPoint2dQ(int transid,TPZVec<REAL> &in,TPZVec<REAL> &out);
                                
#ifdef _AUTODIFF
                                
                                static void TransformPoint2dQ(int transid,TPZVec<FADREAL> &in,TPZVec<FADREAL> &out);
                                
#endif

                                static void ProjectPoint2dQuadToRib(int rib, TPZVec<REAL> &in, REAL &out);
                                static int GetTransformId2dQ(TPZVec<int64_t> &id);

                                static void TransformDerivativeFromRibToQuad(int rib,int num,TPZFMatrix<REAL> &dphi);

                                static REAL gTrans2dQ[8][2][2];
                                static REAL gFaceTr2dQ[6][2][3];
                                static REAL gRibTrans2dQ1d[4][2];

                                static int NConnectShapeF(int side, int order);

                                static int NShapeF(TPZVec<int> &order);
        static TPZTransform<REAL> ParametricTransform(int trans_id);
        static void ShapeInternal(int side, TPZVec<REAL> &x, int order, TPZFMatrix<REAL> &phi, TPZFMatrix<REAL> &dphi);
        
                                
                };
                
};

#endif