An SVD-approach to Jacobi-Davidson solution of nonlinear Helmholtz eigenvalue problems

Numerical solution of the Helmholtz equation in an infinite domain often involves restriction of the domain to a bounded computational window where a numerical solution method is applied. On the boundary of the computational window artificial transparent boundary conditions are posed, for example, widely used perfectly matched layers (PMLs) or absorbing boundary conditions (ABCs). Recently proposed transparent-influx boundary conditions (TIBCs) resolve a number of drawbacks typically attributed to PMLs and ABCs, such as introduction of spurious solutions and the inability to have a tight computational window. Unlike the PMLs or ABCs, the TIBCs lead to a nonlinear dependence of the boundary integral operator on the frequency. Thus, a nonlinear Helmholtz eigenvalue problem arises. \ud This paper presents an approach for solving such nonlinear eigenproblems which is based on a truncated singular value decomposition (SVD) polynomial approximation of the nonlinearity and subsequent solution of the obtained approximate polynomial eigenproblem with the Jacobi-Davidson method

A polynomial Jacobi-Davidson solver with support for non-monomial bases and deflation

[EN] Large-scale polynomial eigenvalue problems can be solved by Krylov methods operating on an equivalent linear eigenproblem (linearization) of size d center dot n where d is the polynomial degree and n is the problem size, or by projection methods that keep the computation in the n-dimensional space. Jacobi-Davidson belongs to the latter class of methods, and, since it is a preconditioned eigensolver, it may be competitive in cases where explicitly computing a matrix factorization is exceedingly expensive. However, a fully fledged implementation of polynomial Jacobi-Davidson has to consider several issues, including deflation to compute more than one eigenpair, use of non-monomial bases for the case of large degree polynomials, and handling of complex eigenvalues when computing in real arithmetic. We discuss these aspects and present computational results of a parallel implementation in the SLEPc library.This work was supported by Agencia Estatal de Investigación (AEI) under Grant TIN2016-75985-P, which includes European Commission ERDF funds.Campos, C.; Jose E. Roman (2020). A polynomial Jacobi-Davidson solver with support for non-monomial bases and deflation. BIT Numerical Mathematics. 60(2):295-318. https://doi.org/10.1007/s10543-019-00778-zS295318602Bai, Z., Su, Y.: SOAR: a second-order Arnoldi method for the solution of the quadratic eigenvalue problem. SIAM J. Matrix Anal. Appl. 26(3), 640–659 (2005)Balay, S., Abhyankar, S., Adams, M., Brown, J., Brune, P., Buschelman, K., Dalcin, L., Eijkhout, V., Gropp, W., Karpeyev, D., Kaushik, D., Knepley, M., May, D., McInnes, L.C., Mills, R., Munson, T., Rupp, K., Sanan, P., Smith, B., Zampini, S., Zhang, H., Zhang, H.: PETSc users manual. Technical report ANL-95/11—revision 3.10, Argonne National Laboratory (2018)Betcke, T., Kressner, D.: Perturbation, extraction and refinement of invariant pairs for matrix polynomials. Linear Algebra Appl. 435(3), 514–536 (2011)Betcke, T., Voss, H.: A Jacobi–Davidson-type projection method for nonlinear eigenvalue problems. Future Gen. Comput. Syst. 20(3), 363–372 (2004)Betcke, T., Higham, N.J., Mehrmann, V., Schröder, C., Tisseur, F.: NLEVP: a collection of nonlinear eigenvalue problems. ACM Trans. Math. Softw. 39(2), 7:1–7:28 (2013)Campos, C., Roman, J.E.: Parallel Krylov solvers for the polynomial eigenvalue problem in SLEPc. SIAM J. Sci. Comput. 38(5), S385–S411 (2016)Effenberger, C.: Robust successive computation of eigenpairs for nonlinear eigenvalue problems. SIAM J. Matrix Anal. Appl. 34(3), 1231–1256 (2013)Effenberger, C., Kressner, D.: Chebyshev interpolation for nonlinear eigenvalue problems. BIT 52(4), 933–951 (2012)Fokkema, D.R., Sleijpen, G.L.G., van der Vorst, H.A.: Jacobi–Davidson style QR and QZ algorithms for the reduction of matrix pencils. SIAM J. Sci. Comput. 20(1), 94–125 (1998)Guo, J.S., Lin, W.W., Wang, C.S.: Numerical solutions for large sparse quadratic eigenvalue problems. Linear Algebra Appl. 225, 57–89 (1995)Hernandez, V., Roman, J.E., Vidal, V.: SLEPc: a scalable and flexible toolkit for the solution of eigenvalue problems. ACM Trans. Math. Softw. 31(3), 351–362 (2005)Higham, N.J., Al-Mohy, A.H.: Computing matrix functions. Acta Numer. 19, 159–208 (2010)Higham, N.J., Mackey, D.S., Tisseur, F.: The conditioning of linearizations of matrix polynomials. SIAM J. Matrix Anal. Appl. 28(4), 1005–1028 (2006)Hochbruck, M., Lochel, D.: A multilevel Jacobi–Davidson method for polynomial PDE eigenvalue problems arising in plasma physics. SIAM J. Sci. Comput. 32(6), 3151–3169 (2010)Hochstenbach, M.E., Sleijpen, G.L.G.: Harmonic and refined Rayleigh–Ritz for the polynomial eigenvalue problem. Numer. Linear Algebra Appl. 15(1), 35–54 (2008)Huang, T.M., Hwang, F.N., Lai, S.H., Wang, W., Wei, Z.H.: A parallel polynomial Jacobi–Davidson approach for dissipative acoustic eigenvalue problems. Comput. Fluids 45(1), 207–214 (2011)Hwang, F.N., Wei, Z.H., Huang, T.M., Wang, W.: A parallel additive Schwarz preconditioned Jacobi–Davidson algorithm for polynomial eigenvalue problems in quantum dot simulation. J.Comput. Phys. 229(8), 2932–2947 (2010)Kressner, D.: A block Newton method for nonlinear eigenvalue problems. Numer. Math. 114, 355–372 (2009)Kressner, D., Roman, J.E.: Memory-efficient Arnoldi algorithms for linearizations of matrix polynomials in Chebyshev basis. Numer. Linear Algebra Appl. 21(4), 569–588 (2014)Lancaster, P.: Linearization of regular matrix polynomials. Electron. J. Linear Algebra 17, 21–27 (2008)Matsuo, Y., Guo, H., Arbenz, P.: Experiments on a parallel nonlinear Jacobi–Davidson algorithm. Procedia Comput. Sci. 29, 565–575 (2014)Meerbergen, K.: Locking and restarting quadratic eigenvalue solvers. SIAM J. Sci. Comput. 22(5), 1814–1839 (2001)Roman, J.E., Campos, C., Romero, E., Tomas, A.: SLEPc users manual. Technical report DSIC-II/24/02—Revision 3.10, D. Sistemes Informàtics i Computació, Universitat Politècnica de València (2018)Romero, E., Roman, J.E.: A parallel implementation of Davidson methods for large-scale eigenvalue problems in SLEPc. ACM Trans. Math. Softw. 40(2), 13:1–13:29 (2014)Rommes, J., Martins, N.: Computing transfer function dominant poles of large-scale second-order dynamical systems. SIAM J. Sci. Comput. 30(4), 2137–2157 (2008)Saad, Y.: Iterative Methods for Sparse Linear Systems, 2nd edn. SIAM Publications, Philadelphia (2003)Sensiau, C., Nicoud, F., van Gijzen, M., van Leeuwen, J.W.: A comparison of solvers for quadratic eigenvalue problems from combustion. Int. J. Numer. Methods Fluids 56(8), 1481–1488 (2008)Sleijpen, G.L.G., van der Vorst, H.A.: A Jacobi–Davidson iteration method for linear eigenvalue problems. SIAM J. Matrix Anal. Appl. 17(2), 401–425 (1996)Sleijpen, G.L.G., Booten, A.G.L., Fokkema, D.R., van der Vorst, H.A.: Jacobi–Davidson type methods for generalized eigenproblems and polynomial eigenproblems. BIT 36(3), 595–633 (1996)Sleijpen, G.L.G., van der Vorst, H.A., Meijerink, E.: Efficient expansion of subspaces in the Jacobi–Davidson method for standard and generalized eigenproblems. Electron. Trans. Numer. Anal. 7, 75–89 (1998)Tisseur, F., Meerbergen, K.: The quadratic eigenvalue problem. SIAM Rev. 43(2), 235–286 (2001)van Gijzen, M.B., Raeven, F.: The parallel computation of the smallest eigenpair of an acoustic problem with damping. Int. J. Numer. Methods Eng. 45(6), 765–777 (1999)van Noorden, T., Rommes, J.: Computing a partial generalized real Schur form using the Jacobi–Davidson method. Numer. Linear Algebra Appl. 14(3), 197–215 (2007)Voss, H.: A Jacobi–Davidson method for nonlinear and nonsymmetric eigenproblems. Comput. Struct. 85(17–18), 1284–1292 (2007

