2. 程式人生 > >量子計算 OpenFermion-Cirq 開源python庫

量子計算 OpenFermion-Cirq 開源python庫

阿新 • • 發佈：2019-01-01

#01 基礎資訊
NISQ，是嘈雜中型量子 (Noisy Intermediate-Scale Quantum) 的簡稱。
擁有50-100量子位元、以及高保真量子門 (Quantum Gate) 的計算機，便可稱為NISQ計算機。
谷歌開源了Cirq框架，這是專為NISQ演算法打造的框架。Cirq主要用來短期 (Near-Term) 問題，
幫助研究者瞭解NISQ計算機到底能不能解決實際應用的中的計算問題。這個框架，經Apache 2.0協議許可，可以修改，可以嵌入任何開源/付費的軟體包。

框架安裝好之後，開發者就可以給特定的量子處理器，編寫量子演算法，使用者可以精確控制量子電路 (Quantum Circuits) ；更加充分地利用NISQ架構。

Cirq支援在模擬器上執行演算法，如果將來有了量子計算機，或者更大的模擬器，也很容易通過雲，把裝置和演算法整合起來。

與Cirq框架一同釋出的，還有OpenFermion-Cirq，這是一個基於Cirq的應用示例。

OpenFermion是一個量子演算法開發平臺，專注解決化學問題。

OpenFermion-Cirq則是一個開源庫，把量子模擬演算法編譯成Cirq能用的樣子。比如，這個新庫，可以用來搭建量子變分演算法 (Variational Algorithms) ，模擬分子或者複雜材料的性質。

git clone https://github.com/quantumlib/OpenFermion-Cirq
cd OpenFermion-Cirq
python -m pip install -e .

##02 安裝
我在win10x64下，基於python2.7安裝沒成功。
後來基於python3.5.3安裝成功。

D:\git\Quantum\OpenFermion-Cirq>python -m pip install -e .
Obtaining file:///D:/git/Quantum/OpenFermion-Cirq
Unknown requires Python '>=3.5.3' but the running Python is 2.7.15
You are using pip version 9.0.3, however version 18.0 is available.
You should consider upgrading via the 'python -m pip install --upgrade pip' command.

D:\git\Quantum\OpenFermion-Cirq>set path=C:\Python35;%path%

D:\git\Quantum\OpenFermion-Cirq>python -m pip install -e .
Obtaining file:///D:/git/Quantum/OpenFermion-Cirq
Collecting cirq==0.3.1.* (from openfermioncirq==0.2.1a35)
  Downloading https://files.pythonhosted.org/packages/dc/c9/f17b365c7ac43c6d9576d44cf894c1c42bf43ac54ef6d76f56a99b32870a/cirq-0.3.1.35-py3-none-any.whl (368kB)
    100% |████████████████████████████████| 378kB 416kB/s
Requirement already satisfied: numpy in c:\python35\lib\site-packages (from openfermioncirq==0.2.1a35)
...
...
  Found existing installation: openfermioncirq 0.2.1a27
    Can't uninstall 'openfermioncirq'. No files were found to uninstall.
  Running setup.py develop for openfermioncirq
Successfully installed Send2Trash-1.5.0 cachetools-2.1.0 certifi-2018.8.24 chardet-3.0.4 cirq-0.3.1.35 google-api-python-client-1.7.4 google-auth-1.5.1 google-auth-httplib2-0.0.3 h5py-2.8.0 httplib2-0.11.3 ipywidgets-7.4.1 jupyter-1.0.0 jupyter-client-5.2.3 jupyter-console-5.2.0 kiwisolver-1.0.1 matplotlib-2.2.3 notebook-5.6.0 openfermion-0.8.1 openfermioncirq prometheus-client-0.3.1 protobuf-3.6.1 pubchempy-1.0.4 pyasn1-modules-0.2.2 pywinpty-0.5.4 pyzmq-17.1.2 requests-2.19.1 rsa-3.4.2 sortedcontainers-2.0.5 terminado-0.8.1 uritemplate-3.0.0 urllib3-1.23 widgetsnbextension-3.4.1
You are using pip version 9.0.1, however version 18.0 is available.
You should consider upgrading via the 'python -m pip install --upgrade pip' command.

D:\git\Quantum\OpenFermion-Cirq>

另外需要 pytest

D:\git\Quantum\OpenFermion-Cirq>pip3 install pytest

##03 執行用例

直接瀏覽器開啟，拷貝程式碼到PyCharm。配置中設定為自己的python3.5.3路徑。執行結果如下：
這裡寫圖片描述

print('01 ############################## start')
import openfermion
import numpy

# Set the number of qubits in our example.
n_qubits = 3
simulation_time = 1.
random_seed = 8317

