qr method failed to converge
Brian Gough
bjg@network-theory.co.uk
Wed Dec 19 13:20:00 GMT 2001
jean-michel.rouet@philips.com writes:
> Can anybody explain why such a hardcoded value "30" is used in the
> qr method ? and how it can we get rid of it (using an other
> hardcoded value such "40" as I did is not a good solution!)
The algorithm is from,
R.S. Martin, G. Peters and J.H. Wilkinson, ``The QR Algorithm for Real
Hessenberg Matrices'', Numerische Mathematik, 14 (1970), 219--231.
They do mention that the value of 30 in their code is adhoc. It might
have been used for single precision and so would need to be increased
now. I'd have to go back and study the paper to see what the rate of
convergence is and how that affects the number of iterations.
regards
Brian Gough
More information about the Gsl-discuss
mailing list