-
Notifications
You must be signed in to change notification settings - Fork 7
/
Copy pathmeyer_quantum_classical_game.py
259 lines (232 loc) · 8.86 KB
/
meyer_quantum_classical_game.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
import pygame
import numpy as np
import pyquil.api as api
from pyquil.gates import *
from pyquil.quil import Program
# define colors
blue = (0, 0, 255)
rect_Q_color = (255, 127.5, 127.5)
rect_P_color = (127.5, 127.5, 255)
# boolean controling when game stops
game_over = False
# frames per second
fps = 60
# screen dimensions
size_x = 1000
size_y = 625
size = (size_x, size_y)
# rectangle coordinates
rect_h = 50
rect_Q_x = 0
rect_Q_y = 0
rect_P_x = 0
rect_P_y = size_y - rect_h
# fire coordinates
fire_x = -20
fire_y = -20
fire_width = 5
fire_length = 30
fire_y_change = 0
# track move_numbers, neither player has made a move at beginning
move_Q = 0
move_P = 0
# track strategy parameters
Q_prob1 = None
P_a1 = None
P_b1 = None
Q_prob2 = None
# initialize destruction log's parameters
dest_y = size_y//2
dest_width = 30
# boolean controling pyquil program run
run_quantum_program = False
# generate unitary matrix
def U_(a, b):
return np.array([[a, b], [b, -a]])
# define flip (X) and not flip (I) operators
X_ = np.array([[0, 1], [1, 0]])
I_ = np.array([[1, 0], [0, 1]])
# initialize pygame
pygame.init()
# setup screen
screen = pygame.display.set_mode(size)
pygame.display.set_caption("Meyer penny quantum/classical game demo")
# setup game clock
clock = pygame.time.Clock()
# define font
font = pygame.font.SysFont('Calibri', 25, True)
while not game_over:
# track moves
first_move = ((move_Q == 0) and (move_P == 0))
second_move = ((move_Q == 1) and (move_P == 0))
third_move = ((move_Q == 1) and (move_P == 1))
for event in pygame.event.get():
if event.type == pygame.QUIT:
game_over = True
if event.type == pygame.KEYDOWN:
if event.key == pygame.K_SPACE:
if first_move:
fire_x = (2 * rect_Q_x + rect_h)/2
fire_y = rect_Q_y + rect_h
elif second_move:
fire_x = (2 * rect_P_x + rect_h)/2
fire_y = rect_P_y - fire_length
elif third_move:
fire_x = (2 * rect_Q_x + rect_h)/2
fire_y = rect_Q_y + rect_h
elif event.type == pygame.KEYUP:
if event.key == pygame.K_SPACE:
if first_move:
Q_prob1 = rect_Q_x / (size_x - rect_h)
move_Q += 1
fire_y_change = 15
elif second_move:
P_b1 = np.sqrt(rect_P_x / (size_x - rect_h))
P_a1 = np.sqrt(1 - np.square(P_b1))
move_P += 1
fire_y_change = -15
elif third_move:
Q_prob2 = rect_Q_x / (size_x - rect_h)
move_Q += 1
fire_y_change = 15
# set game color
screen.fill(blue)
# draw rectangles
pygame.draw.rect(screen, rect_Q_color, (rect_Q_x, rect_Q_y, rect_h, rect_h))
pygame.draw.rect(screen, rect_P_color, (rect_P_x, rect_P_y, rect_h, rect_h))
# draw fire
fire_y += fire_y_change
pygame.draw.rect(screen, (255, 125, 125), [fire_x, fire_y, fire_width, fire_length])
## motions
pressed = pygame.key.get_pressed()
# rectangle motion
if pressed[pygame.K_LEFT]:
if first_move:
rect_Q_x -= 5
elif second_move:
rect_P_x -= 5
elif third_move:
rect_Q_x -= 5
if pressed[pygame.K_RIGHT]:
if first_move:
rect_Q_x += 5
elif second_move:
rect_P_x += 5
elif third_move:
rect_Q_x += 5
# keep rectangle within screen
if rect_Q_x > size_x - rect_h:
rect_Q_x = size_x - rect_h
elif rect_Q_x < 0:
rect_Q_x = 0
if rect_P_x > size_x - rect_h:
rect_P_x = size_x - rect_h
elif rect_P_x < 0:
rect_P_x = 0
if first_move:
# display Q's potential 1st choice
text_rect_Q_x = font.render("Prob(flip): " + str(rect_Q_x / (size_x - rect_h)), True, (0, 0, 0))
screen.blit(text_rect_Q_x, [150, 50])
text_size_ = font.render("Prob(not flip): " + str(1 - rect_Q_x / (size_x - rect_h)), True, (0, 0, 0))
screen.blit(text_size_, [150, 20])
elif second_move:
# display Q's 1st choice
text_rect_Q_x = font.render("Prob(flip): " + str(Q_prob1), True, (0, 0, 0))
screen.blit(text_rect_Q_x, [150, 50])
text_size_ = font.render("Prob(not flip): " + str(1 - Q_prob1), True, (0, 0, 0))
screen.blit(text_size_, [150, 20])
# display P's 1st choice
text_rect_P_x = font.render("|b|^2: " + str(rect_P_x / (size_x - rect_h)), True, (0, 0, 0))
screen.blit(text_rect_P_x, [150, size_y - 50])
text_size_ = font.render("|a|^2: " + str(1 - rect_P_x / (size_x - rect_h)), True, (0, 0, 0))
screen.blit(text_size_, [150, size_y - 100])
elif third_move:
# display Q's 1st choice
text_rect_Q_x = font.render("Prob(flip): " + str(Q_prob1), True, (0, 0, 0))
screen.blit(text_rect_Q_x, [150, 50])
text_size_ = font.render("Prob(not flip): " + str(1 - Q_prob1), True, (0, 0, 0))
screen.blit(text_size_, [150, 20])
# display P's 1st choice
text_rect_P_x = font.render("|b|^2: " + str(np.square(P_b1)), True, (0, 0, 0))
screen.blit(text_rect_P_x, [150, size_y - 50])
text_size_ = font.render("|a|^2: " + str(np.square(P_a1)), True, (0, 0, 0))
screen.blit(text_size_, [150, size_y - 100])
# display Q's potential 2nd choice
text_rect_Q_x = font.render("Prob(flip): " + str(rect_Q_x / (size_x - rect_h)), True, (0, 0, 0))
screen.blit(text_rect_Q_x, [size_x - 400, 50])
text_size_ = font.render("Prob(not flip): " + str(1 - rect_Q_x / (size_x - rect_h)), True, (0, 0, 0))
screen.blit(text_size_, [size_x - 400, 20])
else:
### Display all choices made
# display Q's 1st choice
text_rect_Q_x = font.render("Prob(flip): " + str(Q_prob1), True, (0, 0, 0))
screen.blit(text_rect_Q_x, [150, 50])
text_size_ = font.render("Prob(not flip): " + str(1 - Q_prob1), True, (0, 0, 0))
screen.blit(text_size_, [150, 20])
# display P's 1st choice
text_rect_P_x = font.render("|b|^2: " + str(np.square(P_b1)), True, (0, 0, 0))
screen.blit(text_rect_P_x, [150, size_y - 50])
text_size_ = font.render("|a|^2: " + str(np.square(P_a1)), True, (0, 0, 0))
screen.blit(text_size_, [150, size_y - 100])
# display Q's 2nd choice
text_rect_Q_x = font.render("Prob(flip): " + str(Q_prob2), True, (0, 0, 0))
screen.blit(text_rect_Q_x, [size_x - 400, 50])
text_size_ = font.render("Prob(not flip): " + str(1 - Q_prob2), True, (0, 0, 0))
screen.blit(text_size_, [size_x - 400, 20])
# program is now ready to be run
run_quantum_program = True
if run_quantum_program:
# draw DESTRUCTION! log
pygame.draw.rect(screen, (0, 0, 0), (0, dest_y, size_x, dest_width))
### Carry out program
# initialize program
qvm = api.QVMConnection()
p = Program()
##### Q's 1st stochastic move
# set the parameters for the random move
prob = Q_prob1
# create a dummy gate
a = np.random.rand()
b = np.sqrt(1 - np.square(a))
Uf_ = U_(a, b)
# perform the random move
p.defgate("Uf_Q1", Uf_)
p.inst(("Uf_Q1", 0))
p.define_noisy_gate("Uf_Q1", [0], [np.sqrt(prob) * X_, np.sqrt(1 - prob) * I_])
##### Picard's quantum move
p.defgate("P1", U_(P_a1, P_b1))
p.inst(("P1", 0))
##### Q's 2nd stochastic move
# set the parameters for the random move
prob = Q_prob2
# create a dummy gate
a = np.random.rand()
b = np.sqrt(1 - np.square(a))
Uf_ = U_(a, b)
# perform the random move
p.defgate("Uf_Q2", Uf_)
p.inst(("Uf_Q2", 0))
p.define_noisy_gate("Uf_Q2", [0], [np.sqrt(prob) * X_, np.sqrt(1 - prob) * I_])
##### Measurement
p.inst(MEASURE(0, [0]))
result = qvm.run(p, [0], trials=100)
result_1s = (len([i for i in result if i == [1]]))
result_0s = (len([i for i in result if i == [0]]))
Picard_score = result_1s
Q_score = result_0s
# move DESTRUCTION! log according to scores
delta_y = size_y/1000
dest_y -= int(delta_y * Picard_score)
dest_y += int(delta_y * Q_score)
if dest_y > size_y - dest_width:
dest_y = size_y - dest_width
elif dest_y < 0:
dest_y = 0
text_1s = font.render("Picard's score: " + str(result_1s), True, (0, 0, 0))
screen.blit(text_1s, [size_x/2, size_y/2])
text_0s = font.render("Q's score: " + str(result_0s), True, (0, 0, 0))
screen.blit(text_0s, [size_x/2, size_y/2 + 50])
# avoid repeating program
run_quantum_program = False
pygame.display.flip()
clock.tick(fps)