- Rank = matrix(0, nrow=Size, ncol=k)
- for(r in 1:k)
- {
- # On veut le tableau de toutes les combinaisons de rangs possibles
- # Dans la première colonne : on répète (rangmax-rangmin)^(k-1) chaque chiffre :
- # ça remplit la colonne
- # Dans la deuxieme : on répète (rangmax-rangmin)^(k-2) chaque chiffre,
- # et on fait ça (rangmax-rangmin)^2 fois
- # ...
- # Dans la dernière, on répète chaque chiffre une fois,
- # et on fait ça (rangmin-rangmax)^(k-1) fois.
- Rank[,r] = rangmin + rep(0:(deltaRank-1), deltaRank^(r-1), each=deltaRank^(k-r))
+ RankLambda = matrix(0, nrow=Size*L, ncol=k+1)
+ for (r in 1:k)
+ {
+ # On veut le tableau de toutes les combinaisons de rangs possibles, et des lambdas
+ # Dans la première colonne : on répète (rank.max-rank.min)^(k-1) chaque chiffre :
+ # ça remplit la colonne
+ # Dans la deuxieme : on répète (rank.max-rank.min)^(k-2) chaque chiffre,
+ # et on fait ça (rank.max-rank.min)^2 fois
+ # ...
+ # Dans la dernière, on répète chaque chiffre une fois,
+ # et on fait ça (rank.min-rank.max)^(k-1) fois.
+ RankLambda[,r] = rep(rank.min + rep(0:(deltaRank-1), deltaRank^(r-1), each=deltaRank^(k-r)), each = L)
+ }
+ RankLambda[,k+1] = rep(1:L, times = Size)
+
+ if (ncores > 1)
+ {
+ cl = parallel::makeCluster(ncores, outfile='')
+ parallel::clusterExport( cl, envir=environment(),
+ varlist=c("A1","Size","Pi","Rho","mini","maxi","X","Y","eps",
+ "Rank","m","phi","ncores","verbose") )