[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Help-gsl] gsl_linalg_LU_invert

From: Matthias Kramm
Subject: [Help-gsl] gsl_linalg_LU_invert
Date: Tue, 10 Apr 2007 11:38:21 +0200
User-agent: Mutt/1.5.13 (2006-08-11)

Hi all,

the online documentation gives the impression that
gsl_linalg_LU_invert() will invert a matrix previously 
postprocessed with gsl_linalg_LU_decomp().

I don't quite see how that can be possible, as 
gsl_linalg_LU_decomp() will destroy the original matrix
(transform it into an upper triagonal), and doesn't
save enough information to "reconstruct" it.
(only the permutation of the matrix rows is saved, but not
 the lower triagonal part (U))

Same goes to the function gsl_linalg_LU_solve().

Is this a bug, or am I overlooking something? In any case,
I can get neither of these functions to work correctly.
(A small example C file which should output a 
 unit matrix (but doesn't) is attached)



Attachment: gsl_inverse_test.c
Description: Text Data

reply via email to

[Prev in Thread] Current Thread [Next in Thread]