Skip to content
Licensed Unlicensed Requires Authentication Published by De Gruyter March 12, 2021

Smoothed-Adaptive Perturbed Inverse Iteration for Elliptic Eigenvalue Problems

Stefano Giani ORCID logo, Luka Grubišić ORCID logo, Luca Heltai ORCID logo and Ornela Mulita


We present a perturbed subspace iteration algorithm to approximate the lowermost eigenvalue cluster of an elliptic eigenvalue problem. As a prototype, we consider the Laplace eigenvalue problem posed in a polygonal domain. The algorithm is motivated by the analysis of inexact (perturbed) inverse iteration algorithms in numerical linear algebra. We couple the perturbed inverse iteration approach with mesh refinement strategy based on residual estimators. We demonstrate our approach on model problems in two and three dimensions.

MSC 2010: 65N25; 65N30; 65N50

Funding source: Hrvatska Zaklada za Znanost

Award Identifier / Grant number: IP-2019-04-6268

Funding statement: Luka Grubišić was supported by the Croatian Science Foundation grant HRZZ IP-2019-04-6268. Luka Grubišić is also thankful to the hospitality of the research visit to Scuola Internazionale Superiore di Studi Avanzati where the work has started. Ornela Mulita is thankful to the University of Zagreb for the hospitality during her collaborative research visit there. Luca Heltai was partially supported by the National Research Projects (PRIN 2017) “Numerical Analysis for Full and Reduced Order Methods for the efficient and accurate solution of complex systems governed by Partial Differential Equations”, funded by the Italian Ministry of Education, University, and Research.


[1] M. Arioli, A stopping criterion for the conjugate gradient algorithms in a finite element method framework, Numer. Math. 97 (2004), no. 1, 1–24. 10.1007/s00211-003-0500-ySearch in Google Scholar

[2] M. Arioli, E. H. Georgoulis and D. Loghin, Stopping criteria for adaptive finite element solvers, SIAM J. Sci. Comput. 35 (2013), no. 3, A1537–A1559. 10.1137/120867421Search in Google Scholar

[3] M. Arioli, J. Liesen, A. Miȩdlar and Z. Strakoš, Interplay between discretization and algebraic computation in adaptive numerical solution of elliptic PDE problems, GAMM-Mitt. 36 (2013), no. 1, 102–129.10.1002/gamm.201310006Search in Google Scholar

[4] D. Arndt, W. Bangerth, B. Blais, T. C. Clevenger, M. Fehling, A. V. Grayver, T. Heister, L. Heltai, M. Kronbichler, M. Maier, P. Munch, J.-P. Pelteret, R. Rastak, I. Tomas, B. Turcksin, Z. Wang and D. Wells, The deal.II library, version 9.2, J. Numer. Math. 28 (2020), no. 3, 131–146. 10.1515/jnma-2020-0043Search in Google Scholar

[5] D. Arndt, W. Bangerth, D. Davydov, T. Heister, L. Heltai, M. Kronbichler, M. Maier, J.-P. Pelteret, B. Turcksin and D. Wells, The deal.II finite element library: Design, features, and insights, Comput. Math. Appl. 81 (2021), 407–422. 10.1016/j.camwa.2020.02.022Search in Google Scholar

[6] I. Babuška and J. Osborn, Eigenvalue problems, Handbook of Numerical Analysis. Vol. II, Handb. Numer. Anal. II, North-Holland, Amsterdam (1991), 641–787. 10.1016/S1570-8659(05)80042-0Search in Google Scholar

[7] R. E. Bank, L. Grubišić and J. S. Ovall, A framework for robust eigenvalue and eigenvector error estimation and Ritz value convergence enhancement, Appl. Numer. Math. 66 (2013), 1–29. 10.1016/j.apnum.2012.11.004Search in Google Scholar

[8] R. Becker, C. Johnson and R. Rannacher, Adaptive error control for multigrid finite element methods, Computing 55 (1995), no. 4, 271–288. 10.1007/BF02238483Search in Google Scholar

