Matrix: JGD_Trefethen/Trefethen_20000

Description: Diagonal matrices with primes, Nick Trefethen, Oxford Univ.

JGD_Trefethen/Trefethen_20000 graph
(undirected graph drawing)


JGD_Trefethen/Trefethen_20000

  • Home page of the UF Sparse Matrix Collection
  • Matrix group: JGD_Trefethen
  • Click here for a description of the JGD_Trefethen group.
  • Click here for a list of all matrices
  • Click here for a list of all matrix groups
  • download as a MATLAB mat-file, file size: 1017 KB. Use UFget(2213) or UFget('JGD_Trefethen/Trefethen_20000') in MATLAB.
  • download in Matrix Market format, file size: 883 KB.
  • download in Rutherford/Boeing format, file size: 703 KB.

    Matrix properties
    number of rows20,000
    number of columns20,000
    nonzeros554,466
    structural full rank?yes
    structural rank20,000
    # of blocks from dmperm1
    # strongly connected comp.1
    explicit zero entries0
    nonzero pattern symmetrysymmetric
    numeric value symmetrysymmetric
    typeinteger
    structuresymmetric
    Cholesky candidate?yes
    positive definite?yes

    authorN. Trefethen
    editorJ.-G. Dumas
    date2008
    kindcombinatorial problem
    2D/3D problem?no

    Notes:

    Diagonal matrices with primes, Nick Trefethen, Oxford Univ.          
    From Jean-Guillaume Dumas' Sparse Integer Matrix Collection,         
    http://ljk.imag.fr/membres/Jean-Guillaume.Dumas/simc.html            
                                                                         
    Problem 7 of the Hundred-dollar, Hundred-digit Challenge Problems,   
    SIAM News, vol 35, no. 1.                                            
                                                                         
    7. Let A be the 20,000 x 20,000 matrix whose entries are zero        
    everywhere except for the primes 2, 3, 5, 7, . . . , 224737 along the
    main diagonal and the number 1 in all the positions A(i,j) with      
    |i-j| = 1,2,4,8, . . . ,16384.  What is the (1,1) entry of inv(A)?   
                                                                         
    http://www.siam.org/news/news.php?id=388                             
                                                                         
    Filename in JGD collection: Trefethen/trefethen_20000.sms            
    

    Ordering statistics:result
    nnz(chol(P*(A+A'+s*I)*P')) with AMD86,761,641
    Cholesky flop count7.4e+11
    nnz(L+U), no partial pivoting, with AMD173,503,282
    nnz(V) for QR, upper bound nnz(L) for LU, with COLAMD171,632,676
    nnz(R) for QR, upper bound nnz(U) for LU, with COLAMD187,325,442

    SVD-based statistics:
    norm(A)224737
    min(svd(A))1.12055
    cond(A)200559
    rank(A)20,000
    sprank(A)-rank(A)0
    null space dimension0
    full numerical rank?yes

    singular values (MAT file):click here
    SVD method used:s = svd (full (R)) ; where [~,R,E] = spqr (A) with droptol of zero
    status:ok

    JGD_Trefethen/Trefethen_20000 svd

    For a description of the statistics displayed above, click here.

    Maintained by Tim Davis, last updated 12-Mar-2014.
    Matrix pictures by cspy, a MATLAB function in the CSparse package.
    Matrix graphs by Yifan Hu, AT&T Labs Visualization Group.