量子ウォーク（アダマールウォーク）の実装について

Question

あまり、よろしくない質問の仕方だとは思いますが。。。
(http://d.hatena.ne.jp/ryamada/20151221)のサイトにある、
R言語での実装をpythonで行いたいのですが、以下のR言語（R？Ruby？）をpythonに翻訳したいです。

U <- 1/sqrt(2) * matrix(c(1,1,1,-1),byrow=TRUE,2,2) + 1i*0

my.PQ <- function(U){
	P <- Q <- matrix(0,2,2)
	P[1,] <- U[1,]
	Q[2,] <- U[2,]
	return(list(P=P,Q=Q))
}


Xinit <- 1/sqrt(2) * c(1,1i)

my.HadamardWalk <- function(U,Xinit,n){
	PQ <- my.PQ(U)
	P <- PQ[[1]]
	Q <- PQ[[2]]
	ret <- list()
	ret[[1]] <- matrix(Xinit,nrow=1)
	
	ret2 <- list()
	ret2[[1]] <- apply(Mod(ret[[1]])^2,1,sum)
	ret3 <- matrix(0,n,n)
	ret3[1,1] <- ret2[[1]]
	for(i in 2:n){
		ret[[i]] <- matrix(0,i,2)
		tmp1 <- matrix(ret[[i-1]][1:(i-1),],ncol=2)
		tmp2 <- matrix(ret[[i-1]][1:(i-1),],ncol=2)
		p <- P %*% t(tmp1)
		q <- Q %*% t(tmp2)
		ret[[i]][1:(i-1),] <- ret[[i]][1:(i-1),] + t(p)
		ret[[i]][2:i,] <- ret[[i]][2:i,] + t(q)
		ret2[[i]] <- apply(Mod(ret[[i]])^2,1,sum)
		ret3[i,1:i] <- ret2[[i]]
	}
	return(list(Q=ret,prob=ret2,probMat=ret3))
}

hout <- my.HadamardWalk(U,Xinit,50)
image(hout[[3]])
persp(hout[[3]])

](7d272e84fb9cdc17866c140daf4a4614.jpeg)
一応、私なりにですが、自分なりに努力している所を挙げます。

sympy.physics.quantum.qubit import Qubit
import sympy
import numpy as np
x = Qubit("00")
y = Qubit("01")

def hadamardcoin():
    hadamardcoin = np.array([[1 / sqrt(2), 1 / sqrt(2), 0, 0], [1 / sqrt(2), -1 / sqrt(2), 0, 0]])
    return hadamardcoin
def initQuanStat(x0,x1,y0,y1):
    initQuanStat = np.zeros([4,4],complex) #0+0jの4,4行列
    initQuanStat[0][0] = x0
    initQuanStat[1][1] = x1
    initQuanStat[2][2] = y0
    initQuanStat[3][3] = y1
    return initQuanStat
def initPositionMap(steps):
    positionMap = np.zeros([2*steps + 1, 2*steps + 1, 4, 4],complex)
    return positionMap
def coinOperator(positionMap,coin):
    dimension = shape(positonMap)[0]
    for i in range(dimension):
        for j in range(dimension):
            positonMap[i][j] = dot(positonMap[i][j],coin)
        return positionMap
def shiftOperator(coinMap,step):
    newPositionMap = initPositionMap(step)
＃まだ途中です

Githubなど使って頑張ってはいますが、初心者にはなかなかハードなものだと痛感しており、力を貸していただきたい気持ちです。（私の上記のコードを本題とは少しズレておりますので、あまり気にしないでください）
宜しくお願い致します。

Comments

[Fallout_18]

低評価のオンパレードや～
確かに、すいません。丸投げでした、自分で頑張ります。。

[KojiDoi]

まず、RとRubyの違いを理解するところから始めましょうか。

[Fallout_18]

その通りです

Answer 1

自分でしっかり書き直しました笑

Answer 2

Rで正しく動くコードがあるのであれば、pythonからRを起動するコードを研究したほうがいいのではないでしょうか？
Rのコードを移植したいのであれば、Rのコードを正しく理解するところから始めたほうがいいと思います。その上で、Rで実装できているがpythonだとどの関数を使えばいいのかわからない部分に限定して質問したほうが、回答を得やすくなると思います。

Comments

[Fallout_18]

アドバイスありがとうございます。
pythonとRを併用するやり方があることさえ知らなかったです。。
本当にありがとうございます、頑張ってみます！

Answer 3

目指していることはR言語とpythonを書けるようになってR言語のコードをpythonを移植できるようになることでしょうか。
それとも、物理・数式を理解して、それをpythonのコードに落とし込むことでしょうか。

努力していると言われましても、載っているコードにコメントは特になく、どこまで理解しているのかわかりません。
せめてどう思って、何を試して、どううまくいっていないのかを整理してみてください。

http://www.kisc.meiji.ac.jp/~tz14040/quantumwalk/english/
はじめのうちは順番に結果を確認しながらステップバイステップで実装して行くと、どこまでわかってどこからがわからないかがわかると思います。

Comments

[Fallout_18]

そうですね、相手に明確に伝えることなく投げやりな形になって申し訳なかったです。
物理・数式を理解してコードに落とし込むことを行えることが理想です。
一つ一つ理解しながら、やっていきます。
いつもありがとうございます。