diff options
Diffstat (limited to 'vignettes/web_only/compiled_models.html')
| -rw-r--r-- | vignettes/web_only/compiled_models.html | 81 |
1 files changed, 38 insertions, 43 deletions
diff --git a/vignettes/web_only/compiled_models.html b/vignettes/web_only/compiled_models.html index 1d50ba7c..31d062bb 100644 --- a/vignettes/web_only/compiled_models.html +++ b/vignettes/web_only/compiled_models.html @@ -11,7 +11,7 @@ <meta name="author" content="Johannes Ranke" /> -<meta name="date" content="2020-04-02" /> +<meta name="date" content="2020-05-12" /> <title>Performance benefit by using compiled model definitions in mkin</title> @@ -1583,7 +1583,7 @@ div.tocify { <h1 class="title toc-ignore">Performance benefit by using compiled model definitions in mkin</h1> <h4 class="author">Johannes Ranke</h4> -<h4 class="date">2020-04-02</h4> +<h4 class="date">2020-05-12</h4> </div> @@ -1593,8 +1593,6 @@ div.tocify { <p>When using an mkin version equal to or greater than 0.9-36 and a C compiler is available, you will see a message that the model is being compiled from autogenerated C code when defining a model using mkinmod. Starting from version 0.9.49.9, the <code>mkinmod()</code> function checks for presence of a compiler using</p> <pre class="r"><code>pkgbuild::has_compiler()</code></pre> <p>In previous versions, it used <code>Sys.which("gcc")</code> for this check.</p> -<div id="platform-specific-notes" class="section level3"> -<h3>Platform specific notes</h3> <p>On Linux, you need to have the essential build tools like make and gcc or clang installed. On Debian based linux distributions, these will be pulled in by installing the build-essential package.</p> <p>On MacOS, which I do not use personally, I have had reports that a compiler is available by default.</p> <p>On Windows, you need to install Rtools and have the path to its bin directory in your PATH variable. You do not need to modify the PATH variable when installing Rtools. Instead, I would recommend to put the line</p> @@ -1602,55 +1600,55 @@ div.tocify { <p>into your .Rprofile startup file. This is just a text file with some R code that is executed when your R session starts. It has to be named .Rprofile and has to be located in your home directory, which will generally be your Documents folder. You can check the location of the home directory used by R by issuing</p> <pre class="r"><code>Sys.getenv("HOME")</code></pre> </div> -</div> -<div id="comparison-with-eigenvalue-based-solutions" class="section level2"> -<h2>Comparison with Eigenvalue based solutions</h2> -<p>First, we build a simple degradation model for a parent compound with one metabolite.</p> +<div id="comparison-with-other-solution-methods" class="section level2"> +<h2>Comparison with other solution methods</h2> +<p>First, we build a simple degradation model for a parent compound with one metabolite, and we remove zero values from the dataset.</p> <pre class="r"><code>library("mkin", quietly = TRUE) SFO_SFO <- mkinmod( parent = mkinsub("SFO", "m1"), m1 = mkinsub("SFO"))</code></pre> <pre><code>## Successfully compiled differential equation model from auto-generated C code.</code></pre> -<p>We can compare the performance of the Eigenvalue based solution against the compiled version and the R implementation of the differential equations using the benchmark package. In the output of below code, the warnings about zero being removed from the FOCUS D dataset are suppressed.</p> +<pre class="r"><code>FOCUS_D <- subset(FOCUS_2006_D, value != 0)</code></pre> +<p>We can compare the performance of the Eigenvalue based solution against the compiled version and the R implementation of the differential equations using the benchmark package. In the output of below code, the warnings about zero being removed from the FOCUS D dataset are suppressed. Since mkin version 0.9.49.11, an analytical solution is also implemented, which is included in the tests below.</p> <pre class="r"><code>if (require(rbenchmark)) { b.1 <- benchmark( - "deSolve, not compiled" = mkinfit(SFO_SFO, FOCUS_2006_D, - solution_type = "deSolve", - use_compiled = FALSE, quiet = TRUE), - "Eigenvalue based" = mkinfit(SFO_SFO, FOCUS_2006_D, - solution_type = "eigen", quiet = TRUE), - "deSolve, compiled" = mkinfit(SFO_SFO, FOCUS_2006_D, - solution_type = "deSolve", quiet = TRUE), - replications = 3) + "deSolve, not compiled" = mkinfit(SFO_SFO, FOCUS_D, + solution_type = "deSolve", + use_compiled = FALSE, quiet = TRUE), + "Eigenvalue based" = mkinfit(SFO_SFO, FOCUS_D, + solution_type = "eigen", quiet = TRUE), + "deSolve, compiled" = mkinfit(SFO_SFO, FOCUS_D, + solution_type = "deSolve", quiet = TRUE), + "analytical" = mkinfit(SFO_SFO, FOCUS_D, + solution_type = "analytical", + use_compiled = FALSE, quiet = TRUE), + replications = 1, order = "relative", + columns = c("test", "replications", "relative", "elapsed")) print(b.1) - factor_SFO_SFO <- round(b.1["1", "relative"]) } else { - factor_SFO_SFO <- NA print("R package rbenchmark is not available") }</code></pre> -<pre><code>## test replications elapsed relative user.self sys.self -## 3 deSolve, compiled 3 3.148 1.000 3.146 0.000 -## 1 deSolve, not compiled 3 28.920 9.187 28.904 0.001 -## 2 Eigenvalue based 3 4.442 1.411 4.439 0.000 -## user.child sys.child -## 3 0 0 -## 1 0 0 -## 2 0 0</code></pre> -<p>We see that using the compiled model is by a factor of around 9 faster than using the R version with the default ode solver, and it is even faster than the Eigenvalue based solution implemented in R which does not need iterative solution of the ODEs.</p> +<pre><code>## test replications relative elapsed +## 4 analytical 1 1.000 0.186 +## 3 deSolve, compiled 1 1.769 0.329 +## 2 Eigenvalue based 1 2.371 0.441 +## 1 deSolve, not compiled 1 72.183 13.426</code></pre> +<p>We see that using the compiled model is by more than a factor of 10 faster than using deSolve without compiled code.</p> </div> -<div id="model-that-can-not-be-solved-with-eigenvalues" class="section level2"> -<h2>Model that can not be solved with Eigenvalues</h2> -<p>This evaluation is also taken from the example section of mkinfit.</p> +<div id="model-without-analytical-solution" class="section level2"> +<h2>Model without analytical solution</h2> +<p>This evaluation is also taken from the example section of mkinfit. No analytical solution is available for this system, and now Eigenvalue based solution is possible, so only deSolve using with or without compiled code is available.</p> <pre class="r"><code>if (require(rbenchmark)) { FOMC_SFO <- mkinmod( parent = mkinsub("FOMC", "m1"), m1 = mkinsub( "SFO")) b.2 <- benchmark( - "deSolve, not compiled" = mkinfit(FOMC_SFO, FOCUS_2006_D, + "deSolve, not compiled" = mkinfit(FOMC_SFO, FOCUS_D, use_compiled = FALSE, quiet = TRUE), - "deSolve, compiled" = mkinfit(FOMC_SFO, FOCUS_2006_D, quiet = TRUE), - replications = 3) + "deSolve, compiled" = mkinfit(FOMC_SFO, FOCUS_D, quiet = TRUE), + replications = 1, order = "relative", + columns = c("test", "replications", "relative", "elapsed")) print(b.2) factor_FOMC_SFO <- round(b.2["1", "relative"]) } else { @@ -1658,15 +1656,12 @@ SFO_SFO <- mkinmod( print("R package benchmark is not available") }</code></pre> <pre><code>## Successfully compiled differential equation model from auto-generated C code.</code></pre> -<pre><code>## test replications elapsed relative user.self sys.self -## 2 deSolve, compiled 3 4.879 1.000 4.877 0 -## 1 deSolve, not compiled 3 53.551 10.976 53.525 0 -## user.child sys.child -## 2 0 0 -## 1 0 0</code></pre> -<p>Here we get a performance benefit of a factor of 11 using the version of the differential equation model compiled from C code!</p> -<p>This vignette was built with mkin 0.9.49.9 on</p> -<pre><code>## R version 3.6.3 (2020-02-29) +<pre><code>## test replications relative elapsed +## 2 deSolve, compiled 1 1.00 0.459 +## 1 deSolve, not compiled 1 51.76 23.758</code></pre> +<p>Here we get a performance benefit of a factor of 52 using the version of the differential equation model compiled from C code!</p> +<p>This vignette was built with mkin 0.9.50.2 on</p> +<pre><code>## R version 4.0.0 (2020-04-24) ## Platform: x86_64-pc-linux-gnu (64-bit) ## Running under: Debian GNU/Linux 10 (buster)</code></pre> <pre><code>## CPU model: AMD Ryzen 7 1700 Eight-Core Processor</code></pre> |