[9] R. Becker and R. Rannacher, An optimal control approach to a posteriori error estimation in finite element methods, Acta Numer. 10 (2001), 1–102. 10.1017/S0962492901000010Search in Google Scholar

[10] P. Binev, W. Dahmen and R. DeVore, Adaptive finite element methods with convergence rates, Numer. Math. 97 (2004), no. 2, 219–268. 10.21236/ADA640658Search in Google Scholar

[11] N. Bosner and Z. Drmač, Subspace gap residuals for Rayleigh–Ritz approximations, SIAM J. Matrix Anal. Appl. 31 (2009), no. 1, 54–67. 10.1137/070689425Search in Google Scholar

[12] C. Carstensen, Some remarks on the history and future of averaging techniques in a posteriori finite element error analysis, ZAMM Z. Angew. Math. Mech. 84 (2004), no. 1, 3–21. 10.1002/zamm.200410101Search in Google Scholar

[13] C. Carstensen and J. Gedicke, An oscillation-free adaptive FEM for symmetric eigenvalue problems, Numer. Math. 118 (2011), no. 3, 401–427. 10.1007/s00211-011-0367-2Search in Google Scholar

[14] C. Carstensen and J. Gedicke, An adaptive finite element eigenvalue solver of asymptotic quasi-optimal computational complexity, SIAM J. Numer. Anal. 50 (2012), no. 3, 1029–1057. 10.1137/090769430Search in Google Scholar

[15] C. Carstensen and J. Gedicke, Guaranteed lower bounds for eigenvalues, Math. Comp. 83 (2014), no. 290, 2605–2629. 10.1090/S0025-5718-2014-02833-0Search in Google Scholar

[16] C. Carstensen, J. Gedicke, V. Mehrmann and A. Miȩdlar, An adaptive homotopy approach for non-selfadjoint eigenvalue problems, Numer. Math. 119 (2011), no. 3, 557–583. 10.1007/s00211-011-0388-xSearch in Google Scholar

[17] C. Carstensen, J. Gedicke, V. Mehrmann and A. Miȩedlar, An adaptive finite element method with asymptotic saturation for eigenvalue problems, Numer. Math. 128 (2014), no. 4, 615–634. 10.1007/s00211-014-0624-2Search in Google Scholar

[18] J. M. Cascon, C. Kreuzer, R. H. Nochetto and K. G. Siebert, Quasi-optimal convergence rate for an adaptive finite element method, SIAM J. Numer. Anal. 46 (2008), no. 5, 2524–2550. 10.1137/07069047XSearch in Google Scholar

[19] P. G. Ciarlet, The Finite Element Method for Elliptic Problems, Class. Appl. Math. 40, Society for Industrial and Applied Mathematics, Philadelphia, 2002. 10.1137/1.9780898719208Search in Google Scholar

[20] T. C. Clevenger and T. Heister, The deal.II tutorial step-50: Geometric Multigrid on adaptive meshes distributed in parallel, 2020. Search in Google Scholar

[21] T. C. Clevenger, T. Heister, G. Kanschat and M. Kronbichler, A flexible, parallel, adaptive geometric multigrid method for FEM, ACM Trans. Math. Software 47 (2021), no. 1, 1–27. 10.1145/3425193Search in Google Scholar

[22] W. Dahmen, T. Rohwedder, R. Schneider and A. Zeiser, Adaptive eigenvalue computation: Complexity estimates, Numer. Math. 110 (2008), no. 3, 277–312. 10.1007/s00211-008-0159-5Search in Google Scholar

[23] X. Dai, J. Xu and A. Zhou, Convergence and optimal complexity of adaptive finite element eigenvalue computations, Numer. Math. 110 (2008), no. 3, 313–355. 10.1007/s00211-008-0169-3Search in Google Scholar

[24] P. Daniel, A. Ern and M. Vohralík, An adaptive hp-refinement strategy with inexact solvers and computable guaranteed bound on the error reduction factor, Comput. Methods Appl. Mech. Engrg. 359 (2020), Article ID 112607. 10.1016/j.cma.2019.112607Search in Google Scholar

