tests/testthat/test_a_begin.R

In CharlotteJana/pdmpsim: Simulate Piecewise Deterministic Markov Processes

#-------- simulate pdmps and save the results for testing -----------

if(!dir.exists(paste0(tempdir(), "/pdmpsimtest")))
  dir.create(paste0(tempdir(), "/pdmpsimtest"))
unlink(paste0(tempdir(), "/pdmpsimtest/*"))

simplePdmp <- pdmpModel(
  descr = "a simple pdmp",
  init = c(f = 0, d = 0),
  times = c(from = 0, to = 10, by = 0.1),
  discStates = list(d = -1:1),
  dynfunc = function(t, x, parms) c(x["d"], 0),
  ratefunc = function(t, x, parms) c(1+x["d"], 1-x["d"]),
  jumpfunc = function(t, x, parms, jtype){
    c(0, switch(jtype, x["d"]-1, x["d"]+1))
  }
)

suppressMessages(try({
  multSimCsv(simplePdmp, seeds = 1:5, 
             prefix = paste0(tempdir(), "/pdmpsimtest/","test"))
  saveRDS(multSim(simplePdmp, seeds = 1:5), 
          file = paste0(tempdir(), "/pdmpsimtest/", "test_MultSim.rda"))
}, silent = TRUE))