Log In Sign Up

Quantum algorithms for matrix scaling and matrix balancing

by   Joran van Apeldoorn, et al.

Matrix scaling and matrix balancing are two basic linear-algebraic problems with a wide variety of applications, such as approximating the permanent, and pre-conditioning linear systems to make them more numerically stable. We study the power and limitations of quantum algorithms for these problems. We provide quantum implementations of two classical (in both senses of the word) methods: Sinkhorn's algorithm for matrix scaling and Osborne's algorithm for matrix balancing. Using amplitude estimation as our main tool, our quantum implementations both run in time Õ(√(mn)/ε^4) for scaling or balancing an n × n matrix (given by an oracle) with m non-zero entries to within ℓ_1-error ε. Their classical analogs use time Õ(m/ε^2), and every classical algorithm for scaling or balancing with small constant ε requires Ω(m) queries to the entries of the input matrix. We thus achieve a polynomial speed-up in terms of n, at the expense of a worse polynomial dependence on the obtained ℓ_1-error ε. We emphasize that even for constant ε these problems are already non-trivial (and relevant in applications). Along the way, we extend the classical analysis of Sinkhorn's and Osborne's algorithm to allow for errors in the computation of marginals. We also adapt an improved analysis of Sinkhorn's algorithm for entrywise-positive matrices to the ℓ_1-setting, leading to an Õ(n^1.5/ε^3)-time quantum algorithm for ε-ℓ_1-scaling in this case. We also prove a lower bound, showing that our quantum algorithm for matrix scaling is essentially optimal for constant ε: every quantum algorithm for matrix scaling that achieves a constant ℓ_1-error with respect to uniform marginals needs to make at least Ω(√(mn)) queries.


page 1

page 2

page 3

page 4


Improved quantum lower and upper bounds for matrix scaling

Matrix scaling is a simple to state, yet widely applicable linear-algebr...

Information geometry of operator scaling

Matrix scaling is a classical problem with a wide range of applications....

DeepScaleTool : A Tool for the Accurate Estimation of Technology Scaling in the Deep-Submicron Era

The estimation of classical CMOS "constant-field" or "Dennard" scaling m...

Interior-point methods for unconstrained geometric programming and scaling problems

We provide a condition-based analysis of two interior-point methods for ...

Spectral analysis of matrix scaling and operator scaling

We present a spectral analysis for matrix scaling and operator scaling. ...

Incomplete selected inversion for linear-scaling electronic structure calculations

Pole Expansion and Selected Inversion (PEXSI) is an efficient scheme for...

The Vector Balancing Constant for Zonotopes

The vector balancing constant vb(K,Q) of two symmetric convex bodies K,Q...