Matrix: Quaglino/viscoplastic2
Description: FEM discretization of a viscoplastic collision problem, Alessio Quaglino
![]() |
![]() |
| (bipartite graph drawing) | (graph drawing of A+A') |
![]() |
| Matrix properties | |
| number of rows | 32,769 |
| number of columns | 32,769 |
| nonzeros | 381,326 |
| structural full rank? | yes |
| structural rank | 32,769 |
| # of blocks from dmperm | 1 |
| # strongly connected comp. | 1 |
| explicit zero entries | 0 |
| nonzero pattern symmetry | 57% |
| numeric value symmetry | 0% |
| type | real |
| structure | unsymmetric |
| Cholesky candidate? | no |
| positive definite? | no |
| author | A. Quaglino |
| editor | T. Davis |
| date | 2007 |
| kind | materials problem |
| 2D/3D problem? | yes |
| Additional fields | size and type |
| b | full 32769-by-1 |
| C | cell 7-by-1 |
Notes:
The matrix is in the form [A11 A12 ; A21 A22] where A11 and A22 are diagonal.
Originally, the matrices in this set were poorly scaled, but this was resolved
by a scale factor of the form [A11 A12*e ; A21/e A4] where the scalar e is
of magnitude 1e2 but can be 1e6 or 1e7 for a stiff material. The Problem.A
matrix is the properly scaled problem. The Problem.aux.C{1:7} matrices have
been "unscaled" with a factor e = 10.^(-(1:7)), to give a sequence of matrices
that are well scaled to poorly scaled, and thus well conditioned (C{1}) to
poorly conditioned (C{7}). This mimics the original poorly scaled and ill-
conditioned problem, and may be of interest for those developing algorithms
for automatic scaling. From a FEM discretization of a viscoplastic collision
problem, Alessio Quaglino.
| Ordering statistics: | result |
| nnz(chol(P*(A+A'+s*I)*P')) with AMD | 789,046 |
| Cholesky flop count | 8.7e+07 |
| nnz(L+U), no partial pivoting, with AMD | 1,545,323 |
| nnz(V) for QR, upper bound nnz(L) for LU, with COLAMD | 24,894,893 |
| nnz(R) for QR, upper bound nnz(U) for LU, with COLAMD | 49,751,382 |
| SVD-based statistics: | |
| norm(A) | 20.8475 |
| min(svd(A)) | 9.2434e-05 |
| cond(A) | 225540 |
| rank(A) | 32,769 |
| sprank(A)-rank(A) | 0 |
| null space dimension | 0 |
| full numerical rank? | yes |
| singular values (MAT file): | click here |
| SVD method used: | s = svd (full (A)) |
| status: | ok |

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.