[25] P. Daniel and M. Vohralík, Guaranteed contraction of adaptive inexact h p -refinement strategies with realistic stopping criteria, preprint (2020), Search in Google Scholar

[26] W. Dörfler, A convergent adaptive algorithm for Poisson’s equation, SIAM J. Numer. Anal. 33 (1996), no. 3, 1106–1124. 10.1137/0733054Search in Google Scholar

[27] R. G. Durán, C. Padra and R. Rodríguez, A posteriori error estimates for the finite element approximation of eigenvalue problems, Math. Models Methods Appl. Sci. 13 (2003), no. 8, 1219–1229. 10.1142/S0218202503002878Search in Google Scholar

[28] A. Ern and M. Vohralík, Adaptive inexact Newton methods with a posteriori stopping criteria for nonlinear diffusion PDEs, SIAM J. Sci. Comput. 35 (2013), no. 4, A1761–A1791. 10.1137/120896918Search in Google Scholar

[29] E. M. Garau and P. Morin, Convergence and quasi-optimality of adaptive FEM for Steklov eigenvalue problems, IMA J. Numer. Anal. 31 (2011), no. 3, 914–946. 10.1093/imanum/drp055Search in Google Scholar

[30] E. M. Garau, P. Morin and C. Zuppa, Convergence of adaptive finite element methods for eigenvalue problems, Math. Models Methods Appl. Sci. 19 (2009), no. 5, 721–747. 10.1142/S0218202509003590Search in Google Scholar

[31] J. Gedicke and C. Carstensen, A posteriori error estimators for convection-diffusion eigenvalue problems, Comput. Methods Appl. Mech. Engrg. 268 (2014), 160–177. 10.1016/j.cma.2012.09.018Search in Google Scholar

[32] S. Giani, hp-adaptive composite discontinuous Galerkin methods for elliptic eigenvalue problems on complicated domains, Appl. Math. Comput. 267 (2015), 604–617. 10.1016/j.amc.2015.01.031Search in Google Scholar

[33] S. Giani and I. G. Graham, A convergent adaptive method for elliptic eigenvalue problems, SIAM J. Numer. Anal. 47 (2009), no. 2, 1067–1091.10.1137/070697264Search in Google Scholar

[34] S. Giani, L. Grubišić, A. Miȩdlar and J. S. Ovall, Robust error estimates for approximations of non-self-adjoint eigenvalue problems, Numer. Math. 133 (2016), no. 3, 471–495. 10.1007/s00211-015-0752-3Search in Google Scholar

[35] S. Giani, L. Grubišić and J. S. Ovall, Benchmark results for testing adaptive finite element eigenvalue procedures, Appl. Numer. Math. 62 (2012), no. 2, 121–140. 10.1016/j.apnum.2011.10.007Search in Google Scholar

[36] M. S. Gockenbach, Understanding and Implementing the Finite Element Method, Society for Industrial and Applied Mathematics, Philadelphia, 2006. 10.1137/1.9780898717846Search in Google Scholar

[37] D. S. Grebenkov and B.-T. Nguyen, Geometrical structure of Laplacian eigenfunctions, SIAM Rev. 55 (2013), no. 4, 601–667. 10.1137/120880173Search in Google Scholar

[38] L. Grubišić, A posteriori estimates for eigenvalue/vector approximations, PAMM Proc. Appl. Math. Mech. 6 (2006), no. 1, 59–62. 10.1002/pamm.200610016Search in Google Scholar

[39] L. Grubišić and J. S. Ovall, On estimators for eigenvalue/eigenvector approximations, Math. Comp. 78 (2009), no. 266, 739–770. 10.1090/S0025-5718-08-02181-9Search in Google Scholar

[40] W. Hackbusch, Multi-Grid Methods and Applications. Vol. 4, Springer, Berlin, 2013. Search in Google Scholar

[41] M. R. Hestenes and E. Stiefel, Methods of conjugate gradients for solving linear systems, J. Research Nat. Bur. Standards 49 (1952), 4099–436. 10.1090/psapm/006/0084178Search in Google Scholar

