Matrix: Castrillon/denormal
Description: denormals in A+alpha*I as alpha > 0 varies, J. Castrillon, Teledyne
![]() |
| (undirected graph drawing) |
![]() |
| Matrix properties | |
| number of rows | 89,400 |
| number of columns | 89,400 |
| nonzeros | 1,156,224 |
| structural full rank? | yes |
| structural rank | 89,400 |
| # of blocks from dmperm | 1 |
| # strongly connected comp. | 1 |
| explicit zero entries | 0 |
| nonzero pattern symmetry | symmetric |
| numeric value symmetry | symmetric |
| type | real |
| structure | symmetric |
| Cholesky candidate? | yes |
| positive definite? | yes |
| author | J. Castrillon |
| editor | T. Davis |
| date | 2008 |
| kind | counter-example problem |
| 2D/3D problem? | no |
Notes:
The problem is x=(A+alpha*I)\b as alpha varies from 0 to 100. However, as this occurs, performance suffers because of the many denormals that appear below the diagonal in L.
| Ordering statistics: | result |
| nnz(chol(P*(A+A'+s*I)*P')) with AMD | 8,231,641 |
| Cholesky flop count | 2.5e+09 |
| nnz(L+U), no partial pivoting, with AMD | 16,373,882 |
| nnz(V) for QR, upper bound nnz(L) for LU, with COLAMD | 16,686,436 |
| nnz(R) for QR, upper bound nnz(U) for LU, with COLAMD | 30,789,903 |
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.