def test_def_permutation_gate():
perm_gate = DefPermutationGate("CCNOT", [0, 1, 2, 3, 4, 5, 7, 6])
perm_gate_str = "DEFGATE CCNOT AS PERMUTATION:\n 0, 1, 2, 3, 4, 5, 7, 6".strip(
)
parse_equals(perm_gate_str, perm_gate)
def test_parsing_swap_phase():
parse_equals('SWAP-PHASE 0 "rf" 1 "rf"',
SwapPhase(Frame([Qubit(0)], "rf"), Frame([Qubit(1)], "rf")))
parse_equals(
'SWAP-PHASE 0 1 "ff" 1 0 "ff"',
SwapPhase(Frame([Qubit(0), Qubit(1)], "ff"),
Frame([Qubit(1), Qubit(0)], "ff")),
)
def test_parsing_delay():
parse_equals("DELAY 0 1.0", DelayQubits([Qubit(0)], 1.0))
parse_equals("DELAY 0 1", DelayQubits([Qubit(0)], 1))
parse_equals("DELAY 0 1 1e-6", DelayQubits([Qubit(0), Qubit(1)], 1e-6))
parse_equals('DELAY 0 "rf" 1.0', DelayFrames([Frame([Qubit(0)], "rf")], 1.0))
parse_equals(
'DELAY 0 "ro_tx" "ro_rx" 1.0',
DelayFrames([Frame([Qubit(0)], "ro_tx"), Frame([Qubit(0)], "ro_rx")], 1.0),
)
def test_parse_reset_qubit():
reset = """
RESET
""".strip()
parse_equals(reset, Reset())
reset_qubit = """
RESET 5
""".strip()
parse_equals(reset_qubit, ResetQubit(Qubit(5)))
def test_parsing_raw_capture():
parse_equals(
"DECLARE iqs REAL[200000]\n" 'RAW-CAPTURE 0 "ro_rx" 0.001 iqs',
Declare("iqs", "REAL", 200000),
RawCapture(Frame([Qubit(0)], "ro_rx"), 0.001, MemoryReference("iqs")),
)
parse_equals(
"DECLARE iqs REAL[200000]\n" 'NONBLOCKING RAW-CAPTURE 0 "ro_rx" 0.001 iqs',
Declare("iqs", "REAL", 200000),
RawCapture(Frame([Qubit(0)], "ro_rx"), 0.001, MemoryReference("iqs"), nonblocking=True),
)
def test_parsing_capture():
wf = FlatWaveform(duration=1.0, iq=1.0)
parse_equals(
"DECLARE iq REAL[2]\n" 'CAPTURE 0 "ro_rx" flat(duration: 1.0, iq: 1.0) iq',
Declare("iq", "REAL", 2),
Capture(Frame([Qubit(0)], "ro_rx"), wf, MemoryReference("iq")),
)
parse_equals(
"DECLARE iq REAL[2]\n" 'NONBLOCKING CAPTURE 0 "ro_rx" flat(duration: 1.0, iq: 1.0) iq',
Declare("iq", "REAL", 2),
Capture(Frame([Qubit(0)], "ro_rx"), wf, MemoryReference("iq"), nonblocking=True),
)
def test_def_gate_with_variables():
# Note that technically the RX gate includes -i instead of just i but this messes a bit with
# the test since it's not smart enough to figure out that -1*i == -i
theta = Parameter("theta")
rx = np.array([
[quil_cos(theta / 2), 1j * quil_sin(theta / 2)],
[1j * quil_sin(theta / 2),
quil_cos(theta / 2)],
])
defgate = "DEFGATE RX(%theta):\n" " COS(%theta/2), i*SIN(%theta/2)\n" " i*SIN(%theta/2), COS(%theta/2)\n\n"
parse_equals(defgate, DefGate("RX", rx, [theta]))
def test_jumps():
parse_equals("LABEL @test_1", JumpTarget(Label("test_1")))
parse_equals("JUMP @test_1", Jump(Label("test_1")))
parse_equals("JUMP-WHEN @test_1 ro[0]",
JumpWhen(Label("test_1"), MemoryReference("ro", 0)))
parse_equals("JUMP-UNLESS @test_1 ro[1]",
JumpUnless(Label("test_1"), MemoryReference("ro", 1)))
def test_def_gate():
sqrt_x = DefGate("SQRT-X", np.array([[0.5 + 0.5j, 0.5 - 0.5j], [0.5 - 0.5j, 0.5 + 0.5j]]))
hadamard = DefGate("HADAMARD", np.array([[1 / np.sqrt(2), 1 / np.sqrt(2)], [1 / np.sqrt(2), -1 / np.sqrt(2)]]))
defgates = """
DEFGATE SQRT-X:
0.5+0.5i, 0.5-0.5i
0.5-0.5i, 0.5+0.5i
DEFGATE HADAMARD:
1/SQRT(2), 1/SQRT(2)
1/SQRT(2), -1/SQRT(2)
""".strip()
parse_equals(defgates, sqrt_x, hadamard)
def test_parsing_defcal_measure():
parse_equals("DEFCAL MEASURE 0:\n" " NOP\n", DefMeasureCalibration(Qubit(0), None, [NOP]))
wf = FlatWaveform(duration=1.0, iq=1.0)
# TODO: note that in a calibration body, reference to the formal argument addr parses
# as a memoryreference.
parse_equals(
"DEFCAL MEASURE q addr:\n"
' PULSE q "ro_tx" flat(duration: 1.0, iq: 1.0+0.0*i)\n'
' CAPTURE q "ro_rx" flat(duration: 1.0, iq: 1.0+0*i) addr[0]\n',
DefMeasureCalibration(
FormalArgument("q"),
FormalArgument("addr"),
[
Pulse(Frame([FormalArgument("q")], "ro_tx"), wf),
Capture(Frame([FormalArgument("q")], "ro_rx"), wf, MemoryReference("addr")),
],
),
)
def test_parsing_frame_mutations():
ops = [
("SET-PHASE", SetPhase),
("SHIFT-PHASE", ShiftPhase),
("SET-SCALE", SetScale),
("SET-FREQUENCY", SetFrequency),
("SHIFT-FREQUENCY", ShiftFrequency),
]
frames = [
('0 "rf"', Frame([Qubit(0)], "rf")),
('0 1 "ff"', Frame([Qubit(0), Qubit(1)], "ff")),
('1 0 "ff"', Frame([Qubit(1), Qubit(0)], "ff")),
]
values = [("1", 1), ("1.0", 1.0), ("pi/2", np.pi / 2)] # TODO: should we require a float here?
for op_str, op in ops:
for frame_str, frame in frames:
for val_str, val in values:
parse_equals(f"{op_str} {frame_str} {val_str}", op(frame, val))
def test_parsing_defframe_round_trip_with_json():
fdef = DefFrame(
frame=Frame(qubits=[FormalArgument("My-Cool-Qubit")], name="bananas"),
direction="go west",
initial_frequency=123.4,
center_frequency=124.5,
hardware_object='{"key1": 3.1, "key2": "value2"}',
sample_rate=5,
)
fdef_out = fdef.out()
assert (fdef.out() == r"""DEFFRAME My-Cool-Qubit "bananas":
DIRECTION: "go west"
INITIAL-FREQUENCY: 123.4
CENTER-FREQUENCY: 124.5
HARDWARE-OBJECT: "{\"key1\": 3.1, \"key2\": \"value2\"}"
SAMPLE-RATE: 5
""")
parse_equals(fdef_out, fdef)
def test_parsing_defwaveform():
parse_equals("DEFWAVEFORM foo:\n" " 1.0, 1.0, 1.0\n", DefWaveform("foo", [], [1.0, 1.0, 1.0]))
parse_equals(
"DEFWAVEFORM foo:\n" " 1.0+2.0*i, 1.0-2.0*i, 3.0\n",
DefWaveform("foo", [], [1 + 2j, 1 - 2j, 3 + 0j]),
)
parse_equals(
"DEFWAVEFORM foo(%theta):\n" " 1.0+2.0*i, 1.0-2.0*i, 3.0*%theta\n",
DefWaveform("foo", [Parameter("theta")], [1 + 2j, 1 - 2j, Mul(3.0, Parameter("theta"))]),
)
parse_equals(
"DEFWAVEFORM q0_ro_rx/filter:\n 1.0, 1.0, 1.0",
DefWaveform("q0_ro_rx/filter", [], [1.0, 1.0, 1.0]),
)
def test_def_circuit():
defcircuit = """
DEFCIRCUIT bell a b:
H a
CNOT a b
""".strip()
gate = "bell 0 1"
defcircuit_no_qubits = """
DEFCIRCUIT bell:
H 0
CNOT 0 1
""".strip()
gate_no_qubits = "bell"
defcircuit_param = """
DEFCIRCUIT parameterized(%theta, %phi) a:
RX(%theta) a
RZ(%phi) a
""".strip()
gate_param = "parameterized(0.0, 1.0) 0"
parse_equals(defcircuit + "\n" + gate, RawInstr(defcircuit),
Gate("bell", [], [Qubit(0), Qubit(1)]))
parse_equals(
defcircuit_no_qubits + "\n" + gate_no_qubits,
RawInstr(defcircuit_no_qubits),
RawInstr(gate_no_qubits),
)
parse_equals(
defcircuit_param + "\n" + gate_param,
RawInstr(defcircuit_param),
Gate("parameterized", [0.0, 1.0], [Qubit(0)]),
)
def test_parse_pulse():
wf = FlatWaveform(duration=1.0, iq=1.0)
parse_equals('PULSE 0 "rf" flat(duration: 1.0, iq: 1.0)', Pulse(Frame([Qubit(0)], "rf"), wf))
parse_equals('PULSE 0 1 "ff" flat(duration: 1.0, iq: 1.0)', Pulse(Frame([Qubit(0), Qubit(1)], "ff"), wf))
parse_equals(
'NONBLOCKING PULSE 0 "rf" flat(duration: 1.0, iq: 1.0)',
Pulse(Frame([Qubit(0)], "rf"), wf, nonblocking=True),
)
def test_parsing_defframe():
parse_equals('DEFFRAME 0 "rf"', DefFrame(Frame([Qubit(0)], "rf")))
parse_equals('DEFFRAME 1 0 "ff"', DefFrame(Frame([Qubit(1), Qubit(0)], "ff")))
parse_equals(
'DEFFRAME 0 "rf":\n' " SAMPLE-RATE: 2.0\n",
DefFrame(Frame([Qubit(0)], "rf"), sample_rate=2.0),
)
parse_equals(
'DEFFRAME 0 "rf":\n'
" SAMPLE-RATE: 2.0\n"
" INITIAL-FREQUENCY: 10\n", # TODO: should this parse as a float?
DefFrame(Frame([Qubit(0)], "rf"), sample_rate=2.0, initial_frequency=10),
)
with pytest.raises(UnexpectedToken) as excp:
parse('DEFFRAME 0 "rf":\n' " UNSUPPORTED: 2.0\n")
assert excp.value.token == Token("IDENTIFIER", "UNSUPPORTED")
def test_parsing_defcal():
parse_equals("DEFCAL X 0:\n" " NOP\n", DefCalibration("X", [], [Qubit(0)], [NOP]))
parse_equals(
"DEFCAL X q:\n" " NOP\n" " NOP\n",
DefCalibration("X", [], [FormalArgument("q")], [NOP, NOP]),
)
parse_equals(
"DEFCAL RZ(%theta) 0:\n" ' SHIFT-PHASE 0 "rf" %theta/(-2*pi)\n',
DefCalibration(
"RZ",
[Parameter("theta")],
[Qubit(0)],
[ShiftPhase(Frame([Qubit(0)], "rf"), Div(Parameter("theta"), -2 * np.pi))],
),
)
def wf_agrees(wf_str: str, wf: TemplateWaveform):
frame = Frame([Qubit(0)], "rf")
pulse_str = 'PULSE 0 "rf" ' + wf_str
parse_equals(pulse_str, Pulse(frame, wf))
def test_messy_modifiers():
s = "FORKED DAGGER CONTROLLED FORKED RX(0.1,0.2,0.3,0.4) 0 1 2 3"
parse_equals(
s,
RX(0.1, 3).forked(2,
[0.2]).controlled(1).dagger().forked(0, [0.3, 0.4]))
def test_parse_forked():
s = "FORKED RX(0, pi/2) 0 1"
parse_equals(s, RX(0, 1).forked(0, [np.pi / 2]))
def test_parse_controlled():
s = "CONTROLLED X 0 1"
parse_equals(s, X(1).controlled(0))
def test_parse_dagger():
s = "DAGGER X 0"
parse_equals(s, X(0).dagger())
s = "DAGGER DAGGER X 0"
parse_equals(s, X(0).dagger().dagger())
def test_measure():
parse_equals("MEASURE 0", MEASURE(0, None))
parse_equals("MEASURE 0 ro[1]", MEASURE(0, MemoryReference("ro", 1)))
def test_empty_program():
parse_equals("")
def test_pragma():
parse_equals('PRAGMA gate_time H "10 ns"',
Pragma("gate_time", ["H"], "10 ns"))
parse_equals("PRAGMA qubit 0", Pragma("qubit", [0]))
parse_equals("PRAGMA NO-NOISE", Pragma("NO-NOISE"))
def test_classical():
parse_equals("MOVE ro[0] 1", MOVE(MemoryReference("ro", 0), 1))
parse_equals("MOVE ro[0] 0", MOVE(MemoryReference("ro", 0), 0))
parse_equals("NOT ro[0]", NOT(MemoryReference("ro", 0)))
parse_equals("AND ro[0] 1", AND(MemoryReference("ro", 0), 1))
parse_equals("IOR ro[0] 1", IOR(MemoryReference("ro", 0), 1))
parse_equals("MOVE ro[0] 1", MOVE(MemoryReference("ro", 0), 1))
parse_equals("XOR ro[0] 1", XOR(MemoryReference("ro", 0), 1))
parse_equals("ADD mem[0] 1.2", ADD(MemoryReference("mem", 0), 1.2))
parse_equals("SUB mem[0] 1.2", SUB(MemoryReference("mem", 0), 1.2))
parse_equals("MUL mem[0] 1.2", MUL(MemoryReference("mem", 0), 1.2))
parse_equals("DIV mem[0] 1.2", DIV(MemoryReference("mem", 0), 1.2))
parse_equals("ADD mem[0] -1.2", ADD(MemoryReference("mem", 0), -1.2))
parse_equals("SUB mem[0] -1.2", SUB(MemoryReference("mem", 0), -1.2))
parse_equals("MUL mem[0] -1.2", MUL(MemoryReference("mem", 0), -1.2))
parse_equals("DIV mem[0] -1.2", DIV(MemoryReference("mem", 0), -1.2))
parse_equals(
"EQ comp[1] ro[3] ro[2]",
EQ(MemoryReference("comp", 1), MemoryReference("ro", 3),
MemoryReference("ro", 2)),
)
parse_equals(
"LT comp[1] ro[3] ro[2]",
LT(MemoryReference("comp", 1), MemoryReference("ro", 3),
MemoryReference("ro", 2)),
)
parse_equals(
"LE comp[1] ro[3] ro[2]",
LE(MemoryReference("comp", 1), MemoryReference("ro", 3),
MemoryReference("ro", 2)),
)
parse_equals(
"GT comp[1] ro[3] ro[2]",
GT(MemoryReference("comp", 1), MemoryReference("ro", 3),
MemoryReference("ro", 2)),
)
parse_equals(
"GE comp[1] ro[3] ro[2]",
GE(MemoryReference("comp", 1), MemoryReference("ro", 3),
MemoryReference("ro", 2)),
)
parse_equals("EQ comp[1] ro[3] 0",
EQ(MemoryReference("comp", 1), MemoryReference("ro", 3), 0))
parse_equals("LT comp[1] ro[3] 1",
LT(MemoryReference("comp", 1), MemoryReference("ro", 3), 1))
parse_equals("LE comp[1] ro[3] 2",
LE(MemoryReference("comp", 1), MemoryReference("ro", 3), 2))
parse_equals("GT comp[1] ro[3] 3",
GT(MemoryReference("comp", 1), MemoryReference("ro", 3), 3))
parse_equals("GE comp[1] ro[3] 4",
GE(MemoryReference("comp", 1), MemoryReference("ro", 3), 4))
parse_equals("EQ comp[1] ro[3] 0.0",
EQ(MemoryReference("comp", 1), MemoryReference("ro", 3), 0.0))
parse_equals("LT comp[1] ro[3] 1.1",
LT(MemoryReference("comp", 1), MemoryReference("ro", 3), 1.1))
parse_equals("LE comp[1] ro[3] 2.2",
LE(MemoryReference("comp", 1), MemoryReference("ro", 3), 2.2))
parse_equals("GT comp[1] ro[3] 3.3",
GT(MemoryReference("comp", 1), MemoryReference("ro", 3), 3.3))
parse_equals("GE comp[1] ro[3] 4.4",
GE(MemoryReference("comp", 1), MemoryReference("ro", 3), 4.4))
def test_memory_commands():
parse_equals(
"DECLARE mem OCTET[32] SHARING mem2 OFFSET 16 REAL OFFSET 32 REAL",
Declare("mem",
"OCTET",
32,
shared_region="mem2",
offsets=[(16, "REAL"), (32, "REAL")]),
)
parse_equals(
"STORE mem ro[2] ro[0]",
STORE("mem", MemoryReference("ro", 2), MemoryReference("ro", 0)))
parse_equals("STORE mem ro[2] 7", STORE("mem", MemoryReference("ro", 2),
7))
parse_equals(
"LOAD ro[8] mem mem[4]",
LOAD(MemoryReference("ro", 8), "mem", MemoryReference("mem", 4)))
parse_equals("CONVERT ro[1] ro[2]",
CONVERT(MemoryReference("ro", 1), MemoryReference("ro", 2)))
parse_equals("EXCHANGE ro[0] ro[1]",
EXCHANGE(MemoryReference("ro", 0), MemoryReference("ro", 1)))
parse_equals("MOVE mem[2] 4", MOVE(MemoryReference("mem", 2), 4))
parse_equals("MOVE mem[2] -4", MOVE(MemoryReference("mem", 2), -4))
parse_equals("MOVE mem[2] -4.1", MOVE(MemoryReference("mem", 2), -4.1))
def test_parsing_fence():
parse_equals("FENCE 0", Fence([Qubit(0)]))
parse_equals("FENCE 0 1", Fence([Qubit(0), Qubit(1)]))
parse_equals("FENCE q s", Fence([FormalArgument("q"),
FormalArgument("s")]))
parse_equals("FENCE", FenceAll())
def test_others():
parse_equals("RESET", RESET())
parse_equals("WAIT", WAIT)
parse_equals("NOP", NOP)
def _expr(expression, expected):
parse_equals("RX(" + expression + ") 0", RX(expected, 0))