Matrix: YCheng/psse0

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

YCheng/psse0 graph
(bipartite graph drawing)


YCheng/psse0

  • 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: 817 KB. Use UFget(1870) or UFget('YCheng/psse0') in MATLAB.
  • download in Matrix Market format, file size: 1 MB.
  • download in Rutherford/Boeing format, file size: 1 MB.

    Matrix properties
    number of rows26,722
    number of columns11,028
    nonzeros102,432
    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%
    typereal
    structurerectangular
    Cholesky candidate?no
    positive definite?no

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

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

    Notes:

    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 COLAMD2,433,000
    nnz(R) for QR, upper bound nnz(U) for LU, with COLAMD84,587

    SVD-based statistics:
    norm(A)200252
    min(svd(A))0.185943
    cond(A)1.07695e+06
    rank(A)11,028
    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

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