Sparse approximate matrix-matrix multiplication for density matrix purification with error control

We propose a method for strict error control in sparse approximate matrix-matrix multiplication. The method combines an error bound and a parameter sweep to select an appropriate threshold value. The scheme for error control and the sparse approximate multiplication are implemented using the Chunks and Tasks parallel programming model. We demonstrate the performance of the method in parallel linear scaling electronic structure calculations using density matrix purification with rigorous error control.