編集履歴

質問編集履歴

編集

2018/10/20 15:43

投稿

スコア124

test CHANGED Viewed

	@@ -1 +1 @@
1	- この~~論文の結果を再現できる方い~~ら~~っしゃ~~い~~ますか~~
1	+ 数値計算の問題点がわからない

test CHANGED Viewed

@@ -9,14 +9,6 @@
 [該当論文](https://arxiv.org/abs/1010.4705)
 ちなみに理解しているつもりなんですが、全くできません。コードもエラーなしで動くのですが、、、（ということは何かがわかっていないということなのですが）
-要点は初期状態が等しい重ね合わせ状態であることと、メイン計算は欲しいい点の確率振幅にマイナスを毎回かけながらすべての格子点に拡散行列とフリップフロップシフトというシフトの仕方が必要なようです。
-### イメージ図
-![正方格子量子探索](7ae71615029fd007225a2f924e5c1e7c.png)
-上の図は3×③の正方格子点を表しており、赤い点がマーキングされた点です（マーキングされた点は任意です）。欲しい結果にある結果は190番目の格子点にマーキングをした20×20の正方格子の結果です。
 ```

内容変更

2018/10/20 15:43

投稿

Fallout_18

スコア124

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -36,16 +36,6 @@
             y = i[1]
-            #boundary condition
-            x0 = (x+1) % N
-            x1 = (x-1 + N) % N
-            y0 = (y+1) % N
-            y1 = (y-1) % N
             next_psi[x,y] = np.copy(np.array(np.dot(right,psi[x1,y])+np.dot(down,psi[x,y0])+np.dot(left,psi[x0,y])+ np.dot(up,psi[x,y1])))
         psi = np.copy(next_psi)

内容変更

2018/10/05 07:14

投稿

Fallout_18

スコア124

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -7,10 +7,6 @@
 ![20*20の正方格子](7c8308f5057a6e57ebfeea4b0bdf97ee.jpeg)
 [該当論文](https://arxiv.org/abs/1010.4705)
-[参考論文２](https://arxiv.org/abs/1203.3950)
-[参考論文](https://arxiv.org/abs/1005.3676)
 ちなみに理解しているつもりなんですが、全くできません。コードもエラーなしで動くのですが、、、（ということは何かがわかっていないということなのですが）
@@ -22,49 +18,9 @@
 上の図は3×③の正方格子点を表しており、赤い点がマーキングされた点です（マーキングされた点は任意です）。欲しい結果にある結果は190番目の格子点にマーキングをした20×20の正方格子の結果です。
-### 実行してみた
+```
-```python
-import numpy as np
+#main calculation
-import matplotlib.pyplot as plt
-import itertools
-N = 20 #0<=x,y<=N
-d = 4 #フリップフロップシフトの移動方向(4方向（4状態ともいう）)
-itr =100
-#拡散行列
-C = np.array([[-1,1,1,1],[1,-1,1,1],[1,1,-1,1],[1,1,1,-1]])/2
-x_list = [i for i in range(N)]
-y_list = [i for i in range(N)]
-time = [i for i in range(itr)]
-#フリップフロップシフトに場合分け（ここが怪しい？）
-#right =[1,0,0,0],down=[0,1,0,0],left=[0,0,1,0],up=[0,0,0,1]
-right= np.zeros((4,4));right[2,:] = C[2,:] #stage[x,y][0]
-down = np.zeros((4,4));down[3,:] = C[3,:] #stage[x,y][1]
-left = np.zeros((4,4));left[0,:] = C[0,:] #stage[x,y][2]
-up = np.zeros((4,4));up[1,:] = C[1,:] #stage[x,y][3]にそれぞれ作用させる
-#初期状態
-psi = np.ones((N,N,4),dtype="float")/np.sqrt(N*N*4)
-prob =[]
 for t in range(itr):
@@ -74,11 +30,7 @@
     else:
-        next_psi = np.zeros((N,N,4),dtype="float")　#t+1秒後のpsiを収納する
-        psi[2,2] *= -1 #マーキング
-        for i in itertools.product(x_list,y_list):　#すべての格子点で計算
+        for i in itertools.product(x_list,y_list):
             x = i[0]
@@ -100,16 +52,14 @@
     #print(t,psi)
-    prob.append(np.inner(psi[2,2],psi[2,2]))
+    mark_prob.append(np.inner(psi[8,8],psi[8,8]))
-plt.plot(time,prob)
+    prob.append(np.inner(psi[16,16],psi[16,16]))
-plt.show()
 ```
+### 実行してみた
 ![イメージ説明](3b192866ddf979e8d53194625f9b8147.png)
 こんなへんちくりんな結果になります。
-どこが悪いのかわかりません！！　教えてくれる人も周りにいません。。。ご教授のほどよろしくお願い致します。

追記

2018/10/02 15:16

投稿

Fallout_18

スコア124

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -32,11 +32,11 @@
 import itertools
-N = 20
+N = 20 #0<=x,y<=N
-d = 4
+d = 4 #フリップフロップシフトの移動方向(4方向（4状態ともいう）)
-itr =100
+itr =100
 #拡散行列

内容変更

2018/10/02 04:49

投稿

Fallout_18

スコア124

test CHANGED Viewed

	@@ -1 +1 @@
1	- ~~2次元正~~方~~格子上のグローバーアルゴリズムがう~~まくいかない
1	+ この論文の結果を再現できる方いらっしゃいますか

test CHANGED Viewed

@@ -1,70 +1,28 @@
-###　目的
- 20×20の正方格子上のマークされた一つの座標をグローバーアルゴリズム（量子アルゴリズムの一つ）と周期境界条件をつかって時間ごとの存在確率を見る
 ### 問題点
-やり方をわかっているつもりだが、以下の結果が得られない。
+論文の結果が再現できない
-_________________________________________________
 ### 欲しい結果
-![イメージ説明](6d24453c3dbd07f5291b0900a8572767.jpeg)
+![20*20の正方格子](7c8308f5057a6e57ebfeea4b0bdf97ee.jpeg)
-### 3次元立方体の量子ウォークを用いたグローバー探索の場合
+[該当論文](https://arxiv.org/abs/1010.4705)
-![イメージ説明](0407161182184a928d49aa398be7e063.png)
+[参考論文２](https://arxiv.org/abs/1203.3950)
-上記の図の|100>の解をマーキングしてその確率を見たい場合、グローバーアルゴリズムのやり方にのっとると、
+[参考論文](https://arxiv.org/abs/1005.3676)
-```
+ちなみに理解しているつもりなんですが、全くできません。コードもエラーなしで動くのですが、、、（ということは何かがわかっていないということなのですが）
-initial_Phi = 1/np.sqrt(2**3*3)*(|000>(|0>+|1>+|2>) + |001>(|0>+|1>+|2>) + |010>(|0>+|1>+|2>) + |100>(|0>+|1>+|2>) + |011>(|0>+|1>+|2>) + |101>(|0>+|1>+|2>) + |110>(|0>+|1>+|2>) + |111>(|0>+|1>+|2>))
+要点は初期状態が等しい重ね合わせ状態であることと、メイン計算は欲しいい点の確率振幅にマイナスを毎回かけながらすべての格子点に拡散行列とフリップフロップシフトというシフトの仕方が必要なようです。
-となり、欲しい解をマーキング→コイン演算子→シフト演算子をtstep繰り返します。今は2回繰り返すと、
+### イメージ図
-#marking
+![正方格子量子探索](7ae71615029fd007225a2f924e5c1e7c.png)
-Phi(t=1) = 1/np.sqrt(2**3*3)*(|000>(|0>+|1>+|2>) + |001>(|0>+|1>+|2>) + |010>(|0>+|1>+|2>) - |100>(|0>+|1>+|2>) + |011>(|0>+|1>+|2>) + |101>(|0>+|1>+|2>) + |110>(|0>+|1>+|2>) + |111>(|0>+|1>+|2>))
+上の図は3×③の正方格子点を表しており、赤い点がマーキングされた点です（マーキングされた点は任意です）。欲しい結果にある結果は190番目の格子点にマーキングをした20×20の正方格子の結果です。
-#coin
-C = [[-1,2,2],[2,-1,2],[2,2,-1]]/3
-各方向への重み
-|0>=1/3*(-|0>+2|1>+2|2>)
-|1>=1/3*(2|0>-|1>+2|2>)
-|2>=1/3*(2|0>+2|1>-|2>)
-これを上に代入すると、
-Phi(t=1) = 1/np.sqrt(2**3*3)*(|000>(|0>+|1>+|2>) + |001>(|0>+|1>+|2>) + |010>(|0>+|1>+|2>) - |100>(|0>+|1>+|2>) + |011>(|0>+|1>+|2>) + |101>(|0>+|1>+|2>) + |110>(|0>+|1>+|2>) + |111>(|0>+|1>+|2>))
-#Shift
-|000>|0> → |001>|0>
-|000>|1> → |010>|1>
-|000>|2> → |100>|2>
-とすると、
-Phi(t=1) = 1/np.sqrt(2**3*3)*(|000>(|0>+|1>+|2>) + |001>(|0>+|1>+|2>) + |010>(|0>+|1>+|2>) - |100>(|0>+|1>+|2>) + |011>(|0>+|1>-|2>) + |101>(-|0>+|1>+|2>) + |110>(|0>-|1>+|2>) + |111>(|0>+|1>+|2>))
-同様に繰り返すと、
-phi(t=2) = 1/3*np.sqrt(2**3*3)*(|000>(3|0>+3|1>-3|2>) + |001>(-|0>+3|1>-2|2>) + |010>(3|0>-|1>-2|2>) - |100>(5|0>+5|1>+5|2>) + |011>(3|0>+3|1>+3|2>) + |101>(-3|0>+3|1>+3|2>) + |110>(3|0>-3|1>+3|2>) + |111>(-|0>-|1>+3|2>))
-と、マーキングした|100>が増幅されていることがわかります。
-```
-これを2次元格子且つ、4方向の20*20で考えた場合上の欲しい結果と一致しません。
+### 実行してみた
 ```python
@@ -74,61 +32,41 @@
 import itertools
-##################
+N = 20
-n = 10
+d = 4
-N = (2*n)**2
+itr =100
-step = 100
+#拡散行列
-x_list = [i for i in range(0,2*n)]
+C = np.array([[-1,1,1,1],[1,-1,1,1],[1,1,-1,1],[1,1,1,-1]])/2
-y_list = [i for i in range(0,2*n)]
+x_list = [i for i in range(N)]
-t_list = [i for i in range(0,step+1)]
+y_list = [i for i in range(N)]
-p_list = []
+time = [i for i in range(itr)]
-####################
+#フリップフロップシフトに場合分け（ここが怪しい？）
-G = np.array([[-1,1,1,1],[1,-1,1,1],[1,1,-1,1],[1,1,1,-1]]) / 2
+#right =[1,0,0,0],down=[0,1,0,0],left=[0,0,1,0],up=[0,0,0,1]
-P = np.zeros((4,4)); P[0,:] = G[0,:]
+right= np.zeros((4,4));right[2,:] = C[2,:] #stage[x,y][0]
-Q = np.zeros((4,4)); Q[1,:] = G[1,:]
+down = np.zeros((4,4));down[3,:] = C[3,:] #stage[x,y][1]
-R = np.zeros((4,4)); R[2,:] = G[2,:]
+left = np.zeros((4,4));left[0,:] = C[0,:] #stage[x,y][2]
-S = np.zeros((4,4)); S[3,:] = G[3,:]
+up = np.zeros((4,4));up[1,:] = C[1,:] #stage[x,y][3]にそれぞれ作用させる
-####################
+#初期状態
-O = np.array([[1,-1,-1,-1],[-1,1,-1,-1],[-1,-1,1,-1],[-1,-1,-1,1]]) / 2
+psi = np.ones((N,N,4),dtype="float")/np.sqrt(N*N*4)
-p = np.zeros((4,4)); p[0,:] = O[0,:]
+prob =[]
-q = np.zeros((4,4)); q[1,:] = O[1,:]
-r = np.zeros((4,4)); r[2,:] = O[2,:]
-s = np.zeros((4,4)); s[3,:] = O[3,:]
-###############
-###初期設定
-#initphi
-phi_map = np.ones((2*n, 2*n,4),dtype="complex")
-phi_map /= np.sqrt(4*N)
-next_phi_map = np.zeros((2*n, 2*n,4),dtype="complex")
-######################
-for t in range(0,step+1):
+for t in range(itr):
     if t == 0:
@@ -136,74 +74,42 @@
     else:
-        for i in itertools.product(x_list,y_list):
+        next_psi = np.zeros((N,N,4),dtype="float")　#t+1秒後のpsiを収納する
+        psi[2,2] *= -1 #マーキング
+        for i in itertools.product(x_list,y_list):　#すべての格子点で計算
             x = i[0]
             y = i[1]
-            #q
+            #boundary condition
-            x1 = (x-1 + 2*n) % (2*n)
+            x0 = (x+1) % N
-            #p
+            x1 = (x-1 + N) % N
-            x2 = (x+1) % (2*n)
+            y0 = (y+1) % N
-            #s
+            y1 = (y-1) % N
-            y1 = (y-1 + 2*n) % (2*n)
+            next_psi[x,y] = np.copy(np.array(np.dot(right,psi[x1,y])+np.dot(down,psi[x,y0])+np.dot(left,psi[x0,y])+ np.dot(up,psi[x,y1])))
-            #r
+        psi = np.copy(next_psi)
-            y2 = (y+1) % (2*n)
+    #print(t,psi)
+    prob.append(np.inner(psi[2,2],psi[2,2]))
-            if i == (12,12):
-                next_phi_map[i] = np.copy(np.array([np.dot(p, phi_map[x2,y]) + np.dot(q, phi_map[x1,y]) + np.dot(r, phi_map[x,y2])+ np.dot(s, phi_map[x,y1])]))
-            else:
-                next_phi_map[i] = np.copy(np.array([np.dot(P, phi_map[x2,y]) + np.dot(Q, phi_map[x1,y]) + np.dot(R, phi_map[x,y2])+ np.dot(S, phi_map[x,y1])]))
-            #p_map[i] = np.real(np.vdot(next_phi_map[i], next_phi_map[i]))
-        phi_map = np.copy(next_phi_map)
-    print(t, phi_map)
-    p_list.append(np.real(np.vdot(phi_map[12,12] , phi_map[12,12])))
-########
-plt.xlabel("t",fontsize="24")
-plt.ylabel("probability",fontsize="24")
-plt.plot(t_list,p_list,color="red",label="quantum walk",linewidth="1")
-plt.legend(title="markedvertexprobablity",loc="best",fontsize=10)
-plt.tight_layout()
+plt.plot(time,prob)
 plt.show()
 ```
-![イメージ説明](c8c4eec927043ae73aca2f61c3548fe2.png)
+![イメージ説明](3b192866ddf979e8d53194625f9b8147.png)
-どなたか、上記の理想結果を実装できる方、または量子情報に詳しい方がいらっしゃいましたら、是非ご教授下さい。。
+こんなへんちくりんな結果になります。
-個人ごとですが、結構困っています。宜しくお願い致します。。
+どこが悪いのかわかりません！！　教えてくれる人も周りにいません。。。ご教授のほどよろしくお願い致します。
-この質問はstackoverflowでもすることにします。
-これはほぼ諦め状態です。（なんでできないのかわからないです。。）
-![イメージ説明](f6db0c9b740d730c1e799a23a644bd0c.png)

追加

2018/10/02 04:47

投稿

Fallout_18

スコア124

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -203,3 +203,7 @@
 個人ごとですが、結構困っています。宜しくお願い致します。。
 この質問はstackoverflowでもすることにします。
+これはほぼ諦め状態です。（なんでできないのかわからないです。。）
+![イメージ説明](f6db0c9b740d730c1e799a23a644bd0c.png)

追加

2018/08/31 14:28

投稿

Fallout_18

スコア124

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -201,3 +201,5 @@
 どなたか、上記の理想結果を実装できる方、または量子情報に詳しい方がいらっしゃいましたら、是非ご教授下さい。。
 個人ごとですが、結構困っています。宜しくお願い致します。。
+この質問はstackoverflowでもすることにします。

内容変更

2018/08/18 14:44

投稿

Fallout_18

スコア124

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -4,53 +4,67 @@
 ### 問題点
-論文通りの実装が得られない。
+やり方をわかっているつもりだが、以下の結果が得られない。
-###参考にした試料
-もっと詳細に知りたい方は、以下の論文の3章を参考にしていただきますようお願い致します。
-[Spatial search using the discrete time quantum walk](https://www.arxiv-vanity.com/papers/1010.4705/)
-### グローバーアルゴリズムとは
-以下にグローバーアルゴリズムの流れをポストします。
-______________________________________________
-![イメージ説明](cd8f2c0376cd6947f40c8e6f8874964d.png)
-初期状態については、格子内の全ての可能なサイトを均等に重ね合わせて量子歩行者を開始し、すべての方向を等しく重ね合わせる。
 _________________________________________________
-初期値は原点に[1000]をもっていき、それ以外は[0000]とします。
 ### 欲しい結果
-![](487303f840156a0634d15a4cedd44ccd.jpeg)
-これは正確な答えではありません、あくまでイメージです。黒い線のような(maxが0.25として)周期をとります。
-（上の参考文献に正確な欲しい結果が掲載されています。そちらを参考にしていただければ幸いです。）
-### 追記
-どうやら、上のコインを作用させた後に、フリップフロップウォークという歩き方でシフト演算子を作用させなければいけなかったみたいです。フリップフロップウォークとは
-|→>|x,y> → |←>|x+1,y>
-|←>|x,y> → |→>|x+1,y>
-|↑>|x,y> → |↓>|x,y+1>
-|↓>|x,y> → |↑>|x,y-1>
-のように変化前の内部状態と変化後の内部状態が逆になるように歩かせる必要があります。
-以下を踏まえて。。。
-### 実装
+![イメージ説明](6d24453c3dbd07f5291b0900a8572767.jpeg)
+### 3次元立方体の量子ウォークを用いたグローバー探索の場合
+![イメージ説明](0407161182184a928d49aa398be7e063.png)
+上記の図の|100>の解をマーキングしてその確率を見たい場合、グローバーアルゴリズムのやり方にのっとると、
+```
+initial_Phi = 1/np.sqrt(2**3*3)*(|000>(|0>+|1>+|2>) + |001>(|0>+|1>+|2>) + |010>(|0>+|1>+|2>) + |100>(|0>+|1>+|2>) + |011>(|0>+|1>+|2>) + |101>(|0>+|1>+|2>) + |110>(|0>+|1>+|2>) + |111>(|0>+|1>+|2>))
+となり、欲しい解をマーキング→コイン演算子→シフト演算子をtstep繰り返します。今は2回繰り返すと、
+#marking
+Phi(t=1) = 1/np.sqrt(2**3*3)*(|000>(|0>+|1>+|2>) + |001>(|0>+|1>+|2>) + |010>(|0>+|1>+|2>) - |100>(|0>+|1>+|2>) + |011>(|0>+|1>+|2>) + |101>(|0>+|1>+|2>) + |110>(|0>+|1>+|2>) + |111>(|0>+|1>+|2>))
+#coin
+C = [[-1,2,2],[2,-1,2],[2,2,-1]]/3
+各方向への重み
+|0>=1/3*(-|0>+2|1>+2|2>)
+|1>=1/3*(2|0>-|1>+2|2>)
+|2>=1/3*(2|0>+2|1>-|2>)
+これを上に代入すると、
+Phi(t=1) = 1/np.sqrt(2**3*3)*(|000>(|0>+|1>+|2>) + |001>(|0>+|1>+|2>) + |010>(|0>+|1>+|2>) - |100>(|0>+|1>+|2>) + |011>(|0>+|1>+|2>) + |101>(|0>+|1>+|2>) + |110>(|0>+|1>+|2>) + |111>(|0>+|1>+|2>))
+#Shift
+|000>|0> → |001>|0>
+|000>|1> → |010>|1>
+|000>|2> → |100>|2>
+とすると、
+Phi(t=1) = 1/np.sqrt(2**3*3)*(|000>(|0>+|1>+|2>) + |001>(|0>+|1>+|2>) + |010>(|0>+|1>+|2>) - |100>(|0>+|1>+|2>) + |011>(|0>+|1>-|2>) + |101>(-|0>+|1>+|2>) + |110>(|0>-|1>+|2>) + |111>(|0>+|1>+|2>))
+同様に繰り返すと、
+phi(t=2) = 1/3*np.sqrt(2**3*3)*(|000>(3|0>+3|1>-3|2>) + |001>(-|0>+3|1>-2|2>) + |010>(3|0>-|1>-2|2>) - |100>(5|0>+5|1>+5|2>) + |011>(3|0>+3|1>+3|2>) + |101>(-3|0>+3|1>+3|2>) + |110>(3|0>-3|1>+3|2>) + |111>(-|0>-|1>+3|2>))
+と、マーキングした|100>が増幅されていることがわかります。
+```
+これを2次元格子且つ、4方向の20*20で考えた場合上の欲しい結果と一致しません。
 ```python
@@ -64,13 +78,13 @@
 n = 10
-N = 4*(n**2)
+N = (2*n)**2
-step = 30
+step = 100
-x_list = [i for i in range(0,2*n+1)]
+x_list = [i for i in range(0,2*n)]
-y_list = [i for i in range(0,2*n+1)]
+y_list = [i for i in range(0,2*n)]
 t_list = [i for i in range(0,step+1)]
@@ -80,9 +94,9 @@
 G = np.array([[-1,1,1,1],[1,-1,1,1],[1,1,-1,1],[1,1,1,-1]]) / 2
-P = np.zeros((4,4)); P[1,:] = G[1,:]
+P = np.zeros((4,4)); P[0,:] = G[0,:]
-Q = np.zeros((4,4)); Q[0,:] = G[0,:]
+Q = np.zeros((4,4)); Q[1,:] = G[1,:]
 R = np.zeros((4,4)); R[2,:] = G[2,:]
@@ -92,9 +106,9 @@
 O = np.array([[1,-1,-1,-1],[-1,1,-1,-1],[-1,-1,1,-1],[-1,-1,-1,1]]) / 2
-p = np.zeros((4,4)); p[1,:] = O[1,:]
+p = np.zeros((4,4)); p[0,:] = O[0,:]
-q = np.zeros((4,4)); q[0,:] = O[0,:]
+q = np.zeros((4,4)); q[1,:] = O[1,:]
 r = np.zeros((4,4)); r[2,:] = O[2,:]
@@ -102,19 +116,15 @@
 ###############
-###ここの初期設定があやしい
+###初期設定
+#initphi
-phi_map = np.zeros((2*n+1, 2*n+1,4),dtype="complex")
+phi_map = np.ones((2*n, 2*n,4),dtype="complex")
-phi_map[0,0]= np.array([1,0,0,0])
+phi_map /= np.sqrt(4*N)
-#####初期確率
-next_phi_map = np.zeros((2*n+1, 2*n+1,4),dtype="complex")
+next_phi_map = np.zeros((2*n, 2*n,4),dtype="complex")
-p_map = np.zeros([2*n+1,2*n+1])
-p_map[0,0] = 1.0
 ######################
@@ -128,50 +138,46 @@
         for i in itertools.product(x_list,y_list):
             x = i[0]
             y = i[1]
-            #境界条件
+            #q
-            x1 = (x-1 + 2*n+1) % (2*n+1)
+            x1 = (x-1 + 2*n) % (2*n)
+            #p
-            x2 = (x+1) % (2*n+1)
+            x2 = (x+1) % (2*n)
+            #s
-            y1 = (y-1 + 2*n+1) % (2*n+1)
+            y1 = (y-1 + 2*n) % (2*n)
+            #r
-            y2 = (y+1) % (2*n+1)
+            y2 = (y+1) % (2*n)
-            #マーキング
             if i == (12,12):
-                next_phi_map[i] = np.array([np.dot(p, phi_map[x2,y]) + np.dot(q, phi_map[x1,y])
+                next_phi_map[i] = np.copy(np.array([np.dot(p, phi_map[x2,y]) + np.dot(q, phi_map[x1,y]) + np.dot(r, phi_map[x,y2])+ np.dot(s, phi_map[x,y1])]))
-                + np.dot(r, phi_map[x,y2])+ np.dot(s, phi_map[x,y1])
-            #それ以外
             else:
-                next_phi_map[i] = np.array([np.dot(P, phi_map[x2,y]) + np.dot(Q, phi_map[x1,y])
+                next_phi_map[i] = np.copy(np.array([np.dot(P, phi_map[x2,y]) + np.dot(Q, phi_map[x1,y]) + np.dot(R, phi_map[x,y2])+ np.dot(S, phi_map[x,y1])]))
-                + np.dot(R, phi_map[x,y2])+ np.dot(S, phi_map[x,y1])])
-            p_map[i] = np.real(np.inner(next_phi_map[i], np.conj(next_phi_map[i])))
+            #p_map[i] = np.real(np.vdot(next_phi_map[i], next_phi_map[i]))
         phi_map = np.copy(next_phi_map)
+    print(t, phi_map)
     p_list.append(np.real(np.vdot(phi_map[12,12] , phi_map[12,12])))
 ########
 plt.xlabel("t",fontsize="24")
@@ -186,12 +192,12 @@
 plt.show()
 ```
-### 結果
-![イメージ説明](fe190f0f52a7c553b96575f708cc37b9.png)
+![イメージ説明](c8c4eec927043ae73aca2f61c3548fe2.png)
-周期的になりません。。。
-少し専門寄りなのですが、グローバーのアルゴリズム自体は簡単なルールなので、どなたか再現できる方がいらっしゃいましたら、宜しくお願い致します。。
+どなたか、上記の理想結果を実装できる方、または量子情報に詳しい方がいらっしゃいましたら、是非ご教授下さい。。
+個人ごとですが、結構困っています。宜しくお願い致します。。

修正

2018/08/18 13:42

投稿

Fallout_18

スコア124

test CHANGED Viewed

	@@ -1 +1 @@
1	- 量子アルゴリズム~~実装してい~~ま~~す。本当にどこが悪~~いのか~~見当もつか~~ない~~ので助けて欲しいです。~~
1	+ 2次元正方格子上のグローバーアルゴリズムがうまくいかない

test CHANGED Viewed

@@ -4,9 +4,7 @@
 ### 問題点
-本当にどこが悪いのか見当もつかないので助けが欲しいです
+論文通りの実装が得られない。
-（2週間以上進歩せず流石に教えて欲しいのですが、専門寄りでかつ、周りで知っている方がいないという理由で質問させて頂きます）
 ###参考にした試料
@@ -35,6 +33,22 @@
 これは正確な答えではありません、あくまでイメージです。黒い線のような(maxが0.25として)周期をとります。
 （上の参考文献に正確な欲しい結果が掲載されています。そちらを参考にしていただければ幸いです。）
+### 追記
+どうやら、上のコインを作用させた後に、フリップフロップウォークという歩き方でシフト演算子を作用させなければいけなかったみたいです。フリップフロップウォークとは
+|→>|x,y> → |←>|x+1,y>
+|←>|x,y> → |→>|x+1,y>
+|↑>|x,y> → |↓>|x,y+1>
+|↓>|x,y> → |↑>|x,y-1>
+のように変化前の内部状態と変化後の内部状態が逆になるように歩かせる必要があります。
+以下を踏まえて。。。
 ### 実装
@@ -66,9 +80,9 @@
 G = np.array([[-1,1,1,1],[1,-1,1,1],[1,1,-1,1],[1,1,1,-1]]) / 2
-P = np.zeros((4,4)); P[0,:] = G[0,:]
+P = np.zeros((4,4)); P[1,:] = G[1,:]
-Q = np.zeros((4,4)); Q[1,:] = G[1,:]
+Q = np.zeros((4,4)); Q[0,:] = G[0,:]
 R = np.zeros((4,4)); R[2,:] = G[2,:]
@@ -78,9 +92,9 @@
 O = np.array([[1,-1,-1,-1],[-1,1,-1,-1],[-1,-1,1,-1],[-1,-1,-1,1]]) / 2
-p = np.zeros((4,4)); p[0,:] = O[0,:]
+p = np.zeros((4,4)); p[1,:] = O[1,:]
-q = np.zeros((4,4)); q[1,:] = O[1,:]
+q = np.zeros((4,4)); q[0,:] = O[0,:]
 r = np.zeros((4,4)); r[2,:] = O[2,:]
@@ -92,7 +106,7 @@
 phi_map = np.zeros((2*n+1, 2*n+1,4),dtype="complex")
-phi_map[0,0]= np.array([1,1,-1,-1])/2#/(np.sqrt(4*N))
+phi_map[0,0]= np.array([1,0,0,0])
 #####初期確率
@@ -130,7 +144,7 @@
             #マーキング
-            if i == (8,8):
+            if i == (12,12):
                 next_phi_map[i] = np.array([np.dot(p, phi_map[x2,y]) + np.dot(q, phi_map[x1,y])
@@ -154,7 +168,7 @@
-    p_list.append(np.real(np.vdot(phi_map[8,8] , phi_map[8,8])))
+    p_list.append(np.real(np.vdot(phi_map[12,12] , phi_map[12,12])))
@@ -176,36 +190,8 @@
 ### 結果
-![](d798953e40a861905b167b4282fe6829.png)
+![イメージ説明](fe190f0f52a7c553b96575f708cc37b9.png)
-上の結果と違いすぎで何がなんだがもうお手上げ状態です。
+周期的になりません。。。
-上述した初期状態の取り方の意味がよくわかっていません(参考論文に書いてあるものをgoogle翻訳しました)。
-前述しましたが、周りにだれも教えてくれない中一人で奮闘している状況です。。どんなことでも良いので、ご意見頂戴したいです。。
+少し専門寄りなのですが、グローバーのアルゴリズム自体は簡単なルールなので、どなたか再現できる方がいらっしゃいましたら、宜しくお願い致します。。
-何卒宜しくお願い致します。
-追記】
-初期状態において、格子内の全ての可能なサイトに対して等しく重ね合わせると書いてあるので、上のコードの初期状態を
-```python
-phi_map = np.ones((2*n+1, 2*n+1,4),dtype="complex")
-phi_map /= np.linalg.norm(phi_map)
-#p_map[0,0]=1
-```
-と変更してみました。結果は
-![イメージ説明](b78569dbba646cb0b75ac3c831d3dc2d.png)
-。。。。。周期的にならない。。
-本当に、お手上げ状態です。

タグ追加

2018/08/12 04:38

投稿

Fallout_18

スコア124

test CHANGED Viewed

File without changes

test CHANGED Viewed

File without changes

タグ追加

2018/08/06 12:37

投稿

Fallout_18

スコア124

test CHANGED Viewed

File without changes

test CHANGED Viewed

File without changes

追記

2018/08/06 09:53

投稿

Fallout_18

スコア124

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -198,6 +198,8 @@
 phi_map /= np.linalg.norm(phi_map)
+#p_map[0,0]=1
 ```
 と変更してみました。結果は

追記

2018/08/04 16:39

投稿

Fallout_18

スコア124

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -185,3 +185,25 @@
 前述しましたが、周りにだれも教えてくれない中一人で奮闘している状況です。。どんなことでも良いので、ご意見頂戴したいです。。
 何卒宜しくお願い致します。
+追記】
+初期状態において、格子内の全ての可能なサイトに対して等しく重ね合わせると書いてあるので、上のコードの初期状態を
+```python
+phi_map = np.ones((2*n+1, 2*n+1,4),dtype="complex")
+phi_map /= np.linalg.norm(phi_map)
+```
+と変更してみました。結果は
+![イメージ説明](b78569dbba646cb0b75ac3c831d3dc2d.png)
+。。。。。周期的にならない。。
+本当に、お手上げ状態です。

欲しい結果を明記しなおした

2018/08/04 16:35

投稿

Fallout_18

スコア124

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -32,7 +32,9 @@
 ![](487303f840156a0634d15a4cedd44ccd.jpeg)
+これは正確な答えではありません、あくまでイメージです。黒い線のような(maxが0.25として)周期をとります。
-これが絶対ではないのですが、黒い線のような周期をとります。（上の参考文献に正確な欲しい結果が掲載されていますお手数ですが、宜しくお願い致します。）
+（上の参考文献に正確な欲しい結果が掲載されています。そちらを参考にしていただければ幸いです。）
 ### 実装

情報追加

2018/08/04 13:08

投稿

Fallout_18

スコア124

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -32,7 +32,7 @@
 ![](487303f840156a0634d15a4cedd44ccd.jpeg)
-これが絶対ではないのですが、黒い線のような周期をとります。
+これが絶対ではないのですが、黒い線のような周期をとります。（上の参考文献に正確な欲しい結果が掲載されていますお手数ですが、宜しくお願い致します。）
 ### 実装