Abstract:A framework for constructing circulant and block circulant preconditioners (C) for a symmetric linear system Ax=b arising from certain signal and image processing applications is presented in this paper. The proposed scheme does not make explicit use of matrix elements of A. It is ideal for applications in which A only exists in the form of a matrix vector multiplication routine, and in which the process of extracting matrix elements of A is costly. The proposed algorithm takes advantage of the fact that for many linear systems arising from signal or image processing applications, eigenvectors of A can be well represented by a small number of Fourier modes. Therefore, the construction of C can be carried out in the frequency domain by carefully choosing its eigenvalues so that the condition number of C{sup T} AC can be reduced significantly. We illustrate how to construct the spectrum of C in a way such that the smallest eigenvalues of C{sup T} AC overlaps with those of A extremely well while the largest eigenvalues of C{sup T} AC are smaller than those of A by several orders of magnitude. Numerical examples are provided to demonstrate the effectiveness of the preconditioner on accelerating the solution of linear systems arising from image reconstruction application.

A fast algorithm for solving circulant tensor systems

Efficient Algorithm for Finding the Inverse and the Group Inverse of FLSr-circulant Matrix

Tensor Methods for Solving Symmetric ℳ -tensor Systems

An efficient numerical algorithm for solving linear systems with cyclic tridiagonal coefficient matrices

A New Algorithm for Computing the Inverse and Generalized Inverse of the Scaled Factor Circulant Matrix

Three Dimensional Strongly Symmetric Circulant Tensors

Algorithms for Finding the Inverses of Factor Block Circulant Matrices

Sparse Fast Clifford Fourier Transform.

Matrix Factorization for Fast DCT Algorithms

A fast algorithm for the spectral radii of weakly reducible nonnegative tensors: A fast algorithm for the spectral radii of nonnegative tensors

A New Simple Algorithm for Deriving the Winograd 9-Point FFT by Using New Identical Equations for 3 × 3 Circulant and Quasi-Circulant Matrices

Fast Hankel tensor–vector product and its application to exponential data fitting

A survey and some extensions of T. Chan’s preconditioner

Solving Multi-linear Systems with ℳ -Tensors

Circulant-based approximate inverse preconditioners for a class of fractional diffusion equations

A Fast Algorithm for Cosine Transform Based Tensor Singular Value Decomposition

Some New Parallel Fast Fourier Transform Algorithms

Matrix-free Constructions of Circulant and Block Circulant Preconditioners

Developing iterative algorithms to solve Sylvester tensor equations

Tensor Methods for Solving Symmetric \({\mathcal {M}}\)-tensor Systems