#10795: Fix and upgrade double dense matrix QR decomposition
-------------------------------+--------------------------------------------
    Reporter:  rbeezer         |         Owner:  jason, was 
        Type:  defect          |        Status:  new        
    Priority:  minor           |     Milestone:  sage-4.8   
   Component:  linear algebra  |    Resolution:             
    Keywords:                  |   Work_issues:             
    Upstream:  N/A             |      Reviewer:  Martin Raum
      Author:  Rob Beezer      |        Merged:             
Dependencies:                  |  
-------------------------------+--------------------------------------------
Changes (by jdemeyer):

  * status:  closed => new
  * resolution:  fixed =>
  * merged:  sage-4.8.alpha5 =>


Comment:

 This sometimes gives a doctest failure:
 {{{
 sage -t -long  -force_lib devel/sage/sage/matrix/matrix_double_dense.pyx
 **********************************************************************
 File
 "/mnt/usb1/scratch/buildbot/sage/sage-1/sage_binary/build/sage-4.8.alpha5/devel
 /sage-main/sage/matrix/matrix_double_dense.pyx", line 2243:
     sage: Q.round(6).zero_at(10^-6)
 Expected:
     [-0.359211  0.569326  0.368048  0.641386]
     [ 0.179605 -0.144591  0.925041 -0.301885]
     [ 0.179605  -0.70488  0.077462  0.681825]
     [ 0.898027  0.397625 -0.053281  0.180566]
 Got:
     [-0.359211  0.569326 -0.631992 -0.383955]
     [ 0.179605 -0.144591   -0.6643  0.711014]
     [ 0.179605  -0.70488  -0.39705 -0.559676]
     [ 0.898027  0.397625 -0.040527 -0.183849]
 **********************************************************************
 }}}

 It's not clear why it happens.  Because on the same machine I have one
 Sage build where it works, one where it doesn't with the exact same
 source, but different gcc version.

 Does this code use any external libraries (ATLAS comes to mind)?

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