Add analytical solution for DFOP-SFO

This is about twice as fast as deSolve compiled in the case of FOCUS D

1 files changed, 8 insertions, 1 deletions

diff --git a/man/create_deg_func.Rd b/man/create_deg_func.Rd
index 178661ff..69058038 100644
--- a/man/create_deg_func.Rd
+++ b/man/create_deg_func.Rd
@@ -30,6 +30,13 @@ if (require(rbenchmark))
   benchmark(
     analytical = mkinfit(SFO_SFO, FOCUS_D, solution_type = "analytical", quiet = TRUE),
     deSolve = mkinfit(SFO_SFO, FOCUS_D, solution_type = "deSolve", quiet = TRUE),
-    replications = 1)
+    replications = 2)
+  DFOP_SFO <- mkinmod(
+    parent = mkinsub("DFOP", "m1"),
+    m1 = mkinsub("SFO"))
+  benchmark(
+    analytical = mkinfit(DFOP_SFO, FOCUS_D, solution_type = "analytical", quiet = TRUE),
+    deSolve = mkinfit(DFOP_SFO, FOCUS_D, solution_type = "deSolve", quiet = TRUE),
+    replications = 2)
 }
 }