alpha=1
beta=2
N=7
ered9=list()

for (n in 4:6){
  ered9[[n]]=list()
  for (k in 3:min(5,(n-1))){
    print(paste("n=",n," k=",k," idő: ", date(),sep = ""))
    ered9[[n]][[k]]=list()
    ered9[[n]][[k]][[1]]=paste("n = ",n,", k = ",k, sep="")
    seged=baloldalisv(n,k,alpha,beta)
    ered9[[n]][[k]][[2]]=seged
    ered9[[n]][[k]][[3]]=seged[rovidtabl1(k),]
    a=permutations(n,n)
    b=permutations(n,k)
    if (k==3){
      ered9[[n]][[k]][[6]]=ered9[[n]][[k]][[3]]
    }
    
    if (k>3){
      no_cores=detectCores()-1
      cl=makeCluster(no_cores)
      cl=makeCluster(no_cores)
      clusterExport(cl, "a")
      clusterExport(cl, "b")
      clusterExport(cl, "s1")
      clusterExport(cl, "s2")
      clusterExport(cl, "s3")
      clusterExport(cl, "s4")
      clusterExport(cl, "s5b")
      clusterExport(cl, "seged3")
      seged3=t(parSapply(cl,c(1:dim(b)[1]),function(i){t(s5b(a,b[i,]))%*%seged3}))
      stopCluster(cl)
      sseged4=oszlopegyez2(seged3,seged)
      ered9[[n]][[(k-1)]][[4]]=sseged4[[1]]
      ered9[[n]][[k]][[5]]=sseged4[[2]]
      ered9[[n]][[k]][[6]]=ered9[[n]][[k]][[3]][,ered9[[n]][[k]][[5]]]
    }
    no_cores=detectCores()-1
    cl=makeCluster(no_cores)
    cl=makeCluster(no_cores)
    clusterExport(cl, "a")
    clusterExport(cl, "b")
    clusterExport(cl, "s1")
    clusterExport(cl, "s2")
    clusterExport(cl, "s3")
    clusterExport(cl, "s4")
    clusterExport(cl, "s5")
    clusterExport(cl, "seged")
    seged3=parSapply(cl,c(1:dim(a)[1]),function(i){s5(a[i,],b)%*%seged})
    stopCluster(cl)
    seged3=t(seged3)
    print(date())
  }
}

save(ered9,file="ered9.RData")