Abstract:In this paper, we propose a new choice of poles to define reliable rational Krylov methods. These methods are used for approximating function of positive definite matrices. In particular, the fractional power and the fractional resolvent are considered because of their importance in the numerical solution of fractional partial differential equations. The results of the numerical experiments we have carried out on some fractional models confirm that the proposed approach is promising.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: how to effectively calculate matrix functions involving fractional - order partial differential equations (FPDEs), especially fractional - power and fractional - inverse matrix functions. Specifically, the article focuses on using Krylov - type methods to approximate functions of positive definite matrices in the numerical solution of fractional - order partial differential equations.
### Problem Background
Fractional Partial Differential Equations (FPDEs) are widely used in fields such as physics, engineering, and finance. The numerical solutions of these equations usually involve the calculation of certain special matrix functions, such as the fractional - power matrix \(A^{-\alpha/2}\) and the fractional - inverse matrix \((I + \nu A^{\alpha/2})^{-1}\), where \(A\) is a positive definite matrix and \(\alpha\in(1, 2]\).
### Limitations of Existing Methods
Traditional polynomial Krylov methods and shift - and - invert Krylov methods have some limitations when dealing with such problems:
- **Slow convergence speed**: For some problems, the polynomial Krylov method requires more iterations to reach the required accuracy.
- **High computational complexity**: Although the shift - and - invert method can accelerate convergence, it needs to solve linear systems at each iteration, resulting in high computational costs.
### New Method Proposed in the Paper
To solve the above problems, this paper proposes a new pole - selection strategy for defining reliable rational Krylov methods. Specifically:
- **Pole selection**: Use the Gauss - Jacobi rule to select poles, ensuring that these poles are real and simple (no multiple roots), thereby improving the stability and efficiency of the calculation.
- **Scope of application**: This method is applicable to matrices generated by various discretization techniques, including the finite - difference method, the finite - volume method, and the finite - element method.
### Main Contributions
1. **New pole - selection strategy**: By analyzing the rational approximation of objective functions (such as \(z^{-\alpha/2}\) and \((1 + \nu z^{\alpha/2})^{-1}\)), a new pole - selection scheme is proposed.
2. **Theoretical analysis and implementation details**: Prove the properties of the selected poles and give specific implementation steps to ensure the reliability and efficiency of the method.
3. **Numerical experiment verification**: Through numerical experiments on different model problems, the effectiveness and superiority of the new method are verified, especially in terms of the number of iterations and computational time, which are better than existing methods.
### Conclusion
The rational Krylov method proposed in this paper significantly improves the efficiency and accuracy of calculating matrix functions involving fractional - order partial differential equations by optimizing pole selection. This improvement is not only of great theoretical significance but also shows good performance in practical applications.
### Formula Summary
The key formulas involved in the paper are as follows:
- Spectral decomposition form of the fractional Laplace operator:
\[
(-\Delta)^{\alpha/2}u=\sum_{p = 1}^{\infty}\mu_p^{\alpha/2}c_p\phi_p,\quad u=\sum_{p = 1}^{\infty}c_p\phi_p
\]
where \(\{\phi_p\}_{p = 1}^{\infty}\) are the eigenfunctions of the Laplace operator, and \(\{\mu_p\}_{p = 1}^{\infty}\) are the corresponding positive real eigenvalues.
- Definition of the rational Krylov subspace:
\[
Q_k(A_n, v)=[q_{k - 1}(A_n)]^{-1}K_k(A_n, v)
\]
where \(q_{k - 1}(z)=\prod_{j = 1}^{k - 1}(\sigma_j - z)\), \(K_k(A_n, v)=\text{Span}\{v, A_nv,\ldots, A_n^{k - 1}v\}\).
- Selection of new poles: