Systematic sparse matrix error control for linear scaling electronic structure calculations
2005 (English)In: Journal of Computational Chemistry, ISSN 0192-8651, E-ISSN 1096-987X, Vol. 26, no 15, 1628-1637 p.Article in journal (Refereed) Published
Efficient truncation criteria used in multiatom blocked sparse matrix operations for ab initio calculations are proposed. As system size increases, so does the need to stay on top of errors and still achieve high performance. A variant of a blocked sparse matrix algebra to achieve strict error control with good performance is proposed. The presented idea is that the condition to drop a certain submatrix should depend not only on the magnitude of that particular submatrix, but also on which other submatrices that are dropped. The decision to remove a certain submatrix is based on the contribution the removal would cause to the error in the chosen norm. We study the effect of an accumulated truncation error in iterative algorithms like trace correcting density matrix purification. One way to reduce the initial exponential growth of this error is presented. The presented error control for a sparse blocked matrix toolbox allows for achieving optimal performance by performing only necessary operations needed to maintain the requested level of accuracy.
Place, publisher, year, edition, pages
2005. Vol. 26, no 15, 1628-1637 p.
matrix error control, linear scaling, density purification
IdentifiersURN: urn:nbn:se:kth:diva-6544DOI: 10.1002/jcc.20315ISI: 000232570300008ScopusID: 2-s2.0-27844477884OAI: oai:DiVA.org:kth-6544DiVA: diva2:11286
QC 201011232006-12-112006-12-112010-11-23Bibliographically approved