On the shift-invert Lanczos method for the buckling eigenvalue problem

10/18/2019
by   Chao-Ping Lin, et al.
0

We consider the problem of extracting a few desired eigenpairs of the buckling eigenvalue problem Kx = λ K_Gx, where K is symmetric positive semi-definite, K_G is symmetric indefinite, and the pencil K - λ K_G is singular, namely, K and K_G share a non-trivial common nullspace. Moreover, in practical buckling analysis of structures, bases for the nullspace of K and the common nullspace of K and K_G are available. There are two open issues for developing an industrial strength shift-invert Lanczos method: (1) the shift-invert operator (K - σ K_G)^-1 does not exist or is extremely ill-conditioned, and (2) the use of the semi-inner product induced by K drives the Lanczos vectors rapidly towards the nullspace of K, which leads to a rapid growth of the Lanczos vectors in norms and cause permanent loss of information and the failure of the method. In this paper, we address these two issues by proposing a generalized buckling spectral transformation of the singular pencil K - λ K_G and a regularization of the inner product via a low-rank updating of the semi-positive definiteness of K. The efficacy of our approach is demonstrated by numerical examples, including one from industrial buckling analysis.

READ FULL TEXT
research
12/02/2022

Randomized low-rank approximation for symmetric indefinite matrices

The Nyström method is a popular choice for finding a low-rank approximat...
research
02/01/2017

On the optimality and sharpness of Laguerre's lower bound on the smallest eigenvalue of a symmetric positive definite matrix

Lower bounds on the smallest eigenvalue of a symmetric positive definite...
research
12/18/2021

Revisiting Memory Efficient Kernel Approximation: An Indefinite Learning Perspective

Matrix approximations are a key element in large-scale algebraic machine...
research
05/29/2020

Arnoldi algorithms with structured orthogonalization

We study a stability preserved Arnoldi algorithm for matrix exponential ...
research
08/16/2019

A Shift Selection Strategy for Parallel Shift-Invert Spectrum Slicing in Symmetric Self-Consistent Eigenvalue Computation

The central importance of large scale eigenvalue problems in scientific ...
research
04/07/2023

Hermitian Preconditioning for a class of Non-Hermitian Linear Systems

This work considers the convergence of GMRES for non-singular problems. ...
research
04/14/2022

A New Refinement-Free Preconditioner for the Symmetric Formulation in Electroencephalography

Widely employed for the accurate solution of the electroencephalography ...

Please sign up or login with your details

Forgot password? Click here to reset