Solving the Parametric Eigenvalue Problem by Taylor Series and Chebyshev Expansion

by   Thomas Mach, et al.

We discuss two approaches to solving the parametric (or stochastic) eigenvalue problem. One of them uses a Taylor expansion and the other a Chebyshev expansion. The parametric eigenvalue problem assumes that the matrix A depends on a parameter μ, where μ might be a random variable. Consequently, the eigenvalues and eigenvectors are also functions of μ. We compute a Taylor approximation of these functions about μ_0 by iteratively computing the Taylor coefficients. The complexity of this approach is O(n^3) for all eigenpairs, if the derivatives of A(μ) at μ_0 are given. The Chebyshev expansion works similarly. We first find an initial approximation iteratively which we then refine with Newton's method. This second method is more expensive but provides a good approximation over the whole interval of the expansion instead around a single point. We present numerical experiments confirming the complexity and demonstrating that the approaches are capable of tracking eigenvalues at intersection points. Further experiments shed light on the limitations of the Taylor expansion approach with respect to the distance from the expansion point μ_0.


page 1

page 2

page 3

page 4


Solving two-parameter eigenvalue problems using an alternating method

We present a new approach to compute selected eigenvalues and eigenvecto...

Computation of Tight Enclosures for Laplacian Eigenvalues

Recently, there has been interest in high-precision approximations of th...

On the Numerical Approximation of the Karhunen-Loève Expansion for Lognormal Random Fields

The Karhunen-Loève (KL) expansion is a popular method for approximating ...

Fast non-Hermitian Toeplitz eigenvalue computations, joining matrix-less algorithms and FDE approximation matrices

The present work is devoted to the eigenvalue asymptotic expansion of th...

Match-based solution of general parametric eigenvalue problems

We describe a novel algorithm for solving general parametric (nonlinear)...

The conditional DPP approach to random matrix distributions

We present the conditional determinantal point process (DPP) approach to...

Dynamical perturbation theory for eigenvalue problems

Many problems in physics, chemistry and other fields are perturbative in...

Please sign up or login with your details

Forgot password? Click here to reset