Iterative computation of the smallest singular value and the corresponding singular vectors of a matrix

AbstractInfluenced by some techniques used for computing singular points of nonlinear equations, a generalized inverse iteration method is proposed for approximating the smallest singular value σn and the associated left and right singular vectors u, v of a matrix A∈Rn×n. In the practically relevant case σn>0 the method is mathematically equivalent to inverse iteration with AAT. However, unlike classic inverse iteration, the new method works with matrices Bk∈R(n+1)×(n+1) obtained by bordering A in such a way that the Bk have uniformly bounded condition numbers. This allows using iterative Krylov-type solvers for large problems. If σn−1>σn the singular vector approximations convergence linearly with factor κ=σn/σn−1