## Analyze CPU times for mix and full integration solvers rm(list = ls()) gc() # .libPaths("~/svn-StanPmetrics/script/lib") ## FIX ME - Need more iterations. Don't think 1000 is enough. ########################################################################## ## One Compartment Model with 1st Order Absorption ### dd regime data_s_dd <- read.csv(file = "msSimple_dd.csv") N_s_dd <- length(head(data_s_dd) - 1) data_s_dd_num <- rep(0, N_s_dd) data_s_dd_mix <- data_s_dd_num for (i in 1:N_s_dd) data_s_dd_num[i] <- data_s_dd[1, i] for (i in 1:N_s_dd) data_s_dd_mix[i] <- data_s_dd[2, i] mean(data_s_dd_num) # 5.200799e-06 mean(data_s_dd_mix) # 4.921079e-06 ### vv regime data_s_vv <- read.csv(file = "msSimple_vv.csv") N_s_vv <- length(head(data_s_vv) - 1) data_s_vv_num <- rep(0, N_s_vv) data_s_vv_mix <- data_s_vv_num for (i in 1:N_s_vv) data_s_vv_num[i] <- data_s_vv[1, i] for (i in 1:N_s_vv) data_s_vv_mix[i] <- data_s_vv[2, i] mean(data_s_vv_num) # 2.696004e-05 mean(data_s_vv_mix) # 0.0001372418 ########################################################################## ## Friberg-Karlsson Model with 1st Order Absorption ## dd regime data_dd <- read.csv(file = "mixSolverResult_dd.csv") N_dd <- length(head(data_dd) - 1) data_dd_num <- rep(0, N_dd) data_dd_mix <- data_dd_num for (i in 1:N_dd) data_dd_num[i] <- data_dd[1, i] for (i in 1:N_dd) data_dd_mix[i] <- data_dd[2, i] mean(data_dd_num) # 1.83966e-05 mean(data_dd_mix) # 2.756044e-05 ## vv regime data_vv <- read.csv(file = "mixSolverResult_vv.csv") N_vv <- length(head(data_vv) - 1) data_vv_num <- rep(0, N_vv) data_vv_mix <- data_vv_num for (i in 1:N_vv) data_vv_num[i] <- data_vv[1, i] for (i in 1:N_vv) data_vv_mix[i] <- data_vv[2, i] mean(data_vv_num) # 0.0005151409 mean(data_vv_mix) # 0.004688307