Matrix: GHS_psdef/s3dkq4m2

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

GHS_psdef/s3dkq4m2 graph
(undirected graph drawing)


GHS_psdef/s3dkq4m2
GHS_psdef/s3dkq4m2 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: 11 MB. Use UFget(1275) or UFget('GHS_psdef/s3dkq4m2') in MATLAB.
  • download in Matrix Market format, file size: 12 MB.
  • download in Rutherford/Boeing format, file size: 3 MB.

    Matrix properties
    number of rows90,449
    number of columns90,449
    nonzeros4,427,725
    structural full rank?yes
    structural rank90,449
    # of blocks from dmperm1
    # strongly connected comp.1
    explicit zero entries393,166
    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  s3dkq4m2.mtx                                                            
    %TITLE Cyl shell R/t=1000 unif 150x100 quad mesh DKQ elem with drill rot       
    %KEY   s3dkq4m2                                                                
    %                                                                              
    %                                                                              
    %CONTRIBUTOR Reijo Kouhia (reijo.kouhia@hut.fi)                                
    %                                                                              
    %REFERENCE   M. Benzi, R. Kouhia, M.Tuma: An assesment of some                 
    %            preconditioning techniques in shell problems                      
    %            Technical Report LA-UR-97-3892, Los Alamos National Laboratory    
    %                                                                              
    %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 quadrilateral mesh             
    % element:                                                                     
    % facet-type shell element where the bending part is formulated                
    % using the discrete Kirchhoff quadrilateral (Lyons-Crisfield version),        
    % the membrane part includes drilling rotations using                          
    % the Hughes-Brezzi formulation with (regularizing parameter = G/1000,         
    % where G is the shear modulus)                                                
    % full 2x2 Gauss-Legendre 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 AMD27,399,094
    Cholesky flop count1.8e+10
    nnz(L+U), no partial pivoting, with AMD54,707,739
    nnz(V) for QR, upper bound nnz(L) for LU, with COLAMD50,277,848
    nnz(R) for QR, upper bound nnz(U) for LU, with COLAMD98,614,571

    Note that all matrix statistics (except nonzero pattern symmetry) exclude the 393166 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.