Matrix: YCheng/psse2

Description: Power system state simulation matrix, Yunzhi Cheng, UT Arlington

YCheng/psse2 graph
(bipartite graph drawing)


  • Home page of the UF Sparse Matrix Collection
  • Matrix group: YCheng
  • Click here for a description of the YCheng 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: 894 KB. Use UFget(1872) or UFget('YCheng/psse2') in MATLAB.
  • download in Matrix Market format, file size: 1 MB.
  • download in Rutherford/Boeing format, file size: 1 MB.

    Matrix properties
    number of rows28,634
    number of columns11,028
    structural full rank?yes
    structural rank11,028
    # of blocks from dmperm1
    # strongly connected comp.1
    explicit zero entries0
    nonzero pattern symmetry 0%
    numeric value symmetry 0%
    Cholesky candidate?no
    positive definite?no

    authorY. Cheng
    editorT. Davis
    kindpower network problem
    2D/3D problem?no

    Additional fieldssize and type
    bfull 28634-by-1
    guessfull 11028-by-1


    Power system state simulation matrix from Yunzhi Cheng, UT Arlington.
    In MATLAB, the solution to x=A\b is desired, but this can be slow in 
    MATLAB 7.3 because of the speed of sparse QR as compared to sparse   
    Cholesky.  Using x=(A'*A)\(A'*b) is faster, but of course yields     
    slightly less accurate (but still acceptable) results.  Note that an 
    initial guess to the solution is provided, for use by an iterative   
    method.  However, sparse Cholesky with an AMD ordering is very fast  
    for this matrix and thus iterative methods are unlikely to be        
    competitive.  In MATLAB 7.3 on a 3.2 Ghz Pentium 4 desktop,          
    x=(A'*A)\(A'*b) takes 0.07 seconds.                                  

    Ordering statistics:result
    nnz(V) for QR, upper bound nnz(L) for LU, with COLAMD3,138,002
    nnz(R) for QR, upper bound nnz(U) for LU, with COLAMD109,583

    SVD-based statistics:
    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

    YCheng/psse2 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.