An optimal diagonalization-based preconditioner for parabolic optimal control problems (2410.22686v2)

Abstract: In this work, we propose a novel diagonalization-based preconditioner for the all-at-once linear system arising from the optimal control problem of parabolic equations. The proposed preconditioner is constructed based on an $\epsilon$-circulant modification to the rotated block diagonal (RBD) preconditioning technique and can be efficiently diagonalized by fast Fourier transforms in a parallel-in-time fashion. To our knowledge, this marks the first application of the $\epsilon$-circulant modification to RBD preconditioning. Before our work, the studies of parallel-in-time preconditioning techniques for the optimal control problem are mainly focused on $\epsilon$-circulant modification to Schur complement based preconditioners, which involves multiplication of forward and backward evolutionary processes and thus square the condition number. Compared with those Schur complement based preconditioning techniques in the literature, the advantage of the proposed $\epsilon$-circulant modified RBD preconditioning is that it does not involve the multiplication of forward and backward evolutionary processes. When the generalized minimal residual method is deployed on the preconditioned system, we prove that when choosing $\epsilon=\mathcal{O}(\sqrt{\tau})$ with $\tau$ being the temporal step-size, the convergence rate of the preconditioned GMRES solver is independent of the matrix size and the regularization parameter. Numerical results are provided to demonstrate the effectiveness of our proposed solvers.