Bug???

[Przemyslaw Sliwa <sliwa at euv-frankfurt-o dot de>]

I guess I have found a bug in the gsl:

Look at definition:

Function: int gsl_blas_dtrmv (CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t TransA,
CBLAS_DIAG_t Diag, const gsl_matrix * A, gsl_vector * x).

And additionaly the manual:

These functions compute the matrix-vector product x =\alpha op(A) x for
the triangular matrix A, where op(A) = A, A^T, A^H for TransA =
CblasNoTrans, CblasTrans, CblasConjTrans. When Uplo is CblasUpper then the
upper triangle of A is used, and when Uplo is CblasLower then the lower
triangle of A is used. If Diag is CblasNonUnit then the diagonal of the
matrix is used, but if Diag is CblasUnit then the diagonal elements of the
matrix A are taken as unity and are not referenced.

I do not see the \alpha factor in the definition of this funcion.

Do I miss something?

Cheers,

PShem