tpzquadraticprism.h

Go to the documentation of this file.

#ifndef TPZQUADRATICPRISM_H
#define TPZQUADRATICPRISM_H


#include "pzgeoprism.h"
#include "pzgeoel.h"
#include "pznoderep.h"

#include <iostream>

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

    //virtual void ParametricDomainNodeCoord(int node, TPZVec<REAL> &nodeCoord);
                

                TPZQuadraticPrism(TPZVec<int64_t> &nodeindexes) : TPZRegisterClassId(&TPZQuadraticPrism::ClassId),
    pzgeom::TPZNodeRep<NNodes,pztopology::TPZPrism>(nodeindexes)
                {
                }
                
                TPZQuadraticPrism() : TPZRegisterClassId(&TPZQuadraticPrism::ClassId),
    pzgeom::TPZNodeRep<NNodes,pztopology::TPZPrism>()
                {
                }
                
                TPZQuadraticPrism(const TPZQuadraticPrism &cp,std::map<int64_t,int64_t> & gl2lcNdMap) : TPZRegisterClassId(&TPZQuadraticPrism::ClassId), pzgeom::TPZNodeRep<NNodes,pztopology::TPZPrism>(cp,gl2lcNdMap)
                {
                }
                
                TPZQuadraticPrism(const TPZQuadraticPrism &cp) : TPZRegisterClassId(&TPZQuadraticPrism::ClassId),pzgeom::TPZNodeRep<NNodes,pztopology::TPZPrism>(cp)
                {
                }
                
                TPZQuadraticPrism(const TPZQuadraticPrism &cp, TPZGeoMesh &) : TPZRegisterClassId(&TPZQuadraticPrism::ClassId),pzgeom::TPZNodeRep<NNodes,pztopology::TPZPrism>(cp)
                {
                }
                
                static std::string TypeName() { return "Prism";}
    
    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> &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> &par, 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