2025-11-10T02:59:56.765376

GUPPY: Pythonic Quantum-Classical Programming

Koch, Lawrence, Singhal et al.

We present ongoing work on Guppy, a domain-specific language embedded in Python that allows users to write high-level hybrid quantum programs with complex control flow in Pythonic syntax, aiming to run them on actual quantum hardware.

academic

GUPPY: Pythonic Quantum-Classical Programming

基本信息

论文ID: 2510.12582
标题: GUPPY: Pythonic Quantum-Classical Programming
作者: Mark Koch, Alan Lawrence, Kartik Singhal, Seyon Sivarajah, Ross Duncan (Quantinuum, United Kingdom)
分类: cs.PL (Programming Languages), cs.SE (Software Engineering), quant-ph (Quantum Physics)
发表时间: 2025年10月14日 (arXiv预印本)
论文链接: https://arxiv.org/abs/2510.12582v1

摘要

本文介绍了Guppy，一个嵌入在Python中的领域特定语言(DSL)，允许用户使用Pythonic语法编写具有复杂控制流的高级混合量子程序，目标是在实际量子硬件上运行这些程序。

研究背景与动机

问题定义

Python在量子计算领域的主导地位：根据最新调研，超过90%的量子计算研究人员和从业者使用Python，但现有框架受限于电路级抽象，对高级控制流支持有限。
实时量子-经典计算需求：重复直到成功(repeat-until-success)协议等算法需要基于量子测量结果的实时经典数据处理，这种需求随着下一代量子设备的发展将进一步增长。
现有方法局限性：
- 大多数Python框架通过跟踪Python解释器来构建量子程序表示，无法捕获Python原生控制流语法
- 条件语句和循环必须通过高阶组合子或其他语法结构表达
- 缺乏线性类型系统来防止量子编程错误

研究动机

开发一种既保持Python友好语法和生态系统优势，又能支持复杂量子-经典混合计算的编程语言。

核心贡献

提出了Guppy DSL：一个嵌入Python的领域特定语言，支持Pythonic语法编写量子-经典混合程序
实现了线性类型系统：通过线性Qubit类型在编译时捕获量子编程错误（如违反无克隆定理）
支持复杂控制流：原生支持基于测量结果的条件语句、循环等控制流结构
静态编译到Hugr IR：编译到新型量子中间表示Hugr，支持量子-经典程序的表达和优化
Python互操作性：通过py(...)表达式实现与Python生态系统的无缝集成

方法详解

任务定义

设计一个编程语言，能够：

使用熟悉的Python语法
支持量子-经典混合计算
提供编译时安全保证
生成可在量子硬件上运行的程序

语言设计架构

1. 核心语法特性

Pythonic控制流：支持if、for、while、break、continue、return语句
强类型系统：要求类型注解，变量必须在所有控制流路径中被赋值且具有唯一静态类型
基础类型：支持bool、int、float、None及相应操作

2. 量子操作设计

@guppy
def rx(q: Qubit, a: float) -> Qubit:
    # 线性类型：Qubit不能被复制或隐式丢弃
    return h(rz(h(q), a))

线性类型系统特点：

Qubit值被线性处理，不能复制或隐式丢弃
量子操作签名如：h: Qubit -> Qubit

编译时检测编程错误：

q = Qubit()
return cx(q, q)  # 错误：多次使用同一qubit

q = Qubit()
h(q)  # 错误：返回值未使用

3. 集合类型

线性列表和元组：用于表示qubit寄存器
特殊操作：
- get函数：返回指定索引的qubit和移除该qubit后的新列表
- apply方法：在指定索引上应用函数，如qs.apply(cx, (i, j))
列表推导式：支持[h(q) for q in qs]语法

4. 高阶函数支持

函数是一等值，使用Callable类型构造器
支持嵌套函数定义和变量捕获
支持递归和相互递归

5. Python互操作性

import networkx as nx
g = nx.erdos_renyi_graph(20, 0.25)

@guppy
def apply_graph(qs: list[Qubit]) -> list[Qubit]:
    for i, j in py(g.edges()):  # py()表达式调用Python代码
        qs = qs.apply(zz, (i, j))
    return qs

编译架构

编译目标：Hugr IR

分层统一图表示(Hierarchical Unified Graph Representation)
数据流图：节点表示纯函数调用，边表示数据依赖
结构化控制流：通过Conditional和TailLoop节点表示分支和迭代
线性类型支持：约束节点输出只能有一条数据流边

编译流程

Guppy源码解析
类型检查和静态分析
生成Hugr图表示
优化和目标代码生成

实验设置

示例程序验证

1. 量子隐态传输协议

@guppy
def teleport(src: Qubit, tgt: Qubit) -> Qubit:
    tmp, tgt = cx(h(Qubit()), tgt)
    src, tmp = cx(src, tmp)
    
    if measure(h(src)):
        tgt = z(tgt)
    if measure(tmp):
        tgt = x(tgt)
    return tgt

2. 重复直到成功协议

@guppy
def rus(q: Qubit, tries: int) -> Qubit:
    for _ in range(tries):
        a, b = h(Qubit()), h(Qubit())
        b, a = cx(b, tdg(a))
        
        if not measure(t(a)):
            discard(b)
            continue
            
        q, b = cx(z(t(q)), b)
        if measure(t(b)):
            break
        q = x(q)
    return q