This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: Matrix inversion.
- From: Brian Gough <bjg at network-theory dot co dot uk>
- To: Negin Nejati <negin at Glue dot umd dot edu>
- Cc: gsl-discuss at sourceware dot cygnus dot com
- Date: Wed, 30 Jan 2002 19:10:17 +0000 (GMT)
- Subject: Re: Matrix inversion.
- References: <Pine.GSO.4.21.0201291706550.12978-100000@y.glue.umd.edu>
Negin Nejati writes:
> Hi, I am a new user of this library and I have a silly question:
> For inverting a matrix shouldn't I first use "gsl_linalg_LU_decomp"
> function to decompose my matrix and then use "gsl_linalg_LU_invert"
> passing the results of previous function to it, to invert the
> matrix? If it is correct how come the multiplication of the
> original matrix and it's reverse is not I? I'll appreciate any
> help, NN.
What you've described should work. If it doesn't please send your
program as a bug report.
Thanks,
--
Brian Gough
----------------------------------------------------------------------
Network Theory Ltd Phone: 0117 3179309 (+44 117 3179309)
15 Royal Park WWW: http://www.network-theory.co.uk/
Bristol BS8 3AL Email: bjg@network-theory.co.uk
United Kingdom
----------------------------------------------------------------------