#17124: Matrix constructor broken for symmetric function rings and 1-column
matrices
------------------------------------------------+--------------------------
       Reporter:  darij                         |        Owner:
           Type:  defect                        |       Status:
       Priority:  major                         |  needs_review
      Component:  linear algebra                |    Milestone:  sage-6.4
       Keywords:  matrix, symmetric functions,  |   Resolution:
        Authors:                                |    Merged in:
Report Upstream:  N/A                           |    Reviewers:
         Branch:                                |  Work issues:
   Dependencies:                                |       Commit:
                                                |     Stopgaps:
------------------------------------------------+--------------------------

Comment (by darij):

 Fixed in branch public/linalg/uglyduck. I am not happy with the way it
 went, but at least I can't reproduce the bug anymore. There is now some
 precision issues:

 {{{
 sage -t src/sage/matrix/matrix2.pyx
 **********************************************************************
 File "src/sage/matrix/matrix2.pyx", line 5209, in
 sage.matrix.matrix2.Matrix.diagonal.eigenspaces_right
 Failed example:
     eigenvalues = em[0]; eigenvalues.dense_matrix().zero_at(1e-15)
 Expected:
     [ 13.348469228349...                0.0                 0.0]
     [                0.0 -1.348469228349...                 0.0]
     [                0.0                0.0                 0.0]
 Got:
     [13.348469228349522                0.0                0.0]
     [               0.0 -1.348469228349534                0.0]
     [               0.0                0.0                0.0]
 **********************************************************************
 File "src/sage/matrix/matrix2.pyx", line 5689, in
 sage.matrix.matrix2.Matrix.diagonal.eigenmatrix_right
 Failed example:
     evalues = em[0]; evalues.dense_matrix().zero_at(2e-15)
 Expected:
     [ 13.348469228349...                0.0                 0.0]
     [                0.0 -1.348469228349...                 0.0]
     [                0.0                0.0                 0.0]
 Got:
     [13.348469228349522                0.0                0.0]
     [               0.0 -1.348469228349534                0.0]
     [               0.0                0.0                0.0]
 **********************************************************************
 File "src/sage/matrix/matrix2.pyx", line 12207, in
 sage.matrix.matrix2.Matrix.diagonal.norm
 Failed example:
     A.norm()
 Expected:
     15.0
 Got:
     14.999999999999998
 **********************************************************************
 }}}

--
Ticket URL: <http://trac.sagemath.org/ticket/17124#comment:3>
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 unsubscribe from this group and stop receiving emails from it, send an email 
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/sage-trac.
For more options, visit https://groups.google.com/d/optout.

Reply via email to