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