Boost logo

Ublas :

From: Gunter Winkler (guwi17_at_[hidden])
Date: 2006-12-08 10:20:50

On Friday 08 December 2006 15:59, Seabook wrote:
> I will try that but I still don't know how to use axpy_prod.

1) modify the line in operation.hpp marked with +/i in my diff
2) use axpy_prod(A,X,Y,true) for Y = prod(A,X)
   and axpy_prod(A,X,Y,false) for Y += prod(A,X)

> Could you paste your source code with axpy_prod operation?

I should have atached it already :(

BTW. ublas' product of dense matrices is a lot slower than any optimized blas3
routine. However, you can use the atlas-bindings to use blas instead of prod.