# Generate the random one-body operator.
T = openfermion.random_hermitian_matrix(n_qubits, seed=random_seed)

# Diagonalize T and obtain basis transformation matrix (aka "u").
eigenvalues, eigenvectors = numpy.linalg.eigh(T)
basis_transformation_matrix = eigenvectors.transpose()

# Print out familiar OpenFermion "FermionOperator" form of H.
H = openfermion.FermionOperator()
for p in range(n_qubits):
    for q in range(n_qubits):
        term = ((p, 1), (q, 0))
        H += openfermion.FermionOperator(term, T[p, q])
print(H)
print('01 ############################## end')

print('02 ############################## start')
import openfermioncirq
import cirq

# Initialize the qubit register.
qubits = cirq.LineQubit.range(n_qubits)

# Start circuit with the inverse basis rotation, print out this step.
inverse_basis_rotation = cirq.inverse(openfermioncirq.bogoliubov_transform(qubits, basis_transformation_matrix))
circuit = cirq.Circuit.from_ops(inverse_basis_rotation)
print(circuit)
print('02 ############################## end')

print('03 ############################## start')
# Add diagonal phase rotations to circuit.
for k, eigenvalue in enumerate(eigenvalues):
    phase = -eigenvalue * simulation_time
    circuit.append(cirq.RotZGate(rads=phase).on(qubits[k]))

# Finally, restore basis.
basis_rotation = openfermioncirq.bogoliubov_transform(qubits, basis_transformation_matrix)
circuit.append(basis_rotation)


# Initialize a random initial state.
initial_state = openfermion.haar_random_vector(
    2 ** n_qubits, random_seed).astype(numpy.complex64)

# Numerically compute the correct circuit output.
import scipy
hamiltonian_sparse = openfermion.get_sparse_operator(H)
exact_state = scipy.sparse.linalg.expm_multiply(
    -1j * simulation_time * hamiltonian_sparse, initial_state)

# Use Cirq simulator to apply circuit.
simulator = cirq.google.XmonSimulator()
result = simulator.simulate(circuit, qubit_order=qubits,
                            initial_state=initial_state)
simulated_state = result.final_state

# Print final fidelity.
fidelity = abs(numpy.dot(simulated_state, numpy.conjugate(exact_state)))**2
print(fidelity)


xmon_circuit = cirq.google.optimized_for_xmon(circuit)
print(xmon_circuit.to_qasm())
print('03 ############################## end')

print('04 ############################## start')
xmon_circuit = cirq.google.optimized_for_xmon(circuit)
print(xmon_circuit.to_qasm())
print('04 ############################## end')

執行結果:

C:\Python35\python.exe D:/git/Quantum/OpenFermion-Cirq/examples/test01.py
01 ############################## start
(0.5367212624097257+0j) [0^ 0] +
(-0.26033703159240107+3.3259173741375454j) [0^ 1] +
(1.3433603748462144+1.544987250567917j) [0^ 2] +
(-0.26033703159240107-3.3259173741375454j) [1^ 0] +
(-2.9143303700812435+0j) [1^ 1] +
(-1.52843836446248+1.3527486791390022j) [1^ 2] +
(1.3433603748462144-1.544987250567917j) [2^ 0] +
(-1.52843836446248-1.3527486791390022j) [2^ 1] +
(2.261633626116526+0j) [2^ 2]
01 ############################## end
02 ############################## start
0: ───────────────────────────────────YXXY────────────────────Z^0.0───
                                      │
1: ─────────────YXXY────────Z^0.522───#2^-0.656───YXXY────────Z^0.0───
                │                                 │
2: ───Z^0.762───#2^-0.249─────────────Z^-0.519────#2^-0.479───Z───────
02 ############################## end
03 ############################## start
1.00000014821
// Generated from Cirq

OPENQASM 2.0;
include "qelib1.inc";


// Qubits: [0, 1, 2]
qreg q[3];


// Gate: W(0.238)^-0.5
u3(pi*0.5,pi*0.7381494354,pi*1.2618505646) q[2];

// Gate: Y^0.5
u3(pi*0.5,pi*2.0,0) q[1];

// Gate: Y^0.5
u3(pi*0.5,pi*2.0,0) q[0];

// Gate: Exp11Gate(half_turns=0.24858996064554528)
u3(pi*0.5,pi*1.0,pi*0.75) q[1];
u3(pi*0.5,pi*1.0,pi*0.25) q[2];
rx(pi*0.5) q[1];
cx q[1],q[2];
rx(pi*0.3757050197) q[1];
ry(pi*0.5) q[2];
cx q[2],q[1];
rx(pi*-0.5) q[2];
rz(pi*0.5) q[2];
cx q[1],q[2];
u3(pi*0.5,pi*0.3742949803,0) q[1];
u3(pi*0.5,pi*0.8742949803,0) q[2];

