April 26, 2011
2:25 p.m.
On Tue, 26 Apr 2011 10:33:07 -0500, Jason Grout wrote: [clip]
I was just looking up the documentation for ZGESDD and noticed that the value we have for rwork in the numpy call [1] does not match the Lapack docs. This was changed in Lapack 3.2.2 [2]. I've submitted a pull request:
https://github.com/numpy/numpy/pull/78
but I have not tested the change.
I doubt this fix will fix the problem on this thread, but it makes sense to make the change anyway.
I think it is worth testing if this change fixes it. Note that current LAPACK docs have the number 7 there instead of 5, so there may really be a bug here. Pauli