Matrix: Cylshell/s1rmq4m1
Description: FEM, cylindrical shell, 30x30 quad. mesh, stabilized MITC4 elements, R/t=10
| Matrix properties | |
| number of rows | 5,489 |
| number of columns | 5,489 |
| nonzeros | 262,411 |
| structural full rank? | yes |
| structural rank | 5,489 |
| # of blocks from dmperm | 1 |
| # strongly connected comp. | 1 |
| entries not in dmperm blocks | 0 |
| explicit zero entries | 18,700 |
| 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 |
Notes:
% %FILE s1rmq4m1.mtx %TITLE Cyl shell R/t = 10 unif 30x30 quad mesh stab MITC4 elem with drill rot %KEY s1rmq4m1 % % %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 = 10 % Length to radius ratio R/L = 1 % One octant discretized with uniform 30 x 30 quadrilateral 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 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: | AMD | METIS |
| nnz(chol(P*(A+A'+s*I)*P')) | 846,206 | 651,100 |
| Cholesky flop count | 1.9e+08 | 1.0e+08 |
| nnz(L+U), no partial pivoting | 1,686,923 | 1,296,711 |
| nnz(V) for QR, upper bound nnz(L) for LU | 875,290 | 867,059 |
| nnz(R) for QR, upper bound nnz(U) for LU | 1,735,331 | 1,621,049 |
Note that all matrix statistics (except nonzero pattern symmetry) exclude the 18700 explicit zero entries.
Maintained by Tim Davis, last updated 30-Sep-2008.
Matrix pictures by cspy, a MATLAB function in the CSparse package.
Matrix graphs by Yifan Hu, AT&T Labs Visualization Group.