SVD bugs in gsl version .9.1?

[Brian Gough]

Devin Balkcom writes:
 > I'm using gsl to take singular value decompositions.  I have had
 > problems with all three of the
 > implementations.   I have written a test program.
 > 
 > - gsl_linalg_SV_decomp() does not return for the test matrix.
 > - gsl_linalg_SV_decomp_mod() also does not return.
 > - gsl_linalg_SV_decomp_jacobi() does return, although I have had it fail
 > in other circumstances,
 >     with an incorrect result.  (I will send test code for this case
 > later.)

Thanks for the bug report.  

It is a bug in the SVD.... which is not handling the 2x2 case [0,X;0,0]
correctly, so it never converges for that matrix. I'll fix that.

regards
Brian Gough