Mercurial > repos > siwaa > redelac_stics_s

comparison redelacSticsSimulator.R @ 3:f262814e48a3 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://forgemia.inra.fr/redelac/redelac-toolbox/-/tree/acafb1f877e923efa13655229681753dfa1928c0/tools/REDELACSticsSimulator commit acafb1f877e923efa13655229681753dfa1928c0-dirty"
author siwaa
date Mon, 26 Jun 2023 17:25:19 +0000
parents ea8e91473dc8
children 70ae7516e925
comparison
equal deleted inserted replaced
2:ea8e91473dc8 3:f262814e48a3
38 38
39 # List of successive USMs 39 # List of successive USMs
40 successive_usms_file <- ("successionPlan.csv") 40 successive_usms_file <- ("successionPlan.csv")
41 successive_usms_tab <- read.csv(file = file.path(workspace, successive_usms_file), sep =";") 41 successive_usms_tab <- read.csv(file = file.path(workspace, successive_usms_file), sep =";")
42 42
43 no_cores <- detectCores() - 1 43 #no_cores <- detectCores() - 1
44 #no_cores <- 10 44 #no_cores <- 10
45 cl <- makeCluster(no_cores) 45 #cl <- makeCluster(no_cores)
46 registerDoParallel(cl) 46 #registerDoParallel(cl)
47 47
48 result <- foreach(i=1:nrow(successive_usms_tab), 48 result <- list()
49 .packages = c("SticsRFiles", "SticsOnR") ) %dopar% {
50 49
50 #result <- foreach(i=1:nrow(successive_usms_tab),
51 # .packages = c("SticsRFiles", "SticsOnR") ) %dopar% {
52 for(i in 1:nrow(successive_usms_tab)) {
51 53
52 print(paste("coucou ", i)) 54 print(paste("coucou ", i))
53 nusms <- as.numeric(successive_usms_tab[i,2]) 55 nusms <- as.numeric(successive_usms_tab[i,2])
54 print(paste0("coco ", nusms)) 56 print(paste0("coco ", nusms))
55 57
62 workspace = txt_path, 64 workspace = txt_path,
63 verbose = TRUE, 65 verbose = TRUE,
64 successive = list_successive_usms) 66 successive = list_successive_usms)
65 list_usms <- vec 67 list_usms <- vec
66 68
67 stics_wrapper(model_options = sim_options, situation = list_usms) 69 result = append(result, stics_wrapper(model_options = sim_options, situation = list_usms))
68 } 70 }
69 71
70 stopCluster(cl) 72 #stopCluster(cl)
71 73
72 for (i in 1:nrow(successive_usms_tab)) { #i=1 74 for (i in 1:nrow(successive_usms_tab)) { #i=1
73 if (i == 1) { 75 if (i == 1) {
74 simPCP <- result[[i]] 76 simPCP <- result[[i]]
75 } else { 77 } else {