Taslaman, Leo (2014) An algorithm for quadratic eigenproblems with low rank damping. [MIMS Preprint]
PDF
revisedversion.pdf Download (609kB) 
Abstract
We consider quadratic eigenproblems $\left(M\lambda^2+D\lambda+K\right)x=0$, where all coefficient matrices are real and positive semidefinite, $(M,K)$ is regular and $D$ is of low rank. Matrix polynomials of this form appear in the analysis of vibrating structures with discrete dampers. We develop an algorithm for such problems, which first solves the undamped problem $\left(M\lambda^2+K\right)x=0$ and then accommodates for the low rank term $D\lambda$. For the first part, we develop a new algorithm based on a method proposed by Wang and Zhao [SIAM J. Matrix Anal. Appl. 124 (1991), pp. 654660], which can compute all eigenvalues of definite generalized eigenvalue problems with semidefinite coefficient matrices in a backward stable and symmetry preserving manner. We use this new algorithm to compute the solution to the undamped problem, and then use this solution in order to compute all eigenvalues of the original problem, and the associated eigenvectors if requested. To this end, we use an EhrlichAberth iteration that works exclusively with vectors and tall skinny matrices and contributes only lower order terms to the overall flop count. Numerical experiments show that the proposed algorithm is both fast and accurate. Finally we discuss the application to large scale quadratics and the possibility of generalizations to other problems.
Item Type:  MIMS Preprint 

Uncontrolled Keywords:  quadratic eigenvalue problem, eigenvalue algorithm, matrix polynomial, discrete damper, vibrating system 
Subjects:  MSC 2010, the AMS's Mathematics Subject Classification > 15 Linear and multilinear algebra; matrix theory MSC 2010, the AMS's Mathematics Subject Classification > 65 Numerical analysis MSC 2010, the AMS's Mathematics Subject Classification > 70 Mechanics of particles and systems 
Depositing User:  Leo Taslaman 
Date Deposited:  07 Jan 2015 
Last Modified:  08 Nov 2017 18:18 
URI:  https://eprints.maths.manchester.ac.uk/id/eprint/2233 
Available Versions of this Item

An algorithm for quadratic eigenproblems with low rank damping. (deposited 14 May 2014)
 An algorithm for quadratic eigenproblems with low rank damping. (deposited 07 Jan 2015) [Currently Displayed]
Actions (login required)
View Item 