Ticket #10791: trac_10791-gram-schmidt-qqbar-doctests.patch

File trac_10791-gram-schmidt-qqbar-doctests.patch, 3.3 KB (added by rbeezer, 10 years ago)
  • sage/matrix/matrix2.pyx

    # HG changeset patch
    # User Rob Beezer <beezer@ups.edu>
    # Date 1323061269 28800
    # Node ID 7b10e57e554dd7eac3f96575811c173af236e143
    # Parent  860dc1e11808d23e8e2468587b5ed3eb7db1d3bd
    10791: fix Gram-Schmidt doctests over QQbar
    
    diff --git a/sage/matrix/matrix2.pyx b/sage/matrix/matrix2.pyx
    a b  
    80698069            [0 0 0]
    80708070            [0 0 0]
    80718071            [0 0 0]
    8072             sage: G*G.transpose()
    8073             [      1 0.?e-37 0.?e-35]
    8074             [0.?e-37       1 0.?e-35]
    8075             [0.?e-35 0.?e-35       1]
     8072            sage: (G*G.transpose()-identity_matrix(3)).norm() < 10^-10
     8073            True
    80768074            sage: G.row_space() == A.row_space()
    80778075            True
    80788076
     
    80848082            ...                      [-4*I, -2*I + 17,       0,  9*I + 1],
    80858083            ...                      [   1,  -2*I - 6, -I + 11, -5*I + 1]])
    80868084            sage: G, M = A.gram_schmidt(orthonormal=True)
    8087             sage: G
    8088             [                         -0.3849001794597505?  -0.1924500897298753? - 0.1924500897298753?*I   0.3849001794597505? + 0.7698003589195010?*I                          -0.1924500897298753?]
    8089             [-0.06165497274852388? - 0.1387236886841787?*I  0.8188551068163327? - 0.10018933071635130?*I  0.2003786614327026? + 0.05394810115495839?*I  0.02119389688230508? + 0.5028733714801478?*I]
    8090             [  0.3842387256410419? - 0.5694103019142261?*I  0.1416892863096208? - 0.06139779741542298?*I  0.4633778333528464? - 0.01285039016180503?*I  0.02658516588219101? - 0.5373044261814995?*I]
    8091             sage: M
    8092             [                        5.196152422706632?                                          0                                          0]
    8093             [-3.079201435678004? + 3.464101615137755?*I             19.22286447225071? + 0.?e-37*I                                          0]
    8094             [  4.426352063787131? - 8.66025403784439?*I -5.117362738127481? - 3.502773139275513?*I             7.480012456446966? + 0.?e-35*I]
    8095             sage: M*G-A
    8096             [            0.?e-37 0.?e-37 + 0.?e-37*I 0.?e-37 + 0.?e-37*I             0.?e-37]
    8097             [0.?e-36 + 0.?e-36*I 0.?e-35 + 0.?e-36*I 0.?e-36 + 0.?e-36*I 0.?e-36 + 0.?e-36*I]
    8098             [0.?e-35 + 0.?e-35*I 0.?e-35 + 0.?e-35*I 0.?e-35 + 0.?e-35*I 0.?e-35 + 0.?e-35*I]
    8099             sage: G*G.conjugate().transpose()
    8100             [1.000000000000000? + 0.?e-37*I            0.?e-37 + 0.?e-37*I            0.?e-36 + 0.?e-36*I]
    8101             [           0.?e-37 + 0.?e-37*I 1.000000000000000? + 0.?e-37*I            0.?e-36 + 0.?e-36*I]
    8102             [           0.?e-36 + 0.?e-36*I            0.?e-36 + 0.?e-36*I 1.000000000000000? + 0.?e-35*I]
     8085            sage: (M*G-A).norm() < 10^-10
     8086            True
     8087            sage: id3 = G*G.conjugate().transpose()
     8088            sage: (id3 - identity_matrix(3)).norm() < 10^-10
     8089            True
    81038090            sage: G.row_space() == A.row_space()
    81048091            True
    81058092
     
    81278114            [0 0 0 0]
    81288115            [0 0 0 0]
    81298116            [0 0 0 0]
    8130             sage: G*G.transpose()
    8131             [      1 0.?e-35]
    8132             [0.?e-35       1]
     8117            sage: (G*G.transpose()-identity_matrix(2)).norm() < 10^-10
     8118            True
    81338119            sage: G.row_space() == A.row_space()
    81348120            True
    81358121