tpzellipse3d.h

Go to the documentation of this file.

#ifndef TPZELLIPSE3D_H
#define TPZELLIPSE3D_H

#include "pzgeoel.h"
#include "pznoderep.h"
#include "tpzline.h"

#include <iostream>

class TPZGeoMesh;

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

        
        virtual void ParametricDomainNodeCoord(int64_t node, TPZVec<REAL> &nodeCoord);
        
                                TPZEllipse3D(const TPZEllipse3D &cp,std::map<int64_t,int64_t> & gl2lcNdMap) : TPZRegisterClassId(&TPZEllipse3D::ClassId),pzgeom::TPZNodeRep<NNodes,pztopology::TPZLine>(cp,gl2lcNdMap),
            fsAxeX(cp.fsAxeX), fsAxeY(cp.fsAxeY), fAxes(cp.fAxes), fOrigin(cp.fOrigin), fAngleIni(cp.fAngleIni), fAngleFinal(cp.fAngleFinal)
        {
                                }
                                TPZEllipse3D() : TPZRegisterClassId(&TPZEllipse3D::ClassId),pzgeom::TPZNodeRep<NNodes,pztopology::TPZLine>(),
        fsAxeX(0.), fsAxeY(0.), fAxes(), fOrigin(),fAngleIni(0.), fAngleFinal(0.)
        {
                                }
                                virtual ~TPZEllipse3D()
                                {
                                }
                                TPZEllipse3D(const TPZEllipse3D &cp) : TPZRegisterClassId(&TPZEllipse3D::ClassId),pzgeom::TPZNodeRep<NNodes,pztopology::TPZLine>(cp),
        fsAxeX(cp.fsAxeX), fsAxeY(cp.fsAxeY), fAxes(cp.fAxes), fOrigin(cp.fOrigin), fAngleIni(cp.fAngleIni), fAngleFinal(cp.fAngleFinal)
        {
                                }
                                TPZEllipse3D(const TPZEllipse3D &cp, TPZGeoMesh &) : TPZRegisterClassId(&TPZEllipse3D::ClassId),pzgeom::TPZNodeRep<NNodes, pztopology::TPZLine>(cp),
        fsAxeX(cp.fsAxeX), fsAxeY(cp.fsAxeY), fAxes(cp.fAxes), fOrigin(cp.fOrigin), fAngleIni(cp.fAngleIni), fAngleFinal(cp.fAngleFinal)
        {
                                }
                                TPZEllipse3D(TPZVec<int64_t> &nodeindexes) : TPZRegisterClassId(&TPZEllipse3D::ClassId),pzgeom::TPZNodeRep<NNodes,pztopology::TPZLine>(nodeindexes),
        fsAxeX(0.), fsAxeY(0.), fAxes(), fOrigin(),fAngleIni(0.), fAngleFinal(0.)
        {
                                }
        
        TPZEllipse3D &operator=(const TPZEllipse3D &cp)
        {
            pzgeom::TPZNodeRep<NNodes,pztopology::TPZLine>::operator=(cp);
            fsAxeX = cp.fsAxeX;
            fsAxeY = cp.fsAxeY;
            fAxes = cp.fAxes;
            fOrigin = cp.fOrigin;
            fAngleIni = cp.fAngleIni;
            fAngleFinal = cp.fAngleFinal;
            return *this;
        }
                                
                                virtual void AdjustNodesCoordinates(TPZGeoMesh &mesh);
                                
                                void SetAxes(TPZVec<REAL> Origin, TPZVec<REAL> SemiAxeX, TPZVec<REAL> SemiAxeY, TPZGeoMesh &gmesh);
                                
                                /* 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(TPZFMatrix<REAL> &cornerco, TPZVec<T> &par, TPZFMatrix<T> &gradx) const;
    
        
        template<class T>
                                void X(TPZFMatrix<REAL> &nodeCoord,TPZVec<T> &qsi,TPZVec<T> &x) const;
                                
                                
        static bool IsLinearMapping(int side)
        {
            return false;
        }
        
                                static std::string TypeName() { return "Linear";}
                                // static TPZGeoEl * CreateBCGeoEl(TPZGeoEl *orig, int side,int bc);
                                
                                // static TPZGeoEl *CreateGeoElement(TPZGeoMesh &mesh, MElementType type,
                                //                                                                                                                                TPZVec<int64_t>& nodeindexes,
                                //                                                                                                                                int matid,
                                //                                                                                                                                int64_t& index);
                                
                                void GetNodesCoords(TPZGeoMesh &mesh, TPZFMatrix<REAL> &nodes);
        
                                void SetNodesCoords(TPZGeoMesh &mesh, TPZFMatrix<REAL> &nodes);
                                
                                double sAxeX()
                                {
                                                return this->fsAxeX;
                                }
                                double sAxeY()
                                {
                                                return this->fsAxeY;
                                }
                                TPZVec<REAL> Origin()
                                {
                                                return this->fOrigin;
                                }
                                
        static void InsertExampleElement(TPZGeoMesh &gmesh, int matid, TPZVec<REAL> &lowercorner, TPZVec<REAL> &size);
        

                private:
                                
                                
        double ComputeAngle(TPZVec<REAL> &co) const;
        
                                
                                double fsAxeX;//norma de fSemiAxeX
                                
                                double fsAxeY;//norma de fSemiAxeY
        
        // for any point on the 3D ellips with coordinate co:
        // fAxes*(co-fOrigin) = coordinates in the local system (the third coordinate should be zero)
        TPZFNMatrix<9,REAL> fAxes;
                                
                                TPZManVector<REAL,3> fOrigin;
        
        double fAngleIni = 0.;
        double fAngleFinal = 0.;
        
                };
                
};

#endif