762 research outputs found
A multi-level preconditioned Krylov method for the efficient solution of algebraic tomographic reconstruction problems
Classical iterative methods for tomographic reconstruction include the class
of Algebraic Reconstruction Techniques (ART). Convergence of these stationary
linear iterative methods is however notably slow. In this paper we propose the
use of Krylov solvers for tomographic linear inversion problems. These advanced
iterative methods feature fast convergence at the expense of a higher
computational cost per iteration, causing them to be generally uncompetitive
without the inclusion of a suitable preconditioner. Combining elements from
standard multigrid (MG) solvers and the theory of wavelets, a novel
wavelet-based multi-level (WMG) preconditioner is introduced, which is shown to
significantly speed-up Krylov convergence. The performance of the
WMG-preconditioned Krylov method is analyzed through a spectral analysis, and
the approach is compared to existing methods like the classical Simultaneous
Iterative Reconstruction Technique (SIRT) and unpreconditioned Krylov methods
on a 2D tomographic benchmark problem. Numerical experiments are promising,
showing the method to be competitive with the classical Algebraic
Reconstruction Techniques in terms of convergence speed and overall performance
(CPU time) as well as precision of the reconstruction.Comment: Journal of Computational and Applied Mathematics (2014), 26 pages, 13
figures, 3 table
Graph partitioning using matrix values for preconditioning symmetric positive definite systems
Prior to the parallel solution of a large linear system, it is required to
perform a partitioning of its equations/unknowns. Standard partitioning
algorithms are designed using the considerations of the efficiency of the
parallel matrix-vector multiplication, and typically disregard the information
on the coefficients of the matrix. This information, however, may have a
significant impact on the quality of the preconditioning procedure used within
the chosen iterative scheme. In the present paper, we suggest a spectral
partitioning algorithm, which takes into account the information on the matrix
coefficients and constructs partitions with respect to the objective of
enhancing the quality of the nonoverlapping additive Schwarz (block Jacobi)
preconditioning for symmetric positive definite linear systems. For a set of
test problems with large variations in magnitudes of matrix coefficients, our
numerical experiments demonstrate a noticeable improvement in the convergence
of the resulting solution scheme when using the new partitioning approach
- …