Boost logo

Ublas :

Subject: [ublas] Bug in Ublas' Inner Product between Matrix and Vector?
From: Gundala Viswanath (gundalav_at_[hidden])
Date: 2009-01-29 22:17:45

Hi all,

I am trying to perform an inner product multiplication
with ublas between this Matrix and Vector:

ublas::matrix <double> inputMat (5,5,0);

ublas::vector <double> inputVec (5);

The result I get is this:

Where as the correct one should be (using R):

Note the significant discrepancy in the 3rd element (*).

The product command I use is this:

ublas::vector <double> Result (5);
Result = ublas::prod(inputMat,inputVec);

Is there anything wrong with my way of doing inner product
or it is a bug?

- Gundala Viswanath
Jakarta - Indonesia