mpi.bcast.cmd(np.mpi.initialize(), caller.execute = TRUE)
mpi.bcast.cmd(options(np.messages = FALSE), caller.execute = TRUE)

.np_demo_src <- Sys.getenv("NP_DEMO_SRC", "")
.np_demo_family <- c(if (nzchar(.np_demo_src)) file.path(.np_demo_src, "..", "inst", "demo_family_npcdens.R"),
                     system.file("demo_family_npcdens.R", package = "npRmpi"))
.np_demo_family <- .np_demo_family[nzchar(.np_demo_family) & file.exists(.np_demo_family)]
source(.np_demo_family[[1L]])

npcdens_demo_run_matrix("profile")

mpi.bcast.cmd(mpi.quit(), caller.execute = TRUE)
