[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] bug in gsl_linalg_QRPT_QRsolve
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] bug in gsl_linalg_QRPT_QRsolve |
Date: |
Thu, 8 May 2003 21:37:04 +0100 |
Mario Pernici writes:
> in linalg/qrpt.c
>
> int
> gsl_linalg_QRPT_QRsolve (const gsl_matrix * Q,...)
> {
> // ...
> else
> {
> /* compute b' = Q^T b */
> // gsl_blas_dgemv (CblasNoTrans, 1.0, Q, b, 0.0, x);
> // becomes
> gsl_blas_dgemv (CblasTrans, 1.0, Q, b, 0.0, x);
> // ...
> }
Thanks for the bug report!
I found the function QR_QRsolve was also affected by the same
error.
I have updated the functions in CVS to use the CblasTrans and
added some more tests to the test program to cover these
routines.
best regards
--
Brian Gough
----------------------------------------------------------------------
Network Theory Ltd Phone: +44 (0)117 3179309
15 Royal Park Fax: +44 (0)117 9048108
Bristol BS8 3AL WWW: http://www.network-theory.co.uk/
United Kingdom Email: address@hidden
----------------------------------------------------------------------
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- Re: [Bug-gsl] bug in gsl_linalg_QRPT_QRsolve,
Brian Gough <=