At Wed, 03 Sep 2008 11:37:51 -0500,
Alejandro Cruz Marcelo wrote:
> When using the funciton gsl_blas_dgemv to perform the operation
> x = Ax
> it returns x equal to 0. I guess we are not supposed to use the same
> vector but that warning is not in the manual.
> However, the operation
> x = Ax + x
> it is performed properly.
Thanks for the test program.
The two vectors x and y in "y = alpha * A * x + beta * y" must be
independent, this true of all BLAS functions in general.
I have added a note to the BLAS chapter to make this clearer.
--
Brian Gough
GNU Scientific Library -
http://www.gnu.org/software/gsl/