[42] B. Janssen and G. Kanschat, Adaptive multilevel methods with local smoothing for H 1 - and H curl -conforming high order finite element methods, SIAM J. Sci. Comput. 33 (2011), no. 4, 2095–2114. 10.1137/090778523Search in Google Scholar

[43] P. Jiránek, Z. Strakoš and M. Vohralík, A posteriori error estimates including algebraic error and stopping criteria for iterative solvers, SIAM J. Sci. Comput. 32 (2010), no. 3, 1567–1590. 10.1137/08073706XSearch in Google Scholar

[44] A. V. Knyazev, Toward the optimal preconditioned eigensolver: locally optimal block preconditioned conjugate gradient method, SIAM J. Sci. Comput. 23 (2001), 517–541. 10.1137/S1064827500366124Search in Google Scholar

[45] A. V. Knyazev and K. Neymeyr, A geometric theory for preconditioned inverse iteration. III. A short and sharp convergence estimate for generalized eigenvalue problems, Linear Algebra Appl. 358 (2003), no. 3–4, 95–114. 10.1016/S0024-3795(01)00461-XSearch in Google Scholar

[46] G. Mallik, M. Vohralík and S. Yousef, Goal-oriented a posteriori error estimation for conforming and nonconforming approximations with inexact solvers, J. Comput. Appl. Math. 366 (2020), Article ID 112367. 10.1016/ in Google Scholar

[47] A. Må lqvist and D. Peterseim, Computation of eigenvalues by numerical upscaling, Numer. Math. 130 (2015), no. 2, 337–361. 10.1007/s00211-014-0665-6Search in Google Scholar

[48] V. Mehrmann and A. Miȩdlar, Adaptive computation of smallest eigenvalues of self-adjoint elliptic partial differential equations, Numer. Linear Algebra Appl. 18 (2011), no. 3, 387–409. 10.1002/nla.733Search in Google Scholar

[49] V. Mehrmann and C. Schröder, Nonlinear eigenvalue and frequency response problems in industrial practice, J. Math. Ind. 1 (2011), Article ID 7. 10.1186/2190-5983-1-7Search in Google Scholar

[50] D. Meidner, R. Rannacher and J. Vihharev, Goal-oriented error control of the iterative solution of finite element equations, J. Numer. Math. 17 (2009), no. 2, 143–172. 10.1515/JNUM.2009.009Search in Google Scholar

[51] A. Miȩdlar, Inexact adaptive finite element methods for elliptic pde eigenvalue problems, Ph.D. thesis, Technische Universität Berlin, 2011. Search in Google Scholar

[52] A. Miȩdlar, A story on adaptive finite element computations for elliptic eigenvalue problems, Numerical Algebra, Matrix Theory, Differential-Algebraic Equations and Control Theory, Springer, Cham (2015), 223–255. 10.1007/978-3-319-15260-8_9Search in Google Scholar

[53] A. Miraçi, J. Papež and M. Vohralík, A multilevel algebraic error estimator and the corresponding iterative solver with p-robust behavior, SIAM J. Numer. Anal. 58 (2020), no. 5, 2856–2884. 10.1137/19M1275929Search in Google Scholar

[54] O. Mulita, Smoothed adaptive finite element methods, Ph.D. thesis, SISSA, 2019. Search in Google Scholar

[55] O. Mulita, S. Giani and L. Heltai, Quasi-optimal mesh sequence construction through smoothed adaptive finite element method, preprint (2020), 10.1137/19M1262097Search in Google Scholar

[56] K. Neymeyr, A geometric theory for preconditioned inverse iteration. I. Extrema of the Rayleigh quotient, Linear Algebra Appl. 322 (2001), no. 1–3, 61–85. 10.1016/S0024-3795(00)00239-1Search in Google Scholar

[57] K. Neymeyr, A geometric theory for preconditioned inverse iteration. II. Convergence estimates, Linear Algebra Appl. 322 (2001), no. 1–3, 87–104. 10.1016/S0024-3795(00)00236-6Search in Google Scholar

