--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/parallelSticsRRotation.R	Tue Feb 21 19:03:38 2023 +0000
@@ -0,0 +1,65 @@
+library(SticsRFiles)
+library(SticsOnR)
+library(parallel)
+library(doParallel)
+library(readxl)
+
+workspace <- paste0(getwd(), "/WS")
+javastics_path <- paste0(getwd(), "/JavaSTICS-1.5.1-STICS-10.0.0")
+
+# copy of the var.mod and rap.mod into the config folder
+file.copy(from = file.path(workspace, "var.mod"), to = "config/var.mod", overwrite = TRUE)
+file.copy(from = file.path(workspace, "rap.mod"), to = "config/rap.mod", overwrite = TRUE)
+
+txt_path <- paste0(getwd(), "/WS/txt_files")
+dir.create(txt_path,recursive = T)
+
+print("==================================")
+
+print(paste0("Cores=", detectCores()))
+print(paste0("javastics=", javastics_path))
+print(paste0("ws=", workspace))
+print(paste0("out_dir=", txt_path))
+
+print("==================================")
+
+
+# List of successive USMs
+successive_usms_file <- ("rotationPlan.xlsx")
+successive_usms_tab <- read_excel(path = file.path(workspace, successive_usms_file), sheet = "usms")
+
+no_cores <- detectCores() - 1  
+cl <- makeCluster(no_cores)  
+registerDoParallel(cl)
+
+result <- foreach(i=1:nrow(successive_usms_tab),
+                  .packages = c("SticsRFiles", "SticsOnR") ) %dopar% {
+  
+  nusms <- as.numeric(successive_usms_tab[i,2])
+  
+  vec <- as.character(successive_usms_tab[i,c(3:(2+nusms))])
+  
+  list_successive_usms <- list(vec)
+  
+  sim_options <- stics_wrapper_options(javastics = javastics_path, 
+                                       stics_exe = paste0(javastics_path,"/bin/stics_modulo"),
+                                       workspace = txt_path, 
+                                       verbose = TRUE, 
+                                       successive = list_successive_usms)
+  list_usms <- vec
+  
+  stics_wrapper(model_options = sim_options, situation = list_usms)
+                  }
+
+stopCluster(cl)  
+
+for (i in 1:nrow(successive_usms_tab)) { #i=1
+  if (i == 1) {
+    simPCP <- result[[i]]
+  } else {
+    simPCP$sim_list <- c(simPCP$sim_list, result[[i]]$sim_list)
+  }
+
+}
+
+save(simPCP, file = "results.RData")
\ No newline at end of file