Adaptive System Identification for GMRES(m) Based on DMDc
Abstract
The Restarted Generalized Minimal Residual Method (GMRES(m)) is a widely used iterative solver for large-scale, sparse, and nonsymmetric linear systems [4]. However, its convergence is not guaranteed [1]. To improve convergence, the restart parameter m must be properly tuned; however, there is currently no general rule for selecting its optimal value. The primary challenge in designing a rule is the internal nonlinear dynamics of GMRES(m), which is well understood in its construction, but less well understood in its behavior. In this work, we design a system identification model (SIM) to address the challenge of capturing the internal nonlinear behavior of GMRES(m). The results show the efficiency of the SIM as an observer of the internal dynamics. Future works will be oriented towards using the SIM to enhance the convergence of the method. [...]
Downloads
References
J. C. Cabral, C. E. Schaerer, and A. Bhaya. “Improving GMRES(m) using an adaptive switching controller”. In: Numerical Linear Algebra with Applications 27.5 (2020), e2305. doi: 10.1002/nla.2305.
S. P. Kolodziej, M. Aznaveh, M. Bullock, J. David, T. A Davis, M. Henderson, Y. Hu, and R. Sandstrom. “The suitesparse matrix collection website interface”. In: Journal of Open Source Software 4.35 (2019), p. 1244. doi: 10.21105/joss.01244.
J. L. Proctor, S. L. Brunton, and J. Nathan Kutz. “Dynamic mode decomposition with control”. In: SIAM Journal on Applied Dynamical Systems 15.1 (2016), pp. 142–161. doi: 10.1137/15M1013857.
Y. Saad. Iterative Methods for Sparse Linear Systems. 2nd ed. Philadelphia: Society for Industrial and Applied Mathematics, 2003. isbn: 9780898715347. url: https://epubs.siam.org/doi/book/10.1137/1.9780898718003.
B. Zhong and R. Morgan. “Complementary cycles of restarted GMRES”. In: Numerical Linear Algebra with Applications 15.6 (2008), pp. 559–571. doi: 10.1002/nla.589.