pzshapelinear.h

Go to the documentation of this file.

#ifndef SHAPELINEARHPP
#define SHAPELINEARHPP

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

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

namespace pzshape {
                
                class TPZShapeLinear : public pztopology::TPZLine{
                                
                public:

                                static void (*fOrthogonal)(REAL x,int num,TPZFMatrix<REAL> & phi,TPZFMatrix<REAL> & dphi);
                                
                                static void Chebyshev(REAL x,int num,TPZFMatrix<REAL> & phi,TPZFMatrix<REAL> & dphi);
                                
        static void Expo(REAL x,int num,TPZFMatrix<REAL> & phi,TPZFMatrix<REAL> & dphi);
        
                                static void Legendre(REAL x,int num,TPZFMatrix<REAL> & phi,TPZFMatrix<REAL> & dphi);
                                
                                static void Jacobi(REAL x,int num,TPZFMatrix<REAL> & phi,TPZFMatrix<REAL> & dphi);
        
        static REAL JacobiA;
        static REAL JacobiB;
        
                                static void Hermite(REAL x,int num,TPZFMatrix<REAL> & phi,TPZFMatrix<REAL> & dphi);
                                
                                static void Legendre(REAL x,int num,TPZFMatrix<REAL> & phi,TPZFMatrix<REAL> & dphi, int nderiv);

#ifdef _AUTODIFF

                                static void (*FADfOrthogonal)(FADREAL& x,int num,TPZVec<FADREAL> &phi);
                                static void Chebyshev(FADREAL & x,int num,TPZVec<FADREAL> &phi);
#endif

public:
                                
                                static void Shape(TPZVec<REAL> &pt, TPZVec<int64_t> &id, TPZVec<int> &order,TPZFMatrix<REAL> &phi,TPZFMatrix<REAL> &dphi);
        
        static void ShapeCorner(TPZVec<REAL> &pt,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 ShapeInternal(TPZVec<REAL> &x,int ord,TPZFMatrix<REAL> &phi,TPZFMatrix<REAL> &dphi,int transformation_index);
        
        static void ShapeInternal(TPZVec<REAL> &x,int ord,TPZFMatrix<REAL> &phi,TPZFMatrix<REAL> &dphi);
                                
                                static void ShapeGenerating(TPZVec<REAL> &pt, TPZFMatrix<REAL> &phi, TPZFMatrix<REAL> &dphi);
                                
#ifdef _AUTODIFF

                                static void ShapeInternal(FADREAL & x,int num,TPZVec<FADREAL> & phi,int transformation_index);
#endif
                                
                                static void TransformPoint1d(int transid,REAL in,REAL &out);
#ifdef _AUTODIFF

                                /*
                                 * The transformation applied to compensate for odd/even brokerage between elements can also
                                 * be viewed by the transformation of a variational parameter.
                                 */
                                static void TransformPoint1d(int transid,FADREAL & in,FADREAL &out);
#endif

                                static void TransformDerivative1d(int transid,int num,TPZFMatrix<REAL> &in);

                                static int GetTransformId1d(TPZVec<int64_t> &id);

                                static int NConnectShapeF(int side, int order);
                                
                                static int NShapeF(TPZVec<int> &order);
        static void TransformPoint1d(int transid,double &val) ;
        static TPZTransform<REAL> ParametricTransform(int transid);
        static void ShapeInternal(int side, TPZVec<REAL> &x, int order,TPZFMatrix<REAL> &phi, TPZFMatrix<REAL> &dphi);
                                
                };
                
};

#endif