[Bug-gsl] gsl_blas_dsyrk error

[kaska]

hello,
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.  so i tried using gsl_dgemm (which
does almost the same operation but doesn't take advantage
of symmetric matrix) and then the result was correct.
i also tried using the Fortran dsyrk function, which also
gave the correct result that matched with the answer from
gsl_dgemm.  i've used many other gsl blas routines and they
work fine, dsyrk is the one that hasn't worked properly.

you can compile the attached code with:
> make
and run it with:
> ./test

i also made up a text file of the outputs that i get when
i run all three functions.  you can see the output in 'testout'
file.

thank you
kaska

dsyrkerror.zip
Description: Zip archive