// Gate: W(0.124)^0.5
u3(pi*0.5,pi*1.6242949803,pi*0.3757050197) q[1];

// Gate: W(0.862)^0.5
u3(pi*0.5,pi*0.3624444158,pi*1.6375555842) q[2];

// Gate: Exp11Gate(half_turns=0.24858996064554528)
u3(pi*0.5,pi*1.0,pi*0.75) q[1];
u3(pi*0.5,pi*1.0,pi*0.25) q[2];
rx(pi*0.5) q[1];
cx q[1],q[2];
rx(pi*0.3757050197) q[1];
ry(pi*0.5) q[2];
cx q[2],q[1];
rx(pi*-0.5) q[2];
rz(pi*0.5) q[2];
cx q[1],q[2];
u3(pi*0.5,pi*0.3742949803,0) q[1];
u3(pi*0.5,pi*0.8742949803,0) q[2];

// Gate: W(0.249)^0.5219350773492257
u3(pi*0.5219350773,pi*1.7485899606,pi*0.2514100394) q[1];

// Gate: W(0.987)^-0.5191093380568828
u3(pi*0.5191093381,pi*1.4867393961,pi*0.5132606039) q[2];

// Gate: Exp11Gate(half_turns=0.65648418204460857)
u3(pi*0.5,0,pi*1.4262081912) q[0];
u3(pi*0.5,0,pi*0.9262081912) q[1];
rx(pi*0.5) q[0];
cx q[0],q[1];
rx(pi*0.171757909) q[0];
ry(pi*0.5) q[1];
cx q[1],q[0];
rx(pi*-0.5) q[1];
rz(pi*0.5) q[1];
cx q[0],q[1];
u3(pi*0.5,pi*1.9020338998,pi*1.0) q[0];
u3(pi*0.5,pi*0.4020338998,pi*1.0) q[1];

// Gate: W(0.328)^-0.5
u3(pi*0.5,pi*0.828242091,pi*1.171757909) q[0];

// Gate: W(0.577)^-0.5
u3(pi*0.5,pi*1.0768320517,pi*0.9231679483) q[1];

// Gate: Exp11Gate(half_turns=0.65648418204460857)
u3(pi*0.5,0,pi*1.4262081912) q[0];
u3(pi*0.5,0,pi*0.9262081912) q[1];
rx(pi*0.5) q[0];
cx q[0],q[1];
rx(pi*0.171757909) q[0];
ry(pi*0.5) q[1];
cx q[1],q[0];
rx(pi*-0.5) q[1];
rz(pi*0.5) q[1];
cx q[0],q[1];
u3(pi*0.5,pi*1.9020338998,pi*1.0) q[0];
u3(pi*0.5,pi*0.4020338998,pi*1.0) q[1];

// Gate: W(0.656)^-0.1172916153092598
u3(pi*0.1172916153,pi*1.156484182,pi*0.843515818) q[0];

// Gate: W(0.905)^0.5
u3(pi*0.5,pi*0.4050741427,pi*1.5949258573) q[1];

// Gate: Exp11Gate(half_turns=0.47878385305425208)
u3(pi*0.5,pi*1.0,pi*1.25) q[1];
u3(pi*0.5,pi*1.0,pi*0.75) q[2];
rx(pi*0.5) q[1];
cx q[1],q[2];
rx(pi*0.2606080735) q[1];
ry(pi*0.5) q[2];
cx q[2],q[1];
rx(pi*-0.5) q[2];
rz(pi*0.5) q[2];
cx q[1],q[2];
u3(pi*0.5,pi*1.9893919265,0) q[1];
u3(pi*0.5,pi*0.4893919265,0) q[2];

// Gate: W(0.144)^-0.5
u3(pi*0.5,pi*0.6444660692,pi*1.3555339308) q[1];

// Gate: W(0.226)^-0.5
u3(pi*0.5,pi*0.7261313226,pi*1.2738686774) q[2];

// Gate: Exp11Gate(half_turns=0.47878385305425208)
u3(pi*0.5,pi*1.0,pi*1.25) q[1];
u3(pi*0.5,pi*1.0,pi*0.75) q[2];
rx(pi*0.5) q[1];
cx q[1],q[2];
rx(pi*0.2606080735) q[1];
ry(pi*0.5) q[2];
cx q[2],q[1];
rx(pi*-0.5) q[2];
rz(pi*0.5) q[2];
cx q[1],q[2];
u3(pi*0.5,pi*1.9893919265,0) q[1];
u3(pi*0.5,pi*0.4893919265,0) q[2];

