pzgeopoint.h

Go to the documentation of this file.

#ifndef TPZGEOPOINTH
#define TPZGEOPOINTH

#include <map>            // for map, operator==
#include <sstream>        // for basic_stringbuf<>::int_type, basic_stringbu...
#include "pzeltype.h"     // for MElementType
#include "pzmatrix.h"     // for TPZFMatrix, TPZMatrix, TPZFNMatrix
#include "pznoderep.h"    // for TPZNodeRep
#include "pznoderep.h.h"  // for TPZNodeRep::TPZNodeRep<N, Topology>
#include "pzreal.h"       // for REAL
#include "tpzpoint.h"     // for TPZPoint
class TPZGeoEl;
class TPZGeoMesh;
template <class T> class TPZVec;


namespace pzgeom {
                
                class TPZGeoPoint : public TPZNodeRep<1, pztopology::TPZPoint> {
                                
                public:
                                enum {NNodes = 1};
                
                int ClassId() const override;
                void Read(TPZStream &buf, void *context) override;
                
                void Write(TPZStream &buf, int withclassid) const override;
        
                                struct TMem {
                                };

        typedef pztopology::TPZPoint Top;
                                
                                TPZGeoPoint(TPZVec<int64_t> &nodeindexes) : TPZRegisterClassId(&TPZGeoPoint::ClassId),
        TPZNodeRep<NNodes, pztopology::TPZPoint>(nodeindexes)
                                {
                                }
                                
                                TPZGeoPoint() : TPZRegisterClassId(&TPZGeoPoint::ClassId),TPZNodeRep<NNodes, pztopology::TPZPoint>()
                                {
                                }
                                
                                TPZGeoPoint(const TPZGeoPoint &cp,
                                                                                std::map<int64_t,int64_t> & gl2lcNdMap) : TPZRegisterClassId(&TPZGeoPoint::ClassId),
        TPZNodeRep<NNodes, pztopology::TPZPoint>(cp,gl2lcNdMap)
                                {
                                }
                                
                                TPZGeoPoint(const TPZGeoPoint &cp) : TPZRegisterClassId(&TPZGeoPoint::ClassId),
        TPZNodeRep<NNodes, pztopology::TPZPoint>(cp)
                                {
                                }
                                
                                TPZGeoPoint(const TPZGeoPoint &cp, TPZGeoMesh &) : TPZRegisterClassId(&TPZGeoPoint::ClassId),
        TPZNodeRep<NNodes, pztopology::TPZPoint>(cp)
                                {
                                }
        
        static bool IsLinearMapping(int side)
        {
            return true;
        }
                                
                                static std::string TypeName() { return "Point";}
        
        template<class T>
                                static void X(const TPZFMatrix<REAL> &nodes,TPZVec<T> &loc,TPZVec<T> &result);
        
        template<class T>
        static void GradX(const TPZFMatrix<REAL> &nodes,TPZVec<T> &loc, TPZFMatrix<T> &gradx);
        
                                // static TPZGeoEl *CreateBCGeoEl(TPZGeoEl *gel, int side,int bc);
                                
                public:
        
        /* @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 TPZGeoPoint::X(const TPZFMatrix<REAL> &coord,TPZVec<T> &loc,TPZVec<T> &result){
        for (int i=0;i<coord.Rows();i++){
            result[i] = coord.GetVal(i,0);
        }
    }

    template<class T>
    inline void TPZGeoPoint::GradX(const TPZFMatrix<REAL> &nodes,TPZVec<T> &loc, TPZFMatrix<T> &gradx){
        gradx.Resize(3,0);
    }


};

#endif