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_Inv... 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,