Method discretize

SRC/element/PFEMElement/PFEMMesher3D.cpp:71–236

PLC

// PLC
int 
PFEMMesher3D::discretize(int startnode, const PointVec& points, const FacetVec facets,
                         const PointVec& holes, double maxvol, int ndf, 
                         const std::vector<int>& fix, const std::vector<double>& mass,
                         Domain* theDomain, int& endnode)
{
    if(theDomain == 0) {
        opserr<<"WARNING: null domain";
        opserr<<" -- PFEMMesher3D::discretize\n";
        return -1;
    }

    // structs
    tetgenio in, out;

    // input points
    in.numberofpoints = points.size();
    if(in.numberofpoints > 0) {
        in.pointlist = new double[in.numberofpoints*3];
    } else {
        in.pointlist = NULL;
    }
    for(int i=0; i<in.numberofpoints; i++) {
        const Point& point = points[i];
        in.pointlist[3*i] = point.x();
        in.pointlist[3*i+1] = point.y();
        in.pointlist[3*i+2] = point.z();
        //opserr<<"point "<<i<<": "<<point.x()<<" "<<point.y()<<" "<<point.z()<<"\n";
    }

    // input facets
    in.numberoffacets = facets.size();
    if(in.numberoffacets > 0) {
        in.facetlist = new tetgenio::facet[in.numberoffacets];
    } else {
        in.facetlist = NULL;
    }
    for(int i=0; i<in.numberoffacets; i++) {
        const Facet& facet = facets[i];
        const PolygonVec& polygons = facet.polygons();
        const PointVec& holes = facet.holes();

        // new facet
        tetgenio::facet& tfacet = in.facetlist[i];

        // polygons
        tfacet.numberofpolygons = polygons.size();
        if(tfacet.numberofpolygons > 0) {
            tfacet.polygonlist = new tetgenio::polygon[tfacet.numberofpolygons];
        } else {
            tfacet.polygonlist = NULL;
        }
        for(int j=0; j<tfacet.numberofpolygons; j++) {
            tetgenio::polygon& tpolygon = tfacet.polygonlist[j];
            const Polygon& polygon = polygons[j];
            tpolygon.numberofvertices = polygon.size();
            if(tpolygon.numberofvertices > 0) {
                tpolygon.vertexlist = new int[tpolygon.numberofvertices];