TPZTriangleTorus.cpp

Go to the documentation of this file.

//
//  TPZTriangleTorus.cpp
//  PZ
//
//  Created by Philippe Devloo on 3/21/14.
//
//

#include "TPZTriangleTorus.h"
#include "tpzgeomid.h"
#include "tpzgeoelmapped.h"
#include "tpzgeoelrefpattern.h"

namespace pzgeom {

// TPZGeoEl *TPZTriangleTorus::CreateBCGeoEl(TPZGeoEl *orig, int side,int bc)
// {
    
//         int ns = orig->NSideNodes(side);
//         TPZManVector<int64_t> nodeindices(ns);
//         int in;
//         for(in=0; in<ns; in++)
//         {
//             nodeindices[in] = orig->SideNodeIndex(side,in);
//         }
//         int64_t index;
        
//         TPZGeoMesh *mesh = orig->Mesh();
//         MElementType type = orig->Type(side);
        
//         TPZGeoEl *newel = mesh->CreateGeoBlendElement(type, nodeindices, bc, index);
//         TPZGeoElSide me(orig,side);
//         TPZGeoElSide newelside(newel,newel->NSides()-1);
        
//         newelside.InsertConnectivity(me);
//         newel->Initialize();
        
//         return newel;
// }
    
    // /**
    //  * Creates a geometric element according to the type of the father element
    //  */
    // /** @brief Creates a geometric element according to the type of the father element */
    // TPZGeoEl *TPZTriangleTorus::CreateGeoElement(TPZGeoMesh &mesh, MElementType type,
    //                                   TPZVec<int64_t>& nodeindexes,
    //                                   int matid,
    //                                   int64_t& index)

    // {
    //     return CreateGeoElementMapped(mesh,type,nodeindexes,matid,index);
    // }

    int TPZTriangleTorus::ClassId() const{
        return Hash("TPZTriangleTorus") ^ TPZGeoTriangle::ClassId() << 1;
    }

    void pzgeom::TPZTriangleTorus::InsertExampleElement(TPZGeoMesh &gmesh, int matid, TPZVec<REAL> &lowercorner, TPZVec<REAL> &size)
    {
        REAL R = 1., r = 0.8;
        size[0] = 2.;
        size[1] = 2.;
        size[2] = 1.;
        TPZManVector<REAL,3> origin(lowercorner);
        origin[0] += 1.;
        TPZTriangleTorus torus;
        torus.SetOrigin(origin);
        torus.SetDataRadius(R, r);
        TPZFNMatrix<12,REAL> phitheta(2,3,0.);
        phitheta(0,0) = 0.;
        phitheta(0,1) = M_PI/5.;
        phitheta(0,2) = M_PI;
        
        phitheta(1,0) = M_PI;
        phitheta(1,1) = 2*M_PI;
        phitheta(1,2) = 2*M_PI-M_PI/5.;
        torus.SetDataPhiTheta(phitheta);
        TPZManVector<int64_t,4> indexes(3);
        REAL coords[3][2] = {
            {0,0},{1,0},{0,1}
        };
        for (int i=0; i<3; i++) {
            indexes[i] = gmesh.NodeVec().AllocateNewElement();
            TPZManVector<REAL,3> xco(3), loc(2);
            loc[0] = coords[i][0];
            loc[1] = coords[i][1];
            
            torus.X(phitheta, loc, xco);
            gmesh.NodeVec()[indexes[i]].Initialize(xco, gmesh);
        }
        TPZGeoElRefPattern<pzgeom::TPZTriangleTorus> *gel = new TPZGeoElRefPattern<pzgeom::TPZTriangleTorus>(indexes,matid,gmesh);
        gel->Geom().SetOrigin(origin);
        gel->Geom().SetDataRadius(R, r);
        gel->Geom().SetDataPhiTheta(phitheta);
    }

}

template class TPZRestoreClass< TPZGeoElRefPattern<pzgeom::TPZTriangleTorus>>;


template class TPZGeoElRefLess<pzgeom::TPZTriangleTorus>;