forwarded message from David Clifford

[Brian Gough <bjg at network-theory dot co dot uk>]

--- Begin Message ---

• From: David Clifford <clifford at galton dot uchicago dot edu>
• To: Brian Gough <bjg at network-theory dot co dot uk>
• Date: Thu, 25 Apr 2002 18:28:13 -0500 (CDT)
• Subject: Re: problems with gsl_linalg_SV_decomp

Brian,

thanks for the response.  I have attached 3 files

matrix.txt
matrix2.txt
results.txt
svd.c

the first 2 files each contain a 3 by 3 matrix.  The first matrix is 
singular, the second is non-singular and the SVD decompostion is given in 
results.txt and were calculated using Splus.

the program svd.c reads in a matrix, calculates the svd and outputs the 
various pieces.  For matrix2 is works fine, for matrix1 it hangs at the 
svd step.

I compiled the program using a make file with the lines

svd:    svd.c
        g++ -g -o svd svd.c -lgsl -lgslcblas

I am running the program on red hat linux release 7.2 (enigma)

To run the program type

svd 

or

svd matrix2.txt

matrix.txt is the default file from which it will read in a matrix.  The 
dimensions (3 by 3) are specified in the code itself.

Attachment: matrix.txt
Description: Text document

Attachment: matrix2.txt
Description: Text document

Attachment: results.txt
Description: Text document

Attachment: svd.c
Description: Text document

--- End Message ---