We investigate the properties of sparse-matrix ensembles with particular regard for the spectral ergodicity hypothesis, which claims the identity of ensemble and spectral averages of spectral correlators. An apparent violation of the spectral ergodicity is observed. This effect is studied with the a
Assembling sparse matrices in MATLAB
✍ Scribed by Sergio Zlotnik; Pedro Díez
- Publisher
- Wiley (John Wiley & Sons)
- Year
- 2010
- Tongue
- English
- Weight
- 323 KB
- Volume
- 26
- Category
- Article
- ISSN
- 2040-7939
- DOI
- 10.1002/cnm.1174
No coin nor oath required. For personal study only.
✦ Synopsis
Abstract
The assembly of sparse matrices is a key operation in finite element methods. In this study we analyze several factors that may have an influence on the efficiency of the assembly procedure.
Different insertion strategies are compared using two metrics: a Cost function (the number of memory movements) and actual computing time. An improved algorithm implemented in MATLAB is proposed. It reduces both memory operations and computing time for all tested cases.
The efficiency of the assembly process is found to be highly dependent on node and element numbering. The effect of the classic reverse Cuthill–McKee algorithm is, in most cases, positive and reduces computation costs.
Finally, the case where a sparse matrix has to be re‐assembled at each time step is studied. The efficiency of the assembly is improved if the matrix pattern is entirely or partially inherited from previous steps. Copyright © 2008 John Wiley & Sons, Ltd.
📜 SIMILAR VOLUMES
The extracellular matrix (ECM) of the brain contains hyaluronan and proteoglycans, as does the ECM of cartilage. Aggrecan, the major proteoglycan of cartilage, forms large aggregates with hyaluronan, which then associate with the chondrocyte cell surface through an interaction with surface hyalurona
Investigations of scaling and equilibration of general matrices have been traditionally aimed at the effects on the stability and accuracy of LU factorizations-the so-called scaling problem. Notably, Skeel (1979) concludes that no systematic scaling procedure can be concocted for general matrices ex