Abstract:We propose fast and practical quantum-inspired classical algorithms for solving linear systems. Specifically, given sampling and query access to a matrix $A\in\mathbb{R}^{m\times n}$ and a vector $b\in\mathbb{R}^m$, we propose classical algorithms that produce a data structure for the solution $x\in\mathbb{R}^{n}$ of the linear system $Ax=b$ with the ability to sample and query its entries. The resulting $x$ satisfies $\|x-A^{+}b\|\leq\epsilon\|A^{+}b\|$, where $\|\cdot\|$ is the spectral norm and $A^+$ is the Moore-Penrose inverse of $A$. Our algorithm has time complexity $\widetilde{O}(\kappa_F^4/\kappa\epsilon^2)$ in the general case, where $\kappa_{F} =\|A\|_F\|A^+\|$ and $\kappa=\|A\|\|A^+\|$ are condition numbers. Compared to the prior state-of-the-art result [Shao and Montanaro, arXiv:2103.10309v2], our algorithm achieves a polynomial speedup in condition numbers. When $A$ is $s$-sparse, our algorithm has complexity $\widetilde{O}(s \kappa\log(1/\epsilon))$, matching the quantum lower bound for solving linear systems in $\kappa$ and $1/\epsilon$ up to poly-logarithmic factors [Harrow and Kothari]. When $A$ is $s$-sparse and symmetric positive-definite, our algorithm has complexity $\widetilde{O}(s\sqrt{\kappa}\log(1/\epsilon))$. Technically, our main contribution is the application of the heavy ball momentum method to quantum-inspired classical algorithms for solving linear systems, where we propose two new methods with speedups: quantum-inspired Kaczmarz method with momentum and quantum-inspired coordinate descent method with momentum. Their analysis exploits careful decomposition of the momentum transition matrix and the application of novel spectral norm concentration bounds for independent random matrices. Finally, we also conduct numerical experiments for our algorithms on both synthetic and real-world datasets, and the experimental results support our theoretical claims.

Optimised Hybrid Classical-Quantum Algorithm for Accelerated Solution of Sparse Linear Systems

Solving Linear Systems on Quantum Hardware with Hybrid HHL++

Integrating Quantum Algorithms Into Classical Frameworks: A Predictor-corrector Approach Using HHL

Fast and Practical Quantum-Inspired Classical Algorithms for Solving Linear Systems.

The HHL algorithm: Implementation and research directions

Enhancing the Harrow-Hassidim-Lloyd (HHL) algorithm in systems with large condition numbers

Hybrid quantum-classical and quantum-inspired classical algorithms for solving banded circulant linear systems

A Hybrid Quantum-Classical CFD Methodology with Benchmark HHL Solutions

An Enhanced Hybrid HHL Algorithm

A near-term quantum algorithm for solving linear systems of equations based on the Woodbury identity

Quantum Krylov-Subspace Method Based Linear Solver

Adapting the HHL algorithm to quantum many-body theory

Quantum Algorithm to Cubic Spline Interpolation

An Exponential Separation Between Quantum and Quantum-Inspired Classical Algorithms for Machine Learning

Scalable Quantum Computation of Highly Excited Eigenstates with Spectral Transforms

Quantum linear systems algorithms: a primer

An Efficient Quantum Algorithm for Linear System Problem in Tensor Format

A Hybrid Quantum Algorithm for Load Flow

Distributed hybrid quantum-classical performance prediction for hyperparameter optimization

A Catalyst Framework for the Quantum Linear System Problem via the Proximal Point Algorithm

Quantum computing and preconditioners for hydrological linear systems