-
- model = list()
- tableauRecap = array(0, dim=c(1000,4))
- cpt = 0
- print("main loop: over all k and all lambda")
-
- for (k in kmin:kmax)
+
+ tableauRecap = list()
+ if (verbose)
+ print("main loop: over all k and all lambda")
+
+ if (ncores_outer > 1)
+ {
+ cl = parallel::makeCluster(ncores_outer)
+ parallel::clusterExport( cl=cl, envir=environment(), varlist=c("X","Y","procedure",
+ "selecMod","gamma","mini","maxi","eps","kmin","kmax","rang.min","rang.max",
+ "ncores_outer","ncores_inner","verbose","p","m","k","tableauRecap") )
+ }
+
+ # Compute model with k components
+ computeModel <- function(k)