An efficient parallel implementation of the MSPAI preconditioner |
| |
Authors: | T. Huckle A. Kallischko A. Roy M. Sedlacek T. Weinzierl |
| |
Affiliation: | 1. Department of Statistics, Applied Mathematics and Computing, São Paulo State University – UNESP, Rio Claro/SP, 13506-900, Brazil;2. School of Computer Science, Bangor University, Bangor, LL57 1UT, United Kingdom;3. Instituto de Ciência e Tecnologia, Universidade Federal de São Paulo – UNIFESP, São José dos Campos/SP, 12247-014, Brazil |
| |
Abstract: | We present an efficient implementation of the Modified SParse Approximate Inverse (MSPAI) preconditioner. MSPAI generalizes the class of preconditioners based on Frobenius norm minimizations, the class of modified preconditioners such as MILU, as well as interface probing techniques in domain decomposition: it adds probing constraints to the basic SPAI formulation, and one can thus optimize the preconditioner relative to certain subspaces. We demonstrate MSPAI’s qualities for iterative regularization problems arising from image deblurring.Such applications demand for a fast and parallel preconditioner realization. We present such an implementation introducing two new optimization techniques: First, we avoid redundant calculations using a dictionary. Second, our implementation reduces the runtime spent on the most demanding numerical parts as the code switches to sparse QR decomposition methods wherever profitable. The optimized code runs in parallel with a dynamic load balancing. |
| |
Keywords: | |
本文献已被 ScienceDirect 等数据库收录! |
|