Bug???

Bug???

[Przemyslaw Sliwa]

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

Re: Bug???

[Brian Gough]

Przemyslaw Sliwa writes:
 > 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 alpha in the documentation is spurious.  Cut and paste error from
dgemv I should think.  I'll fix it.

Thanks,

-- 
Brian Gough