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.