Method identify

SRC/analysis/analysis/PFEMAnalysis.cpp:144–256

}

int
PFEMAnalysis::identify()
{
    Domain* theDomain = this->getDomainPtr();

    if (theDomain->getNumPCs() == 0) return 0;

    // disconnect pcs from structural elements
    Pressure_ConstraintIter& thePCs = theDomain->getPCs();
    Pressure_Constraint* thePC = 0;
    while((thePC = thePCs()) != 0) {
        thePC->disconnect();
    }

    // to connect pc to all elements
    Element *elePtr;
    ElementIter &theElemIter = theDomain->getElements();    
    while((elePtr = theElemIter()) != 0) {

        const ID& nodeids = elePtr->getExternalNodes();
        for(int i=0; i<nodeids.Size(); i++) {
            Pressure_Constraint* thePC = theDomain->getPressure_Constraint(nodeids(i));
            if(thePC != 0) {
                thePC->connect(elePtr->getTag(), false);
            }
        }
    }

    // get all SPs and MPs
    SP_ConstraintIter& theSPs = theDomain->getDomainAndLoadPatternSPs();
    SP_Constraint* theSP = 0;
    std::map<int,ID> sps;
    while((theSP = theSPs()) != 0) {
        sps[theSP->getNodeTag()].insert(theSP->getDOF_Number());
    }

    // get number of parameters
    int numParams = theDomain->getNumParameters();

    // move the isolated nodes
    Pressure_ConstraintIter& thePCs1 = theDomain->getPCs();
    while((thePC = thePCs1()) != 0) {

        if(thePC->isIsolated()) {
            int tag = thePC->getTag();
            Node* node = theDomain->getNode(tag);
            if(node == 0) {
                opserr<<"WARNING: node "<<tag;
                opserr<<" does not exist -- PFEMAnalysis::identify\n";
		return -1;
            }
            std::map<int,ID>::iterator it = sps.find(tag);

            // responses
            const Vector& disp = node->getDisp();
            const Vector& vel = node->getVel();
            const Vector& accel = node->getAccel();
            Vector ndisp = disp;