Experiments on a Parallel Nonlinear Jacobi–Davidson Algorithm

AbstractThe Jacobi–Davidson (JD) algorithm is very well suited for the computation of a few eigen-pairs of large sparse complex symmetric nonlinear eigenvalue problems. The performance of JD crucially depends on the treatment of the so-called correction equation, in particular the preconditioner, and the initial vector. Depending on the choice of the spectral shift and the accuracy of the solution, the convergence of JD can vary from linear to cubic. We investigate parallel preconditioners for the Krylov space method used to solve the correction equation.We apply our nonlinear Jacobi–Davidson (NLJD) method to quadratic eigenvalue problems that originate from the time-harmonic Maxwell equation for the modeling and simulation of resonating electromagnetic structures

On the Convergence of Ritz Pairs and Refined Ritz Vectors for Quadratic Eigenvalue Problems

For a given subspace, the Rayleigh-Ritz method projects the large quadratic eigenvalue problem (QEP) onto it and produces a small sized dense QEP. Similar to the Rayleigh-Ritz method for the linear eigenvalue problem, the Rayleigh-Ritz method defines the Ritz values and the Ritz vectors of the QEP with respect to the projection subspace. We analyze the convergence of the method when the angle between the subspace and the desired eigenvector converges to zero. We prove that there is a Ritz value that converges to the desired eigenvalue unconditionally but the Ritz vector converges conditionally and may fail to converge. To remedy the drawback of possible non-convergence of the Ritz vector, we propose a refined Ritz vector that is mathematically different from the Ritz vector and is proved to converge unconditionally. We construct examples to illustrate our theory.Comment: 20 page