[Bug-gsl] GSL _ gsl_blas_zgemv BUG REPORT

[Luca Simone Ronga]

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