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
_______________________________________________
Bug-gsl mailing list
[email protected]
https://lists.gnu.org/mailman/listinfo/bug-gsl
