Matrix: GHS_psdef/s3dkt3m2

Description: FEM, cylindrical shell, 150x100 tri. mesh, R/t=1000

GHS_psdef/s3dkt3m2 graph
(undirected graph drawing)


GHS_psdef/s3dkt3m2
GHS_psdef/s3dkt3m2 graph

  • Home page of the UF Sparse Matrix Collection
  • Matrix group: GHS_psdef
  • Click here for a description of the GHS_psdef 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: 8 MB. Use UFget(1276) or UFget('GHS_psdef/s3dkt3m2') in MATLAB.
  • download in Matrix Market format, file size: 9 MB.
  • download in Rutherford/Boeing format, file size: 3 MB.

    Matrix properties
    number of rows90,449
    number of columns90,449
    nonzeros3,686,223
    structural full rank?yes
    structural rank90,449
    # of blocks from dmperm1
    # strongly connected comp.1
    explicit zero entries67,238
    nonzero pattern symmetrysymmetric
    numeric value symmetrysymmetric
    typereal
    structuresymmetric
    Cholesky candidate?yes
    positive definite?yes

    authorR. Kouhia
    editorR. Boisvert, R. Pozo, K. Remington, B. Miller, R. Lipman, R. Barrett, J. Dongarra
    date1997
    kindstructural problem
    2D/3D problem?yes

    Additional fieldssize and type
    coordfull 90449-by-3

    Notes:

    %                                                                              
    %FILE  s3dkt3m2.mtx                                                            
    %TITLE Cyl shell R/t=1000 unif 150x100 triang mesh DKT elem with drill rot     
    %KEY   s3dkt3m2                                                                
    %                                                                              
    %                                                                              
    %CONTRIBUTOR Reijo Kouhia (reijo.kouhia@hut.fi)                                
    %                                                                              
    %BEGIN DESCRIPTION                                                             
    % Matrix from a static analysis of a cylindrical shell                         
    % Radius to thickness ratio R/t = 1000                                         
    % Length to radius ratio    R/L = 1                                            
    % One octant discretized with uniform 150 x 100 triangular mesh                
    % element:                                                                     
    % facet-type shell element where the bending part is formulated                
    % using the stabilized MITC theory (stabilization paramater 0.4)               
    % the membrane part includes drilling rotations using                          
    % the Hughes-Brezzi formulation with (regularizing parameter = G/1000,         
    % where G is the shear modulus)                                                
    % full 3-point integration                                                     
    % --------------------------------------------------------------------------   
    % Note:                                                                        
    % The sparsity pattern of the matrix is determined from the element            
    % connectivity data assuming that the element matrix is full.                  
    % Since this case the  material model is linear isotropically elastic          
    % and the FE mesh is  uniform there exist some zeros.                          
    % Since the removal of those zero elements is trivial                          
    % but the reconstruction of the current sparsity                               
    % pattern is impossible from the sparsified structure without any further      
    % knowledge of the element connectivity, the zeros are retained in this file.  
    % ---------------------------------------------------------------------------  
    %END DESCRIPTION                                                               
    %                                                                              
    %                                                                              
    

    Ordering statistics:result
    nnz(chol(P*(A+A'+s*I)*P')) with AMD18,632,125
    Cholesky flop count9.7e+09
    nnz(L+U), no partial pivoting, with AMD37,173,801
    nnz(V) for QR, upper bound nnz(L) for LU, with COLAMD35,265,334
    nnz(R) for QR, upper bound nnz(U) for LU, with COLAMD69,484,595

    Note that all matrix statistics (except nonzero pattern symmetry) exclude the 67238 explicit zero entries.

    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.