[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

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

From: Matthias Sitte
Subject: [Bug-gsl] Bugs in I/O for gsl_matrix_complex
Date: Fri, 14 Oct 2011 14:54:12 +0200
User-agent: Mozilla/5.0 (X11; U; Linux i686; en-US; rv: Gecko/20110831 Thunderbird/3.1.13


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!

I've written a small program to demonstrate the bug with some small
complex-valued matrix that is stored and read immediately afterwards
from disk.

I've also included a patch proposal for "matrix/file_source.c".

I'd be happy if you could verify/falsify my findings.



Attachment: matrix_fwrite_bug.c
Description: Text Data

Attachment: file_source.patch
Description: Text Data

Attachment: signature.asc
Description: OpenPGP digital signature

reply via email to

[Prev in Thread] Current Thread [Next in Thread]