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/ _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
