fix recompiling models with PM_load

R/PM_fit.R

@@ -430,15 +430,18 @@ PM_fit <- R6::R6Class(
     binary_path = NULL,
     setup_rust_execution = function() {
 
-      if (!file.exists(self$model$binary_path)) {
-        self$model$compile()
-        if (!file.exists(self$model$binary_path)) {
-          cli::cli_abort(c(
-            "x" = "Rust has not been built.",
-            "i" = "Make sure the model was built correcly when calling {.fn PM_model$new}"
-          ))
-        }
-      }
+      # check if compiled and if not, do so
+      self$model$compile()
+
+      # if (is.null(self$model$binary_path)) {
+      #   self$model$compile()
+      #   if (!file.exists(self$model$binary_path)) {
+      #     cli::cli_abort(c(
+      #       "x" = "Rust has not been built.",
+      #       "i" = "Make sure the model was built correcly when calling {.fn PM_model$new}"
+      #     ))
+      #   }
+      # }
       # cwd <- getwd()
       # self$model$write_rust()
 

R/PM_model.R

@@ -106,9 +106,9 @@ PM_model$new <- function(model, ...) {
   } else {
     cli::cli_abort(c("x" = "Non supported model type: {typeof(model)}"))
   }
-  if (getPMoptions()$backend == "rust") {
-    model$compile()
-  }
+  # if (getPMoptions()$backend == "rust") {
+  #   model$compile()
+  # }
   return(model)
 }
 
@@ -789,7 +789,6 @@ PM_model_list <- R6::R6Class("PM_model_list",
                                  if (!is.null(self$binary_path) && file.exists(self$binary_path)) {
                                    return()
                                  }
-
                                  temp_model <- file.path(tempdir(), "temp_model.txt")
                                  self$write_rust(temp_model)
                                  model_path <- tempfile(pattern = "model_", fileext = ".pmx")

R/PM_sim.R

@@ -1017,7 +1017,6 @@ PM_sim <- R6::R6Class(
       if (length(postToUse) > 0) {
         # simulating from posteriors, each posterior matched to a subject
         # need to set theta as the posterior mean or median for each subject
-        # code here...
 
         ans <- NULL
         data_list <- list()
@@ -1306,6 +1305,7 @@ PM_sim <- R6::R6Class(
     # call simulator and process results
     getSim = function(thisPrior, template, mod, noise2){
       thetas <- thisPrior$thetas %>% select(-prob) %>% as.matrix()
+      mod$compile() # check if compiled and if not, do so
       sim_res <- mod$simulate_all(template, thetas)
       sim_res$.id <- template$standard_data$id[match(sim_res$id, template$standard_data$id)]
       sim_res <- sim_res %>% rename(comp = state_index, nsim = spp_index, amt = state) %>%