library(SticsRFiles)
library(SticsOnR)
library(parallel)
library(doParallel)
library(readxl)

workspace <- paste0(getwd(), "/WS")
javastics_path <- getwd()
default_config_path <- paste0(javastics_path, "/config")

# copy of the var.mod and rap.mod into the config folder
# generalizing maybe should be inside a config folder

#file.copy(from = file.path(workspace, "var.mod"),
#          to = file.path(default_config_path, "var.mod"),
#	  overwrite = TRUE)
#file.copy(from = file.path(workspace, "rap.mod"),
#          to = file.path(default_config_path, "rap.mod"),
#	  overwrite = TRUE)

txt_path <- paste0(getwd(), "/WS/txt_files")
#dir.create(txt_path,recursive = T)

#gen_usms_xml2txt(javastics = javastics_path,
#                 workspace = workspace,
#		 out_dir = txt_path,
#		 check = FALSE)

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 <- ("successionPlan.csv")
successive_usms_tab <- read.csv(file = file.path(workspace, successive_usms_file))

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