Matrix: JGD_Trefethen/Trefethen_200b

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

JGD_Trefethen/Trefethen_200b graph
(undirected graph drawing)


  • 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: 6 KB. Use UFget(2206) or UFget('JGD_Trefethen/Trefethen_200b') in MATLAB.
  • download in Matrix Market format, file size: 5 KB.
  • download in Rutherford/Boeing format, file size: 4 KB.

    Matrix properties
    number of rows199
    number of columns199
    structural full rank?yes
    structural rank199
    # of blocks from dmperm1
    # strongly connected comp.1
    explicit zero entries0
    nonzero pattern symmetrysymmetric
    numeric value symmetrysymmetric
    Cholesky candidate?yes
    positive definite?yes

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


    Diagonal matrices with primes, Nick Trefethen, Oxford Univ.          
    From Jean-Guillaume Dumas' Sparse Integer Matrix Collection,                
    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)?   
    Filename in JGD collection: Trefethen/trefethen_200__199_minor.sms   

    Ordering statistics:result
    nnz(chol(P*(A+A'+s*I)*P')) with AMD8,720
    Cholesky flop count6.0e+05
    nnz(L+U), no partial pivoting, with AMD17,241
    nnz(V) for QR, upper bound nnz(L) for LU, with COLAMD13,770
    nnz(R) for QR, upper bound nnz(U) for LU, with COLAMD17,885

    SVD-based statistics:
    null space dimension0
    full numerical rank?yes

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

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