[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] gsl_blas_dgemv return a 0 vector when using a vector twice
From: |
Alejandro Cruz-Marcelo |
Subject: |
[Bug-gsl] gsl_blas_dgemv return a 0 vector when using a vector twice |
Date: |
Fri, 29 Aug 2008 09:44:30 -0500 |
User-agent: |
Thunderbird 2.0.0.16 (Windows/20080708) |
Hello
When using the function
gsl_blas_dgemv(CblasNoTrans,
1.0, A, x,
1.0, y)
It seems that when the SAME vector is used in the place of both x and y,
the program's output is a vector of 0's. That is, if the call of the
function is
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.
I hope this helps,
Thank you,
Alejandro
- [Bug-gsl] gsl_blas_dgemv return a 0 vector when using a vector twice,
Alejandro Cruz-Marcelo <=