Matrix: Quaglino/viscoplastic2

Description: FEM discretization of a viscoplastic collision problem, Alessio Quaglino

Quaglino/viscoplastic2 graph Quaglino/viscoplastic2 graph
(bipartite graph drawing) (graph drawing of A+A')


Quaglino/viscoplastic2

  • Home page of the UF Sparse Matrix Collection
  • Matrix group: Quaglino
  • Click here for a description of the Quaglino 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: 12 MB. Use UFget(1869) or UFget('Quaglino/viscoplastic2') in MATLAB.
  • download in Matrix Market format, file size: 13 MB.
  • download in Rutherford/Boeing format, file size: 12 MB.

    Matrix properties
    number of rows32,769
    number of columns32,769
    nonzeros381,326
    structural full rank?yes
    structural rank32,769
    # of blocks from dmperm1
    # strongly connected comp.1
    explicit zero entries0
    nonzero pattern symmetry 57%
    numeric value symmetry 0%
    typereal
    structureunsymmetric
    Cholesky candidate?no
    positive definite?no

    authorA. Quaglino
    editorT. Davis
    date2007
    kindmaterials problem
    2D/3D problem?yes

    Additional fieldssize and type
    bfull 32769-by-1
    Ccell 7-by-1

    Notes:

    The matrix is in the form [A11 A12 ; A21 A22] where A11 and A22 are diagonal. 
    Originally, the matrices in this set were poorly scaled, but this was resolved
    by a scale factor of the form [A11 A12*e ; A21/e A4] where the scalar e is    
    of magnitude 1e2 but can be 1e6 or 1e7 for a stiff material.  The Problem.A   
    matrix is the properly scaled problem.  The Problem.aux.C{1:7} matrices have  
    been "unscaled" with a factor e = 10.^(-(1:7)), to give a sequence of matrices
    that are well scaled to poorly scaled, and thus well conditioned (C{1}) to    
    poorly conditioned (C{7}).  This mimics the original poorly scaled and ill-   
    conditioned problem, and may be of interest for those developing algorithms   
    for automatic scaling.  From a FEM discretization of a viscoplastic collision 
    problem, Alessio Quaglino.                                                    
    

    Ordering statistics:result
    nnz(chol(P*(A+A'+s*I)*P')) with AMD789,046
    Cholesky flop count8.7e+07
    nnz(L+U), no partial pivoting, with AMD1,545,323
    nnz(V) for QR, upper bound nnz(L) for LU, with COLAMD24,894,893
    nnz(R) for QR, upper bound nnz(U) for LU, with COLAMD49,751,382

    SVD-based statistics:
    norm(A)20.8475
    min(svd(A))9.2434e-05
    cond(A)225540
    rank(A)32,769
    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

    Quaglino/viscoplastic2 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.