Bug#525113: [Pkg-octave-devel] Bug#525113: Bug#525113: Inconsistant complex matrix multiplication
On Thu, Apr 23, 2009 at 01:28:36PM +0200, Laurent Mazet wrote: > Hi Thomas, > > you was right, my problem is inthe atlas library: Thanks for the follow-up. > > $ LD_LIBRARY_PATH=/usr/lib /usr/bin/octave -q --eval '[1 2; 3 4] * [1; 1i]' > ans = >1 + 2i >3 + 4i > > $ lmt-linux ~ $ /usr/bin/octave -q --eval '[1 2; 3 4] * [1; 1i]' > ans = >1. + 0.i >3. + 0.i > > $ ldd /usr/bin/octave | grep -i lapack > liblapack.so.3gf => /usr/lib/sse2/atlas/liblapack.so.3gf (0xb653f000) > $ dpkg -S /usr/lib/sse2/atlas/liblapack.so.3gf > libatlas3gf-sse2: /usr/lib/sse2/atlas/liblapack.so.3gf > $ dpkg -l libatlas3gf-sse2 > ... > +++-==-==- > ii libatlas3gf-ss 3.6.0-24 Automatically Tuned Linear Algebra Software, > > Does I need to make a follow-up to the atlas package bug-report system? No, we will handle the reassignment. I want to test some things myself before assigning this to ATLAS. Thomas -- To UNSUBSCRIBE, email to debian-bugs-dist-requ...@lists.debian.org with a subject of "unsubscribe". Trouble? Contact listmas...@lists.debian.org
Bug#525113: [Pkg-octave-devel] Bug#525113: Bug#525113: Inconsistant complex matrix multiplication
Hi Thomas, you was right, my problem is inthe atlas library: $ LD_LIBRARY_PATH=/usr/lib /usr/bin/octave -q --eval '[1 2; 3 4] * [1; 1i]' ans = 1 + 2i 3 + 4i $ lmt-linux ~ $ /usr/bin/octave -q --eval '[1 2; 3 4] * [1; 1i]' ans = 1. + 0.i 3. + 0.i $ ldd /usr/bin/octave | grep -i lapack liblapack.so.3gf => /usr/lib/sse2/atlas/liblapack.so.3gf (0xb653f000) $ dpkg -S /usr/lib/sse2/atlas/liblapack.so.3gf libatlas3gf-sse2: /usr/lib/sse2/atlas/liblapack.so.3gf $ dpkg -l libatlas3gf-sse2 ... +++-==-==- ii libatlas3gf-ss 3.6.0-24 Automatically Tuned Linear Algebra Software, Does I need to make a follow-up to the atlas package bug-report system? Regards, Laurent Quoting Thomas Weber : On Wed, Apr 22, 2009 at 11:04:14PM +0200, Thomas Weber wrote: On Wed, Apr 22, 2009 at 12:01:19PM +0200, Laurent Mazet wrote: > Package: octave3.0 > Version: 1:3.0.1-7 > Arch: i386 > Severity: grave > > Hi, > > I've just realized that I can multiply a real 2x2 matrix by a complex vector. Uh, yes. Why shouldn't this work? Or in other words, how do you distinguish the real matrix from a complex matrix with its complex coefficients being zero [ 1, 2; 3,4] is the same as [1+0i, 2+0i; 3+0i, 4+0i], isn't it? Sorry, I just realized that the problem was the result, not the act of multiplication. Anyway, which BLAS/ATLAS libraries are installed on your system? Thomas -- Dr. Laurent Mazet -=- "Use the source, Luke" -=- ma...@softndesign.org -- To UNSUBSCRIBE, email to debian-bugs-dist-requ...@lists.debian.org with a subject of "unsubscribe". Trouble? Contact listmas...@lists.debian.org
Bug#525113: [Pkg-octave-devel] Bug#525113: Bug#525113: Inconsistant complex matrix multiplication
On Wed, Apr 22, 2009 at 11:04:14PM +0200, Thomas Weber wrote: > On Wed, Apr 22, 2009 at 12:01:19PM +0200, Laurent Mazet wrote: > > Package: octave3.0 > > Version: 1:3.0.1-7 > > Arch: i386 > > Severity: grave > > > > Hi, > > > > I've just realized that I can multiply a real 2x2 matrix by a complex > > vector. > > Uh, yes. Why shouldn't this work? Or in other words, how do you > distinguish the real matrix from a complex matrix with its complex > coefficients being zero > > [ 1, 2; 3,4] is the same as [1+0i, 2+0i; 3+0i, 4+0i], isn't it? Sorry, I just realized that the problem was the result, not the act of multiplication. Anyway, which BLAS/ATLAS libraries are installed on your system? Thomas -- To UNSUBSCRIBE, email to debian-bugs-dist-requ...@lists.debian.org with a subject of "unsubscribe". Trouble? Contact listmas...@lists.debian.org
Bug#525113: [Pkg-octave-devel] Bug#525113: Bug#525113: Bug#525113: Inconsistant complex matrix multiplication
On Thu, Apr 23, 2009 at 01:35:01AM +0200, Rafael Laboissiere wrote: > * Thomas Weber [2009-04-22 23:04]: > > > On Wed, Apr 22, 2009 at 12:01:19PM +0200, Laurent Mazet wrote: > > > Package: octave3.0 > > > Version: 1:3.0.1-7 > > > Arch: i386 > > > Severity: grave > > > > > > Hi, > > > > > > I've just realized that I can multiply a real 2x2 matrix by a complex > > > vector. > > > > Uh, yes. Why shouldn't this work? Or in other words, how do you > > distinguish the real matrix from a complex matrix with its complex > > coefficients being zero > > > > [ 1, 2; 3,4] is the same as [1+0i, 2+0i; 3+0i, 4+0i], isn't it? > > I think Laurent meant "I've just realized that I CANNOT multiply [...]" Ah, sorry. I only briefly glanzed over the actual number. On my system, the results were right, though. Thomas -- To UNSUBSCRIBE, email to debian-bugs-dist-requ...@lists.debian.org with a subject of "unsubscribe". Trouble? Contact listmas...@lists.debian.org
Bug#525113: [Pkg-octave-devel] Bug#525113: Bug#525113: Bug#525113: Inconsistant complex matrix multiplication
On 23-Apr-2009, Rafael Laboissiere wrote: | * Thomas Weber [2009-04-22 23:04]: | | > On Wed, Apr 22, 2009 at 12:01:19PM +0200, Laurent Mazet wrote: | > > Package: octave3.0 | > > Version: 1:3.0.1-7 | > > Arch: i386 | > > Severity: grave | > > | > > Hi, | > > | > > I've just realized that I can multiply a real 2x2 matrix by a complex vector. | > | > Uh, yes. Why shouldn't this work? Or in other words, how do you | > distinguish the real matrix from a complex matrix with its complex | > coefficients being zero | > | > [ 1, 2; 3,4] is the same as [1+0i, 2+0i; 3+0i, 4+0i], isn't it? | | I think Laurent meant "I've just realized that I CANNOT multiply [...]" And the multiplication appeared to work, but gave the wrong answer. I don't remember this bug, but I can't duplicate it now on my system. If it wasn't a bug that was fixed in Octave, then I'd guess it was a BLAS bug? Octave just uses BLAS to perform these mutliplications. jwe -- To UNSUBSCRIBE, email to debian-bugs-dist-requ...@lists.debian.org with a subject of "unsubscribe". Trouble? Contact listmas...@lists.debian.org
Bug#525113: [Pkg-octave-devel] Bug#525113: Bug#525113: Inconsistant complex matrix multiplication
* Thomas Weber [2009-04-22 23:04]: > On Wed, Apr 22, 2009 at 12:01:19PM +0200, Laurent Mazet wrote: > > Package: octave3.0 > > Version: 1:3.0.1-7 > > Arch: i386 > > Severity: grave > > > > Hi, > > > > I've just realized that I can multiply a real 2x2 matrix by a complex > > vector. > > Uh, yes. Why shouldn't this work? Or in other words, how do you > distinguish the real matrix from a complex matrix with its complex > coefficients being zero > > [ 1, 2; 3,4] is the same as [1+0i, 2+0i; 3+0i, 4+0i], isn't it? I think Laurent meant "I've just realized that I CANNOT multiply [...]" -- Rafael -- To UNSUBSCRIBE, email to debian-bugs-dist-requ...@lists.debian.org with a subject of "unsubscribe". Trouble? Contact listmas...@lists.debian.org
Bug#525113: [Pkg-octave-devel] Bug#525113: Inconsistant complex matrix multiplication
On Wed, Apr 22, 2009 at 12:01:19PM +0200, Laurent Mazet wrote: > Package: octave3.0 > Version: 1:3.0.1-7 > Arch: i386 > Severity: grave > > Hi, > > I've just realized that I can multiply a real 2x2 matrix by a complex vector. Uh, yes. Why shouldn't this work? Or in other words, how do you distinguish the real matrix from a complex matrix with its complex coefficients being zero [ 1, 2; 3,4] is the same as [1+0i, 2+0i; 3+0i, 4+0i], isn't it? Thomas -- To UNSUBSCRIBE, email to debian-bugs-dist-requ...@lists.debian.org with a subject of "unsubscribe". Trouble? Contact listmas...@lists.debian.org
Bug#525113: Inconsistant complex matrix multiplication
Package: octave3.0 Version: 1:3.0.1-7 Arch: i386 Severity: grave Hi, I've just realized that I can multiply a real 2x2 matrix by a complex vector. Regards, Laurent octave:1> [1, 0.1; 0.68358, 0.72985] * [1; 1i] ans = 1.0 + 0.0i 0.68358 + 0.0i octave:2> [1, 1i] * [1, 0.1; 0.68358, 0.72985]' ans = 1.0 + 0.1i 0.68358 + 0.72985i -- Dr. Laurent Mazet -=- "Use the source, Luke" -=- ma...@softndesign.org -- To UNSUBSCRIBE, email to debian-bugs-dist-requ...@lists.debian.org with a subject of "unsubscribe". Trouble? Contact listmas...@lists.debian.org