library(PerMallows) #elk?sz?ti a a kappa vektorokat m=1,2,..., N-re kappasorozat=function(N){ kappasorozat=rep(0/0,N^2) dim(kappasorozat)=c(N,N) kappasorozat[1,1]=1 for (m in 2:N){ kappasorozat[m,1:(m-1)]=kappasorozat[(m-1),1:(m-1)] v=runif.permutation(n = 1, m)[1] kappasorozat[m,v]=1 if (vseged[(i-1),(i-1)])})) } print(date()) proba=sapply(c(1:1000),function(i){nov(3000)}) print(date()) print(date()) proba2=sapply(c(1:2000),function(i){nov(10000)}) print(date()) no_cores=detectCores()-1 cl=makeCluster(no_cores) clusterExport(cl, "nov") clusterExport(cl, "kappasorozat") clusterExport(cl, "runif.permutation") print(date()) proba2=parSapply(cl,c(1:1000),function(i){nov(10000)}) print(date()) stopCluster(cl) par(mfrow=c(2,1)) hist(proba) hist(proba2)