[Boost-users] matrix inversion

From: Umut Tabak (umut.tabak_at_[hidden])
Date: 2007-01-31 08:33:36

hi list,

I have been experimenting to find the incerse aof a matrix with ublas 
but there is a problem:

I have been looking at the code on

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

But there is sth wrong with that code because assign function template 
does not take a ublas::matrix<> as a direct parameter, it should ne 
converted to a matrix_expression but how may I do that.

I did the inversion by using the assignment operator but wondered how it 
could be done thorough assign function template.

The code on the viki works if you change the line with the assign to

inverse=ublas::identity_matrix<T>(A.size1());

Maybe it can be of help.

My best,