Stepsize control

Stepsize control is an important issue in these algorithms. Too small stepsizes lead to unnecessary work being done, while too big stepsizes can lead to losing details of the curve. An easily implementable and proven to be reliable method is convergence-dependent control.

Consider the computation of a next point using step size $h_i$. If the computation converged, let $n$ denote the number of Newton iterations needed. Then the new step size $h_{i+1}$ will be selected as follows:

\begin{displaymath}
h_{i+1} = \left\{
\begin{array}{ll}
h_i\cdot h_{dec} &...
... n_{thr} \\
h_i & \mbox{otherwise}
\end{array}
\right.
\end{displaymath} (22)

where $h_{dec}<1$, $h_{inc}>1$ and $n_{thr}$ are constants which are experimentally determined.