  Conjugated gradients on the normal equation (CGNE) is a popular method to
regularise linear inverse problems. The idea of the method can be summarised as
minimising the residuum over a suitable Krylov subspace. It is shown that using
the same idea for the shift-and-invert rational Krylov subspace yields an
order-optimal regularisation scheme.