// Gate: W(0.384)^-0.3104017289871259
u3(pi*0.310401729,pi*0.8838579957,pi*1.1161420043) q[1];

// Gate: W(0.466)^-0.843806028656013
u3(pi*0.8438060287,pi*0.9655232491,pi*1.0344767509) q[2];

// Gate: Exp11Gate(half_turns=0.47878385305425208)
u3(pi*0.5,pi*1.0,pi*1.25) q[1];
u3(pi*0.5,pi*1.0,pi*0.75) q[2];
rx(pi*0.5) q[1];
cx q[1],q[2];
rx(pi*0.2606080735) q[1];
ry(pi*0.5) q[2];
cx q[2],q[1];
rx(pi*-0.5) q[2];
rz(pi*0.5) q[2];
cx q[1],q[2];
u3(pi*0.5,pi*1.9893919265,0) q[1];
u3(pi*0.5,pi*0.4893919265,0) q[2];

// Gate: W(0.623)^-0.5
u3(pi*0.5,pi*1.1232499223,pi*0.8767500777) q[1];

// Gate: W(0.705)^-0.5
u3(pi*0.5,pi*1.2049151757,pi*0.7950848243) q[2];

// Gate: Exp11Gate(half_turns=0.47878385305425208)
u3(pi*0.5,pi*1.0,pi*1.25) q[1];
u3(pi*0.5,pi*1.0,pi*0.75) q[2];
rx(pi*0.5) q[1];
cx q[1],q[2];
rx(pi*0.2606080735) q[1];
ry(pi*0.5) q[2];
cx q[2],q[1];
rx(pi*-0.5) q[2];
rz(pi*0.5) q[2];
cx q[1],q[2];
u3(pi*0.5,pi*1.9893919265,0) q[1];
u3(pi*0.5,pi*0.4893919265,0) q[2];

// Gate: W(0.863)^-0.5
u3(pi*0.5,pi*1.3626418488,pi*0.6373581512) q[1];

// Gate: W(0.944)^-0.5191093380568828
u3(pi*0.5191093381,pi*1.4443071022,pi*0.5556928978) q[2];

// Gate: Exp11Gate(half_turns=0.65648418204460857)
u3(pi*0.5,0,pi*1.4262081912) q[0];
u3(pi*0.5,0,pi*0.9262081912) q[1];
rx(pi*0.5) q[0];
cx q[0],q[1];
rx(pi*0.171757909) q[0];
ry(pi*0.5) q[1];
cx q[1],q[0];
rx(pi*-0.5) q[1];
rz(pi*0.5) q[1];
cx q[0],q[1];
u3(pi*0.5,pi*1.9020338998,pi*1.0) q[0];
u3(pi*0.5,pi*0.4020338998,pi*1.0) q[1];

// Gate: W(0.985)^0.5
u3(pi*0.5,pi*0.4847262731,pi*1.5152737269) q[0];

// Gate: W(0.191)^0.5
u3(pi*0.5,pi*1.6908839398,pi*0.3091160602) q[1];

// Gate: Exp11Gate(half_turns=0.65648418204460857)
u3(pi*0.5,0,pi*1.4262081912) q[0];
u3(pi*0.5,0,pi*0.9262081912) q[1];
rx(pi*0.5) q[0];
cx q[0],q[1];
rx(pi*0.171757909) q[0];
ry(pi*0.5) q[1];
cx q[1],q[0];
rx(pi*-0.5) q[1];
rz(pi*0.5) q[1];
cx q[0],q[1];
u3(pi*0.5,pi*1.9020338998,pi*1.0) q[0];
u3(pi*0.5,pi*0.4020338998,pi*1.0) q[1];

// Gate: W(0.813)^0.5
u3(pi*0.5,pi*0.3129683641,pi*1.6870316359) q[0];

// Gate: W(0.519)^0.47806492265077405
u3(pi*0.4780649227,pi*0.0191260308,pi*1.9808739692) q[1];

rz(pi*0.6870316359) q[0];

// Gate: Exp11Gate(half_turns=0.24858996064554528)
u3(pi*0.5,pi*1.0,pi*0.75) q[1];
u3(pi*0.5,pi*1.0,pi*0.25) q[2];
rx(pi*0.5) q[1];
cx q[1],q[2];
rx(pi*0.3757050197) q[1];
ry(pi*0.5) q[2];
cx q[2],q[1];
rx(pi*-0.5) q[2];
rz(pi*0.5) q[2];
cx q[1],q[2];
u3(pi*0.5,pi*0.3742949803,0) q[1];
u3(pi*0.5,pi*0.8742949803,0) q[2];

