Abstract:In this paper, we study the improved block splitting (IBS) iteration method and its accelerated variant, the accelerated improved block splitting (AIBS) iteration method, for solving linear systems of equations stemming from the discretization of the complex Helmholtz equation. We conduct a comprehensive convergence analysis and derive optimal iteration parameters aimed at minimizing the spectral radius of the iteration matrix. Through numerical experiments, we validate the efficiency of both iteration methods.
What problem does this paper attempt to address?
The problems that this paper attempts to solve are: propose and study two block - splitting iterative methods for solving the complex symmetric linear systems generated after the discretization of the complex Helmholtz equation, namely the improved block - splitting (IBS) iterative method and its accelerated variant - the accelerated improved block - splitting (AIBS) iterative method.
Specifically, the objectives of the paper include:
1. **Construct and analyze the IBS iterative method**: Improve the existing block - splitting iterative method by introducing the parameter acceleration technique and conduct the convergence analysis.
2. **Determine the optimal iterative parameters**: Derive the optimal parameters that minimize the spectral radius of the iterative matrix.
3. **Verify the effectiveness of the methods**: Verify the efficiency of the IBS and AIBS iterative methods through numerical experiments.
### Background problems
The complex Helmholtz equation describes the two - dimensional field variable \(u(x, y)\), and its form is:
\[
-\nabla^{2}u+\sigma u = f,\quad\nabla^{2}=\frac{\partial^{2}}{\partial x^{2}}+\frac{\partial^{2}}{\partial y^{2}},\quad\sigma\in\mathbb{C}
\]
where \(\sigma\) is a complex constant. This equation is widely used in periodic forced systems, such as electrochemical impedance spectroscopy, unsteady slow viscous flow, etc.
The discretized complex Helmholtz equation will generate a large - scale, sparse, possibly complex - valued and highly indefinite linear system:
\[
Ax\equiv (W + iT)x=b,
\]
where \(W, T\in\mathbb{R}^{n\times n}\) are symmetric positive semi - definite matrices, and at least one of them is positive definite, \(x, b\in\mathbb{C}^{n}\), \(i = \sqrt{-1}\) is the imaginary unit.
### Solutions
To solve the above complex linear system, the author proposes two block - splitting iterative methods:
- **IBS iterative method**: Improve the block - splitting method and conduct the convergence analysis.
- **AIBS iterative method**: Further accelerate the convergence speed based on the pre - processed version of the AGSOR iterative method.
### Numerical experiments
Through a series of numerical experiments, the author verifies the effectiveness and superiority of these two methods, especially their performance in solving complex symmetric linear systems.
In summary, this paper aims to provide efficient iterative methods for solving the linear systems generated after the discretization of the complex Helmholtz equation, and ensure the feasibility and superiority of these methods through theoretical analysis and numerical experiments.