Re: [ublas] Problems with identity_matrix<> / lu_substitute()

Subject: Re: [ublas] Problems with identity_matrix<> / lu_substitute()
From: Marco Guazzone (marco.guazzone_at_[hidden])
Date: 2009-10-09 07:52:39

On Fri, Oct 9, 2009 at 12:58 PM,  <Pietro.Mele_at_[hidden]> wrote:
> Marco,
>
> with the vector instead of the matrix it works correctly.
>
> Thank you,
> Pietro
>

There is a case in which you have to pass a matrix to lu_substitute
(instead of a vector); this is for matrix inversion.

If you need to compute the matrix inverse, take a look at

http://www.crystalclearsoftware.com/cgi-bin/boost_wiki/wiki.pl?LU_Matrix_Inversion

or in your boost distribution at:

<boost-home>/libs/numeric/ublas/test/test_lu.cpp

Ciao!!

-- Marco