#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):
I posted to the lapack forum about this problem that we're seeing in the
Accelerate framework and they asked for more details. Can people that see
this problem with the SVD provide the following after checking:
1. Download the text to this example program to `test_zgesdd.c`:
http://software.intel.com/sites/products/documentation/hpc/mkl/lapack/mkl_lapack_examples/zgesdd_ex.c.htm
2. Apply the patch found at: http://pastebin.com/cv6ChT4V (raw patch here:
http://pastebin.com/raw.php?i=cv6ChT4V )
3. Compile and run:
{{{
gcc test_zgesdd.c -o test_zgesdd -framework Accelerate
./test_zgesdd
}}}
And post the output here, along with the OSX version number, gcc version
number, maybe the output of `uname -a`,
{{{
sw_vers
xcodebuild -version
gcc --version
uname -a
}}}
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/11248#comment:34>
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.
