Re: [Bug-gsl] gsl_blas_dsyrk error

[Brian Gough]

address@hidden writes:
 > i'm sending you some small code with a makefile that
 > i believe shows an error in the gsl_blas_dsyrk function.
 > i had to use that function to build a cholesky factorization
 > with matrix-matrix operations and i was getting incorrect
 > results with gsl_dsyrk.

Thanks for the bug report.  

I have found the problem and will fix it in the next release (see
below for a temporary patch).

-- 
best regards,

Brian Gough

Network Theory Ltd,
Commercial support for GSL --- http://www.network-theory.co.uk/gsl/

--- blas.c      26 Jul 2003 13:44:31 -0000      1.28
+++ blas.c      8 Dec 2004 17:26:15 -0000
@@ -1650,7 +1650,7 @@
 {
   const size_t M = C->size1;
   const size_t N = C->size2;
-  const size_t K = (Trans == CblasNoTrans) ? A->size1 : A->size2;
+  const size_t K = (Trans == CblasNoTrans) ? A->size2 : A->size1;
 
   if (M != N)
     {