Method solveCurrentStep

SRC/analysis/algorithm/equiSolnAlgo/ExpressNewton.cpp:106–154

}

int 
ExpressNewton::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 *theAnalysisModel = this->getAnalysisModelPtr(); 
    LinearSOE  *theSOE = this->getLinearSOEptr();
    IncrementalIntegrator *theIntegrator = this->getIncrementalIntegratorPtr();

    if ((theAnalysisModel == 0) || (theIntegrator ==0 ) || (theSOE == 0)){
	opserr << "WARNING ExpressNewton::solveCurrentStep() -";
	opserr << "setLinks() has not been called.\n";
	return -5;
    }

	if (factorOnce != 2) {
		if (theIntegrator->formTangent(HALL_TANGENT, kMultiplier1, kMultiplier2) < 0) {
		  opserr << "WARNING ExpressNewton::solveCurrentStep() -";
		  opserr << "the Integrator failed in formTangent()\n";
		  return -1;
		}
		if (factorOnce == 1)
			factorOnce = 2;
    }

    for (int iter = 0; iter <nIter; ++iter)
    {
    if (theIntegrator->formUnbalance() < 0) {
	opserr << "WARNING ExpressNewton::solveCurrentStep() -";
	opserr << "the Integrator failed in formUnbalance()\n";	
	return -2;
    }

    if (theSOE->solve() < 0) {
	opserr << "WARNING ExpressNewton::solveCurrentStep() -";
	opserr << "the LinearSOE failed in solve()\n";	
	return -3;
    }

    if (theIntegrator->update(theSOE->getX()) < 0) {
	opserr << "WARNING ExpressNewton::solveCurrentStep() -";
	opserr << "the Integrator failed in update()\n";	
	return -4;
    }
    }

    return 0;
}

int
ExpressNewton::setConvergenceTest(ConvergenceTest *theNewTest)