Matrix: GHS_psdef/s3dkq4m2
Description: FEM, cylindrical shell, 150x100 quad. mesh, R/t=1000
(undirected graph drawing) |
Matrix properties | |
number of rows | 90,449 |
number of columns | 90,449 |
nonzeros | 4,427,725 |
structural full rank? | yes |
structural rank | 90,449 |
# of blocks from dmperm | 1 |
# strongly connected comp. | 1 |
explicit zero entries | 393,166 |
nonzero pattern symmetry | symmetric |
numeric value symmetry | symmetric |
type | real |
structure | symmetric |
Cholesky candidate? | yes |
positive definite? | yes |
author | R. Kouhia |
editor | R. Boisvert, R. Pozo, K. Remington, B. Miller, R. Lipman, R. Barrett, J. Dongarra |
date | 1997 |
kind | structural problem |
2D/3D problem? | yes |
Additional fields | size and type |
coord | full 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 AMD | 27,399,094 |
Cholesky flop count | 1.8e+10 |
nnz(L+U), no partial pivoting, with AMD | 54,707,739 |
nnz(V) for QR, upper bound nnz(L) for LU, with COLAMD | 50,277,848 |
nnz(R) for QR, upper bound nnz(U) for LU, with COLAMD | 98,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.