Re: [Bug-gsl] bug in gsl_linalg_QRPT_QRsolve

[Brian Gough]

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     
----------------------------------------------------------------------