14 papers · ranked by Valyu relevance
Brendan Gavin, Agnieszka Międlar, Eric Polizzi
The linear FEAST algorithm is a method for solving linear eigenvalue problems. It uses complex contour integration to calculate the eigenvectors whose eigenvalues that are located inside some user-defined region in the complex plane. This makes it possible to parallelize the process of solving eigenvalue problems by…
Benyamin Ghojogh, Fakhri Karray, Mark Crowley
This paper is a tutorial for eigenvalue and generalized eigenvalue problems. We first introduce eigenvalue problem, eigen-decomposition (spectral decomposition), and generalized eigenvalue problem. Then, we mention the optimization problems which yield to the eigenvalue and generalized eigenvalue problems. We also…
Marko Huhtanen, Vesa Kotila, Pauliina Uusitalo
eigenvalue problems Authors: ['Marko Huhtanen' 'Vesa Kotila' 'Pauliina Uusitalo'] Abstract. For the generalized eigenvalue problem, a quotient function is devised for estimating eigenvalues in terms of an approximate eigenvector. This gives rise to an infinite family of quotients, all entirely arguable to be used in…
Henrik Eisenmann, Yuji Nakatsukasa
We present a new approach to compute selected eigenvalues and eigenvectors of the twoparameter eigenvalue problem. Our method requires computing generalized eigenvalue problems of the same size as the matrices of the initial two-parameter eigenvalue problem. The method is applicable for right definite problems…
Michiel E. Hochstenbach, Christian Mehl, Bor Plestenjak
Generalized eigenvalue problems involving a singular pencil are very challenging to solve, both with respect to accuracy and efficiency. The existing package Guptri is very elegant but may be time-demanding, even for small and medium-sized matrices. We propose a simple method to compute the eigenvalues of singular…
Victor Janssens, Karl Meerbergen, Wim Michiels
Nonlinear eigenvalue problems with eigenvector nonlinearities (NEPv) are algebraic eigenvalue problems whose matrix depends on the eigenvector. Applications range from computational quantum mechanics to machine learning. Due to its nonlinear behavior, existing methods almost exclusively rely on fixed-point iterations…
James Demmel, Ioana Dumitriu, Ryan Schneider
We present a randomized, inverse-free algorithm for producing an approximate diagonalization of any n × n matrix pencil (A, B). The bulk of the algorithm rests on a randomized divide-and-conquer eigensolver for the generalized eigenvalue problem originally proposed by Ballard, Demmel, and Dumitriu [Technical Report…
Pavel Osinenko, Grigory Devadze, Stefan Streif
The eigenvalue problem plays a central role in linear algebra and its applications in control and optimization methods. In particular, many matrix decompositions rely upon computation of eigenvalue-eigenvector pairs, such as diagonal or Jordan normal forms. Unfortunately, numerical algorithms computing eigenvectors are…
de Leeuw, Jan
We give formulae for first and second derivatives of generalized eigenvalues/eigenvectors of symmetric matrices and generalized singular values/singular vectors of rectangular matrices when the matrices are linear or nonlinear functions of a vector of parameters. In addition we provide functions in R to compute these…
Qiang Ye
This paper is concerned with computations of a few smallest eigenvalues (in absolute value) of a large extremely ill-conditioned matrix. It is shown that a few smallest eigenvalues can be accurately computed for a diagonally dominant matrix or a product of diagonally dominant matrices by combining a standard iterative…
P. Cheema, M. M. Alamdari, G. A. Vio
Least Eigenvalue of an Unknown Mass Matrix Authors: ['P. Cheema' 'M. M. Alamdari' 'G. A. Vio'] In the field of structural engineering analysis, a common requirement is to calculate the modal frequencies of a structure that has undergone an update, either naturally (such as from material degradation), or due to manmade…
Susanne Bradley, Chen Greif
We derive bounds on the eigenvalues of a generic form of double saddlepoint matrices. The bounds are expressed in terms of extremal eigenvalues and singular values of the associated block matrices. Inertia and algebraic multiplicity of eigenvalues are considered as well. The analysis includes bounds for preconditioned…
Udita N. Katugampola
Diagonalizing a matrix A, that is finding two matrices P and D such that A = P DP −1 with D being a diagonal matrix needs two steps: first find the eigenvalues and then find the corresponding eigenvectors. We show that we do not need the second step when diagonalizing matrices with a spectrum, |σ(A)| ≤ 2 since those…
Juán Tolosa
Starting from a mistake done by a student, we discover an unexpected method of finding both eigenvectors for a 2×2 matrix with distinct eigenvalues in a single computation. We discuss a connection with the Cayley-Hamilton theorem, and show the corresponding generalization for a 3 × 3 matrix. The arguments should be…