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

Jason Grout jason-sage at creativetrax.com
Thu Apr 28 07:04:52 EDT 2011

On 4/26/11 3:18 PM, Jason Grout wrote:
> On 4/26/11 11:49 AM, Pauli Virtanen wrote:
>> But apparently either there's a bug, or the LAPACK man page needs to
>> be understood as you say.
> I've posted a question to the Lapack forum:
> http://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=2&t=2402

We've got a reply now:


Pauli, you're right that the statement should be parsed as:

"If (JOBZ = 'A') or (JOBZ = 'O' and M >= N), VT contains the N-by-N 
unitary matrix V**H;"

so you're right that the function is written correctly in the case we 
are bringing up (but still has an error in assuming too much when JOBZ 
is not 'A').

We'll work on providing an example with the Accelerate framework and the 
zero row.



More information about the NumPy-Discussion mailing list