2012.113: A spatially adaptive iterative method for a class of nonlinear operator eigenproblems
2012.113: Elias Jarlebring and Stefan Güttel (2012) A spatially adaptive iterative method for a class of nonlinear operator eigenproblems.
Full text available as:
|PDF - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader|
We present a new algorithm for the iterative solution of nonlinear operator eigenvalue problems arising from partial differential equations. This algorithm combines automatic spatial resolution of linear operators with the infinite Arnoldi method for nonlinear matrix eigenproblems proposed in [E. JARLEBRING, W. MICHIELS, AND K. MEERBERGEN, A linear eigenvalue algorithm for the nonlinear eigenvalue problem, Numer. Math., 122 (2012)]. The iterates in this infinite Arnoldi method are functions, and each iteration requires the solution of an inhomogeneous differential equation. This formulation is independent of the spatial representation of the functions, which allows us to employ a dynamic representation with an accuracy of about the level of machine precision at each iteration, similar to what is done in the Chebfun system [Z. BATTLES AND L. N. TREFETHEN, An extension of MATLAB to continuous functions and operators, SIAM J. Sci. Comput., 25 (2004)] with its chebop functionality [T. A. DRISCOLL, F. BORNEMANN, AND L. N. TREFETHEN, The chebop system for automatic solution of differential equations, BIT, 48 (2008)], although our function representation is entirely based on coefficients instead of function values. Our approach also allows for nonlinearities in the boundary conditions of the PDE. The algorithm is illustrated with several examples, e.g., the study of eigenvalues of a vibrating string with delayed boundary feedback control.
|Item Type:||MIMS Preprint|
|Subjects:||MSC 2000 > 65 Numerical analysis|
|Deposited By:||Stefan Güttel|
|Deposited On:||28 November 2012|