#-------- 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))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.