The correction equation in the Jacobi-Davidson method is effective in a subspace orthogonal to the current eigenvector approximation, whereas for the continuation of the process only vectors orthogonal to the search subspace are of importance. Such a vector is obtained by orthogonalizing the (approx
A refined jacobi-davidson method and its correction equation
✍ Scribed by Shaoqiang Feng; Zhongxiao Jia
- Publisher
- Elsevier Science
- Year
- 2005
- Tongue
- English
- Weight
- 602 KB
- Volume
- 49
- Category
- Article
- ISSN
- 0898-1221
No coin nor oath required. For personal study only.
✦ Synopsis
A central problem in the Jacobi-Davidson method is to expand a projection subspace by solving a certain correction equation. It has been commonly accepted that the correction equation always has a solution. However, it is proved in this paper that this is not true. Conditions are given to decide when it has a unique solution or many solutions or no solution. A refined Jacobi-Davidson method is proposed to overcome the possible nonconvergence of Ritz vectors by computing certain refined approximation eigenvectors from the subspace. A corresponding correction equation is derived for the refined method. Numerical experiments are conducted and efficiency of the refined method is confirmed. (~) 2005 Elsevier Ltd. All rights reserved.
📜 SIMILAR VOLUMES
Jacobi methods for computing the eigendecomposition of a class of so-called low-rank-plus-shift symmetric matrices are investigated. An order-of-magnitude reduction in the computational complexity can be achieved for this special class of matrices by terminating the Jacobi sweep early in a cyclic or
## Abstract In this paper, a new variant of the Jacobi–Davidson (JD) method is presented that is specifically designed for __real unsymmetric__ matrix pencils. Whenever a pencil has a complex conjugate pair of eigenvalues, the method computes the two‐dimensional real invariant subspace spanned by t
It is proved that the function deÿned by the inÿmum-based Lax formula (for viscosity solutions) provides a solution almost everywhere in x for each ÿxed t ¿ 0 to the Hamilton-Jacobi, Cauchy problem ut + 1 2 ∇u 2 = 0; u(x; 0 + ) = v(x); where the Cauchy data function v is lower semicontinuous on rea