tpzquadraticcube.h

Go to the documentation of this file.

#ifndef TPZQUADRATICCUBE_H
#define TPZQUADRATICCUBE_H


#include "TPZGeoCube.h"
#include "pzgeoel.h"
#include "pznoderep.h"

#include <iostream>

namespace pzgeom {

    class TPZQuadraticCube : public pzgeom::TPZNodeRep<20,pztopology::TPZCube> {
            
    public:
        typedef pztopology::TPZCube Top;
        enum {NNodes = 20};
        
        public:
int ClassId() const override;

        
        //virtual void ParametricDomainNodeCoord(int node, TPZVec<REAL> &nodeCoord);
            
        TPZQuadraticCube(TPZVec<int64_t> &nodeindexes) : 
        TPZRegisterClassId(&TPZQuadraticCube::ClassId),
        pzgeom::TPZNodeRep<NNodes,pztopology::TPZCube>(nodeindexes)
        {
        }
        TPZQuadraticCube() : TPZRegisterClassId(&TPZQuadraticCube::ClassId),
        pzgeom::TPZNodeRep<NNodes,pztopology::TPZCube>()
        {
        }
        TPZQuadraticCube(const TPZQuadraticCube &cp,std::map<int64_t,int64_t> & gl2lcNdMap) : TPZRegisterClassId(&TPZQuadraticCube::ClassId),
        pzgeom::TPZNodeRep<NNodes,pztopology::TPZCube>(cp,gl2lcNdMap)
        {
        }
        TPZQuadraticCube(const TPZQuadraticCube &cp) : TPZRegisterClassId(&TPZQuadraticCube::ClassId),
        pzgeom::TPZNodeRep<NNodes,pztopology::TPZCube>(cp)
        {
        }
        TPZQuadraticCube(const TPZQuadraticCube &cp, TPZGeoMesh &) : TPZRegisterClassId(&TPZQuadraticCube::ClassId),
        pzgeom::TPZNodeRep<NNodes,pztopology::TPZCube>(cp)
        {
        }
        
        static std::string TypeName() { return "Hexa";}
        
        static bool IsLinearMapping(int side)
        {
            return false;
        }
        
        static void Shape(TPZVec<REAL> &loc,TPZFMatrix<REAL> &phi,TPZFMatrix<REAL> &dphi){
            TShape(loc, phi, dphi);
        }
        
        
        template<class T>
        static void TShape(const TPZVec<T> &x,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);

            // 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);

            // TPZGeoEl *CreateBCGeoEl(TPZGeoEl *orig,int side,int bc);         
    };
    
};

#endif