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

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

txt_path <- paste0(getwd(), "/WS/txt_files")

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), sep =";")

#no_cores <- detectCores() - 1  
#no_cores <- 10  
#cl <- makeCluster(no_cores)  
#registerDoParallel(cl)

result <- list()

#result <- foreach(i=1:nrow(successive_usms_tab),
#                  .packages = c("SticsRFiles", "SticsOnR") ) %dopar% {
 for(i in 1:nrow(successive_usms_tab)) {

  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

  singleResult = stics_wrapper(model_options = sim_options, situation = list_usms)

  result = append(result, list(singleResult))
                  }

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