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")