量子程序匹配拓扑结构

量子计算设备存在量子比特之间的有限连接,使得只能在有限的量子位对上应用两个量子位门。 量子程序应用到目标设备时,必须转换原始的量子程序以适应硬件限制,让双量子比特门中的两个量子比特能够满足物理拓扑结构,从而让双量子位门正常作用; 当前解决方案中多数需要在无法相互作用的两个量子比特间插入额外的SWAP操作,以便将逻辑量子位“移动”到它们可以相互作用的位置。 我们称这种解决方法为量子程序匹配拓扑结构。

接口说明

当前QPanda2版本中存在两种思路的匹配拓扑方法:

  • 接口 QProg topology_match(QProg, QVec &, QuantumMachine *, SwapQubitsMethod, ArchType)

    通过采用线路分层以及A*搜索算法,在匹配过程中,让插入的SWAP操作个数近似达到最少,使得算法的整体近似消耗达到最少。 该接口需要传入5个参数,其中分别为 构建的量子程序、使用的量子比特位集合、初始化的虚拟机指针、使用的SWAP操作的方式、拓扑结构的类型; 并且返回映射后的量子程序

  • 接口 QProg qcodar_match(QProg, QVec &, QuantumMachine *, QCodarGridDevice, size_t, size_t, size_t)

    通过采用一种上下文敏感和持续时间感知的重映射算法,该算法能够感知门持续时间差和程序上下文,使得它能够从程序中提取更多的并行性, 在不同体系结构的模拟中平均将量子程序的速度提高1.23。并在原始量子噪声模拟器上运行时保持电路的保真度。 该接口需要传入7个参数,其中分别为 构建的量子程序、使用的量子比特位集合、初始化的虚拟机指针、拓扑结构的类型、拓扑结构的宽度、拓扑结构的长度、匹配拓扑运行的次数; 并且返回映射后的量子程序

实例

topology_match 实例

#include "QPanda.h"
using namespace std;
using namespace QPanda;
int  main()
{
    auto qvm = new CPUQVM();
    qvm->init();
    auto q = qvm->qAllocMany(8);
    auto c = qvm->cAllocMany(8);
    auto srcprog = QProg();

    // 构建量子程序
    srcprog << CNOT(q[0], q[3])
        << CNOT(q[0], q[2])
        << CNOT(q[1], q[3])
        << CZ(q[1], q[2])
        << CZ(q[0], q[2])
        << T(q[1])
        << S(q[2])
        << H(q[3]);

    // 对srcprog进行概率测量,得到结果r1
    qvm->directlyRun(srcprog);
    auto r1 = qvm->pMeasureNoIndex(q);

    // 对srcprog进行拓扑匹配,得到匹配ORIGIN_VIRTUAL_ARCH拓扑结构的量子程序outprog
    QProg  outprog = topology_match(srcprog, q, qvm);

    // 对outprog进行概率测量,得到结果r2
    qvm->directlyRun(outprog);
    auto r2 = qvm->pMeasureNoIndex(q);

    // 对比概率测量结果r1和r2
    int size = std::min(r1.size(), r2.size());
    bool result_equal = true;
    for (int i = 0; i < size; i++)
    {
        if ((fabs(r1[i] - r2[i]) > 1e-6))
            result_equal = false;
    }

    // 如果结果相同,打印结果
    if (result_equal == true)
    {
        std::cout << "The probability measurements are the same, prob list:  " << std::endl;
        for (int i = 0; i < size; i++)
        {
            std::cout << r1[i]  << std::endl;
        }
    }

    qvm->finalize();
    delete qvm;
    return 0 ;
}

具体步骤如下:

  • 首先创建量子虚拟机、量子寄存器、经典寄存器

  • 编写量子程序 srcprog ,对该量子程序进行概率测量得到结果 r1

  • 接着调用 topology_match()srcprog 进行符合特定物理结构的线路映射,得到适配特定物理结构的量子程序 outprog

  • 对量子程序 outprog 进行概率测量得到结果 r2

  • 由于量子程序映射只是对原线路增加额外的 SWAP 操作,以适配物理拓扑结构,并不影响线路的结构。所以对比结果 r1r2 ,如果结果一致,则线路映射正确。

运行结果如下:

The probability measurements are the same, prob list:
0.5
0
0
0
0
0
0
0
0.5
0
0
0
0
0
0
0

qcodar_match 实例

#include "Core/Core.h"
using namespace std;
using namespace QPanda;
int main()
{
    auto  qvm = new CPUQVM();
    qvm->init();
    auto q = qvm->qAllocMany(4);
    auto cv = qvm->cAllocMany(4);
    QProg srcprog;

    // 构建量子程序
    srcprog << CNOT(q[1], q[3])
        << RX(q[0], PI / 2)
        << CNOT(q[0], q[2])
        << CNOT(q[1], q[3])
        << RY(q[1], -PI / 4)
        << CNOT(q[2], q[0])
        << CZ(q[1], q[2])
        << CNOT(q[1], q[3])
        << RZ(q[2], PI / 6)
        << CNOT(q[2], q[0])
        << RZ(q[0], -PI / 4)
        << CNOT(q[0], q[2])
        << H(q[0])
        << T(q[1])
        << RX(q[1], -PI/4)
        << Y(q[2])
        << Z(q[1])
        ;

    // 对srcprog进行概率测量,得到结果r1
    qvm->directlyRun(srcprog);
    auto r1 = qvm->pMeasureNoIndex(q);

    // 对srcprog进行拓扑匹配,得到匹配SIMPLE_TYPE拓扑结构的量子程序outprog
    QProg outprog = qcodar_match_by_simple_type(srcprog, q, qvm, 2, 3, 5);

    // 对outprog进行概率测量,得到结果r2
    qvm->directlyRun(outprog);
    auto r2 = qvm->pMeasureNoIndex(q);

    // 对比概率测量结果r1和r2
    int size = std::min(r1.size(), r2.size());
    bool result_equal = true;
    for (int i = 0; i < size; i++)
    {
        if ((fabs(r1[i] - r2[i]) > 1e-6))
            result_equal = false;
    }

    // 如果结果相同,打印结果
    if (result_equal == true)
    {
        std::cout << "The probability measurements are the same, prob list:  " << std::endl;
        for (int i = 0; i < size; i++)
        {
            std::cout << r1[i] << std::endl;
        }
    }

    qvm->finalize();
    delete qvm;
    return 0;
}

具体步骤如下:

  • 首先创建量子虚拟机、量子寄存器、经典寄存器

  • 编写量子程序 srcprog ,对该量子程序进行概率测量得到结果 r1

  • 接着调用 qcodar_match()srcprog 进行符合特定物理结构的线路映射,得到适配特定物理结构的量子程序 outprog

  • 对量子程序 outprog 进行概率测量得到结果 r2

  • 由于量子程序映射只是对原线路增加额外的 SWAP 操作,以适配物理拓扑结构,并不影响线路的结构。所以对比结果 r1r2 ,如果结果一致,则线路映射正确。

运行结果如下:

The probability measurements are the same, prob list:
0
0
0
0
0.269995
0.458558
0.0463238
0.0786762
0
0
0
0
0.0134987
0.00794791
0.0786762
0.0463238

备注

运行量子程序匹配拓扑结构接口时,需要在本地运行目录下配置 QPandaConfig.json 文件,在配置文件中的 QuantumChipArch 字段对拓扑结构进行配置修改。