Method analyze

SRC/analysis/analysis/PFEMAnalysis.cpp:81–142

}

int 
PFEMAnalysis::analyze(bool flush)
{
    Domain* theDomain = this->getDomainPtr();
    double current = theDomain->getCurrentTime();
    if(newstep) {
        next = current + dtmax;
    }
    bool instep = false;
	
    while(true) {

	// identify domain
	if (this->identify() < 0) {
	    opserr<<"WARNING: failed to identify domain -- PFEMAnalysis\n";
	    return -1;
	}

        // analyze
#ifdef _PARALLEL_INTERPRETERS
	int myid;
	MPI_Comm_rank(MPI_COMM_WORLD, &myid);
	if(myid == 0) {
	    opserr<<"\n\nTime = "<<current<<", dt = "<<dt<<"\n\n";
	}
#else
        opserr<<"\n\nTime = "<<current<<", dt = "<<dt<<"\n\n";
#endif

	// analysis
        int converged  = DirectIntegrationAnalysis::analyze(1, dt);

        // if failed
        if(converged < 0) {
            dt *= ratio;
            if(dt < dtmin) {
                return -1;
            }
            instep = true;
            newstep = false;
            continue;
        }

        // if converged
        if(instep) {
            current = theDomain->getCurrentTime();
            dt = next - current;
        } else {
            newstep = true;
            dt = dtmax;
        }
	
	break;
    }

    Domain* the_Domain = this->getDomainPtr();
    if (the_Domain != 0 && flush) {
        the_Domain->flushRecorders();