  Iterative least-squares MR reconstructions typically use the Conjugate
Gradient algorithm, despite known numerical issues. This paper demonstrates
that the more recent LSMR algorithm has favourable numerical properties, and is
to be preferred in situations where Toeplitz embedding cannot be used to
accelerate the Conjugate Gradient method.
