%0 Journal Article %1 citeulike:14130119 %A Hochstenbach, Michiel E. %A van der Vorst, Henk A. %D 2003 %J SIAM Journal on Scientific Computing %K 65f15-numerical-eigenvalues-eigenvectors %N 2 %P 591--603 %R 10.1137/s1064827502406403 %T Alternatives to the Rayleigh Quotient for the Quadratic Eigenvalue Problem %U http://dx.doi.org/10.1137/s1064827502406403 %V 25 %X We consider the quadratic eigenvalue problem \$łambda^2 Ax + Bx + Cx = 0.\$ Suppose that u is an approximation to an eigenvector x (for instance, obtained by a subspace method) and that we want to determine an approximation to the corresponding eigenvalue \$łambda\$. The usual approach is to impose the Galerkin condition \$r(þeta, u) = (þeta^2 A + B + C)u u\$, from which it follows that \$þeta\$ must be one of the two solutions to the quadratic equation \$(u^*Au) þeta^2 + (u^*Bu) + (u^*Cu) = 0\$. An unnatural aspect is that if \$u=x\$, the second solution has in general no meaning. When u is not very accurate, it may not be clear which solution is the best. Moreover, when the discriminant of the equation is small, the solutions may be very sensitive to perturbations in u. In this paper we therefore examine alternative approximations to \$łambda\$. We compare the approaches theoretically and by numerical experiments. The methods are extended to approximations from subspaces and to the polynomial eigenvalue problem.

@article{citeulike:14130119, abstract = {{We consider the quadratic eigenvalue problem \$\lambda^2 Ax + \lambda Bx + Cx = 0.\$ Suppose that u is an approximation to an eigenvector x (for instance, obtained by a subspace method) and that we want to determine an approximation to the corresponding eigenvalue \$\lambda\$. The usual approach is to impose the Galerkin condition \$r(\theta, u) = (\theta^2 A + \theta B + C)u \perp u\$, from which it follows that \$\theta\$ must be one of the two solutions to the quadratic equation \$(u^*Au) \theta^2 + (u^*Bu) \theta + (u^*Cu) = 0\$. An unnatural aspect is that if \$u=x\$, the second solution has in general no meaning. When u is not very accurate, it may not be clear which solution is the best. Moreover, when the discriminant of the equation is small, the solutions may be very sensitive to perturbations in u. In this paper we therefore examine alternative approximations to \$\lambda\$. We compare the approaches theoretically and by numerical experiments. The methods are extended to approximations from subspaces and to the polynomial eigenvalue problem.}}, added-at = {2017-06-29T07:13:07.000+0200}, author = {Hochstenbach, Michiel E. and van der Vorst, Henk A.}, biburl = {https://www.bibsonomy.org/bibtex/23fab5aa7595714ad91027c87a410f434/gdmcbain}, citeulike-article-id = {14130119}, citeulike-attachment-1 = {hochstenbach_03_alternatives_1082579.pdf; /pdf/user/gdmcbain/article/14130119/1082579/hochstenbach_03_alternatives_1082579.pdf; 0757e90853f0fe74ab5dd4575b7a592f812796bb}, citeulike-linkout-0 = {http://dx.doi.org/10.1137/s1064827502406403}, doi = {10.1137/s1064827502406403}, file = {hochstenbach_03_alternatives_1082579.pdf}, interhash = {ff7470a967d53a309260841a6784fc57}, intrahash = {3fab5aa7595714ad91027c87a410f434}, issn = {1064-8275}, journal = {SIAM Journal on Scientific Computing}, keywords = {65f15-numerical-eigenvalues-eigenvectors}, month = jan, number = 2, pages = {591--603}, posted-at = {2016-09-05 06:08:27}, priority = {2}, timestamp = {2017-06-29T07:13:07.000+0200}, title = {Alternatives to the {R}ayleigh Quotient for the Quadratic Eigenvalue Problem}, url = {http://dx.doi.org/10.1137/s1064827502406403}, volume = 25, year = 2003 }