diff options
Diffstat (limited to 'inst/web/vignettes/compiled_models.html')
| -rw-r--r-- | inst/web/vignettes/compiled_models.html | 59 |
1 files changed, 33 insertions, 26 deletions
diff --git a/inst/web/vignettes/compiled_models.html b/inst/web/vignettes/compiled_models.html index 756a2753..50db7e9a 100644 --- a/inst/web/vignettes/compiled_models.html +++ b/inst/web/vignettes/compiled_models.html @@ -10,7 +10,7 @@ <meta name="author" content="Johannes Ranke" /> -<meta name="date" content="2016-03-23" /> +<meta name="date" content="2016-03-24" /> <title>Performance benefit by using compiled model definitions in mkin</title> @@ -65,7 +65,7 @@ img { <div id="header"> <h1 class="title">Performance benefit by using compiled model definitions in mkin</h1> <h4 class="author"><em>Johannes Ranke</em></h4> -<h4 class="date"><em>2016-03-23</em></h4> +<h4 class="date"><em>2016-03-24</em></h4> </div> <div id="TOC"> @@ -89,8 +89,10 @@ SFO_SFO <- mkinmod( <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 microbenchmark package.</p> <pre class="r"><code>library("microbenchmark") -library("ggplot2") -mb.1 <- microbenchmark( +library("ggplot2")</code></pre> +<pre><code>## Need help? Try the ggplot2 mailing list: +## http://groups.google.com/group/ggplot2.</code></pre> +<pre class="r"><code>mb.1 <- microbenchmark( "deSolve, not compiled" = mkinfit(SFO_SFO, FOCUS_2006_D, solution_type = "deSolve", use_compiled = FALSE, quiet = TRUE), @@ -98,27 +100,29 @@ mb.1 <- microbenchmark( solution_type = "eigen", quiet = TRUE), "deSolve, compiled" = mkinfit(SFO_SFO, FOCUS_2006_D, solution_type = "deSolve", quiet = TRUE), - times = 3, control = list(warmup = 1)) - -smb.1 <- summary(mb.1) + times = 3, control = list(warmup = 0))</code></pre> +<pre><code>## Warning in microbenchmark(`deSolve, not compiled` = mkinfit(SFO_SFO, +## FOCUS_2006_D, : Could not measure overhead. Your clock might lack +## precision.</code></pre> +<pre class="r"><code>smb.1 <- summary(mb.1) print(mb.1)</code></pre> <pre><code>## Unit: milliseconds ## expr min lq mean median uq -## deSolve, not compiled 9539.3064 9543.1547 9554.1987 9547.0031 9561.6448 -## Eigenvalue based 927.5569 928.1716 943.8293 928.7864 951.9656 -## deSolve, compiled 734.6125 737.3273 739.0161 740.0420 741.2179 +## deSolve, not compiled 9280.0854 9299.6757 9323.2559 9319.2659 9344.8411 +## Eigenvalue based 885.7475 891.8548 907.2823 897.9621 918.0498 +## deSolve, compiled 713.2624 721.4990 728.2856 729.7357 735.7972 ## max neval cld -## 9576.2865 3 c -## 975.1447 3 b -## 742.3938 3 a</code></pre> +## 9370.4163 3 c +## 938.1374 3 b +## 741.8588 3 a</code></pre> <pre class="r"><code>autoplot(mb.1)</code></pre> -<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAABUAAAAPACAMAAADDuCPrAAAC7lBMVEUAAAABAQECAgIDAwMEBAQFBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQWFhYXFxcZGRkaGhobGxscHBwdHR0eHh4fHx8gICAiIiIjIyMkJCQlJSUmJiYnJycoKCgpKSkqKiorKyssLCwtLS0uLi4vLy8wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6Ojo7Ozs8PDw9PT0+Pj4/Pz9BQUFCQkJDQ0NERERFRUVGRkZHR0dJSUlKSkpMTExNTU1OTk5PT09QUFBRUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVlaWlpbW1tcXFxdXV1eXl5fX19gYGBhYWFiYmJjY2NkZGRlZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1ubm5vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8fHx9fX1+fn5/f3+AgICBgYGCgoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uMjIyNjY2Ojo6Pj4+QkJCRkZGSkpKTk5OUlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2enp6fn5+goKChoaGioqKjo6OkpKSlpaWmpqanp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKzs7O0tLS1tbW2tra3t7e4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXGxsbHx8fIyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY2NjZ2dna2trb29vc3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr6+vs7Ozt7e3u7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+/v7///+MC2DoAAAACXBIWXMAAB2HAAAdhwGP5fFlAAAgAElEQVR4nO3dfaBU9X3n8UmNeTBJE/uwzbbJJrvtdttuuhnt2iZttulukv1drlxuUEAQ5ElRkRoisWokChoTH8iiofiANmJtrFZMtDYqaBSVGk0wRkTwCQVBEwXlQeAy578958zMvefMfYCZr99zPyPv1x+ZOb85v98542TeztyZey1FAICWlIb7BACgXRFQAGgRAQWAFhFQAGgRAQWAFhFQAGgRAQWAFhFQAGgRAQWAFhFQAGgRAQWAFhFQAGgRAQWAFhFQAGgRAT1w215VtjeKhvsUsB9bo9eH+xQwtF9G0d4hd2iIAgE9cAQUNgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfVDQGFDQOURUD8EFDYEVB4B9UNAYUNA5RFQPwQUNgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfVDQGFDQOURUD8EFDYEVB4B9UNAYUNA5RFQPwQUNgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfVDQGFDQOURUD8EFDYEVB4B9UNAYUNA5RFQPwQUNgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfVDQGFDQOURUD8EFDYEVB4B9UNAYUNA5RFQPwQUNgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfVDQGFDQOURUD8EFDYEVB4B9UNAYUNA5RFQPwQUNgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfVDQGFDQOURUD8EFDYEVB4B9UNAYUNA5RFQPwQUNgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfVDQGFDQOURUD8EFDYEVB4B9UNAYUNA5RFQPwQUNgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfVDQGFDQOURUD8EFDYEVB4B9UNAYUNA5RFQPwQUNgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfVDQGFDQOURUD8EFDYEVB4B9UNAYUNA5RFQPwQUNgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfVDQGFDQOURUD8EFDYEVB4B9UNAYUNA5RFQPwQUNgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfVDQGFDQOURUD8EFDYEVB4B9UNAYUNA5RFQPwQUNgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfVDQGFDQOURUD8EFDYEVB4B9UNAYUNA5RFQPwQUNgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfVDQGFDQOURUD8EFDYEVB4B9UNAYUNA5RFQPwQUNgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfVDQGFDQOURUD8EFDYEVB4B9UNAYUNA5RFQPwQUNgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfVDQGFDQOURUD8EFDYEVB4B9UNAYUNA5RFQPwQUNgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfVDQGFDQOURUD8EFDYEVB4B9UNAYUNA5RFQPwQUNgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfVDQGFDQOURUD8EFDYEVB4B9UNAYUNA5RFQPwQUNgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfVDQGFDQOURUD8EFDYEVB4B9UNAYUNA5RFQPwQUNgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfVDQGFDQOURUD8EFDYEVB4B9UNAYUNA5RFQPwQUNgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfXTdECfXf5yyw9k0wioPgIqj4D6aTagG79U/mrrj2SzCKg+AiqPgPppNqD3lMt/vqX1h7JJBFQfAZVHQP00G9Al5XL5idYfyiYRUH0EVB4B9dNsQL8VB/S+1h/KJhFQfQRUHgH102xAz4wDemvrD2WTCKg+AiqPgPppNqAnxwG9ofWHskkEVB8BlUdA/TQb0OPigC5t/aFsEgHVR0DlEVA/zQb0aAKKHAIqj4D6aTagf0NAkUNA5RFQP80G9CgCihwCKo+A+mkyoC+VCShyCKg8AuqnyYCuIaDII6DyCKifJgP6CAFFHgGVR0D9NBnQ5QQUeQRUHgH102RAbyWgyCOg8gion+YCuuUqAoo8AiqPgPppKqAbuspHEFDkEFB5BNRPUwH913KZV6DII6DyCKifpgL6AwKKRgRUHgH1Q0BhQ0DlEVA/BBQ2BFQeAfXTbEA/P5WAIouAyiOgfpoN6ILrCCiyCKg8AuqHgMKGgMojoH4IKGwIqDwC6oeAwoaAyiOgfggobAioPALqh4DChoDKI6B+CChsCKg8AuqHgMKGgMojoH4IKGwIqDwC6oeAwoaAyiOgfggobAioPALqh4DChoDKEw3ophD6D6696JTu42Yv2nKAu/sY/FD9biGgsCGg8tonoJXLQxhxwpgQOm88kN0dhLCJgKJABFRe+wT0h6H7zj1R9NqVITx4ALs7IKAoFgGV1z4BnRZWVK9cHs4/gN0dLFu2nYCiQARUXtsEdEcI26vXVoeJ+9/dDQFFYQioPLmArp53QtesuzdWY7Th4pO7Js57LLm6PYS11T32ba+WdM2F00edevFTydWkXVeEq9LhDWH0W/mpjZ791kldJy1+fYBFeq6f2T3zyh2vXTK9a+r1u6PkXfvGu07rOmH+6ijKv4XPLZ875z4EFDYEVJ5aQG/oCB0TO8KFaYx+2Bm6p3eHcHVyy8ww4Z5dmT1v7Aidk+O974yqVXsiTKwk49eHyxun5t3RGUZO6QzHvdZ/kblzblsyMswad0ZysShKmnleOOeauSM6/jHKBTS3fO6cE09WbX69CXdUA/qPzcwx6Ymiwo6F1rwZvTHcp4ChbY0DOuQODfVxDuiaEK7ZFW2/NCQxer5zzMOVqLJqbHggvmn9cSF0nX3z+kp1z7Wh45/2RLu+G0b9qlq1ysT0JWplSnKRn5qzobPj1p5o29zwzf6LzNsXRbeFcH58sSyMj5JmdiQLPD4qrMsGNLd87pwT+6r/fbjyimbu+X3VgC4z/eMDIM05oGclWYsrOCuJ0SW1j9tXhq8nF1tvOb0jjtTYa5KXjtHZ4bJ0zzOSV4pp1a5NXw+uDTMq/aZmzQ8Lk4s3Oo/Z12+RZ+KNZ3svoqSZF6aTrg1zswHNLZ875wQBBTAg34BWjg1PplfuTmI0M0yZnpgcxtVu3/7QokkhHP9ifHVs7UeiK8LsWtWeCyfE6Vwcbo0GmtprfG3i8+v29ltkb5S9iJJmVl/BPhOOywY0u3z+nNN7cUbV6reacHc1oDc3M8ck/mdV2LHQmj3RnuE+BexH/IJpP7fn+AZ0awhvpFeeSmJ0TKgb2bdL5dFJ4az4BWQIb9b2HFerWmVGeDraN75zazTI1NSO+iGiARdpuAhhfX3WzkxAs8vnzzmLD5Fgw4dI8rQ+ROqt29NJjGaEzNHnzX2pdm1N6Hirr31Ph9H12t0UlkQ/C/OT0dzUnLh3O+rXB1qkMaDVl6ivh7AtE9Ds8vlzziKgsCGg8rQCGo2rvR1ekcTovPDTdOOtLa9H0aTwo9o+29Nijaml7d5wer12m+P38AvCI8lobmpO/I77mfTKuke3D7RIY0CXpzusDmOyb+Fzy+fOOYuAwoaAyhML6LnVD20qs5MY3Vb7laMrwtIoWhgm1d563x8mR8lHNwvSPeeEK3qjNzusGT2hJ7mWm5p3VlicXOwd27FjwEUaAjo7/dR/bjg7G9Dc8rlzziKgsCGg8sQCGr8NvnpXtHNh+pWg3ZPCkt1Rzx0dXZvjl5fd4YQ7X+mpbL19dLglSn7i2HHT3uQbSF2v9Ebv9nBKuC5dJzc1Wrasp+8QT4XOH/REuxaEcwZepCGgYeHOaOeiEJ7KBjS3fO6cswgobAioPLGARktC6JjUEZamMXpsZOicekzoWJnc8ospcaFGjIz/Z0H6ovDaEI6e0hG67umL3tYRIWysrpObGsLuzCFujCdN6wpjXxl4kYaAzk/Pp+OGKPdF+tzyuXPOIKCwIaDy1AJaeXDe5K5Zq3p/YXL80VO//WL1pp7lc2d0H3vaxetqu666YHrXyRenHy3Vfw393DCnvlB2aj6g0ZPzp3bNWLxtkEUaArrxnpkjJ53/eHWZ7K9y9i2fP+c+BBQ2BFSeWkBd7E2/2dmKtJmtIqCwIaDyDoqAvjyq1ZkEFMOIgMo7KAI6f0GrMwkohhEBlXdQBPS+Vt/BE1AMJwIq76AIaOsIKIYRAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP00G9C//XsCiiwCKo+A+mk2oOWjCCiyCKg8Auqn6YCWCSiyCKg8AuqHgMKGgMojoH4IKGwIqDwC6qepgD5IQNGIgMojoH6aCuirN00goMgjoPIIqJ/mAvrqPxFQ5BFQeQTUT5MB/VcCijwCKo+A+mkyoA8QUOQRUHkE1E+TAV1NQJFHQOURUD9NBvQZAoo8AiqPgPppMqBbjiCgyCGg8gionyYD+upfEVDkEFB5BNRPswH9EgFFDgGVR0D9NBvQ0QQUOQRUHgH102xApxBQ5BBQeQTUT7MB/ds4oN9v/aFsEgHVR0DlEVA/zQb0G3FA72z9oWwSAdVHQOURUD/NBvTyOKCrWn8om0RA9RFQeQTUT7MB/ec4oM+1/lA2iYDqI6DyCKifZgO6+ohyaP2RbBYB1UdA5RFQP80G9NXvHvejlh/IphFQfQRUHgH103RAC0VA9RFQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1A8BhQ0BlUdA/RBQ2BBQeQTUDwGFDQGVR0D9EFDYEFB5BNQPAYUNAZVHQP0QUNgQUHkE1M/wB3TL2d2n/miQ2wioPgIqj4D6Gf6A3l+Ozd084G0EVB8BlUdA/Qx/QO9JAlqetWWg2wioPgIqTzOgoc933q41N4VwQEfe9LYdQCKg4cj4NehAtxFQfQRUnmpAjxlTszgyVa3PQRrQ+bfHBf37AW4joPoIqDzVgG4aYrM1B2tAo38ol4+8qTbw0AXHH3dH9SoB1UdA5bVFQJct225f86ANaHRBufxnNySbd45PfiQ6unobAdVHQOW1RUDfFgdvQPfNjrM57ZrvdKcfKZWPrt5GQPURUHltEdDa5vqLpnZ95ZFapzZcfHLXxHmPJVfjkZ6lx3dNv2FXFF0RrkqnbAij34qinnvmjD96wpy79kW1vtUjV1sxs0jfoe6b1TXxG4+kW9n50eZFJ3cfe/Li6pnlZq6ed0LXrLs3qgY02nNmueZL1xHQNkJA5bVRQO/oDJ2TO8LStFM/7Azd07tDuDpKq7hw5i1LR4eLouiJMLGSTLk+XB4lOQ3jkr0WRQMGNLtI76EuDx0TO0K1w9n567vD0dMmhNC9rnHmDR3plAtlAxpVvv/XST5H3rynQkDbCAGVpxrQydNr1kS13L3U2XHjnmhX3LS4U893jnm4ElVWjQ0PpFU8LX7x+VDo2BlVJoa1STGmJBcvhaMfqUR7bwxH7x4ooLlFeo8cFuyIdiwK4fGG+XPCxTviQHeFMxtmrgnhml3R9ktDPaCVE6se3TPcHq4FNIp23fsP//xk/K+WOKBd1dvijeE9OezX3mjvcJ8C9iN+Tu3n9pzCvwf6s3ruLgyXpXX6atKpS8KD6Y4rw9fTKiYF3JXudm36qnBtmBEH4rHT01eIu9Mb+gc0t0jvkb+SvoKdF77WML87PJ1s3LZoScPMs8I30xObVQ/ovtr75RVO/3QO3E96A1oXB3TU8JwLgGF8Cz8hfWkZRcuTTs0MU9LXp5PDuLSKr/bu9lw4IS7g4nBr7+yeVYMENLdI76FWpJdPhGMq+fmnhTk/7akNZWdWjg1PpoN3Swd0z8ZfJhcEFBhGwxbQHSG8mW48lXTqmN5XqCPTKu7rnVWZEb9S3De+c2syUnnqzsVnfjkMEtDcIr2HSl9nRm+EsDU//9lpIYyZe+O6pKvZmVtDeKPvxFLbql775XBbXg/ohjP/vFz+0ne2JAE9unpb/C+D4T057Ne2aOtwnwKG9qso2jvkDg1pG7aAbqsHdF3SqRkh88PZ/CdDN4Ul0c9C2o1XTg0dJ11y25p+Ad2TjuQW6T1U9XVunMWd+flRz78vPnVECKe/nJ+5ox7Qp3U/RPq3z1RfEx8173k+RGojfIgkT/VDpP5v4cfW0vbjpFPnhZ+mG29teb0xoJvj9/ALQvo1pK+H8+J/QUT7+gV0YzqSW6T3UMvTy8fD8Q3zU7senB3mNMwcV3sLv0I2oLccUS5/duroI/keaJshoPLaJ6DnhgXpxllJp24L56cbV4Sl/b7cOTusGT0h/WnlqJB85Sj6RS6g2+KRf0lHcov0Hur09Gef3wgXNsw/7eT0A7Wnw5cbZp6b7Bq/3Z+tGtB74nKe8syrrz5xzlHpl5mqtxFQfQRUXvsEdH1Hx/f3RruvHpF0aveksGR31HNHR9fmfgG9PZwSrksHpqcfJa05IX1fnu61uzNcV4nWj093zS3Se6ja15g6X2qYPzvZOdr2zeQVaG5m/Nb96l3RzoVBNKA//4ty+YzqX7Rbc/H/LZe/Vb2NgOojoPJUA9r715jG7K6X8eaO0Dmlc8yK0B1vPDYydE49JnSsjPoFdGuc2I3pwB0hnPi1yWHeiWHsrbW9rghh+syOG6q7ZhfpPfJV6bfiO+9onL+6I3RPnzIijF7XOHNJCB2T6t/wz9AI6At/Uy6f2Pc3lZ9dU7tCQPURUHmqAe3TG9DoifMndp+z4edhWrKx4eLxR0/99ovJ1YaAxu+p51SXqfx49rHj5i6vrD+5+8baXnv++aSuiTdXen+Vs3eR3iNveuC0rinfXNdvfvTU/Mkjx5x63ZaocWblwXmTu2at0vxd+GlfLJe/vGGA2wioPgIqTzOgQ7mnnkd5EgFNfgF+7UC3EVB9BFRe+wT0qtE3JBeVM8LCYTuH5gx/QJenn7o/MeBtBFQfAZXXPgFdFbofrkQ7F4fw82E7h+YMf0Cf+2L5M3NfGPg2AqqPgMprn4BWFoQwOvk7SVcM2yk0afgD+uqmn28a7CYCqo+AymufgEaVB86c0DX93OWV/e+qQSCgQyCg+giovDYKaNshoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKueO+U/ltgmoHwIKGwKq5tmjyqfnBgioHwIKGwKq5tFyeVxugID6IaCwIaBqVpXLx+YGCKgfAgobAqpmZbncnRsgoCrBaUgAABlTSURBVH4IKGwIqJofl8sjcwME1A8BhQ0BVXNvuTwiN0BA/RBQ2BBQNcvL5ZAbIKB+CChsCKiau8rlL+UGCKgfAgobAqrmRwS0OAQUNgRUTRzQL+QGCKgfAgobAqomDuj/yQ0QUD8EFDYEVE0c0M/nBgioHwIKGwKqJg7oX+cGCKgfAgobAqomDuj/yg0QUD8EFDYEVE0c0M/lBgioHwIKGwKqJg7oX+UGCKgfAgobAqomDuhncwME1A8BhQ0BVRMH9DO5AQLqh4DChoCqiQP657kBAuqHgMKGgKqJA3pUboCA+iGgsCGgauKA/s/cAAH1Q0BhQ0DVxAE9MjdAQP0QUNgQUDVxQI/IDRBQPwQUNgRUTRzQcm6AgPohoLAhoGqSgG7JDhBQPwQUNgRUTRLQzdkBAuqHgMKGgKpJAropO0BA/RBQ2BBQNQS0QAQUNgRUTRLQjdkBAuqHgMKGgKpJAvpidoCA+iGgsCGgapKAbsgOEFA/BBQ2BFRNEtDnswME1A8BhQ0BVZME9LnsAAH1Q0BhQ0DVJAF9NjtAQP0QUNgQUDVJQJ/JDhBQPwQUNgRUTRLQ9dkBAuqHgMKGgKpJArouO0BA/RBQ2BBQNUlA12YHCKgfAgobAqomCehT2QEC6oeAwoaAqkkCuiY7QED9EFDYEFA1SUCfzA4QUD8EFDYEVE0S0CeyAwTUDwGFDQFVQ0ALREBhQ0DVJAF9PDtAQP0QUNgQUDVJQH+WHSCgfggobAioGgJaIAIKGwKqJgnoY9kBAuqHgMKGgKpJAvqT7AAB9UNAYUNA1SQBfSQ7QED9EFDYEFA1BLRABBQ2BFRNEtBV2QEC6oeAwoaAqkkC+nB2gID6IaCwIaBqkoA+mB0goH4IKGwIqJokoCuzAwTUDwGFDQFVkwT0B9kBAuqHgMKGgKpJAvrF7AAB9UNAYUNA1SQBPTI7QED9EFDYEFA1BLRABBQ2BFQNAS0QAYUNAVVDQAtEQGFDQNUQ0AIRUNgQUDUEtEAEFDYEVA0BLRABhQ0BVUNAC0RAYUNA1RDQAhFQ2BBQNQS0QAQUNgRUDQEtEAGFDQFVQ0ALREBhQ0DVENACEVDYEFA1BLRABBQ2BFQNAS0QAYUNAVVDQAtEQGFDQNUQ0AIRUNgQUDUEtEAEFDYEVA0BLRABhQ0BVUNAC0RAYUNA1RDQAhFQ2BBQNQS0QAQUNgRUDQEtEAGFDQFVQ0ALREBhQ0DVENACEVDYEFA1BLRABBQ2BFQNAS0QAYUNAVVDQAtEQGFDQNUQ0AIRUNgQUDUEtEAEFDYEVA0BLRABhQ0BVUNAC0RAYUNA1RDQAhFQ2BBQNQS0QAQUNgRUDQEtEAGFDQFVQ0ALREBhQ0DVENACEVDYEFA1BLRABBQ2BFQNAS0QAYUNAVVDQAtEQGFDQNUQ0AIRUNgQUDUEtEAEFDYEVA0BLRABhQ0BVUNAC0RAYUNA1RDQAhFQ2BBQNQS0QAQUNgRUDQEtEAGFDQFVQ0ALREBhQ0DVENACEVDYEFA1BLRABBQ2BFQNAS0QAYUNAVVDQAtEQGFDQNUQ0AIRUNgQUDUEtEAEFDYEVA0BLRABhQ0BVUNAC0RAYUNA1RDQAhFQ2BBQNcMd0E0h9B9ce9Ep3cfNXrTlAHcfXoOfUr9bCChsCKgawYBWLg9hxAljQui88UB2b0UIm962ZQgoCkNA1QgG9Ieh+849UfTalSE8eAC7t4KADoCA6iOgagQDOi2sqF65PJx/ALu34m0K6LJl2wkoCkRA1egFdEcI26vXVoeJ+9+9JW9TQFMEFIUhoGqGL6Cr553QNevujdXIbLj45K6J8x5Lrm4PYW11j33bqyVdc+H0Uade/FS9SVeEq9LhDWH0W/mp/frVs/T4ruk37GpY5eKQyHxE9ey3Tuo6afHrAxys5/qZ3TOv3PHaJdO7pl6/O0rSu/Gu07pOmL86ivJv4XOnkbtvfQgobAiommEL6A0doWNiR7gwjcwPO0P39O4Qrk5umRkm3LMrs+eNHaFzcrz3ndWmRU+EiZVk/PpweePUnHjfhTNvWTo6XNSwyl2Xh3DR5W/07nhHZxg5pTMc91r/g82dc9uSkWHWuDOSi0VpM88L51wzd0THP+YDmjuN3H2LVe6qeuENZfuiaLhPAfuxI9o+3KeAnBVJQLMDb0ZRz5Az3p6Argnhml3R9ktDEpnnO8c8XIkqq8aGB+Kb1h8XQtfZN6+vVPdcGzr+aU+067th1K+qtapMTF+iVqYkF/mpjQE9LQ7xQ6FjZ8Mq+bfwGzo7bu2Jts0N3+x/sHlxVW4L4fz4YlkYn87sSA70+KiwLhvQ3Gnk7ltiX7lqRWv/rACIWpUE1LJAiwE9K8lVXMFZSWQuqX3cvjJ8PbnYesvpHXF8xl6TvCSMzg6XpXuekbwCTGt1bfo6b22YUek3NSveN2ndrjRzuVXyAZ0fFiYXb3Qes6/fwZ6JN57tvUhnXphOujbMzQY0dxq5+5YgoMA70zAFtHJseDK9cncSmZlhyvTE5DCudvv2hxZNCuH4F+OrY2s/El0RZtdq9Vw4IU7n4nBrNNDUXvG+6c8b0szlVskHdHztpufX7e13sL1R9iKdWX2l+0w4LhvQ7Gnk71t6b6+tWrddWfwqe7hPAfuxK9o53KeAnPuSgGYHdqQf3gzhbQno1hCqPwt4KonMMaFuZN8ulUcnhbPiF4YhvFnbc1ytVpUZ4elo3/jOrdEgU3sDuq8e0PwquYDuqJ9KNODBGi5CWF+ftTMT0Oxp5O9bFh8iwYYPkdQM04dIvdV6OonMjJBZdd7cl2rX1oSOt/qa9nQYXa/YTWFJ9LMwPxnNTc2rf4soF9B0lVxA497tqF8f6GCNAa2+RH09hG2ZgGZPI3/fsggobAiomuH6FH5c7W3uiiQy54WfphtvbXk9iiaFH9X22Z6WaEwtWfeG0+sV2xy/h18QHklGc1PzsgHNr5ILaPyO+5n0yrpHtw90sMaALk93WB3GZN/C504jd9+yCChsCKia4QroudUPYyqzk8jcVvuVoyvC0ihaGCbV3lLfHyZHyUcyC9I954QremM2O6wZPaEnuZabmpcLaG6V/M9AzwqLk4u9Yzt2DHiwhoDOTr8dMDecnQ1o7jRy9y2LgMKGgKoZroDGb2+v3hXtXJh+1Wf3pLBkd9RzR0fX5vjlZXc44c5Xeipbbx8dbomSnyR23LQ3+WZR1yu9Mbs9nBKuS9fJTY2WLesZJKC5VeKhp/r2eyp0/qAn2rUgnDPwwRoCGhbujHYuSlfoC2juNHL3jYDi7UNA1QzbF+mXhNAxqSMsTSPz2MjQOfWY0LEyueUXU+LyjBgZ/8+C9MXetSEcPaUjdN3TF7OtI0LYWF0nNzWE3YMENLdK/Cb7uDNf693xxnh4WlcY+8rAB2sI6Pz0vDtuiHJfpM+dRu6+EVC8bQiommELaOXBeZO7Zq3q/UXI8UdP/faL1Zt6ls+d0X3saRevq+266oLpXSdf/FI2iueGOfWFslOHCGh2lejh6UeP/2Xfnk/On9o1Y/G2QQ7WENCN98wcOen8x3tXHuAe5O9bHwIKGwKqZrj/mMjba2/6jU1Ppj9DQkBhQ0DVvLMC+vIo7yMQUAwjAqrmnRXQ+Qu8j0BAMYwIqJp3VkDv834HT0AxnAiomndWQP0RUAwjAqqGgBaIgMKGgKohoAUioLAhoGoIaIEIKGwIqBoCWiACChsCqoaAFoiAwoaAqiGgBSKgsCGgaghogQgobAioGgJaIAIKGwKqhoAWiIDChoCqIaAFIqCwIaBqCGiBCChsCKgaAlogAgobAqqGgBaIgMKGgKohoAUioLAhoGoIaIEIKGwIqBoCWiACChsCqoaAFoiAwoaAqiGgBSKgsCGgaghogQgobAioGgJaIAIKGwKqhoAWiIDChoCqIaAFIqCwIaBqCGiBCChsCKgaAlogAgobAqqGgBaIgMKGgKohoAUioLAhoGoIaIEIKGwIqBoCWiACChsCqoaAFoiAwoaAqiGgBSKgsCGgaghogQgobAioGgJaIAIKGwKqhoAWiIDChoCqIaAFIqCwIaBqCGiBCChsCKgaAlogAgobAqqGgBaIgMKGgKohoAUioLAhoGoIaIEIKGwIqBoCWiACChsCqoaAFoiAwoaAqiGgBSKgsCGgaghogQgobAioGgJaIAIKGwKqhoAWiIDChoCqIaAFIqCwIaBqCGiBCChsCKgaAlogAgobAqqGgBaIgMKGgKohoAUioLAhoGoIaIEIKGwIqBoCWiACChsCqoaAFoiAwoaAqkkC+vnsAAH1Q0BhQ0DVJAG9MTtAQP0QUNgQUDVJQFdmBwioHwIKGwKqJgnog9kBAuqHgMKGgKpJAvpwdoCA+iGgsCGgapKArsoOEFA/BBQ2BFRNEtB/zw4QUD8EFDYEVE0S0EeyAwTUDwGFDQFVkwT0J9kBAuqHgMKGgKpJAvpYdoCA+iGgsCGgapKA/jQ7QED9EFDYEFA1SUB/lh0goH4IKGwIqJokoI9nBwioHwIKGwKqJgnoE9kBAuqHgMKGgKohoAUioLAhoGqSgP4iO0BA/RBQ2BBQNUlA12QHCKgfAgobAqqGgBaIgMKGgKpJAro2O0BA/RBQ2BBQNUlA12UHCKgfAgobAqqGgBaIgMKGgKpJAvpMdoCA+iGgsCGgapKAPpsdIKB+CChsCKiaJKDPZQcIqB8CChsCqiYJ6PPZAQLqh4DChoCqSQK6ITtAQP0QUNgQUDVJQF/MDhBQPwQUNgRUTRLQl7IDBNQPAYUNAVWTBHRTdoCA+iGgsCGgaghogQgobAiomiSgm7MDBNQPAYUNAVWTBHRLdoCA+iGgsCGgapKAvpIdIKB+CChsCKiaOKBH5AYIqB8CChsCqiYO6JG5AQLqh4DChoCqiQP6Z7kBAuqHgMKGgKqJA3pUboCA+iGgsCGgaghogQgobAiomjigf5EbIKB+CChsCKiaOKCfzQ0QUD8EFDYEVE0c0L/MDRBQPwQUNgRUTRzQz+UGCKgfAgobAqomDuj/yg0QUD8EFDYEVE0c0L/ODRBQPwQUNgRUTRzQz+cGCKgfAgobAqomDuj/zg0QUD8EFDYEVE0c0C/kBgioHwIKGwKqJg7oF3MDBNQPAYUNAVVzV7n8pdwAAfVDQGFDQNUsL5dDboCA+iGgsCGgau4tl0fkBgioHwIKGwKq5r5yeWRugID6IaCwIaBqVpbLo3IDBNQPAYUNAVWzqlw+JjdAQP0QUNgQUDU/KZePyw0QUD8EFDYEVM36I8un5AYIqB8CChsCKucfvvpobpuA+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D6IaCwIaDyCKgfAgobAiqPgPohoLAhoPIIqB8CChsCKo+A+iGgsCGg8gioHwIKGwIqj4D60Q7oT1aufGW4zwFDe23Pa8N9Chja5pUrHxtyh4YoENB3itHl8t7hPgegzW0rlyc2sz8BfacgoIAZAT1YEVDAjIAerAgoYEZAD1YEFDAjoAcrAgqYEdCDFQEFzAgoABSEgAJAiwgoALSIgAJAiwgoALSIgAJAiwgoALSIgLanfwmbatc2Xzahc8KlL+9vA0Cfpp9AgzybCGhb2n1i/fFf/+UQJoTQvW7oDQB9mn4CDfZsIqBtqGfN10Pt8e+ZHs5+PXrtrDCtZ6gNAL2afwIN+mwioO3nus4Q6o//yjBuV3yxc2x4cKgNAHUtPIEGfTYR0PZz93e+8536439ZuDy9XBgWDLUBoK6FJ9CgzyYC2p7qj/+ccH96+eNw5lAbALKafAIN+mwioO2p/vhPC4+nl6vDiUNtAMhq8gk06LOJgLan+uPfHZ5LL58Lo4faAJDV5BNo0GcTAW1PfY//8+nls6FrqA0AWU0+gQZ9NhHQ9tT/HcjUoTYAZDX5BBr02URA21Pfz8AfSC/vD18bagNAVpNPoEGfTQS0PdUf/0vDovRyUbhsqA0AWU0+gQZ9NhHQ9lR//B8IE/bEF3snJF/wHXwDQFaTT6BBn00EtD31/ibatPCtvdHei8L0nqE2AGQ1+QQa9NlEQNtT/fGP1nWH7tNHhdHPDL0BIKPZJ9BgzyYC2p56H//o5UvHd46/bMv+NgD0afoJNMiziYACQIsIKAC0iIACQIsIKAC0iIACQIsIKAC0iIACQIsIKAC0iIACQIsIKAC0iIACQIsIKAC0iICivXy+lPfx1pdK55cPYLc/yFw4G+goH0jOs4iDo2kEFO1loIC2mJf9BTRfzmI6SkDbCwFFe7ntu6nfLpWqV66PDAH9xIsvbh7qdo2AvvTiiwRUFAFFW/qDUub/ui0HdOhpKgH1PypaRUDRlnIBffPN7a2ssb8q1ZbNl7PFYxnPiYCKIqBoS7mAtugAqzTsHyIVdnA0jYCiLeUCWttIKnPz5z78no8dvy56ecrvHPrxjodrO6w9+Q8/8Ot/cuJP82vUqzTYvMyyfRf1A2/4yl989L2f/Jtr9wxyiKdP+cMPfviIU+v/4YiHjv/E+w/75MSH+468dtLvvec/fvaKHdWRu7780UN/Z8zq+lEaVyegoggo2tLAAf392aXf/eyHS6XDvvfRd/2P/1IqveveZLxyySG1j+z/dm92jb6ADjxvqIBe8e7akn+0daBDVL5d2z5sWbK596T6lwZO3FNb698+8O5Pf/rQUulT2+LtPcdXb33vsupRGlYnoLIIKNrSwAF91/uXRtHrn4vL85urospXSqXPJONXlkofOeOG7512WKk0LbtGX0AHnjdEQH/0rtL7Tlp626X/tVQaM9AhLi+VDpm65MaZ7yt9OPnvOJ4eb0679rrpcRdPr651+K+Pisu5/o9Kpa/E26fEpzD+ykWjSh9Mj9K4OgGVRUDRlgYOaOmi5PLh+MrV8eWbh5Q+FF+89sHSp15Nxl/8z6XSvZk1+gI64LyhAhq/Yrw12Xr9P5Q+tK//IV5+X+kD6ZFWHpoE8slS6QP3JZv3x4FcUz3iX1aSgXtKpT+LosdLpffckWxeWf3CZ8PqEQGVRUDRlgYJ6JvJ5Y5S6dd2JVc+lo4vKJXWVvd7KP8SNBPQgeYNFdByqfTLdO5NF1zwRv9DzC2VvlrdHl06IopOK5W+Vt08o/oSND7i8qh2yHjVk0ulU6u3fyHdblg9IqCyCCja0sAB/Vh1u/4LntXxEaU/re237/DSf8us0RfQAecNFdBjS6WRP+9bqfEQny2Vnqtuv7Dix1H0p6XSC9XN50pJUJMfGuzKnMOflErrq5u3pdsNq0cEVBYBRVsa7FP4KHulOv6p7G9+/nZmjeyn8APMGyqgT384XutTs7+/oTqx8RC/VXrPvsyBPlJ6b21z36Gl30zX+t3sOfx67+7r0u2G1SMCKouAoi01EdBPZuv23swahoBGm076jXS9P/7GmwMc4pD6a9qqQ/r+5MnHSodkj1i91rf7zuot+dUjAiqLgKItNRHQT5e+OPAaloBGUc+/f7Pro3Hkfu/l/of4cOmwSmYz8wr0vaWPRP0CenjpPT3VzWfqt2RXjwioLAKKttREQL9c+r2B17AFNFFZOaJUmtz/EP+9VKp9g37D9763PfMz0Bfit+ZRv4AeUSqtq27+a7aU9dUjAiqLgKItNRHQK0ql+6rjD3384/8vs4YhoGP/aGJ1yq/SP4jXeIjTSqVzq9tfLR2+L5pZKv1ddXNO9fP2hoCe3vsp/Ih0u2H1iIDKIqBoS00EdOsHS7//YrL52qdLv/ZCZg1DQCeWDn0y3fX7pdIX+h9izbtKH3oo2V73kdLIKHqiVPrgymTzvsNKpV9E/QL6ZP17oDdUvwfasHpEQGURULSlJgIaLY77NXvpLed9rFQ6O7uGIaDfK5V+6+9uuv3aCe8ulW4e4BBnlkrvPvHam8/5jdJhSQpPSzavu27aIdVfPGoMaDQr+U2kq64cU/r9dLtxdQIqi4CiLTUT0MoF76p+Pv5rs7Mf7VgC2nNs/TP3d18w0CF6vlK7+UPfTzb3TK3vPmNv7oi1a3smVm/9zcfT7cbVCagsAoq21ExAo+jnUz75vsOPnPZ0fg3Th0j3d//x4e/7xF+d/3I08CEenfzJ9x9+1Mz6H7x/YPx/et/7P3F85q8x5a/d1f3RQ3973DP17cbVCagoAoqDVhtVqY1O9eBCQHHQaqMqtdGpHlwIKA5abVSlNjrVgwsBxUGrjarURqd6cCGgOGiln3MP8d+FV8F/F14XAcVBi4DCioACQIsIKAC0iIACQIsIKAC0iIACQIsIKAC0iIACQIsIKAC0iIACQIsIKAC0iIACQIsIKAC0iIACQIsIKAC06P8DWD3v3ImGWlEAAAAASUVORK5CYII=" title alt width="672" /></p> -<p>We see that using the compiled model is by a factor of 12.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> +<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAABUAAAAPACAMAAADDuCPrAAAC7lBMVEUAAAABAQECAgIDAwMEBAQFBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQWFhYXFxcZGRkaGhobGxscHBwdHR0eHh4fHx8gICAiIiIjIyMkJCQlJSUmJiYnJycoKCgpKSkqKiorKyssLCwtLS0uLi4vLy8wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6Ojo7Ozs8PDw9PT0+Pj4/Pz9AQEBBQUFCQkJDQ0NERERFRUVGRkZHR0dISEhJSUlMTExNTU1OTk5PT09QUFBRUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVlaWlpbW1tcXFxdXV1eXl5fX19gYGBhYWFiYmJjY2NkZGRlZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1ubm5vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8fHx9fX1+fn5/f3+AgICBgYGCgoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uMjIyNjY2Ojo6Pj4+QkJCRkZGSkpKTk5OUlJSVlZWWlpaXl5eYmJiampqbm5ucnJydnZ2enp6fn5+goKChoaGioqKjo6OkpKSlpaWmpqanp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKzs7O0tLS1tbW2tra3t7e4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXGxsbHx8fIyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY2NjZ2dna2trb29vc3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr6+vs7Ozt7e3u7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+/v7///+2USVOAAAACXBIWXMAAB2HAAAdhwGP5fFlAAAgAElEQVR4nO3de5zV9X3n8ZMaczFJE3vZZttkk2673babbn6aza3bJL1kt99hYJyggCDITVGBKhGNGomCt6CSAAWj4EasifGCUWu8EaN4qVGDMSKCN1QEkygqF4Hh/P7b8/udMzPnd+bHMHPmM795f+T1fDzamTO/8/udQ3nPq2eYAUsxAKAppaF+AgDgFQEFgCYRUABoEgEFgCYRUABoEgEFgCYRUABoEgEFgCYRUABoEgEFgCYRUABoEgEFgCYRUABoEgEFgCYR0L5749c+bIvjt4b6OaBXv3UzpgPV63H8du6BhigQ0L7zsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2mt283ed9vUbmv+d7D8CKo+AqiOg9prc/LcPi6Lo7AH8XvYXAZVHQNURUHvNbX5plLpqIL+b/UNA5RFQdQTUXlObv/9/RdHC5VH0lY0D+v3sDwIqj4CqI6D2mtn8prYoOieOp0XRkoH9hvYDAZVHQNURUHvNbP6iKGrdEcdromjUwH5D+4GAyiOg6giovSY2/++HR4c/Ujl1bxSNGODvaN8RUHkEVB0BtdfPzd+34LvnfS6KLklOJaCoR0DVEVB7/dv8OdXvvk/Zk5xKQFGPgKojoPb6tfl70nwe9o0d6akEFPUIqDoCaq9fm/9xFM24484XaqcSUNQjoOoIqL3+BnRB96kEFPUIqDoCao+AwggBVUdA7RFQGCGg6gioPQIKIwRUHQG1R0BhhICqI6D2CCiMEFB1BNQeAYURAqqOgNojoDBCQNURUHsEFEYIqDoCao+AwggBVUdA7RFQGCGg6gioPQIKIwRUHQG1R0BhhICqEwvophB6fnDdhSe2Hz1ryZY+3n1w7PuhehwhoDBCQNXpB7S8KIRhx44KofWavtx9EISwiYBiKBBQdfoBvTm037Y7jl+7LIT7+3D3QUBAMUQIqDr9gE4Jq6rvLArn9uHug2Dlym0EFEOBgKqTD+j2ELZV31sTxu//7oOGgKJ4BFSdTEDXzD22beadL1djtHH+CW3j5z6avLsthHW1uGyrlnTt+VOPOGn+U8m7SbsWh++lH94YRr6dPbXRsxcd33b80tdzLtJx1fT26Zdtf+3iqW2Tr9oVJ1+1v3zHjLZj562J4+yX8JnLZ55zNwIKIwRUnUpAr24JLeNbwvlpjG5uDe1T20O4PDkyPYy7a2fdPa9pCa0TK/e+La5W7Ykwvpx8/KqwqPHUrFtbw4hJreHo13peZM7sm5aNCDPHnJq8WRInzTwnnHXFnGEt/xZnApq5fOY5J56s2vx6P9zaI6D9OXtAdsTxjsIeDM14I9421E8BvXozjnflHmiozyAHdG0IV+yMt10Skhg93zrqwXJcfmh0uK9yaMPRIbSded2GcvWe60LLD3bHO/81HPHbatXK49OXqOVJyZvsqRkbW1tu7IjfmBMu6HmRuXvj+KYQzq28WRnGxkkzW5ILPH5EWF8f0MzlM885sbf6HyiOVvXnV35PY0Dbm/m/HwBlgxzQM5KsVSo4M4nRxbVvt68O30zebL3+5JZKpEZfkbx0jM8Ml6b3PDV5pZhWbXn6enBdmFbucWq9eWFh8ubN1iP39rjIM5Ubz3a9iZNmnp+etDzMqQ9o5vKZ55wgoAByDW5Ay0eFJ9N37kxiND1MmpqYGMbUjm97YMmEEI55sfLu6Nofia4Ks2pVey4cW0nn0nBjnHdql7G1E59fv6fHRfbE9W/ipJnVV7DPhKPrA1p/+exzTn8Vp1atebsf7mwM6BH9OXtAKr/cPYU9GJqxK9491E8BvdoVxx25BxryM7gB3RrCm+k7TyUxOjJ0GtF9l/IjE8IZlReQIbxVu+eYWtXK08LT8d6xrVvjfZya2t75EHHuRRrehLCh86wddQGtv3z2Odfjm0gwwjeR1Gl8E6mrbk8nMZoW6h597pyXau+tDS1vd7fv6TCys3bXhmXxL8K85KOZUzMqvdve+X7eRRoDWn2J+noIb9QFtP7y2edcj4DCCAFVpxHQeEzty+FVSYzOCY+lN97eUnl2E8LttftsS4s1qpa2n4aTO2u3ufI1/ILwcPLRzKkZla+4n0nfWf/ItryLNAb07vQOa8Ko+i/hM5fPPOd6BBRGCKg6kYCeXf2mTXlWEqOban/laHFYEccLw4Tal973holx8q2btDfl2WFxV/RmhbUjx3Uk72VOzTojLE3e7Bndsj33Ig0BnZV+139OOLM+oJnLZ55zPQIKIwRUnUhAK18GX74z3rEw/ZGgXRPCsl1xx60tbZsrLy/bw7G3vdpR3nrLyHB9nPyJY8u1e5KfQGp7tSt6t4QTw5XpdTKnxitXdnQ/xFOh9ccd8c4F4az8izQENCzcEe9YEsJT9QHNXD7znOsRUBghoOpEAhovC6FlQktYkcbo0RGhdfKRoWV1cuRXkyqFGjai8r8WpC8Kl4cwfFJLaLsr7ore1mEhvFy9TubUEHbVPcQ1lZOmtIXRr+ZfpCGg89Ln03J1nPlB+szlM8+5DgGFEQKqTiWg5fvnTmyb+VDXX5gcO3zyt1+sHuq4e8609qNmzF9fu+tD501tO2F++q2lzr+GfnaY3Xmh+lOzAY2fnDe5bdrSN/ZxkYaAvnzX9BETzn28epn6v8rZffnsc+5GQGGEgKpTCeig2JP+ZGcz0mY2i4DCCAFV944O6CtHNHsmAYUCAqruHR3QeQv2f598BBQKCKi6d3RA72n2K3gCCgkEVN07OqDNI6BQQEDVEVB7BBRGCKg6AmqPgMIIAVVHQO0RUBghoOoIqD0CCiMEVB0BtUdAYYSAqiOg9ggojBBQdQTUHgGFEQKqjoDaI6AwQkDVEVB7/dr8zVH0lfa2M6v/Wj4BRQYBVUdA7fVr8+u/nP6nkD+3Kj2VgKIeAVVHQO31b/O/mPyZtKDpf2CJgKIeAVVHQO31f/PPTY2i4TtjAoosAqqOgNprYvMvDo+i+ZVT3yKgqENA1RFQe81sfvVnosMfieP/F0X/MrDf0H4goPIIqDoCaq+pzV8URf9n86tfjqJ7BvT72R8EVB4BVUdA7TW1+U0jo+irw6Jo2oB+O/uFgMojoOoIqL3mNv/EPyXfjP/KkwP53ewfAiqPgKojoPaa3PwTx0RR28PN/1b2GwGVR0DVEVB7TW/+iUea/n1sBgGVR0DVEVB7XjZPQOURUHUE1J6XzRNQeQRUHQG152XzBFQeAVVHQO152TwBlUdA1RFQe142T0DlEVB1BNSel80TUHkEVB0Btedl8wRUHgFVR0Dtedk8AZVHQNURUHteNk9A5RFQdQTUnpfNE1B5BFQdAbXnZfMEVB4BVUdA7XnZPAGVR0DVEVB7XjZPQOURUHUE1J6XzRNQeQRUHQG152XzBFQeAVVHQO152TwBlUdA1RFQe142T0DlEVB1BNSel80TUHkEVB0Btedl8wRUHgFVR0Dtedk8AZVHQNURUHteNk9A5RFQdQTUnpfNE1B5BFQdAbXnZfMEVB4BVUdA7XnZPAGVR0DVEVB7XjZPQOURUHUE1J6XzRNQeQRUHQG152XzBFQeAVVHQO152TwBlUdA1RFQe142T0DlEVB1BNSel80TUHkEVB0Btedl8wRUHgFVR0Dtedk8AZVHQNURUHteNk9A5RFQdQTUnpfNE1B5BFQdAbXnZfMEVB4BVUdA7XnZPAGVR0DVEVB7XjZPQOURUHUE1J6XzRNQeQRUHQG152XzBFQeAVVHQO152TwBlUdA1RFQe142T0DlEVB1BNSel80TUHkEVB0Btedl8wRUHgFVR0Dtedk8AZVHQNURUHteNk9A5RFQdQTUnpfNE1B5BFQdAbXnZfMEVB4BVUdA7XnZPAGVR0DVEVB7XjZPQOURUHUE1J6XzRNQeQRUHQG152XzBFQeAVVHQO152TwBlUdA1RFQe142T0DlEVB1BNSel80TUHkEVB0Btedl8wRUHgFVR0Dtedk8AZVHQNURUHteNk9A5RFQdQTUnpfNE1B5BFQdAbXnZfMEVB4BVUdA7XnZPAGVR0DVEVB7XjZPQOURUHUE1J6XzRNQeQRUHQG152XzBFQeAVVHQO152TwBlUdA1RFQe142T0DlEVB1BNSel80TUHkEVB0Btedl8wRUHgFVR0Dtedk8AZVHQNURUHteNk9A5RFQdQTUnpfNE1B5BFQdAbXnZfMEVB4BVUdA7XnZPAGVR0DVEVB7XjZPQOURUHUE1J6XzRNQeQRUHQG152XzBFQeAVVHQO152TwBlUdA1RFQe142T0DlEVB1BNSel80TUHkEVB0Btedl8wRUHgFVR0Dtedk8AZVHQNURUHteNk9A5RFQdQTUnpfNE1B5BFQdAbXnZfMEVB4BVUdA7XnZPAGVR0DVEVB7XjZPQOURUHUE1J6XzRNQeQRUHQG152XzBFQeAVVHQO152TwBlUdA1RFQe142T0DlEVB1BNSel80TUHkEVB0Btedl8wRUHgFVR0Dtedk8AZVHQNURUHteNk9A5RFQdQTUnpfNE1B5BFQdAbXnZfMEVB4BVUdA7XnZPAGVR0DVEVB7XjZPQOURUHUE1J6XzRNQeQRUHQG152XzBFQeAVVHQO152TwBlUdA1RFQe142T0DlEVB1BNSel80TUHkEVB0Btedl8wRUHgFVR0Dtedk8AZVHQNURUHteNk9A5RFQdQTUnpfNE1B5BFQdAbXnZfMEVB4BVUdA7XnZPAGVR0DVEVB7XjZPQOURUHUE1J6XzRNQeQRUHQG152XzBFQeAVVHQO152TwBlUdA1RFQe142T0DlEVB1BNSel80TUHkEVB0Btedl8wRUHgFVR0Dtedk8AZVHQNURUHteNk9A5RFQdQTUnpfNE1B5BFQdAbXnZfMEVB4BVUdA7XnZPAGVR0DVEVB7XjZPQOURUHUE1J6XzRNQeQRUHQG152XzBFQeAVVHQO152TwBlUdA1RFQe142T0DlEVB1BNSel80TUHkEVB0Btedl8wRUHgFVR0Dtedk8AZVHQNURUHteNk9A5RFQdQTUnpfNE1B5BFQdAbXnZfMEVB4BVUdA7XnZPAGVR0DVEVB7XjZPQOURUHUE1J6XzRNQeQRUHQG152XzBFQeAVVHQO152TwBlUdA1RFQe142T0DlEVB1BNSel80TUHkEVB0Btedl8wRUHgFVR0Dtedk8AZVHQNURUHteNk9A5RFQdQTUnpfNE1B5BFQdAbXnZfMEVB4BVUdA7XnZPAGVR0DVEVB7XjZPQOURUHUE1J6XzRNQeQRUHQG152XzBFQeAVVHQO152TwBlUdA1RFQe142T0DlEVB1BNSel80TUHkEVB0Btedl8wRUHgFVR0Dtedk8AZVHQNURUHteNk9A5RFQdQTUnpfNE1B5BFQdAbXnZfMEVB4BVUdA7XnZPAGVR0DVEVB7XjZPQOURUHUE1J6XzRNQeQRUHQG152XzBFQeAVVHQO152TwBlUdA1RFQe142T0DlEVB1BNSel80TUHkEVB0Btedl8wRUHgFVR0Dtedk8AZVHQNURUHteNk9A5RFQdQTUnpfNE1B5BFQdAbXnZfMEVB4BVUdA7XnZPAGVR0DVEVB7XjZPQOURUHUE1J6XzRNQeQRUHQG152XzBFQeAVVHQO152TwBlUdA1RFQe142T0DlEVB1BNSel80TUHkEVB0Btedl8wRUHgFVR0Dtedk8AZVHQNURUHteNk9A5RFQdQTUnpfNE1B5BFQdAbXnZfMEVB4BVUdA7XnZPAGVR0DVEVB7XjZPQOURUHUE1J6XzRNQeQRUHQG152XzBFQeAVVHQO152TwBlUdA1RFQe142T0DlEVB1BNSel80TUHkEVB0Btedl8wRUHgFVR0Dtedk8AZVHQNURUHteNk9A5RFQdQTUnpfNE1B5BFQdAbXnZfMEVB4BVUdA7XnZPAGVR0DVEVB7XjZPQOURUHUE1J6XzRNQeQRUHQG152XzBFQeAVVHQO152TwBlUdA1RFQe142T0DlEVB1BNSel80TUHkEVB0Btedl8wRUHgFVR0Dtedk8AZVHQNURUHteNk9A5RFQdQTUnpfNE1B5BFQdAbXnZfMEVB4BVUdA7XnZPAGVR0DVEVB7XjZPQOURUHUE1J6XzRNQeQRUHQG152XzBFQeAVVHQO152TwBlUdA1RFQe142T0DlEVB1BNTeEG9+42Ob+nZHAiqPgKojoPaGdPNbLvps9KUr+nRXAiqPgKojoPaGdPPfiBLf7ctdCag8AqqOgNobys3fGEWHH1P5n9v7cF8CKo+AqiOg9oZy8+dF0RXxxVH0z334c1ACKo+AqiOg9oZy8/Oi6N54z1FRtHj/9yWg8gioOgJqb8gDGj8WRV96br/3JaDyCKg6Ampv6AMaz4yiS/Z7XwIqj4CqI6D2BAK67rDo7/b7EpSAyiOg6gioPYGAxqf04SUoAZVHQNURUHsKAX2q8hL02f3cl4DKI6DqCKg9hYAmL0Ev6vroltz7ElB5BFQdAbUnEdANh0efezz90BOnfymauDnnvgRUHgFVR0DtSQQ0PjeKjqlk89XLP5/83c6f59yXgMojoOoIqD2NgL7+j1F02qZfTk3/bnz0UM59Cag8AqqOgNrTCGh87+FR9PnK/0SnH0dAnSKg6gioPZGAxj9K6hl96Zb4NALqFAFVR0DtqQQ0fuyEv2+9+DcxAXWLgKojoPZkAtqJgHpFQNVpBTR0+47VNTeF0KdH3mT2AAQURgioOrWAHjmqZmk8oKp1I6AE1C0Cqk4toJt6udkcAkpA3SKg6qQDunLltoFfk4ASULcIqDrpgJogoATULQKqTjqgtZsbLpzcdsrDtU5tnH9C2/i5jybvVj7SseKYtqlX74zjxeF76Skbw8i347jjrtljh4+bfcfeuNa3zsjVrlh3ke6Humdm2/hvPZzeqj8/3rzkhPajTlhafWaZM9fMPbZt5p0vE1AMEgKqzkFAb20NrRNbwoq0Uze3hvap7SFcHqdVXDj9+hUjw4Vx/EQYX05OuSosipOchjHJvZbEuQGtv0jXQy0KLeNbQrXD9edvaA/Dp4wLoX1945lXt6SnnE9AMUgIqDq1gE6cWrM2ruXupdaWa3bHOytNq3Tq+dZRD5bj8kOjw31pFWdUXnw+EFp2xOXxYV3ljPKk5M1LYfjD5XjPNWH4rryAZi7S9chhwfZ4+5IQHm84f3aYv70S6LZwesOZa0O4Yme87ZLQGdDycVWP7B468/cR0F/l3Lej8lK78CeIfon3DPUzQK/2xPHe3AMNn4OF/xzoL+Ja7s4PlyaHyl9POnVxuD+94+rwzbSKSQF3pndbnr4qXBemVV6IPnpy+gpxV3qgZ0AzF+l65FPSV7Bzw2kN57eHp5MbNy1Z1nDmGeGC9InN7Azo3uq/2xGtGrz/A+3XJfsI6LoheC4AaobuS/hxofrJf3fSqelhUvr6dGIYk1bx1113ey4cWyng0nBj19kdD+0joJmLdD3UqvTtE+HIcvb8GWH2Yx21D9WfWT4qPJl+8E4CCqBXQxbQ7SG8ld54KunUkV2vUEekVdzbdVZ5WuWV4t6xrVuTj5Sfum3p6V8L+who5iJdD5W+zozfDGFr9vxnp4Qwas4165Ou1p+5NYQ3u59Y6o2q134zdPb5Z6A5990Wx9sKf4Loj9fiN4f6KaBXlVq8nXug4XNwyAL6RmdA1yedmhbq/nA2+52ha8Oy+BdhXvKBV08KLcdffNPaHgHdnX4kc5Guh6q+RqtkcUf2/LjjP5aeNCyEk1/Jnrm9M6BP800kDBK+iaRO7ZtIPb+EH11L28+STp0THktvvL3l9caAbq58Db8gpD+G9M1wzm8rb/b2COjL6UcyF+l6qLvTt4+HYxrOT+28f1aY3XDmmNqX8KsIKAYJAVWnH9Czw4L0xhlJp24K56Y3FocVPX64c1ZYO3Jc+qeVR4TkR47iX2UC+kblIzekH8lcpOuhTk7/7PNb4fyG82eckH5D7enwtYYzz07uWvlyfxYBxSAhoOr0A7qhpeWHe+Jdlw9LOrVrQli2K+64taVtc4+A3hJODFemH5iafitp7bHp1+XpvXa1hivL8Yax6V0zF+l6qNqPMbW+1HD+rOTO8RsXJK9AM2dWvnS/fGe8Y2EgoBgkBFSdWkC7/jWmUbs6y3hdS2id1DpqVWiv3Hh0RGidfGRoWR33COjWSmJfTj9wawjHnTYxzD0ujL6xdq/FIUyd3nJ19a71F+l65O+lPxXfemvj+WtaQvvUScPCyPWNZy4LoWVC50/41yGgMEJA1akFtFtXQOMnzh3fftbGX4YpyY2N88cOn/ztF5N3GwJa+Zp6dvUy5Z/NOmrMnLvLG05ov6Z2r90/Or5t/HXlrr/K2XWRrkfedN+MtkkXrO9xfvzUvIkjRp105Za48czy/XMnts18iL8Lj8FCQNVpBbQ3d3XmUR4BhRECqk4/oN8beXXypnxqWDhkz6F/ZAK640dzL01+bICAekVA1ekH9KHQ/mA53rE0hF8O2XPoH5WAbgjJ34s6cTMBdYuAqtMPaHlBCCOTfydp8ZA9hX4SCeimf6z+zdIv305AvSKg6vQDGpfvO31c29Sz7y7v/64aNAK6d3wUDf/3y75USejnCKhTBFSdg4C6oxHQH0TRPzz161+vm5i+DCWgLhFQdQTUnkRA3/xKFP04+cirV/5DFH3p+Zz7ElB5BFQdAbUnEdBFUTSt9rHNd1//ZN59Cag8AqqOgNpTCOjrfxsd/vP93JeAyiOg6gioPYWAfjeK/mV/9yWg8gioOgJqTyCgW/vwApSA6iOg6gioPYGA9uUFKAHVR0DVEVB7Qx/QV7/YhxegBFQfAVVHQO0NfUDPjKLZ+78vAZVHQNURUHtDHNCfxPccFn1h7f7vS0DlEVB1BNTeUG7+u1H01Sv/NooW9eG+BFQeAVVHQO0N5eZfakv/5ubELX24LwGVR0DVEVB7Q7r5te2Vfk5+oS93JaDyCKg6AmpvaDf/yg2Lbnu1T/ckoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgKqE0EsAACAASURBVNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqIZNCy5/Nf8IAbXnZfMEVB4B1fCdKLoh/wgBtedl8wRUHgHVcFwUzcs/QkDtedk8AZVHQDW0R9GM/CME1J6XzRNQeQRUw1ej6Jj8IwTUnpfNE1B5BFTDF6OoLf8IAbXnZfMEVB4BlbDlsCj6av4hAmrPy+YJqDwCKmFjFEVfzD9EQO152TwBlUdAJTxdCehhW3IPEVB7XjZPQOURUAm/rAQ0ei73EAG152XzBFQeAZXw8ySgv8o9REDtedk8AZVHQCWsTgL689xDBNSel80TUHkEVMKqJKD35h4ioPa8bJ6AyiOgEm5PAnp77iECas/L5gmoPAIq4eYkoCtzDxFQe142T0DlEVAJNyYBvSb3EAG152XzBFQeAZXwwySgy3IPEVB7XjZPQOURUAlXJwFdnHuIgNrzsnkCKo+ASliRBPSi3EME1J6XzRNQeQRUQhrQb+UeIqD2vGyegMojoBLSgM7OPURA7XnZPAGVR0AlpAGdlnuIgNrzsnkCKo+ASkgDOjb3EAG152XzBFQeAZWQBrQl9xABtedl8wRUHgGVkAb0C7n/ZXgCas/L5gmoPAIqIQ1o/t/lJKD2vGyegMojoBIqAR2+j2/DE1B7XjZPQOURUAmVgM6PoqPzDhFQe142T0DlEVAJlYBecXj+d5EIqD0vmyeg8giohEpAl385+kreIQJqz8vmCag8AiohCWiIPpt3iIDa87J5AiqPgEpIAjo8OjzvEAG152XzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlUBAC+Vl8wRUHgGVQEAL5WXzBFQeAZVAQAvlZfMEVB4BlTD0Ad0UQs8PrrvwxPajZy3Z0se7D619P6UeR7xsnoDKI6ASJANaXhTCsGNHhdB6TV/u3owQNpldhoCieARUgmRAbw7tt+2O49cuC+H+Pty9GQS0VwRUHgGVIBnQKWFV9Z1F4dw+3L0ZRgFduXIbAcVQIKASFAO6PYRt1ffWhPH7v3tTjAKaIqAoHgGVMJQBXTP32LaZd75cjczG+Se0jZ/7aPLuthDWVe+xd1u1pGvPn3rESfOf6mzS4vC99MMbw8i3s6f26FfHimPapl69s+Eq80Oi7ltUz150fNvxS1/PebCOq6a3T79s+2sXT22bfNWuOEnvy3fMaDt23po4zn4Jn3kamV9bNy+bJ6DyCKiEIQzo1S2hZXxLOD+NzM2toX1qewiXJ0emh3F37ay75zUtoXVi5d63VZsWPxHGl5OPXxUWNZ6aUbnvwunXrxgZLmy4yh2LQrhw0Ztdd7y1NYyY1BqOfq3ng82ZfdOyEWHmmFOTN0vSZp4TzrpizrCWf8sGNPM0Mr+2ivIdVS+86UPl/y/tHOrngF69Fe8Y6qeAN9/8QTWgeYe2x/Hu3HNsAro2hCt2xtsuCUlknm8d9WA5Lj80OtxXObTh6BDazrxuQ7l6z3Wh5Qe7453/Go74bbVW5fHpS9TypORN9tTGgM6ohPiB0LKj4SrZL+E3trbc2BG/MSdc0PPB5u6N45tCOLfyZmUYm57ZkjzQ40eE9fUBzTyNzK8tsTeqWtXc/60ASFpZDehALtFkQM9IclWp4MwkMhfXvt2+OnwzebP1+pNbKvEZfUXykjA+M1ya3vPU5BVgWqvl6eu8dWFaucep9Sr3TVq3M81c5irZgM4LC5M3b7YeubfHgz1TufFs15v0zPPTk5aHOfUBzTyNzK8tQUCBd6IhC2j5qPBk+s6dSWSmh0lTExPDmNrxbQ8smRDCMS9W3h1d+yPRVWFWrVbPhWMr6VwabozzTu1SuW/65w1p5jJXyQZ0bO3Q8+v39HiwPXH9m/TM6ivdZ8LR9QGtfxrZX1v6q11etX6bD7vi+O2hfg7o1XZ+hxT8sBrQvEOVL3v35J5jEtCtIVT/LOCpJDJHhk4juu9SfmRCOKPywjCEt2r3HFOrVXlaeDreO7Z1a7yPU7sCurczoNmrZAK6vfOpxLkP1vAmhA2dZ+2oC2j908j+2up5+XN/vokkj28iSRiybyJ1VevpJDLTQt1V5855qfbe2tDydnfTng4jOyt2bVgW/yLMSz6aOTWr86eIMgFNr5IJaKV32zvfz3uwxoBWX6K+HsIbdQGtfxrZX1s9L5snoPIIqISh+y78mNqXuauSyJwTHktvvL2l8qgTwu21+2xLSzSqlqyfhpM7K7a58jX8gvBw8tHMqVn1Ac1eJRPQylfcz6TvrH9kW96DNQb07vQOa8Ko+i/hM08j82ur52XzBFQeAZUwdAE9u/rNmPKsJDI31f7K0eKwIo4Xhgm1L6nvDRPj5FsyC9J7zg6Lu2I2K6wdOa4jeS9zalYmoJmrZP8M9IywNHmzZ3TL9twHawjorPSnA+aEM+sDmnkamV9bPS+bJ6DyCKiEoQto5cvby3fGOxamP+qza0JYtivuuLWlbXPl5WV7OPa2VzvKW28ZGa6Pkz9JbLl2T/KTRW2vdsXslnBiuDK9TubUeOXKjn0ENHOVyoee6r7fU6H1xx3xzgXhrPwHawhoWLgj3rEkvUJ3QDNPI/NrI6AYBARUwhD+IP2yEFomtIQVaWQeHRFaJx8ZWlYnR341qVKeYSMq/2tB+mJveQjDJ7WEtru6Y7Z1WAgvV6+TOTWEXfsIaOYqlS+yjz79ta47XlP58JS2MPrV/AdrCOi89Hm3XB1nfpA+8zQyvzYCCnsEVMIQBrR8/9yJbTMf6vqLkGOHT/72i9VDHXfPmdZ+1Iz562t3fei8qW0nzH+pPopnh9mdF6o/tZeA1l8lfnDq8LG/6b7nk/Mmt01b+sY+HqwhoC/fNX3EhHMf77pyzq8g+2vr5mXzBFQeAZUw9P+YiK096U9sDqYB/TMkXjZPQOURUAnvtIC+csRgPwIBhQICKuGdFtB5Cwb7EQgoFBBQCe+0gN4z2F/BE1BIIKAS3mkBHXwEFAoIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8gioBAJaKC+bJ6DyCKgEAlooL5snoPIIqAQCWigvmyeg8giohCSgIfps3iECas/L5gmoPAIqIQnoV6Iv5x0ioPa8bJ6AyiOgEioBXfaZ6J/zDhFQe142T0DlEVAJlYBeGkVj8g4RUHteNk9A5RFQCZWAjo6iU/IOEVB7XjZPQOURUAmVgFb8W94hAmrPy+YJqDwCKiEN6Ge35B0ioPa8bJ6AyiOgEtKA5n4PiYAOAi+bJ6DyCKiENKBH5x4ioPa8bJ6AyiOgEtKAHp97iIDa87J5AiqPgEpIA/r13EME1J6XzRNQeQRUQhrQObmHCKg9L5snoPIIqIQ0oBflHiKg9rxsnoDKI6ASrk4CujD3EAG152XzBFQeAZXwgySgl+ceIqD2vGyegMojoBJu2OdfRCKgg8DL5gmoPAIq4cdJQG/IPURA7XnZPAGVR0Al/CQJ6E9yDxFQe142T0DlEVAJdycBvSf3EAG152XzBFQeAZWwOgnow7mHCKg9L5snoPIIqISHk4A+kXuIgNrzsnkCKo+ASng8CeizuYcIqD0vmyeg8giohHVJQF/JPURA7XnZPAGVR0AlvFDp5+fzDxFQe142T0DlEVAJWyoB/cf8QwTUnpfNE1B5BFTDF6JoRP4RAmrPy+YJqDwCquGfomhs/hECas/L5gmoPAKqoS2KpucfIaD2vGyegMojoBomR9G38o8QUHteNk9A5RFQDRdF0TX5RwioPS+bJ6DyCKiGZ8+8IP/HQAnoIPCyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vmyeg8gioOgJqz8vmCag8AqqOgNrzsnkCKo+AqiOg9rxsnoDKI6DqCKg9L5snoPIIqDoCas/L5gmoPAKqjoDa87J5AiqPgKojoPa8bJ6AyiOg6gioPS+bJ6DyCKg6AmrPy+YJqDwCqo6A2vOyeQIqj4CqI6D2vGyegMojoOoIqD0vm1+/evUzQ/0c0Kvf7t461E8BvXpp9eoncg80RIGAvvOsiKIfDfVzAFx7IorO6sv9COg7DwEFBoiAHrgIKDBABPTARUCBASKgBy4CCgwQAT1wEVBggAjogYuAAgNEQAFgkBFQAGgSAQWAJhFQAGgSAQWAJhFQAGgSAQWAJhFQ724Im2rvbb50XOu4S17Z3w0Anfr96dPwuURAndt1XOcCNnwthHEhtK/v/QaATv3+9Gn8XCKgrnWs/WaoLaBjajjz9fi1M8KUjt5uAKjp/6dPj88lAurZla0hdC5gdRizs/Jmx+hwf283AFQ18enT43OJgHp253e+853OBVwaFqVvF4YFvd0AUNXEp0+PzyUC6l3nAmaHe9O3Pwun93YDQLd+fvr0+FwioN51LmBKeDx9uyYc19sNAN36+enT43OJgHrXuYD28Fz69rkwsrcbALr189Onx+cSAfWuewHPp2+fDW293QDQrZ+fPj0+lwiodz2/Bpnc2w0A3fr56dPjc4mAetf9p+D3pW/vDaf1dgNAt35++vT4XCKg3nUu4JKwJH27JFza2w0A3fr56dPjc4mAete5gPvCuN2VN3vGJT/iu+8bALr189Onx+cSAfWu6++iTQkX7Yn3XBimdvR2A0C3fn769PhcIqDedS4gXt8e2k8+Iox8pvcbALr099On8XOJgHrXtYD4lUvGto69dMv+bgDo1O9Pn4bPJQIKAE0ioADQJAIKAE0ioADQJAIKAE0ioADQJAIKAE0ioADQJAIKAE0ioADQJAIKAE0ioADQJAIKX/6+lPXx5i+Vnh/14W5/XvdmkOU9ygeS51nEg6PfCCh8yQtok3nZX0Cz5SymowTUFwIKX27619QflkrVd66KBxDQT7z44ubejmsE9KUXXySgoggoXPrzUt10mw5o76epBHTwHxXNIqBwKRPQt97a1sw19lel2mWz5WzysQb4nAioKAIKlzIBbVIfqzTk30Qq7MHRbwQULmUCWruRVOa6v/vwez52zPr4lUl/dPDHWx6s3WHdCX/xgd/96+Mey16js0r7Oq/ust1vOh944ymf/+h7P/kPy3fv4yGePvEvPvjhw07q/A9GPHDMJ95/yCfHP9j9yOsm/Ml7/vMXF2+vfuSOr3304D8atabzURqvTkBFEVC4lB/QP5tV+uMvfrhUOuT7H33X//yvpdK7fpp8vHzxQbVv2f/LnvprdAc0/7zeArr43bVL/uXWvIcof7t2+5CVyc09x3f+0MBxu2vX+skH3v3pTx9cKn3qjcrt3cdUj753ZfVRGq5OQGURULiUH9B3vX9FHL/+d5Xy/P5DcfmUUukLyccvK5U+curV359xSKk0pf4a3QHNP6+XgN7+rtL7jl9x0yX/rVQalfcQi0qlgyYvu2b6+0ofTv77jSdXbk5ZfuXUShdPrl7r0N89olLODX9ZKp1SuX1i5SmMvWzJEaUPpo/SeHUCKouAwqX8gJYuTN4+WHnn8srbtw4qfajy5rUPlj716+TjL/5pqfTTumt0BzT3vN4CWnnFeGNy6/X/VPrQ3p4P8cr7Sh9IH2n1wUkgnyyVPnBPcvPeSiDXVh/xb8vJB+4qlT4Tx4+XSu+5Nbl5WfUHPhuuHhNQWQQULu0joG8lb7eXSr+zM3nnY+nHF5RK66r3eyD7ErQuoHnn9RbQqFT6TXruteed92bPh5hTKn29entk6bA4nlEqnVa9eWr1JWjlEe+Oaw9ZueoJpdJJ1eNfTW83XD0moLIIKFzKD+jHqrc7/4Jn9ePDSn9Tu9/eQ0v/ve4a3QHNPa+3gB5VKo34ZfeVGh/ii6XSc9XbL6z6WRz/Tan0QvXmc6UkqMkfGuysew5/XSptqN68Kb3dcPWYgMoioHBpX9+Fj+vfqX78U/V/8/MP665R/134nPN6C+jTH65c61OzfrixemLjQ/xB6T176x7oI6X31m7uPbj0++m1/rj+Ofxu193Xp7cbrh4TUFkEFC71I6CfrK/be+uuMYCAxpuO/730en/1rbdyHuKgzte0VQd1/5MnHysdVP+I1fe6776jeiR79ZiAyiKgcKkfAf106f/mX2MgAY3jjv+4oO2jlcj9ySs9H+LDpUPKdTfrXoG+t/SRuEdADy29p6N685nOI/VXjwmoLAIKl/oR0K+V/iT/GgMLaKK8elipNLHnQ/yPUqn2E/Qbv//9bXV/BvpC5UvzuEdADyuV1ldv/nt9KTuvHhNQWQQULvUjoItLpXuqH3/g4x//bt01BhDQ0X85vnrKb9N/EK/xIWaUSmdXb3+9dOjeeHqp9I3qzdnV77c3BPTkru/CD0tvN1w9JqCyCChc6kdAt36w9GcvJjdf+3Tpd16ou8YAAjq+dPCT6V1/WCp9tedDrH1X6UMPJLfXf6Q0Io6fKJU+uDq5ec8hpdKv4h4BfbLz50Cvrv4caMPVYwIqi4DCpX4ENF5a6desFdef87FS6cz6awwgoN8vlf7gG9fesnzcu0ul63Ie4vRS6d3HLb/urN8rHZKkcEZy88orpxxU/YtHjQGNZyZ/E+l7l40q/Vl6u/HqBFQWAYVL/Qlo+bx3Vb8//juz6r+1M5CAdhzV+T33d5+X9xAdp9QOf+iHyc3dkzvvPm1P5hFr7+0eXz36+4+ntxuvTkBlEVC41J+AxvEvJ33yfYcePuXp7DUG9E2ke9v/6tD3feJ/n/tKnP8Qj0z85PsP/ez0zn/w/r6x/+V97//EMXX/GlP2vTvaP3rwH455pvN249UJqCgCigOWoyo5eqoHFgKKA5ajKjl6qgcWAooDlqMqOXqqBxYCigOWoyo5eqoHFgKKA1b6fe5e/rvwKvjvwusioDhgEVAMFAEFg8pgrgAAAENJREFUgCYRUABoEgEFgCYRUABoEgEFgCYRUABoEgEFgCYRUABoEgEFgCYRUABoEgEFgCYRUABoEgEFgCYRUABo0v8Hit7bU7SJrB4AAAAASUVORK5CYII=" title alt width="672" /></p> +<p>We see that using the compiled model is by a factor of 12.8 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 class="r"><code>rownames(smb.1) <- smb.1$expr smb.1["median"]/smb.1["deSolve, compiled", "median"]</code></pre> <pre><code>## median -## deSolve, not compiled 12.900624 -## Eigenvalue based 1.255046 +## deSolve, not compiled 12.770742 +## Eigenvalue based 1.230531 ## deSolve, compiled 1.000000</code></pre> </div> <div id="benchmark-for-a-model-that-can-not-be-solved-with-eigenvalues" class="section level2"> @@ -132,23 +136,26 @@ smb.1["median"]/smb.1["deSolve, compiled", "median" "deSolve, not compiled" = mkinfit(FOMC_SFO, FOCUS_2006_D, use_compiled = FALSE, quiet = TRUE), "deSolve, compiled" = mkinfit(FOMC_SFO, FOCUS_2006_D, quiet = TRUE), - times = 3, control = list(warmup = 1)) -smb.2 <- summary(mb.2) + times = 3, control = list(warmup = 0))</code></pre> +<pre><code>## Warning in microbenchmark(`deSolve, not compiled` = mkinfit(FOMC_SFO, +## FOCUS_2006_D, : Could not measure overhead. Your clock might lack +## precision.</code></pre> +<pre class="r"><code>smb.2 <- summary(mb.2) print(mb.2)</code></pre> <pre><code>## Unit: seconds ## expr min lq mean median uq -## deSolve, not compiled 20.728228 20.867978 20.959811 21.007729 21.075602 -## deSolve, compiled 1.343219 1.382365 1.399697 1.421511 1.427936 -## max neval cld -## 21.143476 3 b -## 1.434362 3 a</code></pre> +## deSolve, not compiled 20.543131 20.661195 20.720383 20.779259 20.809008 +## deSolve, compiled 1.314865 1.316439 1.328049 1.318014 1.334642 +## max neval cld +## 20.83876 3 b +## 1.35127 3 a</code></pre> <pre class="r"><code>smb.2["median"]/smb.2["deSolve, compiled", "median"]</code></pre> <pre><code>## median ## 1 NA ## 2 NA</code></pre> <pre class="r"><code>autoplot(mb.2)</code></pre> -<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAABUAAAAPACAMAAADDuCPrAAAC91BMVEUAAAABAQECAgIDAwMEBAQFBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQWFhYXFxcZGRkaGhobGxscHBwdHR0eHh4fHx8gICAiIiIjIyMkJCQlJSUmJiYnJycoKCgpKSkqKiorKyssLCwtLS0uLi4vLy8wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6Ojo7Ozs8PDw9PT0+Pj4/Pz9AQEBBQUFCQkJDQ0NERERFRUVGRkZHR0dISEhJSUlKSkpLS0tMTExNTU1OTk5PT09QUFBRUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVlaWlpbW1tcXFxdXV1eXl5fX19gYGBhYWFiYmJjY2NkZGRlZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1ubm5vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8fHx9fX1+fn5/f3+AgICBgYGCgoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uMjIyNjY2Ojo6Pj4+QkJCRkZGSkpKTk5OUlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2enp6fn5+goKChoaGioqKjo6OkpKSlpaWmpqanp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKzs7O0tLS1tbW2tra3t7e4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXGxsbHx8fIyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY2NjZ2dna2trb29vc3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr6+vs7Ozt7e3u7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+/v7////z8a3rAAAACXBIWXMAAB2HAAAdhwGP5fFlAAAgAElEQVR4nO3dfZxcdX3o8WMRH1Cr1PbW2+rVe9vba3vt9QcKaqsWrS3tb7NkE0MgDyQkIRAwREgBS4AUEhDkIUqQG+RJAbEUNOGhFHkQIQKlqEEkQAgoECCJ8iB5ALLZ88c9Z2Z2M+e7m9vvd2ZzdvZ7Pp8/nDkz5/ubsy7n/ZrZmd0kKRERtVQy0gdARDRaA1AiohYDUCKiFgNQIqIWA1AiohYDUCKiFgNQIqIWA1AiohYDUCKiFgNQIqIWA1AiohYDUCKiFgNQIqIWA1AiohYDUH0vbyRy10tpunWkj2EUJVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH0ASg4DUFMCBQDVB6DkMAA1JVAAUH3tA7r+zgeG4VtINIwBqCmBAoDqaxvQZw4J4axh+B4SDV8AakqgAKD62gb0zJB1y3B8F4mGKwA1JVAAUH3tAvrEx3JAZw3Lt5FomAJQUwIFANXXLqCXhnDcXmHMsHwbiYYpADUlUABQfe0CemgI9wModVgAakqgAKD62gT0lx8Nn9kOoNRhAagpgQKA6msT0OtDODEFUOqwANSUQAFA9bUJ6GkhrABQ6rQA1JRAAUD1tQnoxBCeAVDqtADUlEABQPW1B+i6j4a/SQGUOi0ANSVQAFB97QH6gxDmAyh1XABqSqAAoPraA3RZCJcDKHVcAGpKoACg+toD9NgQ7gdQ6rgA1JRAAUD1tQfomLD3KwBKHReAmhIoAKi+tgB9au9wQAqg1HEBqCmBAoDqawvQH4RwfAqg1HEBqCmBAoDqawvQi0O4NAVQ6rgA1JRAAUD1tQXoCSGsTAGUOi4ANSVQAFB9bQF6YAi/SgGUOi4ANSVQAFB97QD6/D7hr/M1AJQ6LAA1JVAAUH3tAHp/CEfmawAodVgAakqgAKD62gH0OyF8NV8DQKnDAlBTAgUA1dcOoKeHcHO+BoBShwWgpgQKAKqvHUBnhPBEvgaAUocFoKYECgCqrx1APx0+vj1fA0CpwwJQUwIFANXXBqAPhzC1tgaAUocFoKYECgCqrw1AV4SwuLYGgFKHBaCmBAoAqq8NQM8N4braGgBKHRaAmhIoAKi+NgA9PITVtTUAlDosADUlUABQfSZAN3zti1dtGNjaL+z7em0NAKUOC0BNCRQAVJ8J0NtDCMesb2w8FMKU+hoASh0WgJoSKACoPhOg12eAhq81Nv4lhDPrawAodVgAakqgAKD67IDu81B949QQ/rW+BoBShwWgpgQKAKrPCuhHQ5hf3zgwhGfrawAodVgAakqgAKD6rICe9Imw90/y60/uHf6usQaAUocFoKYECgCqzwrokiUhnJBfXx7CyY01AJQ6LAA1JVAAUH1mQH+1b9hndXZ9QQg3NtYAUOqwANSUQAFA9ZkBTReHcNrGjes/G/Z+obEGgFKHBaCmBAoAqs8O6DMfCR97bOMtIczsXwNAqcMCUFMCBQDVZwc0PSmEEzfODeHa/jUAlDosADUlUABQfS0Aum7fsNdZe4VPbu5fA0CpwwJQUwIFANXXAqDpsvzz9GHZwBoASh0WgJoSKACovlYA7T068/OQ1wbWAFDqsADUlEABQPW1Amjae/Xxy7bsWANAqcMCUFMCBQDV1xKgIgClDgtATYkTGkD1ASg5DEBNiRMaQPUBKDkMQE2JExpA9QEoOQxATYkTGkD1ASg5DEBNiRMaQPUBKDkMQE2JExpA9QEoOQxATYkTGkD1ASg5DEBNiRMaQPUBKDkMQE2JExpA9QEoOQxATYkTuiRAn41x8I2PnnnU+MnzL1yv3H3XtPOHGnQPgJLDANSUOKFHDtC+pTGOOfSgGLuv1uy+C4rxWQClqgegpsQJPXKA3hDH3/x6mr5wUYw/Uuy+CwJQIgC1JU7okQP0sHhH/crSeJpi913Q8uWbAJSqHoCaEif0iAG6OcZN9Wur4rT/fPddFoBStQNQU+KE3uWArlp0aM+8W9fVMXrq7CN7pi36cX51U4yP1vfYvqku6eozZo/7wtmP5Fdzuy6I36jd/FSc8GpxVPbEWUf0HLHsxSEW6b1i7vi5F21+4ZzZPbOuyP+qcYzrvn90z6GLV6Vp8SV8YfnCMe8IQMlhAGpKnNC7GtCrumLXtK54Rg2jG7rj+NnjY7w4v2dunHrb1qY9r+6K3TOyvW9O66o9FKf15bdfEZfK0WI3dcexM7vj5BcGL7Lw+BWXjo3zJh2XX1yY5maeGk+6ZOGYrm+nBUALyxeOuSkAJYcBqClxQu9iQFfHeMnWdNO5McfoF90H3duX9t13cLw7u+vxyTH2LLj28b76no/Gru+8nm79ehz367pqfdNqT1H7ZuYXxdFCT3V3fa83fXlh/PLgRRZtT9MVMZ6WXSyPU9LczK58gQfHxTXNgBaWLxxzXt9x9Va9aujWnQF6gGUVol3d62naO9LHMIoSJ/QuBvTEnLWMoHk5Ruc03m5fGU/OL1667piuDKmDL8mfOqYL4nm1PY/LnynWVLus9nzw0Tinb9Boc4vj+fnFb7oP3D5okbXZxhMDF2lu5hm1ocviwmZAC8sXjjlve6h3h+Urv3NngI61rEJEndyuBbRvYny4duXWHKO5cebsvBlxUuP+TfdcOD3GQ57Orh7c+JHoHXF+Q7Un46EZncvi99KhRgea0hj8xZptgxbZljZfpLmZ9Wewa+PkZkCbly8ecx6AEtGQ7VpAX4rxN7Urj+QYHRj7a0Kk74Hp8cTsCWSMrzT2nNRQrW9OfCzdPqX7pXQno7U29z9EOuQi4iLGx/untjQB2rx88ZhrPVNv4wuGbtzpz0AtqxDt6rL/2l8d6WMYRYkTetcCOqDbYzlGc2LTT2AXLXymcW117Hp1h32PxQn92l0TL01/GhfntxZGC2Xebe6/PtQiEtD6U9QXY3y5CdDm5YvH3BxvIpHDeBPJlDihd/HPQCc1Xg7fkWN0avxJbePV9S+m6fR4S2OfTTWxDmrQ9oN4TL92z2ev4ZfE+/NbC6OFslfca2tX1jywaahFJKC313ZYFQ9qfglfWL5wzM0BKDkMQE2JE3oXA3pK/U2bvvk5Risav3J0QbwyTc+P0xsvve+KM9L8rZuaN33HxwsG0JsfV0+Y2ptfK4wWOzEuyy+2Hdy1echFBKDza+/6L4wLmgEtLF845uYAlBwGoKbECb2LAc1eBl+8Nd1yfu0jQa9Nj5e+lvbe1NXzfPb0cnw89OYNvX0v3TghXpfmP3HsumZb/gmkng0D6N0Yj4qX19YpjKbLl/fueIhHYvf1venWJfGkoRcRgMbzt6RbLozxkWZAC8sXjrk5ACWHAagpcULv6g/SXxpj1/SueGUNox+Pjd2zDoxdK/N7fj4zE2rM2Ox/ltSeFF4W4wEzu2LPbekAei+NiXFdfZ3CaIyvNT3E1dnQYT3x4A1DLyIAXVw7nq6r0sIH6QvLF465KQAlhwGoKXFC72pA+360aEbPvPsGfmFyygGzvvJ0/a7e2xfOGT/x6LPXNHa97/TZPUeeXXtrqf/X0E+Jx/cv1DxaBDR9ePGsnjnLXt7JIgLQdbfNHTv9tAfryzT/KueO5YvHvCMAJYcBqClxQo/Kv0i/rfbJzlaqmdlqAEoOA1BT4oQelYA+N67VSQAlKgagpsQJPSoBXTwUTaoAlKgYgJoSJ/SoBPTOVl/BAyiRCEBNiRN6VALaegBKVAxATYkTumKAthWAksMA1JQ4oQFUH4CSwwDUlDihAVQfgJLDANSUOKEBVB+AksMA1JQ4oQFUH4CSwwDUlDihAVRfa4BuffjlpjUAlDosADUlUABQfS0Betenw8eW71gDQKnDAlBTAgUA1dcKoA/vG0LYe9XAGgBKHRaAmhIoAKi+FgDtnRDCfiFMHVgDQKnDAlBTAgUA1dcCoNeF0PXE/iEMPAUFUOqwANSUQAFA9dkB3fb3Idy4cVkIi/rXAFDqsADUlEABQPXZAV0RwqSNG9d+JOy3vbEGgFKHBaCmBAoAqs8M6PZxIdyQXZ8Rwk8bawAodVgAakqgAKD6zIDeHkLPhux69hr+6401AJQ6LAA1JVAAUH1WQM87KISr8+v/EcL0xhoASh0WgJoSKACoPiug3SHE5/PrG/YL+2ytrwGg1GEBqCmBAoDqswKadU19Y24ID9TXAFDqsADUlEABQPXZAZ28ob6xLIRL6msAKHVYAGpKoACg+kyA3pD5+akHGxs/DOGL9TUAlDosADUlUABQfSZA13wm9NzXv/HsPuFz9TUAlDosADUlUABQfSZAN677edPGgSGsr60BoNRhAagpgQKA6rMBWmhBCHfW1gBQ6rAA1JRAAUD1tQHoN0P4v7U1AJQ6LAA1JVAAUH1tALoyhHm1NQCUOiwANSVQAFB9bQD67EfD/rU1AJQ6LAA1JVAAUH1tALpxXAgv5msAKHVYAGpKoACg+toB9LgQ7svXAFDqsADUlEABQPW1A+iyEK7I1wBQ6rAA1JRAAUD1tQPo90NYkK8BoNRhAagpgQKA6msH0LUhTMjXAFDqsADUlEABQPW1A+jGvwkffS0FUOq4ANSUQAFA9bUF6JwQVqcASh0XgJoSKACovrYA/UoIy1MApY4LQE0JFABUX1uAfjeEL6cASh0XgJoSKACovrYA/Vn930UCUOqwANSUQAFA9bUF6IZPhb/YDqDUcQGoKYECgOprC9D8H4d/AkCp4wJQUwIFANXXHqBnhHADgFLHBaCmBAoAqq89QOvvIgEodVgAakqgAKD62gP05yFMBlDquADUlEABQPW1B+jGz4WPbgFQ6rQA1JRAAUD1tQnovPwv2gEodVgAakqgAKD62gT04hCWAih1WgBqSqAAoPraBPTH+Q9BAZQ6LAA1JVAAUH1tArpx/7DXegClDgtATQkUAFRfu4CeGsI3AZQ6LAA1JVAAUH3tAroyhM+GMG5Yvo1EwxSAmhIoAKi+dgHdODlknTsc30Wi4QpATQkUAFRf24Cu3CeEv3tqOL6LRMMVgJoSKACovrYB3XjHF05YPQzfQ6LhC0BNCRQAVF/7gBJ1XABqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYQBqSqAAoPoAlBwGoKYECgCqD0DJYdUGdP1DxgGBAoDqA1ByWKUB3TA9LLBNCBQAVB+AksMqDeiqEPZ5zjQhUABQfQBKDqs0oDeFEB40TQgUAFQfgJLDKg3olRmgPzRNCBQAVB+AksMqDegFGaDXmyYECgCqD0DJYZUG9MwM0G+bJgQKAKoPQMlhlQb05AzQi00TAgUA1Qeg5LBKAzo/A/RrpgmBAoDqA1ByWKUBnZMBepZpQqAAoPoAlBxWaUCnZ4AuNE0IFABUH4CSwyoN6MQM0BNMEwIFANUHoOSwSgM6NgP0i6YJgQKA6gNQclilAd0/A3SOaUKgAKD6AJQcVmlA/yoDdIZpQqAAoPoAlBxWaUA/ngE62TQhUABQfQBKDqsyoBv2zgCdYBoRKACoPgAlh1UZ0GczP8NY04hAAUD1ASg5rMqAPpEDur9pRKAAoPoAlBxWZUAfyQH9rGlEoACg+gCUHFZlQFflgH7KNCJQAFB9AEoOqzKg9+eAfsw0IlAAUH0ASg6rMqArc0D3No0IFABUH4CSw6oM6B0AWmIASg6rMqA3A2iJASg5rMqAXg+gJQag5LAqA3otgJYYgJLDqgzotwG0xACUHFZlQC8H0BIDUHJYlQH9BoCWGICSw6oM6AUAWmIASg6rMqBLALTEAJQcVmVAzwLQEgNQcliVAV0EoCUGoOSwKgO6EEBLDEDJYVUG9EsAWmIASg6rMqDzAbTEAJQcVmVA5wJoiQEoOazKgB4Rwl4AWlYASg6rMqAzQtgXQMsKQMlhVQZ0agifBNCyAlByWJUBnRjCZwG0rACUHFZlQMeF8HcAWlYASg6rMqBdYe9uAC0rACWHVRnQ/cO+BwBoWQEoOazKgO4XPgWgpQWg5LAqA/qX4a8BtLQAlBxWZUD3CRFASwtAyWEVBnTDXmEsgJYWgJLDKgzosyFMBNDSAlByWIUBfTKEqQBaWgBKDqswoI+FMBNASwtAyWEVBvShEI4E0NICUHJYhQH9cQhfBNDSAlByWIUBvS+EEwC0tACUHFZhQO8M4RQALS0AJYdVGNBbQ1gMoKUFoOSwCgN6UwhnA2hpASg5rMKAfjeE8wG0tACUHFZhQL8TwkUAWloASg6rMKDfCuGbAFpaAEoOqzCgF4fwHQAtLQAlh1UY0KUhLAfQ0gJQcliFAT0nhH8D0NICUHJYhQE9PYQ7AbS0AJQcVmFATwnhPgAtLQAlh1UV0A2Xnzs3hFUAWloASg6rKqDLQ9grhMcAtLQAlBxWVUAX1QB9BkBLC0DJYVUF9KSQ92sALS0AJYdVG9CtAFpaAEoOqzSgH0kBtLQAlBxWaUA/DaDlBaDksEoD2gWg5QWg5LBKAzoJQMsLQMlhlQb0SAAtLwAlh1Ua0JMBtLwAlBxWaUAvBNDyAlByWKUBvRtAywtAyWEVBnTKBX0AWl4ASg6rMKD/kZ3VAFpaAEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUND74pDcAABvRSURBVCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEABQPUBKDkMQE1jAgUA1Qeg5DAANY0JFABUH4CSwwDUNCZQAFB9AEoOA1DTmEChTUCfjXHwjY+eedT4yfMvXK/cfWTb+SENugdAyWEAahoTSgw/oH1LYxxz6EExdl+t2b2VYnx22JYBUKp2AGoaE0oMP6A3xPE3v56mL1wU448Uu7cSgBINVwBqGhNKDD+gh8U76leWxtMUu7fSMAG6fPkmAKWqB6CmMaHEsAO6OcZN9Wur4rT/fPeWGiZAawEoVTsANY0JJVoGdNWiQ3vm3bqujsxTZx/ZM23Rj/Orm2J8tL7H9k11SVefMXvcF85+pN+kC+I3ajc/FSe8Whwd5FfvlYf0zL5qq1jl7JjX9BbVE2cd0XPEsheHeLDeK+aOn3vR5hfOmd0z64rX0pzedd8/uufQxavStPgSvnAYha9tRwBKDgNQ05hgqlVAr+qKXdO64hk1ZG7ojuNnj4/x4vyeuXHqbVub9ry6K3bPyPa+uW5a+lCc1pfffkVcKkcLZfueP/e6KyfEM8Uq318a45lLfzOw403dcezM7jj5hcEPtvD4FZeOjfMmHZdfXFgz89R40iULx3R9Oy0AWjiMwtcGoOQ7ADWNDQ+gq2O8ZGu66dyYI/OL7oPu7Uv77js43p3d9fjkGHsWXPt4X33PR2PXd15Pt349jvt1Xau+abWnqH0z84viqAT06Azie2LXFrFK8SX8U91d3+tNX14Yvzz4wRZtT9MVMZ6WXSyPU2qTXfkDPTgurmkGtHAYha8tb3tXvR/1ErkrOzn6RvoYRqJFA4CaxoYH0BNzrrL/4+flyJzTeLt9ZTw5v3jpumO6MnwOviR/SpguiOfV9jwufwZY0+qy2vO8R+OcvkGjAtDcuq015gqrFAFdHM/PL37TfeD2QQ+2Ntt4YuCiNnlGbeiyuLAZ0MJhFL62vO2h3h2t/X9FRJ3X4gFA21mlNUD7JsaHa1duzZGZG2fOzpsRJzXu33TPhdNjPOTp7OrBjR+J3hHnN7R6Mh6a0bksfi8danSgbN/a0+Uac4VVioBOadz1izXbBj3YtrT5ojZZf6a7Nk5uBrT5MIpfWx6AErlrJAF9Kcb6zyAfyZE5MPY3dscufQ9MjydmTwxjfKWx56SGVn1z4mPp9indL6U7Ga2X7bs9v8yZK65SAHRz/6GkQz6YuIjx8f6pLU2ANh9G8Wtrjp+BksP4GahpTDDVGqADaj2WIzMnNq26aOEzjWurY9erO0x7LE7oV+yaeGn607g4v7UwWqz/U0QFQGurFADNvNvcf32oB5OA1p+ivhjjy02ANh9G8WtrDkDJYQBqGhNMtfgz0EmNl7l35MicGn9S23h1/YtpOj3e0thnU02igxpk/SAe06/Y89lr+CXx/vzWwmixZkCLqxQAzV5xr61dWfPApqEeTAJ6e22HVfGg5pfwhcMofG3NASg5DEBNY4KpFgE9pf5mTN/8HJkVjV85uiBemabnx+mNl9R3xRlp/pbMktqex8cLBjCbH1dPmFp7O6swWqwAaGGV4s9AT4zL8ottB3dtHvLBBKDza58OWBgXNANaOIzC19YcgJLDANQ0JphqEdDs5e3FW9Mt59c+6vPa9Hjpa2nvTV09z2dPL8fHQ2/e0Nv30o0T4nVp/pPErmu25Z8s6tkwgNmN8ah4eW2dwmi6fHnThwQKgBZWyW56ZMd+j8Tu63vTrUviSUM/mAA0nr8l3XJhbYUdgBYOo/C1ASg5D0BNY8MDaHppjF3Tu+KVNWR+PDZ2zzowdq3M7/n5zEyeMWOz/1lSe7J3WYwHzOyKPbftwOylMTGuq69TGI3xtZ0AWlgle5E9+UsvDOx4dXbzYT3x4A1DP5gAdHHtuLuuSgsfpC8cRuFrA1DyHYCaxoYJ0L4fLZrRM+++gV+EnHLArK88Xb+r9/aFc8ZPPPrsNY1d7zt9ds+RZz/TjOIp8fj+hZpH/z+ANq+S3jv7gCm/2rHnw4tn9cxZ9vJOHkwAuu62uWOnn/bgwMpDfAXFr21HAEoOA1DTmICwo/4i/bbaJzZ3ZW39GRIAJYcBqGlMoNBRgD43blc/AoASFQNQ05hAoaMAXbxkVz8CgBIVA1DTmEChowC9c1e/ggdQIhGAmsYECh0F6K4PQImKAahpTKBQMUDbCkDJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHAahpTKAAoPoAlBwGoKYxgQKA6gNQchiAmsYECgCqD0DJYQBqGhMoAKg+ACWHVRjQE24B0BIDUHJYhQENYRWAlheAksMqDejFAFpeAEoOqzSg/wSg5QWg5LBKAzoXQMsLQMlhlQZ0KoCWF4CSwyoN6AEAWl4ASg6rNKCfAdDyAlByWKUB3QdAywtAyWGVBjS8BqClBaDksGoD+iKAlhaAksOqCujCGqDrALS0AJQcVlVAr6wB+jiAlhaAksOqCuizC+fNCOFnAFpaAEoOqyqgG2s/Bv13AC0tACWHVRjQxSH8EEBLC0DJYRUG9OwQbgHQ0gJQcliFAf1aCCsAtLQAlBxWYUAvCuEaAC0tACWHVRjQy0O4AkBLC0DJYRUG9OoQLgbQ0gJQcliFAf1uCEsBtLQAlBxWYUBvDOEcAC0tACWHVRjQW0I4A0BLC0DJYRUG9IchLATQ0gJQcliFAb0nhC8BaGkBKDmswoD+RwjHAmhpASg5rMKAPhTCUQBaWgBKDqswoI+GMAtASwtAyWEVBvSJEKYBaGkBKDmswoCuC2EigJYWgJLDKgzo+hDGAWhpASg5rMKAbtwnRAAtLQAlh1UZ0L8InwPQ0gJQcliVAd0vfBpASwtAyWFVBvRvwscBtLQAlBxWZUC7wt4AWloASg6rMqA9Ifw9gJYVgJLDqgzoxBD+GkDLCkDJYVUGdEoInwLQsgJQcliVAZ0RwscAtKwAlBxWZUAPD2EvAC0rACWHVRnQuSELQEsKQMlhVQb0WAAtMQAlh1UZ0BMAtMQAlBxWZUBPBtASA1ByWJUBPQ1ASwxAyWFVBvRMAC0xACWHVRnQcwG0xACUHFZlQJcCaIkBKDmsyoBeBKAlBqDksCoDejmAlhiAksOqDOhVAFpiAEoOqzKg/wKgJQag5LAqA7oCQEsMQMlhVQb0ZgAtMQAlh1UZ0NsBtMQAlBxWZUDvBtASA1ByWJUB/fcc0H1MIwIFANUHoOSwKgP60xzQT5pGBAoAqg9AyWFVBnR1DuhnTCMCBQDVB6DksCoDujYHdH/TiEABQPUBKDmsyoA+kwPabRoRKACoPgAlh1UZ0A17ZYB+3jQiUABQfQBKDqsyoBv3zQCdZJoQKACoPgAlh1Ua0E9ngM4wTQgUAFQfgJLDKg3o32aAHmGaECgAqD4AJYdVGtDuDNB5pgmBAoDqA1ByWKUBnZgBerxpQqAAoPoAlBxWaUAPyQBdaJoQKACoPgAlh1Ua0MMzQL9smhAoAKg+ACWHVRrQL2aALjFNCBQAVB+AksMqDeiCDNCLTBMCBQDVB6DksEoDekYG6JWmCYECgOoDUHJYpQH9Wgbo90wTAgUA1Qeg5LBKA/rNDNDbTRMCBQDVB6DksEoDen0G6E9MEwIFANUHoOSwSgN6fwh7rzNNCBQAVB+AksMqDej6rnCYbUKgAKD6AJQcVmlAN669/hnbgEABQPUBKDms2oCaEygAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1ByGICaEigAqD4AJYcBqCmBAoDqA1By2FMrV/58pI9hFCVQAFCiSvdACItG+hhGbwBKVOkAtJ0AlKjSAWg7AShRpQPQdgJQokoHoO0EoESVDkDbCUCJKh2AthOAEhG1GIASEbUYgBIRtRiAEhG1GIASEbUYgBIRtRiAEhG1GIASVbXvxmcb154/b2r31HOfG9GjGZUBKFFFe+3wfkAf/3yMU2Mcv2ZkD2gUBqBElax39cmxAWjv7LjgxfSFE+NhvSN8UKMuACWqYpd3x9gP6Mo4aWt2seXg+KORPajRF4ASVbFbv/rVr/YDel5cWrs8Py4ZyUMajQEoUVXrB/T4eFft8ofxSyN5OKMxACWqav2AHhYfrF2uioeP5OGMxgCUqKr1Azo+Plm7fDJOGMGjGZUBKFFV2wHoL2qXT8SekTyc0RiAElW1wS/hZ43k4YzGAJSoqu14E+nu2uVd8YSRPJzRGIASVbV+QM+NF9YuL4znjeThjMYAlKiq9QN6d5z6enaxbSofpLcGoERVbeBXOQ+LZ21Lt50ZZ/OrnMYAlKiq9QOarhkfxx8zLk5YO7LHMwoDUKKqNgBo+ty5U7qnnLd+RI9mVAagREQtBqBERC0GoERELQagREQtBqBERC0GoERELQagREQtBqBERC0GoERELQagREQtBqBERC0GoERELQag5KX9kmLva32p2nwYvkMbvP4fD7rtbfljDr6ZOjoAJS8NBWiLJAEo6QJQ8tKKr9f6vSSpX7kibQPQ9z/99PPDfHyF9Qcf1jNPPw2goy4AJWf9cdL0H3XLgO5ayXayPoCOugCUnFUA9JVXNrWyBoCSLgAlZxUAbTEAJV0ASs4qANrYyGW69pPvfNN7D1mTPjfz93d/X9e9jR0ePfJP3vbbf3b4T4prDEj21LEfe8+bP/CZy17fye6PHfUnb3/nXl/o/5cx7jnk/W/d4wPT7t2xyqPT//BN//UTF2yu3/L9z79n998/aFX/+nJ1AB11ASg5a2hA/2h+8gefeGeS7PGt97zh//yPJHnDD/Lb+87ZrfGW/Re3Na/RL9kFb2zc/cGXhtq97yuN7T2W55vbjuj/AMDhrzdW+be3vfHDH949ST70crb9+iH1e9+8vL6+WB1AR2EASs4aGtA3vPXKNH3xk5lW774v7Ts2ST6e335RkrzruKu+dfQeSXJY8xoNyW55Q/KWI65cce7/TJKDhtp9aZLsNuvSq+e+JXln/u9ZHpNtHnbZ5bMzF4+pr7Lnb4/L5Hz8g0lybLZ9VHYYUy66cFzy9tr6cnUAHYUBKDlraECTM/PLe7MrF2eXr+yWvCO7eOHtyYc25rc//d+T5AdNazQky54xfi+/fPG/JO/YPnj3596SvK02tXL3HMiHk+Rtd+abd2VArq4/6l/05TfcliQfSdMHk+RNN+WbF9U/8ClWTwF0FAag5KydAPpKfrk5SX5ra37lvbXblyTJo/X97ik+BW1IFpLkV7Xta04//TeDd1+YJP9Q356Q7JWmRyfJCfXN4+pPQbNHvT1tPGy23pFJ8oX6/Z+rbYvVUwAdhQEoOWtoQN9b3+7/Bc/67WOSP2/st33P5H81rdGQbGKSjP3Zjlvl7p9Ikifr27+844dp+udJ8sv65pNJDmr+g4OtTev9WZI8Xt9cUdsWq6cAOgoDUHLWzt6FT5uv1G//UPNvfv5e0xqN3R57Z3b7h+b/81P1W+Xuv5u8aXvT0LuSNzc2t++evLu2yh80r/fbA7uvqW2L1VMAHYUBKDnLAOgHmkV8c9Ma/bs9e8Tv1O770396ZYjdd+t/Xltvtx1/vuS9yW7Nj1q/tmP3LfV7iqunADoKA1BylgHQDyd/O/QaOyTr/fcv97wnQ+4Pnxu8+zuTPfqaNpuegb45eVc6CNA9kzf11jfX9t/TvHoKoKMwACVnGQD9fPKHQ69RlKxv5ZgkmTF49/+dJI1P0D/1rW9tavoZ6C+zl+bpIED3SpI19c1/bV6/f/VBD0ujIAAlZxkAvSBJ7qzffs/73ve1pjUaux38wWn17V/X/rid3P3oJDmlvv0PyZ7b07lJ8o/1zePr77cLQI8ZeBd+TG1brJ4C6CgMQMlZBkBfenvyR0/nmy98OPmtXzat0dhtWrL7w7Xtf06Szw3effUbknfck2+veVcyNk0fSpK3r8w379wjSX6eDgL04f7PgV5V/xyoWD0F0FEYgJKzDICmyzLz5l953anvTZIFzWs0dvtWkvzuP15z42VT35gk1w6x+5eS5I2HX3btSb+T7JFTeHS+efnlh+1W/8UjCWg6L/9NpG9cdFDyR7VtuTqAjsIAlJxlAbTv9DfU31P/rfnNbwf179Y7sf899zeePtTuvcc27n7HP+ebr8/q333OtsKjNq69Pq1+77sfrG3L1QF0FAag5CwLoGn6s5kfeMueex/2WHGNgf3vGv+ne77l/X952nPp0Ls/MOMDb91zn7n9f7z+7in/7S1vff8h94pVBq59f/x7dv+9SWv7t+XqADrqAlCiQY2QZAA66gJQokEBKOkCUKJBASjpAlCiQQEo6QJQokHV3hvfhf8u/FDx78KPxgCUaFAASroAlIioxQCUiKjFAJSIqMUAlIioxQCUiKjFAJSIqMUAlIioxQCUiKjFAJSIqMUAlIioxQCUiKjFAJSIqMUAlIioxQCUiKjF/h/dRVyxsh6FSgAAAABJRU5ErkJggg==" title alt width="672" /></p> -<p>Here we get a performance benefit of a factor of 14.8 using the version of the differential equation model compiled from C code using the inline package!</p> +<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAABUAAAAPACAMAAADDuCPrAAAC7lBMVEUAAAABAQECAgIDAwMEBAQFBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQWFhYXFxcZGRkaGhobGxscHBwdHR0eHh4fHx8gICAiIiIjIyMkJCQlJSUmJiYnJycoKCgpKSkqKiorKyssLCwtLS0uLi4vLy8wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6Ojo7Ozs8PDw9PT0+Pj4/Pz9AQEBBQUFDQ0NERERFRUVGRkZHR0dISEhJSUlKSkpLS0tNTU1OTk5PT09QUFBRUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVlaWlpbW1tcXFxdXV1eXl5fX19gYGBhYWFiYmJjY2NkZGRlZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1ubm5vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8fHx9fX1+fn5/f3+AgICBgYGCgoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uMjIyNjY2Ojo6Pj4+QkJCRkZGSkpKTk5OUlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2enp6fn5+goKChoaGjo6OkpKSlpaWmpqanp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKzs7O0tLS1tbW2tra3t7e4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXGxsbHx8fIyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY2NjZ2dna2trb29vc3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr6+vs7Ozt7e3u7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+/v7///8pyQqHAAAACXBIWXMAAB2HAAAdhwGP5fFlAAAgAElEQVR4nO3de5xcZZ3n8eMgXvDKOLPjzujqrrOzzqyzHnC8X2cdmafTpImBXEkgJBAhQYyAcssA4SZCFDALIndRhJEoyspNkOs6XsIoAUJAglwzcpGEAEn6/Ld1qqqTql83M99TnnO66/d83n9Qdaqe56nTPPB5VXVVdycZAKAnyXifAAD0KwIKAD0ioADQIwIKAD0ioADQIwIKAD0ioADQIwIKAD0ioADQIwIKAD0ioADQIwIKAD0ioADQIwIKAD0ioLpn1gNOvZBlT473OfQFEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCguoIBfbCaDQQqQEBFJgoEVFcsoF9MF1W0hUDpCKjIRIGA6goF9O5d0vRfqtpEoGQEVGSiQEB1hQJ6dZqml1a1iUDJCKjIRIGA6goF9HuNgJ5a1SYCJSOgIhMFAqorHNDDqtpEoGQEVGSiQEB1hQM6v6pNBEpGQEUmCgRUVzigU6vaRKBkBFRkokBAdYUD+g9VbSJQMgIqMlEgoLrCAX1/VZsIlIyAikwUCKiucEDTR6raRaBcBFRkokBAdcUDek9VuwiUi4CKTBQIqK54QPlRJPQJAioyUSCguuIBvbGqXQTKRUBFJgoEVFc8oD+saheBchFQkYkCAdUVD+h3qtpFoFwEVGSiQEB1xQN6QVW7CJSLgIpMFAiornhAv1bVLgLlIqAiEwUCqiseUH4dE/oEARWZKBBQXfGAHlfVLgLlIqAiEwUCqise0C9UtYtAuQioyESBgOqKB/TgqnYRKBcBFZkoEFBd8YDyC0HRJwioyESBgOqKB3R2VbsIlIuAikwUCKiueECnVLWLQLkIqMhEgYDqigc0VLWLQLkIqMhEgYDqigf041XtIlAuAioyUSCguuIB/UBVuwiUi4CKTBQIqK54QHd5vKptBEpFQEUmCgRUVzyg6bqqthEoFQEVmSgQUF0PAeVveqA/EFCRiQIB1fUQ0F9WtY1AqQioyESBgOp6COjtVW0jUCoCKjJRIKC6HgJ6Q1XbCJSKgIpMFAioroeA8keR0B8IqMhEgYDqegjoFVVtI1AqAioyUSCguqIB3SVNL65qG4FSEVCRiQIB1RUN6PvT9NyqthEoFQEVmSgQUF3RgH4iTc+oahuBUhFQkYkCAdUVDWhI01Oq2kagVARUZKJAQHVFAzo1Tf+pqm0ESkVARSYKBFRXNKBz+Kty6BcEVGSiQEB1RQN6AH9VDv2CgIpMFAiormhAl6Tpgqq2ESgVARWZKBBQXdGAHsVflUO/IKAiEwUCqisa0JPS9NNVbSNQKgIqMlEgoLqiAV2+K39VDn2CgIpMFAiornBAP5x+oqptBEpFQEUmCgRUVzign2r+VbmHfvpERXsJlIWAikwUCKiucED3yP+q3P27pSdWtZlASQioyESBgOoKB3RWmj6w/qw03auqzQRKQkBFJgoEVFc4oPun6a/XzyegmPgIqMhEgYDqCgf0c2n60/WfJKCY+AioyESBgOoKB/ToNL3x/pSAYuIjoCITBQKqKxzQU9L0hzcRUPQBAioyUSCgusIB/VqaXv4tAoo+QEBFJgoEVFc4oBem6QVfJaDoAwRUZKJAQHWFA/rPaXrWkQQUfYCAikwUCKiucECvTdNT5hNQ9AECKjJRIKC6wgG9I02PnkxA0QcIqMhEgYDqCgd0dZp+9oMEFH2AgIpMFAiornBAH07TaSkBRR8goCITBQKqKxzQjWn6PgKKfkBARSYKBFRXOKBZM58EFBMfARWZKBBQXfGAfoqAoj8QUJGJAgHVFQ/oDAKK/kBARSYKBFRXPKCLCSj6AwEVmSgQUF3xgC4joOgPBFRkokBAdcUDen6jnh8ioJj4CKjIRIGA6ooH9IZGQPcmoJj4CKjIRIGA6ooH9PFGQL9MQDHxEVCRiQIB1RUPaHZA+v67CSgmPgIqMlEgoLoeArpx5ZqNBBQTHwEVmSgQUF0PAW0goOgDBFRkokBAdQQUbhFQkYkCAdURULhFQEUmCgRUR0DhFgEVmSgQUB0BhVsEVGSiQEB1BBRuEVCRiQIB1RFQuEVARSYKBFRHQOEWARWZKBBQHQGFWwRUZKJAQHUEFG4RUJGJAgHVEVC4RUBFJgoEVEdA4RYBFZkoEFAdAYVbBFRkokBAdQQUbhFQkYkCAdURULhFQEUmCgRUR0DhFgEVmSgQUB0BhVsEVGSiQEB1BBRuEVCRiUJNAX0khNE33nPygVNmLlnxuDi8Gi/9UKPuIaBwi4CKTCXGL6DDZ4YwaZ9pIQxeqgyvQAiPEFAgR0BFphLjF9DvhylXv9jYtXNCuFUYXgECCrQRUJGpxPgFdH64oXXlzHCcMLwCV165gYACOQIqMpUYt4BuDGFD69qqMOc/Hl4ZAgoQUJmpROUBXXX8PkMHX/twK0brTv3M0Jzjf55f3RDCPa0RWze0Srr6xAV7HHTq3fnVvF1nha83b14Xpj7fPdW6/5QDhg44+6kxFtly0aIpi87Z+OSXFwztd1HjP5HGq/aHr1k8tM+yVVnW/RK+a/muc96OgMItAioy8ak6oJcMhIE5A+HEZoy+PximLJgSwrn5PYvC7Os2dYy8dCAM7tsYfXXWqtqvwpzh/PaLwpl2arcfDIbJ8wbDzCdHL7L0sJXnTQ4Hzzg0v1iR5c08Nhz1jaWTBr6ZdQW0a/muc+5AQOEWARWZ+lQc0NUhfGNTtuG0kMfoN4PTbh/Ohu+YHm5u3HXfzBCGjrz8vuHWyHvCwLdezDZ9Lezxu1bVhuc0n6IOz8svuqd2WTc48N0t2TNLw0mjFzl+a5atDOG4xsWVYVaWN3MgX+DOPcKazoB2Ld91zrnhQ1tWPV/AtR0BnVFkIjAOGv+LvDDe59AXTH4qDugRedYaCTo4j9GX22+33xKOzi+evuKQgUakpn8jf+qYHRlOb8Uqf6bYrNr5zeeD94SFw6OmdloWzsgvfj+459ZRi6xtHNy/7SLLm3lic9L5YWlnQLuW7zrn3Na05YYiX/mNHQGdWWQigL5RbUCH9wp3Na9cm8doUZi3ILdvmNG+f8NtK+aGsPdDjavT298SvSEsaVftgbBPI51nh+9mY03dZlZ74m/WbB61yOas8yLLm9l6Brs2zOwMaOfy3eecI6AAxlRtQJ8O4ffNK3fnMdozjJi8fcjwz+aGIxpPIEN4tj1yRrtqwwvDvdnWWYNPZy8xtWnjyENkYy5iLkK4b2TWcx0B7Vy++5ybftuy/skCruoI6LQiE4Fx8GKWPTPe59AXTH+qDei2ut2bx2hh6PgO7PFLf9u+tjoMPL+9ffeGqSO1uyycl/0yLMtv7ZrapdG7jSPXx1rEBrT1FPWpEJ7pCGjn8t3n3Ik3keAWbyKJTH8q/h7ojPbL4RvyGB0bftE8eP7xp7JsbvhRe8yGZrGmtdP243DISO0ea7yGXx5+mt/aNbVL4xX32uaVNT/bMNYiNqDXNwesCtM6X8J3Ld91zp0IKNwioCLTn4oDekzrTZvhJXmMVrZ/5OiscHGWnRHmtl96/yTsm+Vv3TR7M3xYOGtb9JaE1VNnb8mvdU3tdkQ4O7/YPH1g45iLmIAuab7rvzQc2RnQruW7zrkTAYVbBFRk8lNxQBsvg8/dlD13RvMjQS/MDee9kG35wcDQY42nl1PCPlc/sWX46aumhiuy/DuOA5dtzj+BNPTEtuhdFQ4MFzTX6ZqaXXnllu0PcXcY/N6WbNPycNTYi5iAhjOey55bEcLdnQHtWr7rnDsRULhFQEWmcFV/kP68EAbmDoSLmzH6+eQwuN+eYeCW/J5fz2sUatLkxj+WN58Unh/C7vMGwtB122qXPT0phIdb63RNDeGFjoe4tDFp/lCY/sTYi5iALmuez8AlWdcH6buW7zrnDgQUbhFQkQlc1QEdvvX4fYcOvmPbD0zO2n2/Lz3UumvL9UsXTtlr8alr2kPvOGHB0GdObb61NPJj6MeEw0YW6pzaHdDsrmX7DS08+5mXWMQE9OHrFk2ee9ydrWU6f5Rz+/Ld57wdAYVbBFRkAteXv5F+c/OTnb1oNrNXBBRuEVCRiUJfBvTRPXqdSUCBsRBQkYlCXwZ02fJeZxJQYCwEVGSi0JcBvbHXV/AEFBgTARWZKPRlQHtHQIGxEFCRiUJkAf2DEFC4RUBFJgoEVEdA4RYBFZkoEFAdAYVbBFRkokBAdQQUbhFQkYkCAdURULhFQEUmCgRUR0DhFgEVmSgQUB0BhVsEVGSiQEB1BBRuEVCRiQIB1RFQuEVARSYKBFRHQOEWARWZKBBQHQGFWwRUZKJAQHUEFG4RUJGJAgHVEVC4RUBFJgoEVEdA4RYBFZkoEFBdDwFd+fdLniGgmPgIqMhEgYDqigf0hQ+n6XcIKCY+AioyUSCguuIB/WWapocQUEx8BFRkokBAdcUDenkjoP9IQDHxEVCRiQIB1RUP6PJGQHchoJj4CKjIRIGA6ooH9PC0iYBioiOgIhMFAqorHtB5BBT9gYCKTBQIqK54QCcTUPQHAioyUSCguuIB/RgBRX8goCITBQKqKxzQzbvkbyERUEx8BFRkokBAdYUD+rs03Y2Aoh8QUJGJAgHVFQ7o2jQ9YFcCij5AQEUmCgRUVzigv0jTw/+BgKIPEFCRiQIB1RUO6E1pevwMAoo+QEBFJgoEVFc4oFel6Vc+S0DRBwioyESBgOoKB/RbaXruiQQUfYCAikwUCKiucEDPTdNLzyOg6AMEVGSiQEB1hQO6PE2/90MCij5AQEUmCgRUVzigx6fpdasIKPoAARWZKBBQXeGAHp6mtz32PgKKiY+AikwUCKiucEAPTNNV6/ckoJj4CKjIRIGA6goHdJ80XbP+aAKKiY+AikwUCKiucEAbTz4fWX9Dms6rajOBkhBQkYkCAdUVDuhA+neNa5cc+i8V7SVQFgIqMlEgoLrCAf14+tGqthEoFQEVmSgQUF3hgL4v3a2qbQRKRUBFJgoEVFc0oKem6eSqthEoFQEVmSgQUF3RgB6bpjOq2kagVARUZKJAQHVFA3p4mu5b1TYCpSKgIhMFAqorGtBFaXpQVdsIlIqAikwUCKiuaEDnpennq9pGoFQEVGSiQEB1RQM6I02PrmobgVIRUJGJAgHVFQ3o5DQ9saptBEpFQEUmCgRUVzSgn0rT06vaRqBUBFRkokBAdUUD+uE0/T9VbSNQKgIqMlEgoLqiAd01TS+oahuBUhFQkYkCAdUVDWjDt6vaRqBUBFRkokBAdT0E9HtVbSNQKgIqMlEgoLoeAnpNVdsIlIqAikwUCKiuh4DeXNU2AqUioCITBQKq6yGg/Cpl9AcCKjJRIKC6HgL666q2ESgVARWZKBBQXQ8BXVvVNgKlIqAiEwUCqushoI9UtY1AqQioyESBgOqKB/S9Ve0iUC4CKjJRIKC64gHlb8qhTxBQkYkCAdUVD+inKtpEoGQEVGSiQEB1xQPK35RDnyCgIhMFAqorHtDpVe0iUC4CKjJRIKC64gHlb8qhTxBQkYkCAdUVD+iBVe0iUC4CKjJRIKC64gHlb8qhTxBQkYkCAdUVD+gxVe0iUC4CKjJRIKC64gE9qapdBMpFQEUmCgRUVzygX61qF4FyEVCRiQIB1RUP6LlV7SJQLgIqMlEgoLriAf1mVbsIlIuAikwUCKiueEBXVrWLQLkIqMhEgYDqigf0uqp2ESgXARWZKBBQXfGA3lHVLgLlIqAiEwUCqise0F9VtYtAuQioyESBgOqKB/SBqnYRKBcBFZkoEFBd4YDu+kRVuwiUi4CKTBQIqK5wQD9W1SYCJSOgIhMFAqorHNDBqjYRKBkBFZkoEFBd4YDOrmoTgZIRUJGJAgHVFQ7o4qo2ESgZARWZKBBQXeGAHlvVJgIlI6AiEwUCqisc0HOq2kSgZARUZKJAQHWFAnprI6A/qmoTgZIRUJGJAgHVFQroY59MP7quqk0ESkZARSYKBFRXKKDr/+WU2yraQqB0BFRkokBAdcUCCvQRAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3Io2oA+uLTbeRIGA6ggo3Io1oKs+8t7rC00wUSCgOgIKt2IN6Llp+k+FJpgoEFAdAYVbsQb0lDRdXGiCiQIB1RFQuBVrQI9J0/0KTTBRIKA6Agq3Yg3ooWk6s9AEEwUCqiOgcCvWgC5K008XmmCiQEB1BBRuxRrQ/dN0sNAEEwUCqiOgcCvWgO6bprsVmmCiQEB1BBRuxRrQmWn6iUITTBQIqI6Awq1YA7pXmn600AQTBQKqI6BwK9aADqXp+wtNMFEgoDoCCrdiDeikNH1PoQkmCgRUR0DhVqwB3S1Ndy00wUSBgOoIKNyKNaD/m4DWhoDCrVgD+jECWhsCCrdiDeiHCGhtCCjcijWg7yOgtSGgcCvWgL6HgNaGgMKtSAP6eEpAa0NA4VakAX2EgNaHgMKtSAO6joDWh4DCrUgDej8BrQ8BhVuRBnQNAa0PAYVbkQb0bgJaHwIKtyIN6K8IaH0IKNyKNKCrCGh9CCjcijSgPyeg9SGgcCvSgP6UgNaHgMKtSAN6OwGtDwGFW5EG9BYCWh8CCrciDehNBLQ+BBRuRRrQ6wlofQgo3Io0oNcQ0PoQULgVaUCvJqD1IaBwK9KA/oCA1oeAwq1IA7qSgNaHgMKtSAP6XQJaHwIKtyIN6HcIaH0IKNyKNKDfIqD1IaBwK9KAXkJA60NA4VakAb2IgNaHgMKtSAN6AQGtDwGFW5EG9FwCWh8CCrciDeg5BLQ+BBRuRRrQFQS0PgQUbkUa0DMJaH0IKNyKNKBfJaD1IaBwK9KAnkZA60NA4VakAf0SAa0PAYVbkQb0ZAJaHwIKtyIN6AkEtD4EFG5FGtDjCGh9CCjcijSgSwlofQgo3Io0oEcR0PoQULgVaUC/SEDrQ0DhVqQBPZyA1oeAwq1IA/p5AlofAgq3Ig3oIQS0PgQUbkUa0MUEtD4EFG5FGtCDCGh9CCjcijSgCwlofQgo3Io0oPMJaH0IKNyKNKDzCGh9CCjcijSgcwlofQgo3Io0oLMJaH0IKNyKNKAzCWh9CCjcijSgexHQ+hBQuBVpQD9NQOtDQOFWpAEdIqD1IaBwK9KA7k5A60NA4VakAQ0EtD4EFG5FGtDdCGh9CCjcijSgnySg9SGgcCvSgH6CgNaHgMKtSAP6MQJaHwIKtyIN6IcJaH0IKNyKNKAfIKD1IaBwK9KAvpeA1oeAwq1IA7orAa0PAYVbcQb08ZSA1oeAwq0YA/roZTcR0BoRULgVY0BPTncloDUioHArxoAuTNNdCGh9CCjcijSgPAOtEQGFWwRUZKJAQHUEFG4RUJGJAgHVEVC4RUBFJgoEVEdA4RYBFZkoEFAdAYVbBFRkokBAdQQUbhFQkYkCAdURULhFQEUmCgRUR0DhFgEVmSgQUB0BhVsEVGSiQEB1BBRuEVCRiQIB1RFQuEVARSYKBFRHQOEWARWZKBBQHQGFWwRUZKJAQHUEFG4RUJGJAgHVEVC4RUBFJgoEVEdA4RYBFZkoEFAdAYVbBFRkokBAdQQUbhFQkYkCAdURULhFQEUmCgRUR0DhFgEVmSgQUB0BhVsEVGSiQEB1BBRuEVCRiQIB1RFQuBVtQMJb+/MAABl8SURBVD9AQOtCQOFWrAHd5R8JaF0IKNyKNaCDuxPQuhBQuBVpQBetIqC1IaBwK9KAPpgR0NoQULhFQEUmCgRUR0DhFgEVmSgQUB0BhVsEVGSiQEB1BBRuEVCRiQIB1RFQuEVARSYKBFRHQOEWARWZKBBQHQGFWwRUZKJAQHUEFG4RUJGJAgHVEVC4RUBFJgoEVEdA4RYBFZkoEFAdAYVbBFRkokBAdQQUbhFQkYkCAdURULhFQEUmCgRUR0DhFgEVmSgQUB0BhVsEVGSiQEB1BBRuEVCRiQIB1RFQuEVARSYKBFRHQOEWARWZKBBQHQGFWwRUZKJAQHUEFG4RUJGJAgHVEVC4RUBFJgoEVEdA4RYBFZkoEFAdAYVbBFRkokBAdQQUbhFQkYkCAdURULhFQEUmCgRUR0DhFgEVmSgQUB0BhVsEVGSiQEB1BBRuEVCRiQIB1RFQuEVARSYKBFRHQOEWARWZKBBQHQGFWwRUZKJAQHUEFG4RUJGJAgHVEVC4RUBFJgoEVEdA4RYBFZkoEFAdAYVbBFRkokBAdQQUbhFQkYkCAdURULhFQEUmCgRUR0DhFgEVmSgQUB0BhVsEVGSiQEB1BBRuEVCRiQIB1RFQuEVARSYKBFRHQOEWARWZKBBQHQGFWwRUZKJAQHUEFG4RUJGJAgHVEVC4RUBFJgoEVEdA4RYBFZkoEFAdAYVbBFRkokBAdQQUbhFQkYkCAdURULhFQEUmCgRUR0DhFgEVmSgQUB0BhVsEVGSiQEB1BBRuEVCRiQIB1RFQuEVARSYKBFRHQOEWARWZKBBQHQGFWwRUZKJAQHUEFG4RUJGJAgHVEVC4RUBFJgoEVEdA4RYBFZko/IEBfSSE0Tfec/KBU2YuWfG4OHx8vfQpjbqHgMItAioylSg/oMNnhjBpn2khDF6qDO9FCI+UtgwBBQiozFSi/IB+P0y5+sXGbpwTwq3C8F4QUKBcBFRkKlF+QOeHG1pXzgzHCcN7UVJAr7xyAwEFcgRUZCpRekA3hrChdW1VmPMfD+9JSQFtIqAAAZWZSvQc0FXH7zN08LUPtyKz7tTPDM05/uf51Q0h3NMasXVDq6SrT1ywx0Gn3j3SpLPC15s3rwtTn++eOqpfWy7ee2jBJZvMKqeGXMdbVPefcsDQAWc/NcaDbblo0ZRF52x88ssLhva7qPGfSCO9D1+zeGifZauyrPslfNdpdH1t2xFQuEVARSZTvQb0koEwMGcgnNiMzPcHw5QFU0I4N79nUZh93aaOkZcOhMF9G6OvbjUt+1WYM5zfflE4007t0hh7xqIrLp4aTjarXHNmCCef+fttA38wGCbPGwwznxz9YEsPW3ne5HDwjEPzixXNZh4bjvrG0kkD3+wOaNdpdH1tBBQxIKCicgK6OoRvbMo2nBbyyPxmcNrtw9nwHdPDzY277psZwtCRl9833Bp5Txj41ovZpq+FPX7XqtXwnOZT1OF5+UX3VBvQxY0Q3xYGnjOrdL+EXzc48N0t2TNLw0mjH+z4rVm2MoTjGhdXhlnNmQP5A925R1jTGdCu0+j62nJbB1pu3QI41fjfdbxPoXaL2wEtNKmcgB6R56pRwYPzyHy5/Xb7LeHo/OLpKw4ZaMRn+jfyp4TZkeH05shD82eAzVqd33yed09YODxqqglo3rpNzcx1rdId0GXhjPzi94N7bh31YGsbB/dvu2jOPLE56fywtDOgXafR9bU1A5q23NDbvysAE9HB7YD+IWv0FtDhvcJdzSvX5pFZFOYtyO0bZrTv33Dbirkh7P1Q4+r09rdEbwhL2rV6IOzTSOfZ4bvZWFO3aYxtPl1uZq5rle6Azmrf9Zs1m0c92Oas86I5s/VMd22Y2RnQztPo/tpyBBRwaPwC+nQIre9B3p1HZs8wYvL2IcM/mxuOaDwxDOHZ9sgZ7VoNLwz3ZltnDT6dvcTUlsbYrfllnrnuVboCunHkVLIxH8xchHDfyKznOgLaeRrdX1snvgcKt/geqMhkqreAbqvWvXlkFoaOVY9f+tv2tdVh4PntTbs3TB2p2GXhvOyXYVl+a9fUbiOfIuoKaHOVroA2erdx5PpYD2YD2nqK+lQIz3QEtPM0ur+2TgQUbhFQkclUj98DndF+mXtDHpljwy+aB88//lSWzQ0/ao/Z0CzRtHayfhwOGanYY43X8MvDT/Nbu6Z26wxo9ypdAW284l7bvLLmZxvGejAb0OubA1aFaZ0v4btOo+tr60RA4RYBFZlM9RjQY1pvxgwvySOzsv0jR2eFi7PsjDC3/ZL6J2HfLH9LZnlz5GHhrG0xWxJWT53dfDura2q3roB2rdL9PdAjwtn5xebpAxvHfDAT0CXNTwcsDUd2BrTrNLq+tk4EFG4RUJHJVI8Bbby8PXdT9twZzY/6vDA3nPdCtuUHA0OPNZ5eTgn7XP3EluGnr5oarsjy7yQOXLY5/2TR0BPbYnZVODBc0Fyna2p25ZUdHxLoCmjXKo2b7t4+7u4w+L0t2abl4aixH8wENJzxXPbciuYK2wPadRpdXxsBRRQIqKicgGbnhTAwdyBc3IzMzyeHwf32DAO35Pf8el6jPJMmN/6xvPlk7/wQdp83EIau2x6zpyeF8HBrna6pIbzwEgHtWqXxInvmF57cNvDSxs3zh8L0J8Z+MBPQZc3zHrgk6/ogfddpdH1tBBQxIKCikgI6fOvx+w4dfMe2H4Sctft+X3qoddeW65cunLLX4lPXtIfeccKCoc+c+tvOKB4TDhtZqHPqvxPQzlWy2xfsPuvfto+8a9l+QwvPfuYlHswE9OHrFk2ee9yd21Ye4yvo/tq2I6Bwi4CKTAgn1G+k39z8xGaV/qBfQ0JA4RYBFZkoTKiAPrpH1Y9AQIGxEFCRicKECuiy5VU/AgEFxkJARSYKEyqgN1b9Cp6AAmMioCIThQkV0OoRUGAsBFRkohBZQP8gBBRuEVCRiQIB1RFQuEVARSYKBFRHQOEWARWZKBBQHQGFWwRUZKJAQHUEFG4RUJGJAgHVEVC4RUBFJgoEVEdA4RYBFZkoEFAdAYVbBFRkokBAdQQUbhFQkYkCAdURULhFQEUmCgRUR0DhFgEVmSgQUB0BhVsEVGSiQEB1BBRuEVCRiQIB1RFQuEVARSYKBFRHQOEWARWZKBBQHQGFWwRUZKJAQHUEFG4RUJGJAgHVEVC4RUBFJgoEVEdA4RYBFZkoEFAdAYVbBFRkokBAdQQUbhFQkYkCAdURULhFQEUmCgRUR0DhFgEVmSgQUB0BhVsEVGSiQEB1BBRuEVCRiQIB1RFQuEVARSYKBFRHQOEWARWZKBBQHQGFWwRUZKJAQHUEFG4RUJGJAgHVEVC4RUBFJgoEVEdA4RYBFZkoEFAdAYVbBFRkokBAdQQUbhFQkYkCAdURULhFQEUmCgRUR0DhFgEVmSgQUB0BhVsEVGSiQEB1BBRuEVCRiQIB1RFQuEVARSYKBFRHQOEWARWZKBBQHQGFWwRUZKJAQHUEFG4RUJGJAgHVEVC4RUBFJgoEVEdA4RYBFZkoEFAdAYVbBFRkokBAdQQUbhFQkYkCAdURULhFQEUmCgRUR0DhFgEVmSgQUB0BhVsEVGSiQEB1BBRuEVCRiQIB1RFQuEVARSYKBFRHQOEWARWZKBBQHQGFWwRUZKJAQHUEFG4RUJGJAgHVEVC4RUBFJgoEVEdA4VakAT3pYQJaGwIKtyINaDqDgNaGgMKtWAO6ayCgdSGgcCvWgKYfIqB1IaBwK9qA7kJA60JA4Va0AU0JaF0IKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMKtGAP6xfyvwhPQ2hBQuBVjQNeeeDHPQGtEQOFWjAFdv/5xAlojAgq34gzo+l0JaH0IKNyKNKB/R0DrQ0DhVqQBfT8BrQ8BhVuRBvTDBLQ+BBRuRRrQjxHQ+hBQuBVpQD9OQOtDQOFWpAH9JAGtDwGFW5EGdDcCWh8CCrciDWggoPUhoHAr0oAOEtD6EFC4FWlAhwhofQgo3Io0oFMIaH0IKNyKNKB7EdD6EFC4FWlAZxDQ+hBQuBVpQGcR0PoQULgVaUDnEtD6EFC4FWlA5xHQ+hBQuBVpQOcT0PoQULgVaUAPIKD1IaBwK9KAHkhA60NA4VakAV1EQOtDQOFWpAH9LAGtDwGFW5EG9PMEtD4EFG5FGtDDCGh9CCjcijSgXySg9SGgcCvSgB5FQOtDQOFWpAFdSkDrQ0DhVqQBPY6A1oeAwq1IA3oCAa0PAYVbkQb0JAJaHwIKtyIN6CkEtD4EFG5FGtDTCGh9CCjcijSgXyWg9SGgcCvSgJ5BQOtDQOFWpAH9GgGtDwGFW5EG9GwCWh8CCrciDei5BLQ+BBRuRRrQ8wlofQgo3Io0oBcS0PoQULgVaUAvIaD1IaBwK9KAXkpA60NA4VakAf0OAa0PAYVbkQb0nwlofQgo3Io0oCsJaH0IKNyKNKBXEdD6EFC4FWlAryag9SGgcCvSgP6IgNaHgMKtSAN6PQGtDwGFW5EG9CYCWh8CCrciDejNBLQ+BBRuRRrQ2whofQgo3Io0oD8loPUhoHAr0oD+jIDWh4DCrUgDuoqA1oeAwq1IA/orAlofAgq3Ig3oagJaHwIKtyIN6BoCWh8CCrciDej9BLQ+BBRuRRrQdQS0PgQUbkUa0EcIaH0IKNyKNKCPE9D6EFC4FWlA1+9KQGtDQOFWrAF9LwGtDQGFW7EG9IMEtDYEFG7FGtCPEdDaEFC4FWtA/56A1oaAwq1YA7obAa0NAYVbsQZ0IE3fU2iCiQIB1RFQuBVrQIfS9H2FJpgoEFAdAYVbsQZ0zzT9SKEJJgoEVEdA4VasAZ2Rph8vNMFEgYDqCCjcijWgc9P0U4UmmCgQUB0BhVuxBnRBmk4qNMFEgYDqCCjcijWgB6XplEITTBQIqI6Awq1YA/r5NJ1ZaIKJAgHVEVC4FWtAj0rTeYUmmCgQUB0BhVuxBvTENF1UaIKJAgHVEVC4FWtAV6TpUYUmmCgQUB0BhVuxBvSWXdMrCk0wUSCgOgIKt2IN6Po7flxsvIkCAdURULgVbUCLMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awi0CKjJRIKA6Agq3CKjIRIGA6ggo3CKgIhMFAqojoHCLgIpMFAiojoDCLQIqMlEgoDoCCrcIqMhEgYDqCCjcIqAiEwUCqiOgcIuAikwUCKiOgMItAioyUSCgOgIKtwioyESBgOoIKNwioCITBQKqI6Bwi4CKTBQIqI6Awq1Vt9zyyHifQ18wUSCgALKD0/TB8T6HfkRAARDQHhFQAAS0RwQUAAHtEQEFQEB7REABENAeEVAABLRHBBQAekRAAaBHBBQAekRAAaBHBBQAekRAAaBHBBQAekRAgbj9c3ikfe2x02cPzj7t0XE9mz5DQIGovbD/SEDv+3QIs0OYsmZ8T6ivEFAgYltWHx3aAd2yIBz5VPbkEWH+lnE+qT5CQIF4XTAYwkhAbwkzNjUunpsebh3fk+onBBSI17Vf+cpXRgJ6ejizeXlGWD6ep9RfCCgQt5GAHhZ+0ry8KXxhPE+nvxBQIG4jAZ0f7mxergr7j+fp9BcCCsRtJKBTwgPNywfC1HE8mz5DQIG4bQ/ob5qX94eh8Tyd/kJAgbiNfgm/33ieTn8hoEDctr+JdHPz8ifh8PE8nf5CQIG4jQT0tLCiebkinD6ep9NfCCgQt5GA3hxmv9i42DybD9LrCCgQt20/yjk/nLI523xyWMCPcsoIKBC3kYBma6aEKYfsEaauHd/z6SsEFIjbtoBmj542a3DW6Y+P69n0GQIKAD0ioADQIwIKAD0ioADQIwIKAD0ioADQIwIKAD0ioADQIwIKAD0ioADQIwIKAD0ioADQIwIKLz6edHtr70s156flndro9f9y1G2vyR9z9M2Y0AgovBgroD0miYBCQ0DhxcqvNf1pkrSuXJT9AQF920MPPVby+XWtP/q0fvvQQwS07xBQOPOXScd/1D0HtNqSvcT6BLTvEFA40xXQZ5/d0MsaBBQaAgpnugLaIwIKDQGFM10BbR/kZbr8I294xVv2XpM9Ou/PdnzrwO3tAfd85q9e8/q/2f8X3WtsK9m6z73/za98+yfOf/Elht974F+99g27HDTyNzFu2/ttr97p7XNu377KPXP/4hX/+YNnbWzdcs2n37zjn01bNbK+XZ2A9h0CCmfGDug7liR//sE3JMlOF775Zf/rvyXJy36c3z785R3ab9l/dnPnGiMlO+vl7bvf+fRYw4e/1D7e6cr8cPMBIx8A2P/F9ir/9zUvf/e7d0ySdz3TOH5x79a9r7yytb5ZnYD2IQIKZ8YO6MtefXGWPfWRRq3edEc2/Lkk+UB++zlJ8sZDL7lw8U5JMr9zjXbJfvSy5FUHXLzytP+eJNPGGn5mkuyw33mXLnpV8ob8L1ke0jicf/4FCxpdPKS1ys6v36NRzvvemSSfaxwf2DiNWees2CN5bXN9uzoB7UMEFM6MHdDk5Pzy9saVcxuXz+6QvK5x8eRrk3etz29/6L8myY871miXrPGM8bv55VP/KXnd1tHDH31V8prmrFt2zAN5V5K85sb88CeNQK5uPeqHhvMbrkuS92TZnUnyih/kh+e0PvBpVs8IaB8ioHDmJQL6bH65MUn+aFN+5S3N25cnyT2tcbd1PwVtlyxNkn9rHl92wgm/Hz18aZJ8vnU8NdklyxYnyeGtw0NbT0Ebj3p91n7YxnqfSZKDWvd/snlsVs8IaB8ioHBm7IC+pXU88gOerdsnJX/bHrd15+R/dKzRLtleSTL5X7ffaod/MEkeaB0/eMNNWfa3SfJg6/CBJA9q/o2DTR3r/U2S3Nc6XNk8NqtnBLQPEVA481LvwmedV1q3v6vzJz//tGON9rB739C4/V1Lvr2udasd/ifJK7Z2THpj8sr24dYdkzc1V/nzzvVev234muaxWT0joH2IgMKZAgF9e2cRX9mxxsiwRw744+Z9f/1Pz44xfIeR57UtO2z/9SVvSXbofNTWte3Dn2vd0716RkD7EAGFMwUC+u7kU2Ovsb1kW/7fSUNvbkTuLx4dPfwNyU7DHYcdz0BfmbwxGxXQnZNXbGkdrh25p3P1jID2IQIKZwoE9NPJX4y9RnfJhm+ZlCT7jh7+P5Ok/Qn6dRdeuKHje6APNl6aZ6MCukuSrGkd/rBz/ZHVRz0s+gABhTMFAnpWktzYuv22t771qx1rtIdNf+ec1vHvmr/czg5fnCTHtI4/n+y8NVuUJF9sHR7Wer/dBPSQbe/CT2oem9UzAtqHCCicKRDQp1+bvOOh/PDJdyd/9GDHGu1hc5Id72oefztJPjl6+OqXJa+7LT9e88Zkcpb9Kklee0t+eONOSfLrbFRA7xr5HOglrc+BmtUzAtqHCCicKRDQ7OxG85ZcfMWxb0mSIzvXaA+7MEn+5IuXXXX+7JcnyeVjDP9Ckrx8//MvP+qPk53yFC7ODy+4YP4OrR88sgHNDs5/Eunr50xL3tE8tqsT0D5EQOFMkYAOn/Cy1nvqf7Sk8+2gkWFb9hp5z/3lJ4w1fMvn2ne/7tv54Yv7jQxfuLnrUdvXXpzTuvdNdzaP7eoEtA8RUDhTJKBZ9q/z3v6qnXedf2/3GtvG/2TKX+/8qrd9+LhHs7GH/2zft7965/cuGvnl9TfP+i+vevXb9r7drLLt2jVT3rzjn85YO3JsVyegfYeAAqOMU8kIaN8hoMAoBBQaAgqMQkChIaDAKAQUGgIKjNJ8b7zCvws/Fv4ufD8ioMAoBBQaAgoAPSKgANAjAgoAPSKgANAjAgoAPSKgANAjAgoAPSKgANAjAgoAPSKgANAjAgoAPSKgANAjAgoAPSKgANCj/w8b9FliMAeATgAAAABJRU5ErkJggg==" title alt width="672" /></p> +<p>Here we get a performance benefit of a factor of 15.8 using the version of the differential equation model compiled from C code using the inline package!</p> <p>This vignette was built with mkin 0.9.42 on</p> <pre><code>## R version 3.2.4 Revised (2016-03-16 r70336) ## Platform: x86_64-pc-linux-gnu (64-bit) |