[58] K. Neymeyr, A posteriori error estimation for elliptic eigenproblems, Numer. Linear Algebra Appl. 9 (2002), no. 4, 263–279. 10.1002/nla.272Search in Google Scholar

[59] K. Neymeyr, A geometric theory for preconditioned inverse iteration. IV. On the fastest convergence cases, Linear Algebra Appl. 415 (2006), no. 1, 114–139. 10.1016/j.laa.2005.06.022Search in Google Scholar

[60] S. Oliveira, A convergence proof of an iterative subspace method for eigenvalues problems, Foundations of Computational Mathematics, Springer, Berlin (1997), 316–325. 10.1007/978-3-642-60539-0_25Search in Google Scholar

[61] J. Papež and Z. Strakoš, On a residual-based a posteriori error estimator for the total error, IMA J. Numer. Anal. 38 (2018), no. 3, 1164–1184. 10.1093/imanum/drx037Search in Google Scholar

[62] J. Papež, Z. Strakoš and M. Vohralík, Estimating and localizing the algebraic and total numerical errors using flux reconstructions, Numer. Math. 138 (2018), no. 3, 681–721. 10.1007/s00211-017-0915-5Search in Google Scholar

[63] R. Rannacher, A. Westenberger and W. Wollner, Adaptive finite element solution of eigenvalue problems: Balancing of discretization and iteration error, J. Numer. Math. 18 (2010), no. 4, 303–327. 10.1515/jnum.2010.015Search in Google Scholar

[64] T. Rohwedder, R. Schneider and A. Zeiser, Perturbed preconditioned inverse iteration for operator eigenvalue problems with applications to adaptive wavelet discretization, Adv. Comput. Math. 34 (2011), no. 1, 43–66. 10.1007/s10444-009-9141-8Search in Google Scholar

[65] Y. Saad, Analysis of subspace iteration for eigenvalue problems with evolving matrices, SIAM J. Matrix Anal. Appl. 37 (2016), no. 1, 103–122. 10.1137/141002037Search in Google Scholar

[66] A. Sartori, N. Giuliani, M. Bardelloni and L. Heltai, deal2lkit: A toolkit library for high performance programming in deal.II, SoftwareX 7 (2018), 318–327. 10.1016/j.softx.2018.09.004Search in Google Scholar

[67] R. Stevenson, Optimality of a standard adaptive finite element method, Found. Comput. Math. 7 (2007), no. 2, 245–269. 10.1007/s10208-005-0183-0Search in Google Scholar

[68] L. N. Trefethen and T. Betcke, Computed eigenmodes of planar regions, Recent Advances in Differential Equations and Mathematical Physics, Contemp. Math. 412, American Mathematical Society, Providence (2006), 297–314. 10.1090/conm/412/07783Search in Google Scholar

[69] R. S. Varga, Matrix Iterative Analysis, 2nd ed., Springer Ser. Comput. Math. 27, Springer, Berlin, 2009. Search in Google Scholar

[70] T. Vejchodský, Flux reconstructions in the Lehmann-Goerisch method for lower bounds on eigenvalues, J. Comput. Appl. Math. 340 (2018), 676–690. 10.1016/ in Google Scholar

[71] T. Vejchodský, Three methods for two-sided bounds of eigenvalues—a comparison, Numer. Methods Partial Differential Equations 34 (2018), no. 4, 1188–1208. 10.1002/num.22251Search in Google Scholar

[72] J. Xu and A. Zhou, A two-grid discretization scheme for eigenvalue problems, Math. Comp. 70 (2001), no. 233, 17–25. 10.1090/S0025-5718-99-01180-1Search in Google Scholar

[73] A. Zeiser, On the optimality of the inexact inverse iteration coupled with adaptive finite element methods, Report 57, Philipps-Universität Marburg, 2010, Search in Google Scholar

Received: 2020-03-02
Revised: 2021-02-13
Accepted: 2021-02-16
Published Online: 2021-03-12
Published in Print: 2021-04-01

© 2021 Walter de Gruyter GmbH, Berlin/Boston

Scroll Up Arrow