An algorithm for the complete solution of the quartic eigenvalue problem

DOI10.1145/3494528arXiv1905.07013WikidataQ113309850 ScholiaQ113309850MaRDI QIDQ6318896FDOQ6318896

Authors: Zlatko Drmač, Ivana Šain Glibić 

Publication date: 16 May 2019

Abstract: Quartic eigenvalue problem $(lambd{a}^{4}A+lambd{a}^{3}B+lambd{a}^{2}C+lambdaD+E)x=mathbf0$ naturally arises e.g. when solving the Orr-Sommerfeld equation in the analysis of the stability of the {Poiseuille} flow, in theoretical analysis and experimental design of locally resonant phononic plates, modeling a robot with electric motors in the joints, calibration of catadioptric vision system, or e.g. computation of the guided and leaky modes of a planar waveguide. This paper proposes a new numerical method for the full solution (all eigenvalues and all left and right eigenvectors) that is based on quadratification, i.e. reduction of the quartic problem to a spectraly equivalent quadratic eigenvalue problem, and on a careful preprocessing to identify and deflate zero and infinite eigenvalues before the linearized quadratification is forwarded to the QZ algorithm. Numerical examples and backward error analysis confirm that the proposed algorithm is superior to the available methods.

Mathematics Subject Classification ID

Numerical analysis (65-XX)

This page was built for publication: An algorithm for the complete solution of the quartic eigenvalue problem

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6318896)