Strictly balancing matrices in polynomial time using osborne's iteration

Osborne's iteration is a method for balancing n × n matrices which is widely used in linear algebra packages, as balancing preserves eigenvalues and stabilizes their numeral computation. The iteration can be implemented in any norm over Rⁿ, but it is normally used in the L₂ norm. The choice of norm not only a ects the desired balance condition, but also defines the iterated balancing step itself. In this paper we focus on Osborne's iteration in any L_p norm, where p < ∞. We design a specific implementation of Osborne's iteration in any L_p norm that converges to a strictly - balanced matrix in O(^{− 2}n⁹K) iterations, where K measures, roughly, the number of bits required to represent the entries of the input matrix. This is the first result that proves a variant of Osborne's iteration in the L₂ norm (or any L_p norm, p < ∞) strictly balances matrices in polynomial time. This is a substantial improvement over our recent result (in SODA 2017) that showed weak balancing in L_p norms. Previously, Schulman and Sinclair (STOC 2015) showed strict balancing of another variant of Osborne's iteration in the L_∞ norm. Their result does not imply any bounds on strict balancing in other norms.