#11248: SVD matrix decomposition may return a non-invertible "unitary" matrix
------------------------------+---------------------------------------------
Reporter: rbeezer | Owner: jason, was
Type: defect | Status: new
Priority: major | Milestone: sage-4.7
Component: linear algebra | Keywords:
Work_issues: | Upstream: N/A
Reviewer: | Author:
Merged: | Dependencies:
------------------------------+---------------------------------------------
Comment(by jason):
Can you doing the following:
1. Put the code at
http://software.intel.com/sites/products/documentation/hpc/mkl/lapack/mkl_lapack_examples/zgesdd_ex.c.htm
in a file called `test_zgesdd.c`
2. On the command line, compile this file with `gcc test_zgesdd.c -o
test_zgesdd -framework Accelerate`
3. Run the program: `./test_zgesdd` and post the result. This uses the
system zgesdd to compute a 3x4 SVD.
Numpy by default uses a different option of zgesdd, though. To test that,
change the first argument of zgesdd from "Singular Vectors" to "A". Then
run the program again.
I think those are correct instructions.
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/11248#comment:22>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en.
