notes/indexing/outflow_test.R

In bbolker/McMasterPandemic: Pandemic Model

library(McMasterPandemic)
library(tools)
library(dplyr)

spec_version = "0.1.1"
cpp <- file.path("../../inst/tmb", spec_version, "macpan.cpp")
dll <- file_path_sans_ext(cpp)
options(MP_flex_spec_version = spec_version)
options(MP_flex_spec_dll = basename(dll))

compile(cpp)
dyn.load(dynlib(dll))

params <- read_params("ICU1.csv")
state <- make_state(params = params)
start_date = "2021-05-10"
end_date = "2021-12-10"
model <- make_base_model(
  params, state,
  start_date = start_date, end_date = end_date
)
model$do_hazard = FALSE


report <- tmb_fun(model)$report()

unpack(model$tmb_indices$initialization_mapping)

length(model$state[all_to_infected_idx])
length(model$state[all_to_eigen_idx][eigen_to_infected_idx])