Method save

SRC/element/PFEMElement/PFEMMesher2D.cpp:1174–1241

}

int 
PFEMMesher2D::save(const char* filename, Domain* theDomain, int maxelenodes)
{
    if(theDomain == 0) {
        opserr<<"WARNING: null domain";
        opserr<<" -- PFEMMesher2D::save\n";
        return -1;
    }

    // get nodes
    std::map<int,int> fluidNodes;
    ID regions;
    theDomain->getRegionTags(regions);
    getNodes(regions,fluidNodes,theDomain);

    std::ofstream outfile(std::string(filename).append(".node").c_str());
    outfile<<theDomain->getCurrentTime()<<" ";
    for(int i=0; i<8; i++) outfile<<0<<" ";
    outfile<<"\n";

    // write nodes
    for(std::map<int,int>::iterator it=fluidNodes.begin(); it!=fluidNodes.end(); it++) {
        int tag = it->first;
        int type = it->second;
        Node* node = theDomain->getNode(tag);
        if(node == 0) {
            opserr<<"WARNING: node "<<tag<<" dose not exist -- ";
            opserr<<" -- PFEMMesher2D::save\n";
            return -1;
        }
        const Vector& coord = node->getCrds();
        const Vector& disp = node->getDisp();
        const Vector& vel = node->getVel();
        const Vector& accel = node->getAccel();
        outfile<<tag<<" "<<coord(0)+disp(0)<<" "<<coord(1)+disp(1)<<" ";
        outfile<<type<<" ";
        outfile<<vel(0)<<" "<<vel(1)<<" "<<accel(0)<<" "<<accel(1);
        Pressure_Constraint* thePC = theDomain->getPressure_Constraint(tag);
        if(thePC != 0) {
            outfile<<" "<<thePC->getPressure()<<"\n";
        } else {
            outfile<<" "<<0<<"\n";
        }
    }
    outfile.close();

    // write elements
    outfile.open(std::string(filename).append(".ele").c_str());
    ElementIter& theEles = theDomain->getElements();
    Element* theEle = 0;
    while((theEle = theEles()) != 0) {
        const ID& ntags = theEle->getExternalNodes();
        int num = 0;
        for(int i=0; i<ntags.Size(); i++) {
            if(fluidNodes.find(ntags(i)) != fluidNodes.end()) {
                outfile<<ntags(i)<<" ";
                num++;
            }