[Numpy-discussion] Zero row in SVD's unitary matrix on some Mac's

Pauli Virtanen pav at iki.fi
Tue Apr 26 14:25:22 EDT 2011


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




More information about the NumPy-Discussion mailing list