We prove a Saad's type bound for harmonic Ritz vectors of a Hermitian matrix. The new bound reveals a dependence of the harmonic Rayleigh--Ritz procedure on the condition number of a shifted problem operator. Several practical implications are discussed. In particular, the bound motivates incorporation of preconditioning into the harmonic Rayleigh--Ritz scheme.
In particular, if @math is Hermitian, the result of @cite_4 states that eq:chen.jia (x, v) 1 + ^2 |B^ -1 |^2 sep ( ,G)^2 (x, K), where @math , @math , and @math denotes an orthogonal projector onto @math . The quantity @math describes separation of the targeted eigenvalue @math from the harmonic Ritz values different from the one associated with @math , which are the eigenvalues of @math ; see @cite_4 for the precise definition of @math .
@cite_4
1975625541
Let A be a real m×n matrix with m≧n. It is well known (cf. [4]) that @math (1) where @math The matrix U consists of n orthonormalized eigenvectors associated with the n largest eigenvalues of AA T , and the matrix V consists of the orthonormalized eigenvectors of A T A. The diagonal elements of ∑ are the non-negative square roots of the eigenvalues of A T A; they are called singular values. We shall assume that @math Thus if rank(A)=r, σ r+1 = σ r+2=⋯=σ n = 0. The decomposition (1) is called the singular value decomposition (SVD). We derive a new upper bound on the diameter of a polyhedron @math P = x ? R n : A x ≤ b , where @math A ? Z m × n . The bound is polynomial in @math n and the largest absolute value of a sub-determinant of @math A , denoted by @math Δ . More precisely, we show that the diameter of @math P is bounded by @math O ( Δ 2 n 4 log n Δ ) . If @math P is bounded, then we show that the diameter of @math P is at most @math O ( Δ 2 n 3.5 log n Δ ) . For the special case in which @math A is a totally unimodular matrix, the bounds are @math O ( n 4 log n ) and @math O ( n 3.5 log n ) respectively. This improves over the previous best bound of @math O ( m 16 n 3 ( log m n ) 3 ) due to Dyer and Frieze (Math Program 64:1---16, 1994). Introduction Preliminaries Special linear groups: @math Orthogonal groups: @math Symplectic groups: @math Unitary groups: @math Proofs of Theorems 1.1 and 1.1, and of corollaries 1.2-1.4 Permutation group algorithms Concluding remarks References. The little Grothendieck problem consists of maximizing @math źijCijxixj for a positive semidefinite matrix C, over binary variables @math xiź ?1 . In this paper we focus on a natural generalization of this problem, the little Grothendieck problem over the orthogonal group. Given @math CźRdn?dn a positive semidefinite matrix, the objective is to maximize @math źijtrCijTOiOjT restricting @math Oi to take values in the group of orthogonal matrices @math Od, where @math Cij denotes the (ij)-th @math d?d block of C. We propose an approximation algorithm, which we refer to as Orthogonal-Cut, to solve the little Grothendieck problem over the group of orthogonal matrices @math Od and show a constant approximation ratio. Our method is based on semidefinite programming. For a given @math dź1, we show a constant approximation ratio of @math źR(d)2, where @math źR(d) is the expected average singular value of a @math d?d matrix with random Gaussian @math N0,1d i.i.d. entries. For @math d=1 we recover the known @math źR(1)2=2 ź approximation guarantee for the classical little Grothendieck problem. Our algorithm and analysis naturally extends to the complex valued case also providing a constant approximation ratio for the analogous little Grothendieck problem over the Unitary Group @math Ud. Orthogonal-Cut also serves as an approximation algorithm for several applications, including the Procrustes problem where it improves over the best previously known approximation ratio of @math 122. The little Grothendieck problem falls under the larger class of problems approximated by a recent algorithm proposed in the context of the non-commutative Grothendieck inequality. Nonetheless, our approach is simpler and provides better approximation with matching integrality gaps. Finally, we also provide an improved approximation algorithm for the more general little Grothendieck problem over the orthogonal (or unitary) group with rank constraints, recovering, when @math d=1, the sharp, known ratios.