ダメ出し：ぐるぐる回るために

http://d.hatena.ne.jp/ryamada/20111002　について

Npt <- 10000
n <- 3
X <- matrix(rnorm(Npt*n), Npt, n)
X <- X/sqrt(rowSums(X^2))
Y <- rowSums(X[,1:(n-1)]*X[,2:n])

は

Npt<-10000
n<-3
X<-matrix(rnorm(Npt*n),Npt,n)

X<-X/sqrt(apply(X^2,1,sum))
#apply(X^2,1,sum)

Y<-rep(0,Npt)

for(i in 1:(n-1)){
    Y<-Y+X[,i]*X[,i+1]
}

より，20 倍速い