Convergence order and efficiency

Definition of convergence order and the error equation, their computational estimators COC and ACOC, the efficiency indices and the Kung-Traub conjecture defining optimal methods.

Convergence order and the error equation

Equivalently, writing the error as ek=xkαe_k=x_k-\alpha, the method has order pp if and only if it satisfies the error equation:

ek+1=Cekp+O(ekp+1)e_{k+1}=C\,e_k^p+\mathcal{O}\bigl(e_k^{p+1}\bigr)

The error equation is obtained with Taylor expansions, as in the proof of Newton's order or in the fixed-point theorem.

Measuring the order in practice: COC and ACOC

Given the sequence of iterates, the theoretical order is estimated with the computational order of convergence (COC), which requires knowing α\alpha:

COC=ln(xk+1α/xkα)ln(xkα/xk1α),k=1,2,COC=\frac{\ln\bigl(|x_{k+1}-\alpha|/|x_k-\alpha|\bigr)}{\ln\bigl(|x_k-\alpha|/|x_{k-1}-\alpha|\bigr)},\qquad k=1,2,\dots

Since in practice α\alpha is unknown, the error is replaced by the difference between consecutive iterates, giving the approximate computational order of convergence (ACOC):

ACOC=ln(xk+1xk/xkxk1)ln(xkxk1/xk1xk2),k=2,3,ACOC=\frac{\ln\bigl(|x_{k+1}-x_k|/|x_k-x_{k-1}|\bigr)}{\ln\bigl(|x_k-x_{k-1}|/|x_{k-1}-x_{k-2}|\bigr)},\qquad k=2,3,\dots

It is the same philosophy as order estimation in ODE methods: compare how errors shrink under refinement. The Newton exercise shows the ACOC tending to 2.

Efficiency and optimal methods

High order is not free: each iteration may require several evaluations of ff and its derivatives. To compare methods one defines the efficiency index I=p1/dI=p^{1/d}, where pp is the order and dd the number of distinct functional evaluations per iteration, and the computational efficiency index IC=p1/(d+op)IC=p^{1/(d+op)}, which also counts the opop products and quotients of each iteration.

MethodppddIIOptimal?
Newton2221/21.4142^{1/2}\approx 1.414
Halley3331/31.4423^{1/3}\approx 1.442
Chebyshev3331/31.4423^{1/3}\approx 1.442
Super-Halley3331/31.4423^{1/3}\approx 1.442
Ostrowski4341/31.5874^{1/3}\approx 1.587
Order, evaluations per iteration and efficiency index of several methods.

Newton is optimal (2=2212=2^{2-1}); third-order methods with three evaluations are not (3<231=43<2^{3-1}=4). Multipoint methods like Ostrowski or Jarratt reach order 4 with only three evaluations: they are optimal, which is why they stand out in the numerical comparison.