Matrix: JGD_Trefethen/Trefethen_500

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

JGD_Trefethen/Trefethen_500 graph
(undirected graph drawing)


JGD_Trefethen/Trefethen_500

  • 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: 15 KB. Use UFget(2209) or UFget('JGD_Trefethen/Trefethen_500') in MATLAB.
  • download in Matrix Market format, file size: 13 KB.
  • download in Rutherford/Boeing format, file size: 11 KB.

    Matrix properties
    number of rows500
    number of columns500
    nonzeros8,478
    structural full rank?yes
    structural rank500
    # 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_500.sms              
    

    Ordering statistics:result
    nnz(chol(P*(A+A'+s*I)*P')) with AMD55,480
    Cholesky flop count1.0e+07
    nnz(L+U), no partial pivoting, with AMD110,460
    nnz(V) for QR, upper bound nnz(L) for LU, with COLAMD90,513
    nnz(R) for QR, upper bound nnz(U) for LU, with COLAMD109,344

    SVD-based statistics:
    norm(A)3571.25
    min(svd(A))1.12105
    cond(A)3185.64
    rank(A)500
    sprank(A)-rank(A)0
    null space dimension0
    full numerical rank?yes

    singular values (MAT file):click here
    SVD method used:s = svd (full (A)) ;
    status:ok

    JGD_Trefethen/Trefethen_500 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.