From: Kresimir Fresl (fresl_at_[hidden])
Date: 2008-02-04 01:59:36
Arun Tom Abraham wrote:
> however, what does one need to do if the "b" in the equation Ax = b is
> a general matrix .
> I think the lapack sgesv does allow one to enter a matrix in general .
Yes, it does.
ublas::matrix<float, ublas::column_major> A (n, n); // system matrix
ublas::matrix<float, ublas::column_major> b (n, nrhs);
// right-hand side matrix
// nrhs - number of right-hand side vectors
lapack::gesv (A, b);
I think that in your first post (dated 2/2/08) problem was that you
didn't specify that matrices are column_major:
Hope this helps,