Method analyze

SRC/analysis/analysis/TransientDomainDecompositionAnalysis.cpp:161–232

int 
TransientDomainDecompositionAnalysis::analyze(double dT)
{
  int result = 0;
  Domain *the_Domain = this->getDomainPtr();

  // check for change in Domain since last step. As a change can
  // occur in a commit() in a domaindecomp with load balancing
  // this must now be inside the loop

  //opserr << "TransientDomainDecomp - hasDomainChanged\n";
  int stamp = the_Domain->hasDomainChanged();

  if (stamp != domainStamp) {
    domainStamp = stamp;
    result = this->domainChanged();
    if (result < 0) {
      opserr << "TransientDomainDecompositionAnalysis::analyze() - domainChanged failed";
      return -1;
    }	
  }
  // opserr << "TransientDomainDecomp - newStep\n";
  // result = theAnalysisModel->newStepDomain();
  if (result < 0) {
    opserr << "TransientDomainDecompositionAnalysis::analyze() - the AnalysisModel failed";
    opserr << " with domain at load factor ";
    opserr << the_Domain->getCurrentTime() << endln;
    the_Domain->revertToLastCommit();
    
    return -2;
  }

  // opserr << "TransientDomainDecomp - integrator newStep\n";
  result = theIntegrator->newStep(dT);

  // barrierCheck in PartitionedDomain
 // result = this->checkAllResult(result);

  if (result < 0) {
    opserr << "TransientDomainDecompositionAnalysis::analyze() - the Integrator failed";
    opserr << " with domain at load factor ";
    opserr << the_Domain->getCurrentTime() << endln;
    the_Domain->revertToLastCommit();
    theIntegrator->revertToLastStep();    
    return -2;
  }
 //  opserr << "TransientDomainDecomp - solveCurrentStep\n";
  result = theAlgorithm->solveCurrentStep();
  if (result < 0) {
    opserr << "TransientDomainDecompositionAnalysis::analyze() - the Algorithm failed";
    opserr << " with domain at load factor ";
    opserr << the_Domain->getCurrentTime() << endln;
    the_Domain->revertToLastCommit();	    
    theIntegrator->revertToLastStep();
    
    return -3;
  }