Method solveCurrentStep

SRC/analysis/algorithm/equiSolnAlgo/NewtonHallM.cpp:111–201

int 
NewtonHallM::solveCurrentStep(void)
{
    // set up some pointers and check they are valid
    // NOTE this could be taken away if we set Ptrs as protecetd in superclass
    AnalysisModel   *theAnaModel = this->getAnalysisModelPtr();
    IncrementalIntegrator *theIntegrator = this->getIncrementalIntegratorPtr();
    //IncrementalIntegrator *theIntegratorSens=this->getIncrementalIntegratorPtr();//Abbas
    LinearSOE  *theSOE = this->getLinearSOEptr();

    if ((theAnaModel == 0) || (theIntegrator == 0) || (theSOE == 0)
	|| (theTest == 0)){
	opserr << "WARNING NewtonHallM::solveCurrentStep() - setLinks() has";
	opserr << " not been called - or no ConvergenceTest has been set\n";
	return -5;
    }	

    if (theIntegrator->formUnbalance() < 0) {
      opserr << "WARNING NewtonHallM::solveCurrentStep() -";
      opserr << "the Integrator failed in formUnbalance()\n";	
      return -2;
    }	    

    // set itself as the ConvergenceTest objects EquiSolnAlgo
    theTest->setEquiSolnAlgo(*this);
    if (theTest->start() < 0) {
      opserr << "NewtnRaphson::solveCurrentStep() -";
      opserr << "the ConvergenceTest object failed in start()\n";
      return -3;
    }

    int result = -1;
    numIterations = 0;

    do {

      int tangent = HALL_TANGENT;
      SOLUTION_ALGORITHM_tangentFlag = tangent;

      double iFact, cFact;
      if (method == 0) {
	iFact = iFactor*exp(-alpha*numIterations);
	cFact = 1.0 - iFact;
      } else if (method ==1) {
	double iFact0 = 1.0/(1. + exp(-alpha*c));
	iFact = 1/(1 + exp(alpha*(numIterations-c)));
	iFact = iFactor*iFact/iFact0;
	cFact = 1.0 - iFact;
      } else {
	iFact = iFactor;
	cFact = cFactor;
      }
      
      if (theIntegrator->formTangent(tangent, iFact, cFact) < 0){
	opserr << "WARNING NewtonHallM::solveCurrentStep() -";
	opserr << "the Integrator failed in formTangent()\n";
	return -1;
      }