rref returns wrong results for sparse matrices

John W. Eaton jwe at bevo.che.wisc.edu
Mon Mar 3 19:05:40 CST 2008


On  4-Mar-2008, David Bateman wrote:

| Jaroslav Hajek wrote:
| > I've tracked this problem down to the statement
| > 	A ([pivot, r], c:cols) = A ([r, pivot], c:cols);
| >
| > which apparently confuses Octave when pivot == r. The attached
| > changeset is a quickfix for this; however, according to what John
| > observed there may be bugs buried more deeply in liboctave.
| >
| >   
| I believe the attached is a fix for this issue that addresses it within
| the sparse assignment operators.

I applied it to the main and 3.0.x branches of my archive.

Thanks,

jwe


More information about the Bug-octave mailing list