Abstract:In this work, we propose an absolute value block $\alpha$-circulant preconditioner for the minimal residual (MINRES) method to solve an all-at-once system arising from the discretization of wave equations. Motivated by the absolute value block circulant preconditioner proposed in [E. McDonald, J. Pestana, and A. Wathen. SIAM J. Sci. Comput., 40(2):A1012-A1033, 2018], we propose an absolute value version of the block $\alpha$-circulant preconditioner. Since the original block $\alpha$-circulant preconditioner is non-Hermitian in general, it cannot be directly used as a preconditioner for MINRES. Our proposed preconditioner is the first Hermitian positive definite variant of the block $\alpha$-circulant preconditioner for the concerned wave equations, which fills the gap between block $\alpha$-circulant preconditioning and the field of preconditioned MINRES solver. The matrix-vector multiplication of the preconditioner can be fast implemented via fast Fourier transforms. Theoretically, we show that for a properly chosen $\alpha$ the MINRES solver with the proposed preconditioner achieves a linear convergence rate independent of the matrix size. To the best of our knowledge, this is the first attempt to generalize the original absolute value block circulant preconditioner in the aspects of both theory and performance the concerned problem. Numerical experiments are given to support the effectiveness of our preconditioner, showing that the expected optimal convergence can be achieved.
What problem does this paper attempt to address?
### What problem does this paper attempt to solve?
This paper aims to solve the problem of solving the all - at - once linear systems formed after the discretization of the wave equation, especially how to use the preconditioned minimum residual method (MINRES) to solve these systems efficiently. Specifically:
1. **Background and motivation**:
- The wave equation is an important partial differential equation (PDE) for describing wave phenomena and has wide applications in physics, engineering and other fields.
- Traditional methods for solving the wave equation are usually solved sequentially, but this method is less efficient in parallel computing.
- Parallel - in - time solvers (PinT solvers) can solve in parallel over multiple time steps, improving computational efficiency.
2. **Existing challenges**:
- Direct application of the existing block α - circulant preconditioning methods cannot be directly used for the MINRES solver, because these preconditioning matrices are usually non - Hermitian positive definite (non - Hermitian positive definite, non - HPD), while MINRES requires that the preconditioning matrix must be Hermitian positive definite.
- The existing absolute value block circulant (ABC) preconditioning method, although suitable for parabolic equations, performs poorly for wave equations because its eigenvalue distribution is not ideal, resulting in a slow convergence rate.
3. **Proposed new method**:
- The paper proposes a new absolute value block α - circulant (ABAC) preconditioning method, which is Hermitian positive definite and can be effectively applied to the MINRES solver.
- The new ABAC pre - conditioner, by introducing the parameter α, makes the pre - conditioned matrix have better spectral properties, thereby accelerating the convergence rate of MINRES.
- It is theoretically proven that for an appropriately chosen α, the MINRES solver can achieve a linear convergence rate independent of the matrix size when using the ABAC pre - conditioner.
4. **Theoretical contributions**:
- It proposes the first non - trivial symmetric positive definite version of the block α - circulant matrix preconditioning method, filling the gap in this field.
- New pre - conditioners are constructed through techniques such as matrix square roots, providing new ideas for solving other more complex PDEs.
5. **Numerical experiments**:
- Numerical experiments verify the effectiveness of the newly proposed ABAC pre - conditioner, indicating that it is superior to the existing ABC pre - conditioner in solving the wave equation.
### Summary
The main goal of this paper is to develop a new preconditioning method (ABAC) that can effectively accelerate the solution of the all - at - once linear system of the wave equation and ensure that the convergence rate of the MINRES solver is independent of the matrix size. This not only solves the deficiencies of existing methods in solving the wave equation but also provides new research directions and technical means for the field of parallel computing.