Thu 23 Feb 2012 11:42:32 AM UTC, comment #1:
I forgot to mention that the dgesvd parameters m,3,n=5,lda=3,ldvt=1 is causing a call to dorglq with m=5,n=5,lda=lvdt=1.
LDA violates the requirement for dorglq for LDA>=max(1,m): http://www.netlib.no/netlib/lapack/double/dorglq.f
I don't understand the dgesvd requirements enough to know if the call dgesvd is correct, because it's its internal call to dorglq that fails.
http://www.netlib.no/netlib/lapack/double/dgesvd.f
For the gdb backtrace before, I used "break dorglq_" and then run the test in octave.
|
Thu 23 Feb 2012 11:33:15 AM UTC, original submission:
I'm trying octave built with liblapack 3.4.0 on GNU/Linux x86_64 and qrupdate 1.1.2. qrupdate, lapack, and atlas tests work all fine. I get the same outcome with octave 3.4.3 and 3.6.0.
Some tests of the octave "make check" fail for me, related to lapack. I don't list them all here.
I pick one of the tests as an example, as at least someone got it working fine in his system, but not me. I run this test:
This fails for me this way:
Building lapack and octave with debug symobls, I could check with gdb where all this happens. Backtrace:
jwe at #octave claimed that the test works for his octave.
Can someone reproduce this?
|