Abstract
AbstractWe present a low-rank greedily adapted hp-finite element algorithm for computing an approximation to the solution of the Lyapunov operator equation. We show that there is a hidden regularity in eigenfunctions of the solution of the Lyapunov equation which can be utilized to justify the use of high order finite element spaces. Our numerical experiments indicate that we achieve eight figures of accuracy for computing the trace of the solution of the Lyapunov equation posed in a dumbbell-domain using a finite element space of dimension of only $$10^4$$
10
4
degrees of freedom. Even more surprising is the observation that hp-refinement has an effect of reducing the rank of the approximation of the solution.
Funder
Hrvatska Zaklada za Znanost
Aalto University
Publisher
Springer Science and Business Media LLC
Subject
Applied Mathematics,Computational Mathematics,Computer Networks and Communications,Software
Reference37 articles.
1. Antoulas, A.C., Sorensen, D.C.: Approximation of large-scale dynamical systems: an overview. Int. J. Appl. Math. Comput. Sci. 11(5), 1093–1121 (2001). Numerical analysis and systems theory (Perpignan, 2000)
2. Babuška, I., Guo, B.Q.: Approximation properties of the $$h$$-$$p$$ version of the finite element method. Comput. Methods Appl. Mech. Eng. 133(3–4), 319–346 (1996)
3. Bank, R.E., Ovall, J.S.: Some remarks on interpolation and best approximation. Numer. Math. 137(2), 289–302 (2017)
4. Bonito, A., Lei, W., Pasciak, J.E.: Numerical approximation of space-time fractional parabolic equations. Comput. Methods Appl. Math. 17(4), 679–705 (2017)
5. Bonito, A., Lei, W., Pasciak, J.E.: On sinc quadrature approximations of fractional powers of regularly accretive operators. J. Numer. Math. 27(2), 57–68 (2019)