// Gate: W(0.643)^-0.5
u3(pi*0.5,pi*1.1434210112,pi*0.8565789888) q[1];

// Gate: W(0.0686)^-0.5
u3(pi*0.5,pi*0.5686020825,pi*1.4313979175) q[2];

// Gate: Exp11Gate(half_turns=0.24858996064554528)
u3(pi*0.5,pi*1.0,pi*0.75) q[1];
u3(pi*0.5,pi*1.0,pi*0.25) q[2];
rx(pi*0.5) q[1];
cx q[1],q[2];
rx(pi*0.3757050197) q[1];
ry(pi*0.5) q[2];
cx q[2],q[1];
rx(pi*-0.5) q[2];
rz(pi*0.5) q[2];
cx q[1],q[2];
u3(pi*0.5,pi*0.3742949803,0) q[1];
u3(pi*0.5,pi*0.8742949803,0) q[2];

// Gate: W(0.268)^0.5
u3(pi*0.5,pi*1.7677159915,pi*0.2322840085) q[1];

// Gate: W(0.693)^-0.5
u3(pi*0.5,pi*1.1928970628,pi*0.8071029372) q[2];

rz(pi*-0.4547476274) q[2];
rz(pi*-0.7677159915) q[1];

03 ############################## end
04 ############################## start
// Generated from Cirq

OPENQASM 2.0;
include "qelib1.inc";


// Qubits: [0, 1, 2]
qreg q[3];


// Gate: W(0.238)^-0.5
u3(pi*0.5,pi*0.7381494354,pi*1.2618505646) q[2];

// Gate: Y^0.5
u3(pi*0.5,pi*2.0,0) q[1];

// Gate: Y^0.5
u3(pi*0.5,pi*2.0,0) q[0];

// Gate: Exp11Gate(half_turns=0.24858996064554528)
u3(pi*0.5,pi*1.0,pi*0.75) q[1];
u3(pi*0.5,pi*1.0,pi*0.25) q[2];
rx(pi*0.5) q[1];
cx q[1],q[2];
rx(pi*0.3757050197) q[1];
ry(pi*0.5) q[2];
cx q[2],q[1];
rx(pi*-0.5) q[2];
rz(pi*0.5) q[2];
cx q[1],q[2];
u3(pi*0.5,pi*0.3742949803,0) q[1];
u3(pi*0.5,pi*0.8742949803,0) q[2];

// Gate: W(0.124)^0.5
u3(pi*0.5,pi*1.6242949803,pi*0.3757050197) q[1];

// Gate: W(0.862)^0.5
u3(pi*0.5,pi*0.3624444158,pi*1.6375555842) q[2];

// Gate: Exp11Gate(half_turns=0.24858996064554528)
u3(pi*0.5,pi*1.0,pi*0.75) q[1];
u3(pi*0.5,pi*1.0,pi*0.25) q[2];
rx(pi*0.5) q[1];
cx q[1],q[2];
rx(pi*0.3757050197) q[1];
ry(pi*0.5) q[2];
cx q[2],q[1];
rx(pi*-0.5) q[2];
rz(pi*0.5) q[2];
cx q[1],q[2];
u3(pi*0.5,pi*0.3742949803,0) q[1];
u3(pi*0.5,pi*0.8742949803,0) q[2];

// Gate: W(0.249)^0.5219350773492257
u3(pi*0.5219350773,pi*1.7485899606,pi*0.2514100394) q[1];

// Gate: W(0.987)^-0.5191093380568828
u3(pi*0.5191093381,pi*1.4867393961,pi*0.5132606039) q[2];

// Gate: Exp11Gate(half_turns=0.65648418204460857)
u3(pi*0.5,0,pi*1.4262081912) q[0];
u3(pi*0.5,0,pi*0.9262081912) q[1];
rx(pi*0.5) q[0];
cx q[0],q[1];
rx(pi*0.171757909) q[0];
ry(pi*0.5) q[1];
cx q[1],q[0];
rx(pi*-0.5) q[1];
rz(pi*0.5) q[1];
cx q[0],q[1];
u3(pi*0.5,pi*1.9020338998,pi*1.0) q[0];
u3(pi*0.5,pi*0.4020338998,pi*1.0) q[1];

// Gate: W(0.328)^-0.5
u3(pi*0.5,pi*0.828242091,pi*1.171757909) q[0];

// Gate: W(0.577)^-0.5
u3(pi*0.5,pi*1.0768320517,pi*0.9231679483) q[1];

