Method analyzeVariable

SRC/runtime/runtime/BasicAnalysisBuilder.cpp:571–654

int 
BasicAnalysisBuilder::analyzeVariable(int numSteps, double dT, double dtMin, double dtMax, int Jd)
{

  // set some variables
  int result = 0;  
  double totalTimeIncr = numSteps * dT;
  double currentTimeIncr = 0.0;
  double currentDt = dT;

  // loop until analysis has performed the total time incr requested
  while (currentTimeIncr < totalTimeIncr) {

    if (theAnalysisModel->analysisStep(currentDt) < 0) {
      opserr << "DirectIntegrationAnalysis::analyze() - the AnalysisModel failed in newStepDomain";
      opserr << " at time " << theDomain->getCurrentTime() << "\n";
      theDomain->revertToLastCommit();
      return -2;
    }

    // check if domain has undergone change
    int stamp = theDomain->hasDomainChanged();
    if (stamp != domainStamp) {
      domainStamp = stamp;
      if (this->domainChanged() < 0) {
        opserr << "DirectIntegrationAnalysis::analyze() - domainChanged() failed\n";
        return -1;
      }
    }

    //
    // do newStep(), solveCurrentStep() and commit() as in regular
    // DirectINtegrationAnalysis - difference is we do not return
    // if a failure - we stop the analysis & resize time step if failure
    //

    if (theTransientIntegrator->newStep(currentDt) < 0) {
      result = -2;
    }


    if (result >= 0) {
      result = theAlgorithm->solveCurrentStep();
      if (result < 0) 
        result = -3;
    }    

    if (result >= 0) {
      result = theTransientIntegrator->commit();
      if (result < 0) 
        result = -4;
    }

    // if the time step was successful increment delta T for the analysis
    // otherwise revert the Domain to last committed state & see if can go on

    if (result >= 0) 
      currentTimeIncr += currentDt;