pzshapecube.h

Go to the documentation of this file.

#ifndef SHAPECUBEHPP
#define SHAPECUBEHPP

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

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

namespace pzshape {
                
                class TPZShapeCube : public pztopology::TPZCube {
                                
                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);
        
                                
#ifdef _AUTODIFF

                                static void ShapeCube(TPZVec<REAL> &point, TPZVec<int64_t> &id, TPZVec<int> &order, TPZVec<FADREAL> &phi);

                                static void ShapeCornerCube(TPZVec<FADREAL> &pt, TPZVec<FADREAL> &phi);
                                
                                static void Shape3dCubeInternal(TPZVec<FADREAL> &x, int order,TPZVec<FADREAL> &phi);//,int quad_transformation_index
#endif
                                
                                static void ShapeCorner(TPZVec<REAL> &pt, TPZFMatrix<REAL> &phi, TPZFMatrix<REAL> &dphi);
                                
        static void ShapeGenerating(TPZVec<REAL> &pt, TPZFMatrix<REAL> &phi, TPZFMatrix<REAL> &dphi);
        
                private:
                                
                

                                static void ShapeInternal(TPZVec<REAL> &x, int order,TPZFMatrix<REAL> &phi,
                                                                                                                                  TPZFMatrix<REAL> &dphi);//,int quad_transformation_index
                                static void ProjectPoint3dCubeToRib(int face, TPZVec<REAL> &in, REAL &outval);
                                
#ifdef _AUTODIFF

                                static void ProjectPoint3dCubeToRib(int face, TPZVec<FADREAL> &in, FADREAL &outval);

                                static void ProjectPoint3dCubeToFace(int face, TPZVec<FADREAL> &in, TPZVec<FADREAL> &outval);
#endif
                                
                                static void ProjectPoint3dCubeSide(int face, TPZVec<REAL> &in, REAL &out);
                                
                                static void ProjectPoint3dCubeFace(int face, TPZVec<REAL> &in, TPZVec<REAL> &out);
                                
                                static void ProjectPoint3dCubeToFace(int face, TPZVec<REAL> &in, TPZVec<REAL> &outval);
                                
                                static void TransformDerivativeFromRibToCube(int rib,int num,TPZFMatrix<REAL> &dphi);

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

                                static REAL gFaceTrans3dCube2d[6][2][3];
                                static REAL gRibTrans3dCube1d[12][3];
                                
                public:
                                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);

                };
                
};

#endif