diff --git a/OndselSolver/NewtonRaphson.cpp b/OndselSolver/NewtonRaphson.cpp index f73c55f..2023125 100644 --- a/OndselSolver/NewtonRaphson.cpp +++ b/OndselSolver/NewtonRaphson.cpp @@ -112,7 +112,7 @@ bool NewtonRaphson::isConvergedToNumericalLimit() size_t nDivergenceMax = 3; auto dxNormIterNo = dxNorms->at(iterNo); if (iterNo > 0) { - auto dxNormIterNoOld = dxNorms->at(iterNo); + auto dxNormIterNoOld = dxNorms->at(iterNo - 1); auto farTooLargeError = dxNormIterNo > tooLargeTol; auto worthIterating = dxNormIterNo > (smallEnoughTol * pow(10.0, (iterNo / iterMax) * nDecade)); bool stillConverging;