Quartic eigenvalue problem $(łambda^4 A + łambda^3 B + łambda^2C + łambda
D + E)x = 0$ 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.
Description
An algorithm for the complete solution of the quartic eigenvalue problem
%0 Generic
%1 drmac2019algorithm
%A Drmač, Zlatko
%A Glibić, Ivana Šain
%D 2019
%K 65f15-numerical-eigenvalues-eigenvectors
%T An algorithm for the complete solution of the quartic eigenvalue problem
%U http://arxiv.org/abs/1905.07013
%X Quartic eigenvalue problem $(łambda^4 A + łambda^3 B + łambda^2C + łambda
D + E)x = 0$ 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.
@misc{drmac2019algorithm,
abstract = {Quartic eigenvalue problem $(\lambda^4 A + \lambda^3 B + \lambda^2C + \lambda
D + E)x = \mathbf{0}$ 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.},
added-at = {2020-07-08T07:28:30.000+0200},
author = {Drmač, Zlatko and Glibić, Ivana Šain},
biburl = {https://www.bibsonomy.org/bibtex/2461d919600d69d71ad0300949c247b95/gdmcbain},
description = {An algorithm for the complete solution of the quartic eigenvalue problem},
interhash = {e48d4d7968299cb7f36c61516b3bc6d7},
intrahash = {461d919600d69d71ad0300949c247b95},
keywords = {65f15-numerical-eigenvalues-eigenvectors},
note = {cite arxiv:1905.07013},
timestamp = {2020-07-08T07:28:30.000+0200},
title = {An algorithm for the complete solution of the quartic eigenvalue problem},
url = {http://arxiv.org/abs/1905.07013},
year = 2019
}