feat(speed): set log_() level so only record when verbose - massive speed increase! 1m16s -> 0m33s

Author: amyheather

R/create_asu_trajectory.R

@@ -26,7 +26,7 @@ create_asu_trajectory <- function(env, patient_type, param) {
   # Set up simmer trajectory object...
   trajectory(paste0("ASU_", patient_type, "_path")) |>
 
-    log_("🚶 Arrived at ASU") |>
+    log_("🚶 Arrived at ASU", level = 1) |>
 
     seize("asu_bed", 1L) |>
 
@@ -42,7 +42,7 @@ create_asu_trajectory <- function(env, patient_type, param) {
       dest <- dest_names[dest_index]
       # Create log message
       paste0("🎯 Planned ASU -> ", dest_index, " (", dest, ")")
-    }) |>
+    }, level = 1) |>
 
     set_attribute("asu_los", function() {
       # Retrieve attribute, and use to get post-ASU destination as a string
@@ -75,11 +75,11 @@ create_asu_trajectory <- function(env, patient_type, param) {
     log_(function() {
       paste0("⏳ ASU length of stay: ",
              round(get_attribute(env, "asu_los"), 3L))
-    }) |>
+    }, level = 1) |>
 
     timeout(function() get_attribute(env, "asu_los")) |>
 
-    log_("🏁 ASU stay completed") |>
+    log_("🏁 ASU stay completed", level = 1) |>
 
     release("asu_bed", 1L) |>
 

R/create_rehab_trajectory.R

@@ -26,7 +26,7 @@ create_rehab_trajectory <- function(env, patient_type, param) {
   # Set up simmer trajectory object...
   trajectory(paste0("rehab_", patient_type, "_path")) |>
 
-    log_("🚶 Arrived at rehab") |>
+    log_("🚶 Arrived at rehab", level = 1) |>
 
     seize("rehab_bed", 1L) |>
 
@@ -42,7 +42,7 @@ create_rehab_trajectory <- function(env, patient_type, param) {
       dest <- dest_names[dest_index]
       # Create log message
       paste0("🎯 Planned rehab -> ", dest_index, " (", dest, ")")
-    }) |>
+    }, level = 1) |>
 
     set_attribute("rehab_los", function() {
       # Retrieve attribute, and use to get post-rehab destination as a string
@@ -74,11 +74,11 @@ create_rehab_trajectory <- function(env, patient_type, param) {
     log_(function() {
       paste0("⏳ Rehab length of stay: ",
              round(get_attribute(env, "rehab_los"), 3L))
-    }) |>
+    }, level = 1) |>
 
     timeout(function() get_attribute(env, "rehab_los")) |>
 
-    log_("🏁 Rehab stay completed") |>
+    log_("🏁 Rehab stay completed", level = 1) |>
 
     release("rehab_bed", 1L)
 }

R/get_occupancy_stats.R

@@ -54,7 +54,7 @@ get_occupancy_stats <- function(occupancy) {
     occ_stats[["pct"]] <- occ_stats[["freq"]] / sum(occ_stats[["freq"]])
     occ_stats[["c_pct"]] <- cumsum(occ_stats[["pct"]])
 
-    # Calculate probability of delay
+    # Calculate probability of delay using the Erlang loss formula
     occ_stats[["prob_delay"]] <- occ_stats[["pct"]] / occ_stats[["c_pct"]]
 
     # Calculate 1 in every n patients delayed

R/model.R

@@ -22,15 +22,17 @@ model <- function(run_number, param, set_seed = TRUE) {
   }
 
   # Determine whether to get verbose activity logs
-  verbose <- any(c(param[["log_to_console"]], param[["log_to_file"]]))
+  param[["verbose"]] <- any(c(param[["log_to_console"]],
+                              param[["log_to_file"]]))
 
   # Transform LOS parameters to lognormal scale
   param[["asu_los_lnorm"]] <- transform_to_lnorm(param[["asu_los"]])
   param[["rehab_los_lnorm"]] <- transform_to_lnorm(param[["rehab_los"]])
 
   # Create simmer environment - set verbose to FALSE as using custom logs
   # (but can change to TRUE if want to see default simmer logs as well)
-  env <- simmer("simulation", verbose = FALSE)
+  env <- simmer("simulation", verbose = FALSE,
+                log_level = if (param[["verbose"]]) 1 else 0)
 
   # Add ASU and rehab direct admission patient generators
   for (unit in c("asu", "rehab")) {
@@ -70,7 +72,7 @@ model <- function(run_number, param, set_seed = TRUE) {
   )
 
   # Save and/or display the log
-  if (isTRUE(verbose)) {
+  if (isTRUE(param[["verbose"]])) {
     # Create full log message by adding parameters
     param_string <- paste(names(param), param, sep = "=", collapse = ";\n ")
     full_log <- append(c("Parameters:", param_string, "Log:"), sim_log)