A Preconditioned MINRES Method for Optimal Control of Wave Equations and its Asymptotic Spectral Distribution Theory
In this work, we propose a novel preconditioned Krylov subspace method for solvingan optimal control problem of wave equations, after explicitly identifying the asymptotic spectraldistribution of the involved sequence of linear coefficient matrices from the optimal control prob-lem. Namely, we first...
Uložené v:
| Vydané v: | SIAM journal on matrix analysis and applications Ročník 44; číslo 4; s. 1477 - 1509 |
|---|---|
| Hlavní autori: | , , |
| Médium: | Journal Article |
| Jazyk: | English |
| Vydavateľské údaje: |
01.01.2023
|
| Predmet: | |
| ISSN: | 0895-4798, 1095-7162, 1095-7162 |
| On-line prístup: | Získať plný text |
| Tagy: |
Pridať tag
Žiadne tagy, Buďte prvý, kto otaguje tento záznam!
|
| Shrnutí: | In this work, we propose a novel preconditioned Krylov subspace method for solvingan optimal control problem of wave equations, after explicitly identifying the asymptotic spectraldistribution of the involved sequence of linear coefficient matrices from the optimal control prob-lem. Namely, we first show that the all-at-once system stemming from the wave control problem isassociated to a structured coefficient matrix-sequence possessing an eigenvalue distribution. Then, based on such a spectral distribution of which the symbol is explicitly identified, we develop an ideal preconditioner and two parallel-in-time preconditioners for the saddle point system composed of two block Toeplitz matrices. For the ideal preconditioner, we show that the eigenvalues of the precon-ditioned matrix-sequence all belong to the set (-3/2, -1/2) U (1/2, 3/2) well separated from zero, leadingto mesh-independent convergence when the minimal residual method is employed. The proposed parallel-in-time preconditioners can be implemented efficiently using fast Fourier transforms or dis-crete sine transforms, and their effectiveness is theoretically shown in the sense that the eigenvaluesof the preconditioned matrix-sequences are clustered around ±1, which leads to rapid convergence.When these parallel-in-time preconditioners are not fastly diagonalizable, we further propose modi-fied versions which can be efficiently inverted. Several numerical examples are reported to verify our derived localization and spectral distribution result and to support the effectiveness of our proposed preconditioners and the related advantages with respect to the relevant literature. |
|---|---|
| ISSN: | 0895-4798 1095-7162 1095-7162 |
| DOI: | 10.1137/23M1547251 |