Ticket #11595: trac_11595-fix_noisy_zero_doctest_errors.reviewer.patch

File trac_11595-fix_noisy_zero_doctest_errors.reviewer.patch, 1.9 KB (added by leif, 10 years ago)

Reviewer patch. Fixes doctest errors due to noisy zeroes in eigenvalues. Apply on top of other patches.

  • sage/matrix/matrix2.pyx

    # HG changeset patch
    # User Leif Leonhardy <not.really@online.de>
    # Date 1317011108 14400
    # Node ID 37fcb898921032f04b0e236d042f7b5757e5975e
    # Parent  10f49e34d981dccab46b3c81dc60a14f0ee2d747
    #11595 Reviewer patch: Fix doctest errors due to noisy zeroes in eigenvalues.
    
    diff --git a/sage/matrix/matrix2.pyx b/sage/matrix/matrix2.pyx
    a b  
    46414641            consult numerical or symbolic matrix classes for other options
    46424642
    46434643            sage: em = A.change_ring(RDF).eigenmatrix_left()
    4644             sage: eigenvalues = em[0]; eigenvalues
     4644            sage: eigenvalues = em[0]; eigenvalues.dense_matrix().zero_at(2e-15)
    46454645            [    13.3484692...                 0                 0]
    46464646            [                0    -1.34846922...                 0]
    4647             [                0                 0 -6.2265089...e-16]
     4647            [                0                 0                 0]
    46484648            sage: eigenvectors = em[1]; eigenvectors
    46494649            [ 0.440242867...  0.567868371...  0.695493875...]
    46504650            [ 0.897878732...  0.278434036... -0.341010658...]
     
    49004900            consult numerical or symbolic matrix classes for other options
    49014901
    49024902            sage: em = B.change_ring(RDF).eigenmatrix_right()
    4903             sage: eigenvalues = em[0]; eigenvalues
     4903            sage: eigenvalues = em[0]; eigenvalues.dense_matrix().zero_at(1e-15)
    49044904            [     13.3484692...                  0                  0]
    49054905            [                 0     -1.34846922...                  0]
    4906             [                 0                  0 -8.86256604...e-16]
     4906            [                 0                  0                  0]
    49074907            sage: eigenvectors = em[1]; eigenvectors
    49084908            [ 0.164763817...  0.799699663...  0.408248290...]
    49094909            [ 0.505774475...  0.104205787... -0.816496580...]