#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 kcrisman):

 {{{
 Last login: Thu Apr 21 15:18:29 on ttyp4
 Welcome to Darwin!
 Dasher-03:~ student$ cd Desktop/
 Dasher-03:~/Desktop student$ gcc test_zgesdd.c -o test_zgesdd -framework
 Accelerate
 Dasher-03:~/Desktop student$ ./test_zgesdd
  ZGESDD Example Program Results

  Singular values
   21.76  16.60   3.97

  Left singular vectors (stored columnwise)
  (  0.55,  0.00) (  0.76,  0.00) ( -0.34,  0.00)
  ( -0.04, -0.15) (  0.27, -0.23) (  0.55, -0.74)
  (  0.81,  0.12) ( -0.52, -0.14) (  0.13, -0.11)

  Right singular vectors (stored rowwise)
  (  0.23,  0.21) (  0.37,  0.39) (  0.24,  0.33) ( -0.56, -0.37)
  ( -0.58,  0.40) (  0.11,  0.17) (  0.60, -0.27) (  0.16,  0.06)
  (  0.60,  0.12) ( -0.19,  0.30) (  0.39,  0.20) (  0.45,  0.31)
 }}}
 With the change
 {{{
 Dasher-03:~/Desktop student$ gcc test_zgesdd1.c -o test_zgesdd1 -framework
 Accelerate
 Dasher-03:~/Desktop student$ ./test_zgesdd1
  ZGESDD Example Program Results

  Singular values
   21.76  16.60   3.97

  Left singular vectors (stored columnwise)
  (  0.55,  0.00) (  0.76,  0.00) ( -0.34,  0.00)
  ( -0.04, -0.15) (  0.27, -0.23) (  0.55, -0.74)
  (  0.81,  0.12) ( -0.52, -0.14) (  0.13, -0.11)

  Right singular vectors (stored rowwise)
  (  0.23,  0.21) (  0.37,  0.39) (  0.24,  0.33) ( -0.56, -0.37)
  ( -0.58,  0.40) (  0.11,  0.17) (  0.60, -0.27) (  0.16,  0.06)
  (  0.60,  0.12) ( -0.19,  0.30) (  0.39,  0.20) (  0.45,  0.31)
 }}}

-- 
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/11248#comment:23>
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.

Reply via email to