well I think the gesv functions allows you to enter a matrix for "b" so that you can do multiple Ax=b operations in a single function call. It will do the calculation Ax=b for each column of b where A is an mxn matrix and x and b are n x nrhs matrices.
Jesse
Hi,
Thanks it worked by adding the line :however, what does one need to do if the "b" in the equation Ax = b is
#include <boost/numeric/bindings/traits/ublas_vector2.hpp>
a general matrix .
I think the lapack sgesv does allow one to enter a matrix in general .
Arun_______________________________________________
ublas mailing list
ublas@lists.boost.org
http://lists.boost.org/mailman/listinfo.cgi/ublas