// Gate: Exp11Gate(half_turns=0.65648418204460857)
u3(pi*0.5,0,pi*1.4262081912) q[0];
u3(pi*0.5,0,pi*0.9262081912) q[1];
rx(pi*0.5) q[0];
cx q[0],q[1];
rx(pi*0.171757909) q[0];
ry(pi*0.5) q[1];
cx q[1],q[0];
rx(pi*-0.5) q[1];
rz(pi*0.5) q[1];
cx q[0],q[1];
u3(pi*0.5,pi*1.9020338998,pi*1.0) q[0];
u3(pi*0.5,pi*0.4020338998,pi*1.0) q[1];

// Gate: W(0.656)^-0.1172916153092598
u3(pi*0.1172916153,pi*1.156484182,pi*0.843515818) q[0];

// Gate: W(0.905)^0.5
u3(pi*0.5,pi*0.4050741427,pi*1.5949258573) q[1];

// Gate: Exp11Gate(half_turns=0.47878385305425208)
u3(pi*0.5,pi*1.0,pi*1.25) q[1];
u3(pi*0.5,pi*1.0,pi*0.75) q[2];
rx(pi*0.5) q[1];
cx q[1],q[2];
rx(pi*0.2606080735) q[1];
ry(pi*0.5) q[2];
cx q[2],q[1];
rx(pi*-0.5) q[2];
rz(pi*0.5) q[2];
cx q[1],q[2];
u3(pi*0.5,pi*1.9893919265,0) q[1];
u3(pi*0.5,pi*0.4893919265,0) q[2];

// Gate: W(0.144)^-0.5
u3(pi*0.5,pi*0.6444660692,pi*1.3555339308) q[1];

// Gate: W(0.226)^-0.5
u3(pi*0.5,pi*0.7261313226,pi*1.2738686774) q[2];

// Gate: Exp11Gate(half_turns=0.47878385305425208)
u3(pi*0.5,pi*1.0,pi*1.25) q[1];
u3(pi*0.5,pi*1.0,pi*0.75) q[2];
rx(pi*0.5) q[1];
cx q[1],q[2];
rx(pi*0.2606080735) q[1];
ry(pi*0.5) q[2];
cx q[2],q[1];
rx(pi*-0.5) q[2];
rz(pi*0.5) q[2];
cx q[1],q[2];
u3(pi*0.5,pi*1.9893919265,0) q[1];
u3(pi*0.5,pi*0.4893919265,0) q[2];

// Gate: W(0.384)^-0.3104017289871259
u3(pi*0.310401729,pi*0.8838579957,pi*1.1161420043) q[1];

// Gate: W(0.466)^-0.843806028656013
u3(pi*0.8438060287,pi*0.9655232491,pi*1.0344767509) q[2];

// Gate: Exp11Gate(half_turns=0.47878385305425208)
u3(pi*0.5,pi*1.0,pi*1.25) q[1];
u3(pi*0.5,pi*1.0,pi*0.75) q[2];
rx(pi*0.5) q[1];
cx q[1],q[2];
rx(pi*0.2606080735) q[1];
ry(pi*0.5) q[2];
cx q[2],q[1];
rx(pi*-0.5) q[2];
rz(pi*0.5) q[2];
cx q[1],q[2];
u3(pi*0.5,pi*1.9893919265,0) q[1];
u3(pi*0.5,pi*0.4893919265,0) q[2];

// Gate: W(0.623)^-0.5
u3(pi*0.5,pi*1.1232499223,pi*0.8767500777) q[1];

// Gate: W(0.705)^-0.5
u3(pi*0.5,pi*1.2049151757,pi*0.7950848243) q[2];

// Gate: Exp11Gate(half_turns=0.47878385305425208)
u3(pi*0.5,pi*1.0,pi*1.25) q[1];
u3(pi*0.5,pi*1.0,pi*0.75) q[2];
rx(pi*0.5) q[1];
cx q[1],q[2];
rx(pi*0.2606080735) q[1];
ry(pi*0.5) q[2];
cx q[2],q[1];
rx(pi*-0.5) q[2];
rz(pi*0.5) q[2];
cx q[1],q[2];
u3(pi*0.5,pi*1.9893919265,0) q[1];
u3(pi*0.5,pi*0.4893919265,0) q[2];

// Gate: W(0.863)^-0.5
u3(pi*0.5,pi*1.3626418488,pi*0.6373581512) q[1];

// Gate: W(0.944)^-0.5191093380568828
u3(pi*0.5191093381,pi*1.4443071022,pi*0.5556928978) q[2];

