量子迭代法

经典 CFD 从不直接对大矩阵求逆——内存与代价决定了工业求解器是松弛迭代 + 多重网格的世界。Williams 等（Exeter/PASQAL/Siemens，2024）把这一实践哲学带进量子算法：用量子 Jacobi 与量子 Gauss-Seidel 的块编码—LCU 实现逐步松弛，不需要一次性的大矩阵求逆，在粘性 Burgers 与 Euler 激波管基准上验证收敛Williams et al. 2024 (arXiv)。其资源账对早期容错设备友好，并与量子多重网格构成嵌套愿景。

为什么迭代、不直接求逆

论文开宗明义的比较：现有量子微分方程求解器多依赖（直接或变相的）全矩阵求逆，而这在 CFD 中是非典型做法——受内存与计算约束，经典求解器用基于松弛的迭代法逼近解Williams et al. 2024 (arXiv)。更深一层，QLSA 无法利用初始猜测迭代改进（量子多重网格论文同样强调这一点，见quantum-multigrid），而迭代法天然支持”从上一时间步的解出发暖启动”——这恰是时间步进 CFD 的日常形态。

量子 Jacobi：分裂 + LCU 编码迭代轨迹

对角占优时 $A=D+R$（$D$ 对角、$R$ 余项），Jacobi 迭代写为 $x_{k+1}=x_0-D^{-1}R{x}_k$、$x_0=D^{-1}b$。展开成 $k+1$ 项的有限和：

$$x_k=\sum _{l=0}^k(-D^{-1}R{)}^l{x}_0.$$

量子实现用 LCU 把全部迭代轨迹编码进一个量子态：$a=\lceil {\log }_2(k+1)\rceil$ 个控制辅助比特选择阶数 $l$，非 2 幂的项数由归一化酉 $V$ 处理，$|x_k⟩=\mathcal{U}|\varnothing ⟩$——迭代历史以叠加形式存进态里，末态即第 $k$ 次迭代解Williams et al. 2024 (arXiv)。矩阵件全部块编码化：$D^{-1}$、$R$、向量 $b$、初值 $x_0$ 各有专用块编码（Givens 旋转积实现）；块编码相乘时用延拓寄存器使各块部分重叠，消去直接乘法引入的非对角交叉项。

图:量子迭代方法的完整算法管线——问题矩阵与向量块编码、LCU 轨迹求和、末态测量读出(来源:QC007)

图:量子 Jacobi 的电路组件——(a) 块编码乘法(延拓寄存器消除交叉项)与 LCU 结构,块编码由 Givens 旋转积实现(来源:QC007)

量子 Gauss-Seidel：Woodbury 恒等式绕开串行性

Gauss-Seidel 的逐分量顺序更新（下三角即时回代）本质串行，无法照搬成并行量子算符。Williams 等的改写：把 $R=B+T$ 分裂后经 Woodbury 恒等式把更新算符化为截断几何级数

$$\Omega \stackrel{\mathrm{def}}{=}\sum _{l=0}^L(-D^{-1}B{)}^l,$$

以有限项和近似顺序依赖——串行性被级数截断取代，重新落回 LCU 可编码的形状Williams et al. 2024 (arXiv)。

资源账

论文给出的显式账本（$k$ 为迭代数、$C$ 为所选块编码协议的成本，依赖稀疏度 $s$ 与矩阵元素重复度 $z$）：

• 宽度：${\log }_{2}N+2k+\lceil {\log }_2(k+1)\rceil$ 个量子比特；换用 QSP 乘法后压到 ${\log }_{2}N+4+\lceil {\log }_2(k+1)\rceil$（对 $k$ 只剩对数依赖）；
• 深度：$O(k^{2}C)$——主导项是 $O(k^2)$ 次块编码矩阵乘法；
• 内存：解存为量子态，$O({\log }_{2}N)$ 对经典迭代法 $O(N)$ 的数据操作内存是指数节省；收敛所需的迭代次数与经典实现相同，数值上 $K\sim O(\kappa )$（256×256 三对角测试，$\kappa$ 从 3 到 70 时迭代数线性增长，且单次迭代成本与 $\kappa$ 无关）Williams et al. 2024 (arXiv)。

读出沿用泛函/矩路线：测 $⟨x_k|M|x_k⟩$（统计矩、特定格点值）或 shadow tomography；迭代历史在态叠加里，完整提取仍受测量瓶颈约束（见io-problem）。基准：粘性 Burgers 方程 80 次迭代后的行波解、Euler 方程激波管/声传播问题均得到高质量解；作者特别指出该求解器适合含间断/不稳定的解，并优于需要变分态制备的技术Williams et al. 2024 (arXiv)。

生态位：早期容错设备与多重网格嵌套

与 HHL/QSVT 路线相比，量子迭代法的深度账以 $k^2$（而非 ${\kappa }^2\log N/ϵ$ 型）为主，且无相位估计的长相干链——定位在早期容错设备上先跑起来的工程路线。更大的愿景是作为光滑子的角色嵌入量子多重网格的 V-cycle（见quantum-multigrid：量子 V-cycle 的平滑步正是 Gauss-Seidel/Richardson 型松弛），迭代法量子化与多重网格量子化因此是同一蓝图的两个层次；预条件思想则在块编码框架下另有系统展开（见block-encoding、quantum-multigrid）。

延伸阅读

• quantum-multigrid
• block-encoding
• hhl
• quantum-burgers
• multigrid