PreviousIndexNext

[bug-gsl] Re: gsl_blas_dgemv return a 0 vector when using a vector twice (2008-08-30)

From: Brian Gough <bjg@gnu.org>
Subject: Re: gsl_blas_dgemv return a 0 vector when using a vector twice
Date: Sat, 30 Aug 2008 11:25:19 +0100
To: Alejandro Cruz-Marcelo <ac6@........>
Cc: bug-gsl@gnu.org

At Fri, 29 Aug 2008 09:44:30 -0500,
Alejandro Cruz-Marcelo wrote:
> gsl_blas_dgemv(CblasNoTrans,
> 1.0, A, y,
> 1.0, y)
> then y becomes a vector of 0's and no error is produced, even though the
> operation A*y + y makes perfect sense. I believe that in the
> documentation there is no warning about using the same vector twice in
> the function.

Please can you send an example program that demonstrates the problem,
thanks.

--
Brian Gough