// Gate: Exp11Gate(half_turns=0.65648418204460857)
u3(pi*0.5,0,pi*1.4262081912) q[0];
u3(pi*0.5,0,pi*0.9262081912) q[1];
rx(pi*0.5) q[0];
cx q[0],q[1];
rx(pi*0.171757909) q[0];
ry(pi*0.5) q[1];
cx q[1],q[0];
rx(pi*-0.5) q[1];
rz(pi*0.5) q[1];
cx q[0],q[1];
u3(pi*0.5,pi*1.9020338998,pi*1.0) q[0];
u3(pi*0.5,pi*0.4020338998,pi*1.0) q[1];

// Gate: W(0.985)^0.5
u3(pi*0.5,pi*0.4847262731,pi*1.5152737269) q[0];

// Gate: W(0.191)^0.5
u3(pi*0.5,pi*1.6908839398,pi*0.3091160602) q[1];

// Gate: Exp11Gate(half_turns=0.65648418204460857)
u3(pi*0.5,0,pi*1.4262081912) q[0];
u3(pi*0.5,0,pi*0.9262081912) q[1];
rx(pi*0.5) q[0];
cx q[0],q[1];
rx(pi*0.171757909) q[0];
ry(pi*0.5) q[1];
cx q[1],q[0];
rx(pi*-0.5) q[1];
rz(pi*0.5) q[1];
cx q[0],q[1];
u3(pi*0.5,pi*1.9020338998,pi*1.0) q[0];
u3(pi*0.5,pi*0.4020338998,pi*1.0) q[1];

// Gate: W(0.813)^0.5
u3(pi*0.5,pi*0.3129683641,pi*1.6870316359) q[0];

// Gate: W(0.519)^0.47806492265077405
u3(pi*0.4780649227,pi*0.0191260308,pi*1.9808739692) q[1];

rz(pi*0.6870316359) q[0];

// Gate: Exp11Gate(half_turns=0.24858996064554528)
u3(pi*0.5,pi*1.0,pi*0.75) q[1];
u3(pi*0.5,pi*1.0,pi*0.25) q[2];
rx(pi*0.5) q[1];
cx q[1],q[2];
rx(pi*0.3757050197) q[1];
ry(pi*0.5) q[2];
cx q[2],q[1];
rx(pi*-0.5) q[2];
rz(pi*0.5) q[2];
cx q[1],q[2];
u3(pi*0.5,pi*0.3742949803,0) q[1];
u3(pi*0.5,pi*0.8742949803,0) q[2];

// Gate: W(0.643)^-0.5
u3(pi*0.5,pi*1.1434210112,pi*0.8565789888) q[1];

// Gate: W(0.0686)^-0.5
u3(pi*0.5,pi*0.5686020825,pi*1.4313979175) q[2];

// Gate: Exp11Gate(half_turns=0.24858996064554528)
u3(pi*0.5,pi*1.0,pi*0.75) q[1];
u3(pi*0.5,pi*1.0,pi*0.25) q[2];
rx(pi*0.5) q[1];
cx q[1],q[2];
rx(pi*0.3757050197) q[1];
ry(pi*0.5) q[2];
cx q[2],q[1];
rx(pi*-0.5) q[2];
rz(pi*0.5) q[2];
cx q[1],q[2];
u3(pi*0.5,pi*0.3742949803,0) q[1];
u3(pi*0.5,pi*0.8742949803,0) q[2];

// Gate: W(0.268)^0.5
u3(pi*0.5,pi*1.7677159915,pi*0.2322840085) q[1];

// Gate: W(0.693)^-0.5
u3(pi*0.5,pi*1.1928970628,pi*0.8071029372) q[2];

rz(pi*-0.4547476274) q[2];
rz(pi*-0.7677159915) q[1];

04 ############################## end

相關推薦

量子計算 OpenFermion-Cirq 開源python庫

#01 基礎資訊 NISQ，是嘈雜中型量子 (Noisy Intermediate-Scale Quantum) 的簡稱。 擁有50-100量子位元、以及高保真量子門 (Quantum Gate) 的計算機，便可稱為NISQ計算機。 谷歌開源了Cirq框架，這是

谷歌開源量子計算框架Cirq

AI 前線導讀：量子計算是不少科學家們辛勤耕耘的領域，在他們的努力下，量子計算領域有了令人驚喜的進步。在昨天舉行的第一屆量子軟體和量子機器學習國際研討會（QSML）上，Google AI Quantum 團隊釋出了 Cirq 的公開測試版，這是一款用於 NISQ 計算機的開源

量子計算 Quantum++ c++11開源庫

01 基本資料 QC simulators是一個現代C++ 11量子計算程式庫。程式碼提供了CMakeLists.txt和vs2017的工程檔案。依賴eigen開源庫，並且支援OpenMP。 下載最新的qpp程式碼。 cd d:\git git clone h

