Saved in:
Bibliographic Details
Main Authors: Lefterovici, Andreea-Iulia, Perk, Michael, Ramacciotti, Debora, Rotundo, Antonio F., Skelton, S. E., Steinbach, Martin
Format: Preprint
Published: 2025
Subjects:
Online Access:https://arxiv.org/abs/2503.21420
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866911523736125440
author Lefterovici, Andreea-Iulia
Perk, Michael
Ramacciotti, Debora
Rotundo, Antonio F.
Skelton, S. E.
Steinbach, Martin
author_facet Lefterovici, Andreea-Iulia
Perk, Michael
Ramacciotti, Debora
Rotundo, Antonio F.
Skelton, S. E.
Steinbach, Martin
contents Solving systems of linear equations is a key subroutine in many quantum algorithms. In the last 15 years, many quantum linear solvers (QLS) have been developed, competing to achieve the best asymptotic worst-case complexity. Most QLS assume fault-tolerant quantum computers, so they cannot yet be benchmarked on real hardware. Because an algorithm with better asymptotic scaling can underperform on instances of practical interest, the question of which of these algorithms is the most promising remains open. In this work, we implement a method to partially address this question. We consider four well-known QLS algorithms which directly implement an approximate matrix inversion function: the Harrow-Hassidim-Lloyd algorithm, two algorithms utilizing a linear combination of unitaries, and one utilizing the quantum singular value transformation (QSVT). These methods, known as functional QLS, share nearly identical assumptions about the problem setup and oracle access. Their computational cost is dominated by query calls to a matrix oracle encoding the problem one wants to solve. We provide formulas to count the number of queries needed to solve specific problem instances; these can be used to benchmark the algorithms on real-life instances without access to quantum hardware. We select three data sets: random generated instances that obey the assumptions of functional QLS, linear systems from simplex iterations on MIPLIB, and Poisson equations. Our methods can be easily extended to other data sets and provide a high-level guide to evaluate the performance of a QLS algorithm. In particular, our work shows that HHL underperforms in comparison to the other methods across all data sets, often by orders of magnitude, while the QSVT-based method shows the best performance.
format Preprint
id arxiv_https___arxiv_org_abs_2503_21420
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Beyond asymptotic scaling: Comparing functional quantum linear solvers
Lefterovici, Andreea-Iulia
Perk, Michael
Ramacciotti, Debora
Rotundo, Antonio F.
Skelton, S. E.
Steinbach, Martin
Quantum Physics
Solving systems of linear equations is a key subroutine in many quantum algorithms. In the last 15 years, many quantum linear solvers (QLS) have been developed, competing to achieve the best asymptotic worst-case complexity. Most QLS assume fault-tolerant quantum computers, so they cannot yet be benchmarked on real hardware. Because an algorithm with better asymptotic scaling can underperform on instances of practical interest, the question of which of these algorithms is the most promising remains open. In this work, we implement a method to partially address this question. We consider four well-known QLS algorithms which directly implement an approximate matrix inversion function: the Harrow-Hassidim-Lloyd algorithm, two algorithms utilizing a linear combination of unitaries, and one utilizing the quantum singular value transformation (QSVT). These methods, known as functional QLS, share nearly identical assumptions about the problem setup and oracle access. Their computational cost is dominated by query calls to a matrix oracle encoding the problem one wants to solve. We provide formulas to count the number of queries needed to solve specific problem instances; these can be used to benchmark the algorithms on real-life instances without access to quantum hardware. We select three data sets: random generated instances that obey the assumptions of functional QLS, linear systems from simplex iterations on MIPLIB, and Poisson equations. Our methods can be easily extended to other data sets and provide a high-level guide to evaluate the performance of a QLS algorithm. In particular, our work shows that HHL underperforms in comparison to the other methods across all data sets, often by orders of magnitude, while the QSVT-based method shows the best performance.
title Beyond asymptotic scaling: Comparing functional quantum linear solvers
topic Quantum Physics
url https://arxiv.org/abs/2503.21420