This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: gsl_linalg_SV_decomp goes into infinite loop on reasonable data
- From: Brian Gough <bjg at network-theory dot co dot uk>
- To: tlb at trevorblackwell dot com
- Cc: gsl-discuss at sources dot redhat dot com
- Date: Fri, 14 Jun 2002 21:59:13 +0100 (BST)
- Subject: Re: gsl_linalg_SV_decomp goes into infinite loop on reasonable data
- References: <200206131921.g5DJLe511585@tlb.blackwells.net>
tlb@trevorblackwell.com writes:
> A much simpler matrix for which gsl_linalg_SV_decomp goes into an
> infinite loop is:
>
> 0 0 0
> 0 0 1
> 1 1 0
>
> 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 don't understand the Golub-Reinsch algorithm well enough to figure
> out what the problem is.
Thanks. I see what the problem is (Golub & Van Loan, p454, missing
case of d_n = 0)