『python』科學計算專項_科學繪圖庫matplotlib學習之繪制動畫（待續）

同時 func 動作 .com block save init [] first 示例代碼 簡單調用繪圖 from matplotlib import pyplot as plt import matplotlib.animation as animation impor

OpenAI 開源機器人模擬 Python 庫，並行模擬處理速度提升400%

class 功能 lan pic 速度 tail pytho 實現 模擬 10000da.cnvboyule.cnjiaeidaypt.cn 在過去一年的研究中，OpenAI團隊開源一個使用 MuJoCoengine開發的用於機器人模擬的高性能Python庫。雷

『python』科學計算專項_科學繪圖庫matplotlib學習(下)

時序 nco 字符 color 由於 enc -- angle fig 基本的讀取csv文件並繪制餅圖 由於之前沒有過實際處理的經驗，所以這個程序還是值得一看，涉及了處理表格數據的基本方法： import matplotlib.pyplot as plt import p

利用python庫計算person相關係數

使用numpy庫，可以實現person相關係數的計算，例如對於矩陣a。 a Out[235]: array([[1, 1, 2, 2, 3], [2, 2, 3, 3, 5],

不再神祕的量子計算，用Python就能實現（視訊+程式碼）

python庫學習筆記——分組計算利器：pandas中的groupby技術

        最近處理資料需要分組計算，又用到了groupby函式，溫故而知新。   分組運算的第一階段，pandas 物件（無論是 Series、DataFrame 還是其他的）中的資料會根據你

20個必不可少的Python庫

數據庫 科學 str 腳本 tor 圖形用戶界面 ren 不為 api 轉載：http://www.python123.org/tutorials/58b41f2a28c8f30100bd41dc 讀者們好。今天我將介紹20個屬於我常用工具的Python庫，我相信你看完之

python庫資料匯總

and html scikit sci gen res taf afr pandas numpy 1.11.3 scipy (0.18.1) http://www.scipy-lectures.org/index.html 裏面包含numpy pandas0.19.2

Linux下經常使用的C/C++開源Socket庫

bsd 面向對象 sql數據庫 高速 com telnet ade ftp 版本 1. Linux Socket Programming In C++ : http://tldp.org/LDP/LG/issue74/tougher.html2. A

Python的6大Python庫

python庫Python程序員在編寫程序的時候對某一些庫都會有自己的偏愛，但是有幾個卻是所有python程序員都喜歡的庫，也學你也在用，但是你沒有在意，今天千鋒就簡單的說上幾個，也許你看完之後也會引起你的共鳴。　　在編程時，一些小的挫折可能和大難題一樣讓人痛苦。沒有人希望在自己費勁心思之後，只是做出彈消息的

好用的python庫（轉）

郵件 oot 1.3 內容 一個 你在 登錄 後臺管理 sin 這個專區就是用來分享你在使用 Python 過程中發現的或者自己寫的 Python 庫。有時候一個好用的庫將大大節省一個開發者的時間，也會讓開發者多這個庫的作者感激涕零的。例如我在開發 Pythonzhcn 時

各種 Python 庫/模塊/工具

開發框架 light 動態庫 ase received pyo 協議 environ self. 1 算法 1.1 字符串處理 re 正則表達式的標準庫。 StringIO / cStringIO 以讀寫文件的方式來操作字符串（有點類似於內存文件）。 cStringIO 是

批量更新python庫

ins distrib spa 批量更新 () blog all pre style 1 import pip 2 from subprocess import call 3 4 for dist in pip.get_installed_distributions()

python庫： scrapy （深坑未填）

cfg file pipeline lin () eating tps uci pip3 scrapy　　一個快速高級的屏幕爬取及網頁采集框架 http://scrapy.org/　　官網 https://docs.scrapy.org/en/latest/　　文檔

python庫： arrow （時間）

asi inf 環境 smi docs ace 派生 row -a arrow是個時間日期庫，簡潔易用。支持python3.6 https://arrow.readthedocs.io/en/latest/　　arrow官網api https://github.com/

利用pip批量更新python庫

imp python環境 如果 代碼 logs upgrade 批量 .get roc 如果python庫比較舊，需要更新到最新版本，可以利用pip工具。 DOS命令行下，輸入pip -V查看pip版本，可以先把pip更新到新版本。 查看系統裏過期的python庫 p

如何在Pypi上發表自己的Python庫

record 註意 esc code com script ... rip 大神 背景 最近興趣使然寫了幾個Python庫，也發布到了Pypi上，雖然沒什麽人下載，但自己在其他機器上用著也會很方便。這裏我向大家介紹一下如何在Pypi上發表自己的Python庫。 準備 註冊