[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] GSL _ gsl_blas_zgemv BUG REPORT
From: |
Luca Simone Ronga |
Subject: |
[Bug-gsl] GSL _ gsl_blas_zgemv BUG REPORT |
Date: |
Wed, 08 Sep 2004 18:24:52 +0200 |
GSL v.1.1.1 running on Linux 2.4 SMP i686. Compiler gcc version 3.2
20020903 (Red Hat Linux 8.0 3.2-7).
The following code WORKS:
gsl_blas_zgemv(CblasTrans,
gsl_complex_rect(1.0,0),
transform_mat,
tmp,
gsl_complex_rect(0,0),
tmp1);
but the following code returns an error:
gsl_blas_zgemv(CblasConjTrans,
gsl_complex_rect(1.0,0),
transform_mat,
tmp,
gsl_complex_rect(0,0),
tmp1);
The error is
gsl: blas.c:685: ERROR: invalid length
Since the only difference is a traspose operation instead of a conjugate
transpose (the matrices are correctly allocated) I suspect a BUG in the
gls-to-blas interface.
I did not find this bug corrected in the successive versions of GSL (1.5 is
the most updated at the time of writing).
Thank you.
---
Luca Simone Ronga - Ph.D.
CNIT
Head of Research
University of Florence Unit
email: address@hidden
sms: address@hidden
www: http://lenst.det.unifi.it/~ronga
phone: +39-3488757025
fax:+39-3488700766
- [Bug-gsl] GSL _ gsl_blas_zgemv BUG REPORT,
Luca Simone Ronga <=