tpzquadratictrig.h

Go to the documentation of this file.

#ifndef TPZQUADRATICTRIG_H
#define TPZQUADRATICTRIG_H

#include "pznoderep.h"
#include "tpztriangle.h"

class TPZGeoMesh;

namespace pzgeom {
    
    class TPZQuadraticTrig : public pzgeom::TPZNodeRep<6,pztopology::TPZTriangle> {
                                
                public:
        typedef pztopology::TPZTriangle Top;
                                enum {NNodes = 6};
                
                public:
int ClassId() const override;

        
        //irtual void ParametricDomainNodeCoord(int node, TPZVec<REAL> &nodeCoord);
        
        static bool IsLinearMapping(int side)
        {
            return false;
        }
        
                                TPZQuadraticTrig(TPZVec<int64_t> &nodeindexes) : 
        TPZRegisterClassId(&TPZQuadraticTrig::ClassId),
        pzgeom::TPZNodeRep<NNodes,pztopology::TPZTriangle>(nodeindexes) {
                                }
                                TPZQuadraticTrig() : TPZRegisterClassId(&TPZQuadraticTrig::ClassId),
        pzgeom::TPZNodeRep<NNodes,pztopology::TPZTriangle>() {
                                }
                                TPZQuadraticTrig(const TPZQuadraticTrig &cp,std::map<int64_t,int64_t> & gl2lcNdMap) : TPZRegisterClassId(&TPZQuadraticTrig::ClassId),
        pzgeom::TPZNodeRep<NNodes,pztopology::TPZTriangle>(cp,gl2lcNdMap) {
                                }
                                TPZQuadraticTrig(const TPZQuadraticTrig &cp) : TPZRegisterClassId(&TPZQuadraticTrig::ClassId),
        pzgeom::TPZNodeRep<NNodes,pztopology::TPZTriangle>(cp) {
                                }
                                TPZQuadraticTrig(const TPZQuadraticTrig &cp, TPZGeoMesh &) : TPZRegisterClassId(&TPZQuadraticTrig::ClassId),
        pzgeom::TPZNodeRep<NNodes,pztopology::TPZTriangle>(cp) {
                                }
                                static std::string TypeName() { return "Triangle";}
                                
                                // /**
                                //  * @brief Method which creates a geometric boundary condition 
                                //  * element based on the current geometric element, 
                                //  * a side and a boundary condition number
                                //  */
                                // static  TPZGeoEl * CreateBCGeoEl(TPZGeoEl *orig,int side,int bc);
                                
        static void Shape(TPZVec<REAL> &loc,TPZFMatrix<REAL> &phi,TPZFMatrix<REAL> &dphi){
            TShape(loc, phi, dphi);
        }
        
                                
                                /* brief compute the coordinate of a point given in parameter space */
//        template<class T>
//        void X(const TPZGeoEl &gel,TPZVec<T> &loc,TPZVec<T> &result) const
//        {
//            TPZFNMatrix<3*NNodes> coord(3,NNodes);
//            CornerCoordinates(gel, coord);
//            X(coord,loc,result);
//        }
        
//        template<class T>
//        void GradX(const TPZGeoEl &gel, TPZVec<T> &loc, TPZFMatrix<T> &gradx) const
//        {
//            TPZFNMatrix<3*NNodes> coord(3,NNodes);
//            CornerCoordinates(gel, coord);
//            GradX(coord,loc,gradx);
//        }
                                
        template<class T>
        static void TShape(const TPZVec<T> &param,TPZFMatrix<T> &phi,TPZFMatrix<T> &dphi);
        
        template<class T>
                                static void X(const TPZFMatrix<REAL> &coord, TPZVec<T> &par, TPZVec< T > &result);
        
        template<class T>
        static void GradX(const TPZFMatrix<REAL> &nodes,TPZVec<T> &loc, TPZFMatrix<T> &gradx);
                                
                public:
                                // /** @brief Creates a geometric element according to the type of the father element */
                                // static TPZGeoEl *CreateGeoElement(TPZGeoMesh &mesh, MElementType type,
                                //                                                                                                                                TPZVec<int64_t>& nodeindexes,
                                //                                                                                                                                int matid, int64_t& index);

        static void InsertExampleElement(TPZGeoMesh &gmesh, int matid, TPZVec<REAL> &lowercorner, TPZVec<REAL> &size);

                };

};

#endif