Gundala Viswanath wrote:
Is there a boost library that allow inner product
between sparse matrix and vector?
The uBlas library
(http://www.boost.org/doc/libs/1_37_0/libs/numeric/ublas/doc/index.htm)
might help you here. The following code compiles at least, so there is
good hope that it does more or less what you want. (There is a special
mailing list for ublas, if you need more information/help.)
#include
#include
void test()
{
boost::numeric::ublas::mapped_matrix<double> A(3,3, 3*3);
boost::numeric::ublas::compressed_matrix<double> B(3,3, 3*3);
boost::numeric::ublas::coordinate_matrix<double> C(3,3, 3*3);
boost::numeric::ublas::vector<double> x(3);
boost::numeric::ublas::vector<double> a =
boost::numeric::ublas::prod(A,x);
boost::numeric::ublas::vector<double> b =
boost::numeric::ublas::prod(B,x);
boost::numeric::ublas::vector<double> c =
boost::numeric::ublas::prod(C,x);
boost::numeric::ublas::vector<double> d =
boost::numeric::ublas::prod(x,C);
}
Regards,
Thomas
