A QSVT-Based Quantum Jacobi Algorithm for Linear Systems with Application to the Poisson Equation
Abstract
Many computational fluid dynamics (CFD) algorithms solve partial differential equations by discretization, resulting in large and sparse systems of linear equations. While iterative methods are widely used to solve these systems classically, most existing quantum linear system solvers target the solution through matrix inversion rather than approximating it using an iterative procedure. In this work, we develop a quantum implementation of the Jacobi method based on quantum singular value transformation (QSVT). By reformulating the Jacobi iteration as a polynomial transformation of a block-encoded operator, the algorithm requires only a constant ancilla overhead with respect to the number of iterations while maintaining a circuit depth that scales linearly with the iteration number. We demonstrate the algorithm for one- and two-dimensional Poisson problems, including the pressure Poisson equation arising in Chorin's projection method for the lid-driven cavity flow. The proposed algorithm provides a promising building block for future quantum implementations of multigrid methods and preconditioning techniques, bringing quantum algorithms closer to established CFD solution strategies.
Related
- Application a quantum machine learning algorithm to the applied problem of modeling the thermal regime of a servers
- Gorini-Kossakowski-Sudarshan-Lindblad equation in different bases: application to driven-dissipative two- and multilevel systems
- Quantum-Classical Iterative Refinement for the Pressure-Poisson Equation on NISQ Hardware
- Quantum HHL Algorithm for Solving 1D Linear Advection Equation Using Explicit and Implicit Finite Difference Schemes
- An Optimal Quantum Linear Systems Algorithm
BibTeX
@misc{piskol2026qsvtbased,
title = {{A QSVT-Based Quantum Jacobi Algorithm for Linear Systems with Application to the Poisson Equation}},
author = {Louisa M. Piskol and Thorsten Grahs and Stefan Langer and Oleksandr Kyriienko},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.24266},
eprint = {2609.24266},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.24266},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - A QSVT-Based Quantum Jacobi Algorithm for Linear Systems with Application to the Poisson Equation
AU - Louisa M. Piskol
AU - Thorsten Grahs
AU - Stefan Langer
AU - Oleksandr Kyriienko
PY - 2026
DA - 2026/09/21
PB - arXiv
DO - 10.48550/arxiv.2609.24266
UR - https://arxiv.org/abs/2609.24266
AB - Many computational fluid dynamics (CFD) algorithms solve partial differential equations by discretization, resulting in large and sparse systems of linear equations. While iterative methods are widely used to solve these systems classically, most existing quantum linear system solvers target the solution through matrix inversion rather than approximating it using an iterative procedure. In this work, we develop a quantum implementation of the Jacobi method based on quantum singular value transformation (QSVT). By reformulating the Jacobi iteration as a polynomial transformation of a block-encoded operator, the algorithm requires only a constant ancilla overhead with respect to the number of iterations while maintaining a circuit depth that scales linearly with the iteration number. We demonstrate the algorithm for one- and two-dimensional Poisson problems, including the pressure Poisson equation arising in Chorin's projection method for the lid-driven cavity flow. The proposed algorithm provides a promising building block for future quantum implementations of multigrid methods and preconditioning techniques, bringing quantum algorithms closer to established CFD solution strategies.
ER -