Method doTriangulation

SRC/element/PFEMElement/PFEMMesher2D.cpp:847–907

linear elements

// linear elements
int
PFEMMesher2D::doTriangulation(int starteletag, double alpha, const ID& groups, 
                              const ID& addgroups,Domain* theDomain, 
                              double rho, double mu, double b1, double b2, 
                              double thk, double kappa, int type, int& endele)
{
    if(theDomain == 0) {
        opserr<<"WARNING: null domain";
        opserr<<" -- PFEMMesher2D::doTriangulation\n";
        return -1;
    }

    //Timer timer;
    //timer.start();

    // do triangulation
    ID eles;
    int res = 0;
    res = doTriangulation(alpha,groups,addgroups,theDomain,eles);
    if(res < 0) {
        opserr<<"WARNING: failed to do triangulation --";
        opserr<<"PFEMMesher2D::soTriangulation\n";
        return res;
    }

    // add PFEM elements
    int numeles = eles.Size()/3;
    if(numeles == 0) return 0;
    int etag = starteletag-1;
    for(int i=0; i<numeles; i++) {
        Element* theEle = 0;
        if(type == 1) {
            theEle = new PFEMElement2D(++etag, eles(3*i), eles(3*i+1), eles(3*i+2),rho, mu, b1, b2, thk);
        } else if(type == 3) {
            theEle = new PFEMElement2DCompressible(++etag, eles(3*i), eles(3*i+1), eles(3*i+2),rho, mu, b1, b2, thk, kappa);
        } else if(type == 4) {
            theEle = new PFEMElement2DBubble(++etag, eles(3*i), eles(3*i+1), eles(3*i+2),rho, mu, b1, b2, thk, kappa);
        }
        
        if(theEle == 0) {
            opserr<<"WARNING: no enough memory -- ";
            opserr<<" -- PFEMMesher2D::doTriangulation\n";
            return -1;
        }
        if(theDomain->addElement(theEle) == false) {
            opserr<<"WARNING: failed to add element to domain -- ";
            opserr<<" -- PFEMMesher2D::doTriangulation\n";
            delete theEle;
            return -1;
        }
    }

    // identify
    identify(b2,theDomain);
    //timer.pause();
    //opserr<<"meshing :"<<timer.getCPU()<<"\n";

    endele = etag;