gsl_linalg_SV_decomp goes into infinite loop on reasonable data

Brian Gough bjg@network-theory.co.uk
Tue Dec 31 09:55:00 GMT 2002


tlb@trevorblackwell.com writes:
 > Below is a patch for test.c to try SVDs on various 3x3 matrices. It
 > provokes an infinite loop on about the 15th matrix.

I've fixed the SVD routine to handle this case properly now.  

The updated files are in CVS or you can download them individually
(linalg/svd.c, linalg/svdstep.c) from CVSweb on http://sources.redhat.com/gsl/

regards
Brian Gough



More information about the Gsl-discuss mailing list