TPZGeoLinear.h

Go to the documentation of this file.

#ifndef TPZGEOLINEARH
#define TPZGEOLINEARH

#include "pznoderep.h"

#include "pzvec.h"
#include "pzeltype.h"
#include "pzfmatrix.h"
#include "tpzline.h"

#include <string>

class TPZGeoEl;
class TPZGeoMesh;

namespace pzgeom
{
    
    
    class TPZGeoLinear : public TPZNodeRep<2, pztopology::TPZLine> {
        
    public:
        typedef pztopology::TPZLine Top;
        enum {NNodes = 2};
        
        TPZGeoLinear(TPZVec<int64_t> &nodeindexes) : TPZRegisterClassId(&TPZGeoLinear::ClassId), TPZNodeRep<NNodes, pztopology::TPZLine>(nodeindexes)
        {
        }
        
        TPZGeoLinear() : TPZRegisterClassId(&TPZGeoLinear::ClassId),TPZNodeRep<NNodes, pztopology::TPZLine>()
        {
        }
        
        TPZGeoLinear(const TPZGeoLinear &cp, std::map<int64_t,int64_t> & gl2lcNdMap) : 
        TPZRegisterClassId(&TPZGeoLinear::ClassId), TPZNodeRep<NNodes, pztopology::TPZLine>(cp,gl2lcNdMap)
        {
        }
        
        TPZGeoLinear(const TPZGeoLinear &cp) : TPZRegisterClassId(&TPZGeoLinear::ClassId),
        TPZNodeRep<NNodes, pztopology::TPZLine>(cp)
        {
        }
        
        TPZGeoLinear(const TPZGeoLinear &cp, TPZGeoMesh &) : TPZRegisterClassId(&TPZGeoLinear::ClassId),
        TPZNodeRep<NNodes, pztopology::TPZLine>(cp)
        {
        }
        
        static bool IsLinearMapping(int side)
        {
            return true;
        }
        
        static std::string TypeName() { return "Linear";}
        
        template<class T>
        static void X(const TPZFMatrix<REAL> &nodes,TPZVec<T> &loc,TPZVec<T> &x);
        
        template<class T>
        static void GradX(const TPZFMatrix<REAL> &nodes,TPZVec<T> &loc, TPZFMatrix<T> &gradx);

        
        static void Jacobian(const TPZFMatrix<REAL> &nodes,TPZVec<REAL> &param,TPZFMatrix<REAL> &jacobian,
                             TPZFMatrix<REAL> &axes,REAL &detjac,TPZFMatrix<REAL> &jacinv);
        
        // /**
        //  * @brief Method which creates a geometric boundary condition
        //  * element based on the current geometric element, \n
        //  * a side and a boundary condition number
        //  */
        // static TPZGeoEl *CreateBCGeoEl(TPZGeoEl *gel, int side,int bc);
        
        public:
        int ClassId() const override;
        void Read(TPZStream &buf, void *context) override;
        void Write(TPZStream &buf, int withclassid) const override;

        
        /* @param gmesh mesh in which the element should be inserted
         @param matid material id of the element
         @param lowercorner (in/out) on input lower corner o the cube where the element should be created, on exit position of the next cube
         @param size (in) size of space where the element should be created
         */
        static void InsertExampleElement(TPZGeoMesh &gmesh, int matid, TPZVec<REAL> &lowercorner, TPZVec<REAL> &size);
        
        // static TPZGeoEl *CreateGeoElement(TPZGeoMesh &mesh,
        //                                   MElementType type,
        //                                   TPZVec<int64_t>& nodeindexes,
        //                                   int matid,
        //                                   int64_t& index);
        
    };

    
    template<class T>
    inline void TPZGeoLinear::X(const TPZFMatrix<REAL> &nodes,TPZVec<T> &loc,TPZVec<T> &x){
        
        T xi = loc[0];
        int nrow = nodes.Rows();
        for(int i = 0; i < nrow; i++)
        {
            x[i] = nodes.GetVal(i,0)*(1.-xi)*0.5+nodes.GetVal(i,1)*(1.+xi)*0.5;
        }
        
    }
    
    
    template<class T>
    inline void TPZGeoLinear::GradX(const TPZFMatrix<REAL> &nodes,TPZVec<T> &loc, TPZFMatrix<T> &gradx){
        
        int nrow = nodes.Rows();
        int ncol = nodes.Cols();
        
        gradx.Resize(nrow,1);
        gradx.Zero();
        
        TPZFNMatrix<3,T> phi(2,1);
        TPZFNMatrix<6,T> dphi(2,2);
        TShape(loc,phi,dphi);
        for(int i = 0; i < ncol; i++)
        {
            for(int j = 0; j < nrow; j++)
            {
                gradx(j,0) += nodes.GetVal(j,i)*dphi(0,i);

            }
        }
        
    }
    
};

#endif