Re: [Bug-gsl] Bugs in I/O for gsl_matrix_complex

[Rhys Ulerich]

Hi Matthias,

> I've discovered a bug in the I/O for complex matrices. If tda is larger
> than the actual size of the matrix, the fwrite() function called
> internally writes the wrong blocks of data to disk. In particular, the
> code containing the address of the memory block ("m->data + i * tda")
> does /not/ respect the MULTIPLICITY!

Thank you for the pleasantly thorough bug report and patch.  I added
some unit tests confirming the bug
(http://bzr.savannah.gnu.org/lh/gsl/trunk/revision/4764) and applied
your patch (http://bzr.savannah.gnu.org/lh/gsl/trunk/revision/4765).
It repaired the problem without regressing any of the previously
working real-valued cases.  I missed the NEWS update in 4765 and added
it in http://bzr.savannah.gnu.org/lh/gsl/trunk/revision/4766.

- Rhys