• UF Sparse Matrix Collection
  • Matrix group: Cylshell
  • Click here for a description of the Cylshell group.
  • Click here for a list of all matrices
  • Click here for a list of all matrix groups


  • Matrix: Cylshell/s3rmt3m3
  • Description: FEM, cylindrical shell, graded tri. mesh w/ 1666 triangles. , R/t=1000
  • download as a MATLAB mat-file, file size: 1 MB. Use UFget(1611) or UFget('Cylshell/s3rmt3m3') in MATLAB.
  • download in Matrix Market format
  • download in Rutherford/Boeing format

    Cylshell/s3rmt3m3

    Cylshell/s3rmt3m3 graph

    Matrix properties
    number of rows5,357
    number of columns5,357
    nonzeros207,123
    structural full rank?yes
    structural rank5,357
    # of blocks from dmperm1
    # strongly connected comp.1
    entries not in dmperm blocks0
    explicit zero entries572
    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

    Notes:

    %                                                                              
    %FILE  s3rmt3m3.mtx                                                            
    %TITLE Cyl shell R/t=1000 grad trian mesh 1666 stab MITC3 elem with drill rot  
    %KEY   s3rmt3m3                                                                
    %                                                                              
    %                                                                              
    %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 graded triangular mesh (1666 elements)           
    % 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          
    % 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:AMD METIS
    nnz(chol(P*(A+A'+s*I)*P'))438,361 490,790
    Cholesky flop count5.7e+07 6.9e+07
    nnz(L+U), no partial pivoting871,365 976,223
    nnz(V) for QR, upper bound nnz(L) for LU566,229 644,967
    nnz(R) for QR, upper bound nnz(U) for LU1,115,039 1,199,942

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

    Maintained by Tim Davis, last updated 04-May-2008.
    Matrix pictures by cspy, a MATLAB function in the CSparse package.
    Matrix graphs by Yifan Hu, AT&T Labs Visualization Group.