diff options
Diffstat (limited to 'vignettes/compiled_models.html')
| -rw-r--r-- | vignettes/compiled_models.html | 38 |
1 files changed, 19 insertions, 19 deletions
diff --git a/vignettes/compiled_models.html b/vignettes/compiled_models.html index c7f4fbea..bd39ae2d 100644 --- a/vignettes/compiled_models.html +++ b/vignettes/compiled_models.html @@ -10,7 +10,7 @@ <meta name="author" content="Johannes Ranke" /> -<meta name="date" content="2015-11-13" /> +<meta name="date" content="2015-12-09" /> <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>2015-11-13</em></h4> +<h4 class="date"><em>2015-12-09</em></h4> </div> <div id="TOC"> @@ -104,21 +104,21 @@ smb.1 <- summary(mb.1) print(mb.1)</code></pre> <pre><code>## Unit: milliseconds ## expr min lq mean median uq -## deSolve, not compiled 9538.4007 9570.3211 9605.6503 9602.2416 9639.2752 -## Eigenvalue based 881.9438 885.9337 901.1558 889.9236 910.7618 -## deSolve, compiled 692.0913 695.6109 697.9629 699.1304 700.8987 +## deSolve, not compiled 9508.4631 9522.5843 9634.9196 9536.7055 9698.1479 +## Eigenvalue based 872.6560 877.4544 888.3598 882.2527 896.2117 +## deSolve, compiled 698.8148 700.5031 708.8625 702.1914 713.8864 ## max neval cld -## 9676.3087 3 c -## 931.5999 3 b -## 702.6669 3 a</code></pre> +## 9859.5902 3 b +## 910.1707 3 a +## 725.5815 3 a</code></pre> <pre class="r"><code>autoplot(mb.1)</code></pre> -<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAABUAAAAPACAMAAADDuCPrAAAC+lBMVEUAAAABAQECAgIDAwMEBAQFBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQVFRUWFhYXFxcYGBgZGRkaGhobGxscHBwdHR0eHh4fHx8gICAhISEiIiIjIyMkJCQlJSUmJiYnJycoKCgpKSkqKiorKyssLCwtLS0uLi4wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6Ojo7Ozs8PDw9PT0+Pj4/Pz9AQEBBQUFCQkJDQ0NERERFRUVGRkZHR0dISEhJSUlKSkpLS0tMTExNTU1OTk5PT09QUFBRUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVlbW1tcXFxdXV1eXl5fX19gYGBhYWFiYmJjY2NkZGRlZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1ubm5vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8fHx9fX1+fn5/f3+AgICBgYGCgoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uMjIyNjY2Ojo6Pj4+QkJCRkZGSkpKTk5OUlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2enp6fn5+goKChoaGioqKjo6OkpKSlpaWmpqanp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKzs7O0tLS1tbW2tra3t7e4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXGxsbHx8fIyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY2NjZ2dna2trb29vc3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr6+vs7Ozt7e3u7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+/v7////Lq4d6AAAACXBIWXMAAB2HAAAdhwGP5fFlAAAgAElEQVR4nO3dCZgc9Xnn8cr6jpPYTuJkc9mb3STOvUkZ7NiwiWPn8JHTiZN4N3GcOOsjdu5kvVnnP5qRhAQIAQIkIS4JcVnGWAgTEIfBHOKwIMgSSEgIJCGBJUD3OZrp59muru6ZqZ4eTVWNqN/7dn8/z7Oeru66mu73uz3dPUpUAwCUEqlPAAC8IqAAUBIBBYCSCCgAlERAAaAkAgoAJRFQACiJgAJASQQUAEoioABQEgEFgJIIKACUREABoCQCCgAlEVAAKImAAkBJBBQASiKgAFASAQWAkggoAJREQAGgJAIKACURUAAoiYACQEkEFABKIqAAUBIBBYCSCCgAlERAAaAkAgoAJRFQACiJgAJASQQUAEoioABQEgEFgJIIKACUREABoCQCerIc2O3CkaNH1KfQ8/buUZ9Br9tz9OihydbJN/YE9GR5cbsLQ7Uh9Sn0vD1Onivd67la7chk6+QbewJ6sjgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegRUjYDa42QoCKgeAVUjoPY4GQoCqkdA1QioPU6GgoDqEVA1AmqPk6EgoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9jgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegRUjYDa42QoCKgeAVUjoPY4GQoCqkdA1QioPU6GgoDqEVA1AmqPk6EgoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9jgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegRUjYDa42QoCKgeAVUjoPY4GQoCqkdA1QioPU6GgoDqEVA1AmqPk6EgoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNpTcCi2XvHvZR/9KSGgegRUjYDaU3Ao5sZvv6vswz8VBFSPgKoRUHsKDsXH4nhe2Yd/KgioHgFVI6D2FByKD8bxv5V9+KeCgOoRUDUCak/BoTgtjv+27MM/FQRUj4CqEVB7ig3FM3Ec/2Xpx38KCKgeAVUjoPYUG4p19YD+SenHfwoIqB4BVSOg9hQbim/UA/q7pR//KSCgegRUjYDaU2wovl4P6G+UfvyngIDqEVA1AmpPsaG4tR7QXy79+E8BAdUjoGoE1J5iQ7G8HtB3lH78p4CA6hFQNQJqT7GhuK4e0Hhr6SdAeQRUj4CqEVB7ig3F4iSgG0s/AcojoHoEVI2A2lNsKBYmAV1b+glQHgHVI6BqBNSeYkNxfhLQR0o/AcojoHoEVI2A2lNsKM4ioL2LgKoRUHuKDcV0Atq7CKgaAbWn2FB8gYD2LgKqRkDtKTYU/0xAexcBVSOg9hQbir8loL2LgKoRUHuKDcWnCGjvIqBqBNSeYkPxcQLauwioGgG1p9hQ/E8C2rsIqBoBtafYUHyEgPYuAqpGQO0pNhS/S0B7FwFVI6D2FBuKDxDQ3kVA1QioPcWG4n0EtHcRUDUCak+xoTidgPYuAqpGQO0pNhSnEtDeRUDVCKg9hYZia0xAexcBVSOg9hQaio0EtIcRUDUCak+hoVhHQHsYAVUjoPYUGopHCGgPI6BqBNSeQkOxioD2MAKqRkDtKTQUdxPQHkZA1QioPYWGYiUB7WEEVI2A2lNoKFYQ0B5GQNUIqD2FhmIZAe1hBFSNgNpTaCiWEtAeRkDVCKg9hYbiUgLawwioGgG1p9BQXERAexgBVSOg9hQairMJaA8joGoE1J4CQ7Ft2Sfj+O0EtFcRUDUCak+BoZiXvPw8nYD2KgKqRkDtKTAUn0gC+n4C2qsIqBoBtafAUPx5EtA/JKC9ioCqEVB7igb04wS0VxFQNQJqT8GAnvo5AtqrCKgaAbWnYEB/828IaK8ioGoE1J6CAf1zAtqzCKgaAbWnYEDnEdCeRUDVCKg9xQJ64WEC2rMIqBoBtadYQF+oEdCeRUDVCKg9BBQ5EVA1AmoPAUVOBFSNgNpDQJETAVUjoPYQUOREQNUIqD0EFDkRUDUCag8BRU4EVI2A2kNAkRMBVSOg9hBQ5ERA1QioPQQUORFQNQJqDwFFTgRUjYDaQ0CREwFVI6D2EFDkREDVCKg9BBQ5EVA1AmoPAUVOBFSNgNpDQJETAVUjoPYQUOREQNUIqD0EFDkRUDUCag8BRU4EVI2A2kNAkRMBVSOg9hBQ5ERA1QioPQQUORFQNQJqDwFFTgRUjYDaQ0CREwFVI6D2EFDkREDVCKg9BBQ5EVA1AmoPAUVOBFSNgNpDQJETAVUjoPYQUOREQNUIqD0EFDkRUDUCag8BRU4EVI2A2kNAkRMBVSOg9hBQ5ERA1QioPQQUORFQNQJqDwFFTgRUjYDaQ0CREwFVI6D2EFDkREDVCKg9BBQ5EVA1AmoPAUVOBFSNgNpDQJETAVUjoPYQUOREQNUIqD0EFDkRUDUCag8BRU4EVI2A2kNAkRMBVSOg9hBQ5ERA1QioPQQUORFQNQJqDwFFTgRUjYDaQ0CREwFVI6D2EFDkREDVCKg9BBQ5EVA1AmoPAUVOBFSNgNpDQJETAVUjoPYQUOREQNUIqD0EFDkRUDUCag8BRU4EVI2A2kNAkRMBVfMc0BB2dbj24O1XzJl+wdUPH+9w264QXuZzKnTEzneAgCIvAqrWdQF99IyQmrt5/I1VBjQ9FgHFy4iAqnVbQDeFsHjjvqM7H5obpr807lYCehIRUD0CqtZtAb0sfGm4ceHQuWH5uFurDOjeZcsmOyIBxdQQULUuC+jQjJHrVodzxm3Be6AnEQHVI6BqXRbQPSEcbl584fobxm1BQE8iAqpHQNVcBnTz0lkDC1YPN/sz+OAlZ0yfd0vjHc/BvnBf28ovrjhvYMZFtx1ILjdy9tWwKL1lTZhxNLt5uxdvnNt/zpWbhjvt6Phdc/vPXnFo+NElc6ZfcPexWqOH2649s3/ODTtHjjUS0OxRsndg3FHzP3wEtLcRUDWPAV2Vfsx+a9qf/QvSxRkbktuuC+Ha9cfGrLxuIL159rZaM2dbQ9jTuOmq5F3S7OZZG5rbruy0o+vOufn6vnDx0r6rVy4NIXm1Wz+lvjOuXDo79K+ttQU0e5TsHRiHgCInAqrmMKAv9YdFzxzednlo9Gf48nDu2oODz14TZuyu37hvfv3qgSX3PZ++aqztHAjzNx858PiccPbhZs6Gzw33JzcdnBa2tW+ePdCMcMX2wX23h/BUhx1ddKhWe6h+rGfqa94dpg0mAQ3XHKm/2FwRBl7IBjR7lOwdSB07d9SG47n972ZAt+bf5KSpn7XgqBhraEh9Bj2vPt2TrLE/X9cqC+iXwwX1WtWOz2/056kwu/F6cvi6cHPyc/A/Fvc3vgZ6b+N16A3h/GTl2v6Z4d5Wzu5If4d/OFw4PG7zsVaEBUPpAVd02NGT9YUDIXytcWXjVEKY01hl+JJwYzag2aNk70DqWBi1Lv9/i082A5rz/5MDULWD+VarLKBzwprGz8ca/flauCm9enO4qLnCsY23XlTv0MXJiZ/XXLl2e7i6lbOd6e/wlyUvRDtsPmJeWNv4+dytqzrsqP4CtDYcwpZa+rMR0HvTVdaH87MBzR4lewea50xAga5kLKD11LzYuPBioz/LRsMze8xaL9zUl7xuPN7XXLketXkj70guSNK5J/QfnHDzxNC01ra1iXY08jl6M6Ab01VeCtOGMgHNHKXtDozHe6DIifdA1fy9B7o7hPRDomON/lwRzpzXdGH9Rd1jh1rr3RcGhkZXru0I00e6d1/yO/y94Yu19s07Hyi7MHZH7QHdka5yNHmJOzagmaO03YHxCChyIqBq/gI68lrwpUZ/lqfvQjbNCo+3Lh5KKjb6wnFDOG+ke3v76jfND5tq7ZtnHOtrflqfOWpmR+0BfSJd5fnQN5gJaOYobXdgPAKKnAiomr+A1uY230Jc1+jPPeGa9Oq9G3ck72t+pbXarjB9aMxbl3eGq0a7d0W4f2c4J/mEKLt5+4HWN36+eP1NE+2oLaAr01VWjTY2/d/sUbJ3YDwCipwIqJrDgC4PFyZfoRla2OjP9tD3bHLt8OKwqvG9okfTLzANLw+X1JJPvOe1Pjy/Z7R7q8OiO9MXhdnNs74UFjX2dVuSzM47agvo9MZL1gOzk7dfxwY0e5TsHRiPgCInAqrmMKC7+8MlW45sXzKQ9ufqMPuR/UM7l4XZB+q/HtejdOmjW3bvWLs4hORL6zsHwoKnjxxMvr55aLR7h/vDWX3p9z4zm9eWLTsweqDn+8OS7YOHVvUl7wt03lFbQMNZ39y3b+3ZYeb+ti/SZ47SdgfGIaDIiYCqOQxo7YH08+zVaX8OLmz+kU/yF0K1g1e2Pu4eeLix8tr+dPHM5CXgyF+mXxvCkvRSdvMQxn6d/rFp6W3/PuGO2gJ6R19jlf51tbaAZo+SvQPjEFDkREDVPAY0+VPy/vnrWvUaXLVo5oz5tzRfOg5vWn7ZOQPnLbmt9cftL6w4d+CsxXeM/i184vHQ/I5n2+bZgNaeX75w+tylT028o7aA7tp+7expnf8WPnOS2TvQjoAiJwKq5jKgL5/jYz54L2yCHhZGQJETAVUjoBm7w7HJV5oIAUXFCKgaAc342qVT2JiAomIEVI2AZizr+M+C5kRAUTECqkZATx4CiooRUDUCevIQUFSMgKoRUHsIKHIioGoE1B4CipwIqBoBtYeAIicCqkZA7SGgyImAqhFQewgociKgagTUHgKKnAioGgG1h4AiJwKqRkDtIaDIiYCqEVB7CChyIqBqBNQeAoqcCKgaAbWHgCInAqpGQO0hoMiJgKoRUHsIKHIioGoE1B4CipwIqBoBtYeAIicCqkZA7SGgyImAqhFQewgociKgagTUHgKKnAioGgG1h4AiJwKqRkDtIaDIiYCqEVB7CChyIqBqBNQeAoqcCKgaAbWHgCInAqpGQO0hoMiJgKoRUHsIKHIioGoE1B4CipwIqBoBtYeAIicCqkZA7SGgyImAqhFQewgociKgagTUHgKKnAioGgG1h4AiJwKqRkDtIaDIiYCqEVB7CChyIqBqBNQeAoqcCKgaAbWHgCInAqpGQO0hoMiJgKoRUHsIKHIioGoE1B4CipwIqBoBtYeAIicCqkZA7SGgyImAqhFQewgociKgagTUHgKKnAioGgG1h4AiJwKqRkDtIaDIiYCqEVB7CChyIqBqBNQeAoqcCKgaAbWHgCInAqpGQO0hoMiJgKoRUHsIKHIioGoE1B4CipwIqBoBtadYQJcPE9CeRUDVCKg9xQIaLyGgPYuAqhFQewoG9K8IaM8ioGoE1J6CAf0gAe1ZBFSNgNpTMKC/9DkC2qsIqBoBtadgQONPENBeRUDVCKg9RQP6RwS0VxFQNQJqT4Gh+HgS0A8S0F5FQNUIqD0FhmJGEtBfJqC9ioCqEVB7CgzFU2f9SRyfQkB7FQFVI6D2FBqKs5PXoAS0RxFQNQJqT6GhuICA9jACqkZA7Sk0FJcS0B5GQNUIqD2FhuJKAtrDCKgaAbWn0FAsI6A9jICqEVB7Cg3FjQS0hxFQNQJqT6GhuJWA9jACqkZA7Sk0FHcR0B5GQNUIqD2FhmIVAe1hBFSNgNpTaChWE9AeRkDVCKg9hYZiLQHtYQRUjYDaU2goNhLQHkZA1QioPYWGYgsB7WEEVI2A2lNsKE4hoL2LgKoRUHuKDcVpBLR3EVA1AmpPsaF4LwHtXQRUjYDaU2wofpOA9i4CqkZA7Sk2FL9DQHsXAVUjoPYUG4o/JKC9i4CqEVB7ig3FRwlo7yKgagTUnmJD8XEC2rsIqBoBtafYUHySgPYuAqpGQO0pNhSfI6C9i4CqEVB7ig3FPxLQ3kVA1QioPcWG4l8JaO8ioGoE1J5iQ9FPQHsXAVUjoPYUG4rZBLR3EVA1AmpPsaE4l4D2LgKqRkDtKTYU85OAri39BCiPgOoRUDUCak+xobg8CeiTpZ8A5RFQPQKqRkDtKTYU1yQB3VL6CVAeAdUjoGoE1J5iQ3FDvZ+nln78p4CA6hFQNQJqT7GhuKUe0NNLP/5TQED1CKgaAbWn2FDcVQ/or5d+/KeAgOoRUDUCak+xoXioHtDfLv34TwEB1SOgagTUnmJDsbYe0D8u/fhPAQHVI6BqBNSeYkPxdD2gHy/9+E8BAdUjoGoE1J6CQ/HOOP5s2Yd/KgioHgFVI6D2FByK34zjfy378E8FAdUjoGoE1J6CQ/HROJ5b9uGfCgKqR0DVCKg9BYeiP45vLvvwTwUB1SOgagTUnoJD8eSMRWUf/SkhoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9jgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegRUjYDa42QoCKgeAVUjoPY4GQoCqkdA1QioPU6GgoDqEVA1AmqPk6EgoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9jgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegRUjYDa42QoCKgeAVUjoPY4GQoCqkdA1QioPU6GgoDqEVA1AmqPk6EgoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9jgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegRUjYDa42QoCKgeAVUjoPY4GQoCqkdA1QioPU6GgoDqEVA1AmqPk6EgoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9jgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegRUjYDa42QoCKgeAVUjoPY4GQoCqkdA1QioPU6GgoDqEVA1AmqPk6EgoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9jgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegRUjYDa42QoCKgeAVUjoPY4GQoCqkdA1QioPU6GgoDqEVA1AmqPk6EgoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9jgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegRUjYDa42QoCKgeAVUjoPY4GQoCqkdA1QioPU6GgoDqEVA1AmqPk6EgoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9jgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegRUjYDa42QoCKgeAVUjoPY4GQoCqkdA1QioPU6GgoDqEVA1AmqPk6EgoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9jgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegRUjYDa42QoCKgeAVUjoPY4GQoCqkdA1QioPU6GgoDqEVA1AmqPk6EgoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9miH4uZ/nLUx14oEVI+AqhFQe6RDseyUOP6Dp/OsSUD1CKgaAbVHORTb3hvXnZNnVQKqR0DVCKg9yqF4Io5/7ZT4l5/KsSoB1SOgagTUHnFAP/FPcbw4x6oEVI+AqhFQe9QBXRXHf5FjVQKqR0DVCKg96oAO/Vp86vrJVyWgegRUjYDaow5obUYcL518VQKqR0DVCKg98oB+PY7/bvJVCageAVUjoPbIA3rol+L3bJt0VQKqR0DVCKg98oDWPhXHd0+6KgHVI6BqBNQefUAXx/F5k65KQPUIqBoBtUcf0Mfj+OOTrkpA9QioGgG1Rx/QoV+N3znp38MTUD0CqkZA7dEHtPYvcXzzZKsSUD0CqkZA7TEQ0C/H8YzJViWgegRUjYDaYyCg2+P4w5OtSkD1CKgaAbXHQEBrH47fvqZ51aM3d/63mQioHgFVI6D2WAjo3Di+tHHF1mmnTPCJPAHVI6BqBNQeCwH9Rhz/ZeOKf4nj+H90XJWA6hFQNQJqj4WADr0vPnVdffnymIAaRkDVCKg9FgJamxXHF23fvvq0OH47ATWLgKoRUHtMBHRdHH9g69aPxvG//CYBNYuAqhFQe0wEtPbncTz3/8Xxb+whoHYRUDUCao+NgD6YvPsZn/JAjYDaRUDVCKg9NgJaO7Pez7dfUyOghhFQNQJqj5GADl/955++t0ZALSOgagTUHiMBbSGgdhFQNQJqDwFFTgRUzW9Awxh3n9Q97woh7ynsmtKBJtiegCInAqrWPQHNnb1JEdAMAmoXAVXzHNBn9rccrRHQk4OAukJA1TwHNFufvcuWnaQ9E9AMAmoXAVXrnoCePAQ0g4DaRUDVCOh4BDSDgNpFQNW6J6Ct7D33xTMHFj7RunnwwUvOmD7vlpdaa2y4dOYZFz8yVKt9NSxKt1sTZhyt1bYsO3/6zHk3PDu6p5GMdtjT6Ck8/8Uzpy/4yo50eexOagdvX3DGrIX3HUtvym69eemsgQWrhwkopoaAqnVdQNdOa3wuvzq9ef+C9GP6GRvSNR6advXKxSGsrNW2hrCnsd1VYXmtVl899NX/37T1tQkCmtnTyCms6m9c2fdIspjZyZ6z0vUvG6qN23pVunArAcXUEFC1bgvowZnhyucGn7uqv3Hz8OXh3LUHB5+9JszY3Vhj5ub6OneG/sHa8Lnh/mSzg9PCtuR/vvLi0LEtC8KCWueAZvc0cgp9Z67Zt/ebc8L0PbW2nVwTFj57/Ohj/eGxWvvWL/WHRc8c3nZ5GHMXjp07asNxnRc6BfRXOq5av6nik0O7oSH1GfS8+nRPssb+fEFTfg90b60VvFvCxclLvqFLGnF6KsxuvMwcvi7c3FjjtmRpf+O2O9Lf4R8OFw7X1ofzhpOFLaGv1jmg2T2NnML0Rk0PnhW+XGvbyezwTLJw89zkO6rZrb8cLhisLxyfPzagY+7NupfjP1dOuzsGVHIqQFc4mG+16gM6d15L0vg0eOeFxxu3rmvE6WvhpnTlzeGixhrPJQvpe48709/hL0teiB7dm97JLWPSmQ1odk8jp3BbeuGhMLd9J2eHu4ZHVsxuPSesaSw8RkCBrmc2oB1+hT8+LaS/Yb/UuHnZaJVmN9Y4MmbTBUk694T+9P4NvbT5GzecOWFAs3saOYX16YWdoe9Y206+HsI5KzceTm/PbF1v5YuNa1/kPVBMDe+BqnXZe6C7QxhsLB1r3HxFOLP1IvXCcR8L3Zf8Dn9v+GJyzdarptfjNueyDgEd7rCnkVNofvx+tPFqNrOT4W/OTz5PumJt8jo0s3X9HI+NOcfxCChyIqBqXRbQY33NV6C7GzcvD19rX2N007199erND5vqF5/oC5fds/FAJp2ttQ902NPIKYy8Ah0YattJ3b7/WDEnhFtrbVsf72u+An2JgGJqCKhalwW0Nrf5Huj6xs33hGvSG/du3DH+m51XhPt3hnOSz5wuCHc1btiWDWjjG0ibO+xp5BRWphceDBe27+S5byWXhx5t/HKf3Xpu8z3QdQQUU0NA1botoMvDouR35uHLGjdvD32Nb7UPLw6rxgd0dVh0Z/rScHrY1rjhpjEBfTGERgOv7rCnkVNIP4U/cGa4s30nZ/cdSBb2h3Cgfevl4cLkO0BDCwkopoaAqnkO6Oi/xjT6Kfze6WHpc4PfunZmGqerw+xH9g/tXBZmHxgf0MP94ay+RgEXhi8dqg3vvCGEsLO13tCMcOne2sGvntFhTyOn0H/mmn37vnl2OPto+06WhiU7h47vWBLm1dq33t0fLtlyZPuSAQKKqSGgap4DOuqM2kgev5l8kBMGngyNN0MPLmz+BdC22viA1q4NYUnjwpN9oe+s6aHv7jND3+2t9e5JPjMP0zaka2f2NHIKDw+kVz49bifP9zfXb7z0zG79QLqwmoBiagioWtcFtPb8snNnLN6+r/lx/OCqRTNnzL+l8apxXEAfD2Ftes2WpXP65y7/Vu3p8864t7Xe8JpFZ8y4/JnW2mP3NHIKu7513ez+uTftGb+T2s7rz58+a/7Kvemq2a03L53VP3/duHchmggociKgan4DemI7GlH1iYAiJwKq1m0BffD69B/7WJW+9+gSAUVOBFSt2wJ6b7gi+RR+cF64UX0qpVkJ6G1//4Xkr0oJqF0EVK3bArp7RrjxxaPPXhb6dky+slFGAnpJHMfvuJWAWkZA1botoLV109PPle5Tn0h5NgK69pR6QONfWktADSOgal0X0NreW6+cf83KnCdtko2AfjqOp302jn/nIAG1i4CqdV9A/TMR0Kfi+H1PP/37cTyTgNpFQNUIqD0mAjo3js/Zvv3ed8RvP5WAmkVA1QioPRYCOvT++JQ19eXzkzdCCahVBFSNgNpjIaBr4vjPkuVtnyKghhFQNQJqj4WAXhjHFzWuePqzcfzHHVcloHoEVI2A2mMhoH8Sx6ubV911xbqOqxJQPQKqRkDtMRDQnW+Pf3uyVQmoHgFVI6D2GAjoTXHcN9mqBFSPgKoRUHsMBPQLcbx8slUJqB4BVSOg9ugDOvwb8ambJluVgOoRUDUCao8+oJvi+H9NuioB1SOgagTUHn1Ar4rjsyddlYDqEVA1AmqPPqB/E8d3TLoqAdUjoGoE1BW5ZL8AACAASURBVB55QI+eFp++ddJVCageAVUjoPbIA/pAHH928lUJqB4BVSOg9sgDelYcXz75qgRUj4CqEVB71AEd/mB8ytrJVyWgegRUjYDaow7oI3m+xERALSCgagTUHnVAvxDHF+dYlYDqEVA1AmqPOKAffEf87g05ViWgegRUjYDaoxyKzacn/wT99DyrElA9AqpGQO2RDsWl9X5+cGOeNQmoHgFVI6D2aIdi6Sf+eU2uFQmoHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9jgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegRUjYDa42QoCKgeAVUjoPY4GQoCqkdA1QioPU6GgoDqEVA1AmqPk6EgoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9jgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegRUjYDa42QoCKgeAVUjoPY4GQoCqkdA1QioPU6GgoDqEVA1AmqPk6EgoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9jgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegRUjYDa42QoCKgeAVUjoPY4GQoCqkdA1QioPU6GgoDqEVA1AmqPk6EgoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9jgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegRUjYDa42QoCKgeAVUjoPY4GQoCqkdA1QioPU6GgoDqEVA1AmqPk6EgoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9jgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegRUjYDa42QoCKgeAVUjoPY4GQoCqkdA1QioPU6GgoDqEVA1AmqPk6EgoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9jgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegRUjYDa42QoCKgeAVUjoPY4GQoCqkdA1QioPU6GgoDqEVA1AmqPk6EgoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9jgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegRUjYDa42QoCKgeAVUjoPY4GQoCqkdA1QioPU6GgoDqEVA1AmqPk6EgoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9jgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegRUjYDa42QoCKgeAVUjoPY4GQoCqkdA1QioPU6GgoDqEVA1AmqPk6EgoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9jgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegRUjYDa42QoCKgeAVUjoPY4GQoCqkdA1QioPU6GgoDqEVA1AmqPk6EgoHoEVI2A2uNkKAioHgFVI6D2OBkKAqpHQNUIqD1OhoKA6hFQNQJqj5OhIKB6BFSNgNrjZCgIqB4BVSOg9jgZCgKqR0DVCKg9ToaCgOoRUDUCao+ToSCgegS0ejee/nvfHF0ioPY4GQoCqkdAq/d7cTxtdImA2uNkKAioHgGt3Po4jj88ukhA7XEyFARUj4BW7rZ6QN+xZWSRgNrjZCgIqB4BrdxV9YDGD4wsElB7nAwFAdUjoJW7IAnoDSOLBNQeJ0NBQPUIaOVmJAFdMLJIQO1xMhQEVI+AVu5fkoCOfgxPQO1xMhQEVI+AVu4zSUA/M7JIQO1xMhQEVI+AVu5jSUD/YGSRgNrjZCgIqB4BrdxH4vg98bufbS0SUHucDAUB1SOglftQ/M6/iuNHWosE1B4nQ0FA9Qho5d4b/+qcOL66tUhA7XEyFARUj4BW7rT4g3fF8WdbiwTUHidDQUD1CGjl3hX/1tFfjt/V+mNOAmqPk6EgoHoEtHL1gNb+Lo7vay4SUHucDAUB1SOglUsCen4cL2suElB7nAwFAdUjoJVLAnrd6B9zElB7nAwFAdUjoJVLAnpHHM9qLhJQe5wMBQHVI6CVSwK6Oo7/b3ORgNrjZCgIqB4BrVwS0Cfi+B+biwTUHidDQUD1CGjlCKh5ToaCgOoR0MoRUPOcDAUB1SOglSOg5jkZCgKqR0ArR0DNczIUBFSPgFaOgJrnZCgIqB4BrRwBNc/JUBBQPQJaOQJqnpOhIKB6BLRyBNQ8J0NBQPUIaOUIqHlOhoKA6hHQyhFQ85wMBQHVI6CVI6DmORkKAqpHQCtHQM1zMhQEVI+AVo6AmudkKAioHgGtHAE1z8lQEFA9Alo5Amqek6EgoHoEtHIE1DwnQ0FA9Qho5QioeU6GgoDqEdDKEVDznAwFAdUjoJUjoOY5GQoCqkdAK0dAzXMyFARUj4BWjoCa52QoCKgeAa0cATXPyVAQUD0CWjkCap6ToSCgegS0cgTUPCdDQUD1CGjlCKh5ToaCgOoR0MoRUPOcDAUB1SOglSOg5jkZCgKqR0ArR0DNczIUBFSPgFaOgJrnZCgIqB4BrRwBNc/JUBBQPQJaOQJqnpOhIKB6BLRyBNQ8J0NBQPUIaOUIqHlOhoKA6hHQyhFQ85wMBQHVI6CVI6DmORkKAqpHQCtHQM1zMhQEVI+AVo6AmudkKAioHgGtHAE1z8lQEFA9Alo5Amqek6EgoHoEtHIE1DwnQ0FA9Qho5QioeU6GgoDqEdDKEVDznAwFAdUjoJUjoOY5GQoCqkdAK0dAzXMyFARUj4BWjoCa52QoCKgeAa0cATXPyVAQUD0CWjkCap6ToSCgegS0cgTUPCdDQUD1CGjlCKh5ToaCgOoR0MoRUPOcDAUB1SOglSOg5jkZCgKqR0ArR0DNczIUBFSPgFaOgJrnZCgIqB4BrRwBNc/JUBBQPQJaOQJqnpOhIKB6BLRyBNQ8J0NBQPUIaOUIqHlOhoKA6hHQyhFQ85wMBQHVI6CVI6DmORkKAqpHQCtHQM1zMhQEVI+AVo6AmudkKAioHgGtHAE1z8lQEFA9Alo5Amqek6EgoHoEtHIE1DwnQ0FA9Qho5QioeU6GgoDqEdDKEVDznAwFAdUjoJUjoOY5GQoCqkdAK0dAzXMyFARUj4BWjoCa52QoCKgeAa0cATXPyVAQUD0CWjkCap6ToSCgegS0cgTUPCdDQUD1CGjlCKh5ToaCgOoR0MqZD2gIuzpce/D2K+ZMv+Dqh493uG1XCCfp2CffCc+t810loMiLgFbOZ0AfPSOk5m4ef6PNgKZnRUDxMiKglXMZ0E0hLN647+jOh+aG6S+Nu/UkBnSCpJVBQPGyI6CVcxnQy8KXhhsXDp0blo+71WZA9y5bViOgeFkR0Mp5DOjQjJHrVodzxm1hM6ApAoqXEQGtnMeA7gnhcPPiC9ffMG4LAipBQPUIaOXsBnTz0lkDC1YPN6sy+OAlZ0yfd0vjHc/BvnBfe2RWnDcw46LbDiSXG5H6aliU3rImzDia3TwrWXvDpTPPuPiRofY9LWp8TvXEmKPcOLf/nCs3DXc65PG75vafveLQ8KNL5ky/4O5jtUYPt117Zv+cG3aOnNVIQLPnk72r7ZwMBQHVI6CVMxvQVenH7LemVdm/IF2csSG57boQrl1/bMzK6wbSm2dvqzUjtTWEPY2brkreJc1unlFf+6FpV69cHMLK9j19tX7losVbRlbd0Lxt3IrJTq475+br+8LFS/uuXrk0hOR1cf3k+864cuns0L+21hbQ7Plk7+o4ToaCgOoR0MpZDehL/WHRM4e3XR4aVRm+PJy79uDgs9eEGbvrN+6bX796YMl9z6evBWs7B8L8zUcOPD4nnH24Ganhc8P9yU0Hp4Vt7Ztn1NeemXwX6s7QP9i+p+wv1S/NCFdsH9x3ewhPdTjkRYdqtYfqZ/VMfc27w7TBZONwzZH6i80VYeCFbECz55O9q6lj547acNyF+lmrT6HnDQ2pz6DnvDsN6P9rLdene5It9ucr4BQD+uVwwWAykvMbVXkqzG68nhy+Ltyc/Bz8j8X9ja+B3tt4HXpDOD9ZubZ/Zri3Fak70t/hHw4XDo/bfKz62rc1tm0cKLunbEBXhAVD6amt6HDIJ+sLB0L42ui+QpjTWGX4knBjNqDZ88ne1dSxMGrd1P5LAnjZNAP6hfxbHMy32hQDOiesafx8rFGVr4Wb0qs3h4uaKxzbeOtF9bpcnJzOec2Va7eHq1uR2pn+Dn9Z8kK0w+Yj6ms/l/xM34HM7ikb0HlhbePnc7eu6nDIQ+k+tozuKzQaXLc+nJ8NaPZ8sne1ee9mjXpiyIX6WatPoecND6vPoOe0XoG2luvTP8kWB/IVcGoBrb8Ce7Fx4cVGVZaNvh6bPWatF27qS14NHu9rrlxP1byR9xkXJOncE/oPTrh5Q33tI40LyYHa9pQJ6NC01m21iQ45sn4zoBvTVV4K04YyAc2cT9tdHc/J+1q8B6rHe6CVM/oe6O4Q0g+JjjWqckU4c17ThfWXao8daq13XxgYGl25tiNMH6nZfcnv8PeGL9baN88a+WJRcqC2PWUCOnpbZmHsIdsDuiNd5WjyYnhsQDPn03ZXx3MyFARUj4BWzmhAR17hvdSoyvL0vcWmWeHx1sVDSZtGXw5uCOeN1GxvX/2m+WFTrX3zrExA2/aUCeixvubn+pnzyxyyPaDNb0A9H/oGMwHNnE/bXR3PyVAQUD0CWjmjAa3Nbb4xuK5RlXvCNenVezfuSN7X/EprtV1h+tCYNyTvDFeN1uyKcP/OcE7y3lx286xMQNv2lH0PdG5Y3/j54vU3TXTItoCuTFdZNdrY9H+z55O9q+M5GQoCqkdAK2c1oMvDhckXY4YWNqqyPfQ9m1w7vDisanxb6NH0C0zDy8MlteRz7Hmtj8TvGa3Z6rDozvSlXnbzrGxAs3vKBvRLYVHjqLclyex8yLaATm+8ZD0wO3mjdmxAs+eTvavjORkKAqpHQCtnNaC7+8MlW45sXzKQVuXqMPuR/UM7l4XZB+q/9NZTc+mjW3bvWLs4hOSr6DsHwoKnjxxMvpR5aLRmh/vDWX3p9z4zm9eWLRvzQVg2oNk91TJ/h/R8f1iyffDQqr7kHYTOh2wLaDjrm/v2rT07zNzf9kX6zPm03dVxnAwFAdUjoJWzGtDaA+mn1KvTqhxc2PzTneTvfmoHr2x9iD3wcGPltf3p4pnP1sYk8doQlqSXspuHMObr9NmAZvdUO6fewNE/XnpsWnrbv094yLaA3tHXWKV/Xa0toNnzyd7VcZwMBQHVI6CVMxvQ5A/E++evazVpcNWimTPm39J86Ti8afll5wyct+S21h+3v7Di3IGzFt8x+rfwicdD85ubbZufKKCZPdU2X9R/1lOj6z6/fOH0uUufmviQbQHdtf3a2dM6/y185u5k72o7J0NBQPUIaOXsBvTlc3zMx+kvo5P1Tzk5GQoCqkdAK9eLAd0djk2+0tQRUFSMgFauFwP6tUsrOQwBRcUIaOV6MaDLOvyzoC8DAoqKEdDK9WJAK0JAUTECWjkC+rIhoKgYAa0cATXPyVAQUD0CWjkCap6ToSCgegS0cgTUPCdDQUD1CGjlCKh5ToaCgOoR0MoRUPOcDAUB1SOglSOg5jkZCgKqR0ArR0DNczIUBFSPgFaOgJrnZCgIqB4BrRwBNc/JUBBQPQJaOQJqnpOhIKB6BLRyBNQ8J0NBQPUIaOUIqHlOhoKA6hHQyhFQ85wMBQHVI6CVI6DmORkKAqpHQCtHQM1zMhQEVI+AVo6AmudkKAioHgGtHAE1z8lQEFA9Alo5Amqek6EgoHoEtHIE1DwnQ0FA9Qho5QioeU6GgoDqEdDKEVDznAwFAdUjoJUjoOY5GQoCqkdAK0dAzXMyFARUj4BWjoCa52QoCKgeAa0cATXPyVAQUD0CWjkCap6ToSCgegS0cgTUPCdDQUD1CGjlCKh5ToaCgOoR0MoRUPOcDAUB1SOglSOg5jkZCgKqR0ArR0DNczIUBFSPgFaOgJrnZCgIqB4BrRwBNc/JUBBQPQJaOQJqnpOhIKB6BLRyBNQ8J0NBQPUIaOUIqHlOhoKA6hHQyhFQ85wMBQHVI6CVI6DmORkKAqpHQCtHQM1zMhQEVI+AVo6AmudkKAioHgGtHAE1z8lQEFA9Alo5Amqek6EgoHoEtHIE1DwnQ0FA9Qho5QioeU6GgoDqEdDKEVDznAwFAdUjoJUjoOY5GQoCqkdAK0dAzXMyFARUj4BWjoCa52QoCKgeAa0cATXPyVAQUD0CWjkCap6ToSCgegS0cgTUPCdDQUD1CGjlCKh5ToaCgOoR0MoRUPOcDAUB1SOglSOg5jkZCgKqR0ArR0DNczIUBFSPgFaOgJrnZCgIqB4BrRwBNc/JUBBQPQJaOQJqnpOhIKB6BLRyBNQ8J0NBQPUIaOUIqHlOhoKA6hHQyhFQ85wMBQHVI6CVI6DmORkKAqpHQCtHQM1zMhQEVI+AVo6AmudkKAioHgGtHAE1z8lQEFA9Alo5Amqek6EgoHoEtHIE1DwnQ0FA9Qho5QioeU6GgoDqEdDKEVDznAwFAdUjoJUjoOY5GQoCqkdAK5cEdHUc/2tzkYDa42QoCKgeAa1cEtA74nh2c5GA2uNkKAioHgGtXBLQL8bx/OYiAbXHyVAQUD0CWrkkoPPieFlzkYDa42QoCKgeAa1cEtB/iON7m4sE1B4nQ0FA9Qho5eoBPfqe+J1bmosE1B4nQ0FA9Qho5U6LP/j1OP5Ua5GA2uNkKAioHgGt3HvjXz0njpe2FgmoPU6GgoDqEdDKfSh+5yfjeHVrkYDa42QoCKgeAa3cR+L4V+N3bWstElB7nAwFAdUjoJX7s7ju90cWCag9ToaCgOoR0Mp9Ognop0cWCag9ToaCgOoR0Mr9cxLQaSOLBNQeJ0NBQPUIaOWmJwGdP7JIQO1xMhQEVI+AVm5eEtAbRhYJqD1OhoKA6hHQyi1NArpqZJGA2uNkKAioHgGt3Mp6P099ZmSRgNrjZCgIqB4BrdwTmW8xEVCDnAwFAdUjoNX77TjuG10ioPY4GQoCqkdAq3f9uz7w2OgSAbXHyVAQUD0CKvDMtjELBNQeJ0NBQPUIqBoBtcfJUBBQPQKqRkDtcTIUBFSPgKoRUHucDAUB1SOgagTUHidDQUD1CKgaAbXHyVAQUD0CqkZA7XEyFARUj4CqEVB7nAwFAdUjoGoE1B4nQ0FA9QioGgG1x8lQEFA9AqpGQO1xMhQEVI+AqhFQe5wMBQHVI6BqBNQeJ0NBQPUIqBoBtcfJUBBQPQKqRkDtcTIUBFSPgKoRUHucDAUB1SOgagTUHidDQUD1CKgaAbXHyVAQUD0CqkZA7XEyFARUj4CqEVB7nAwFAdUjoGoE1B4nQ0FA9QioGgG1x8lQEFA9AqpGQO1xMhQEVI+AqhFQe5wMBQHVI6BqBNQeJ0NBQPUIqBoBtcfJUBBQPQKqRkDtcTIUBFSPgKoRUHucDAUB1SOgagTUHidDQUD1CKgaAbXHyVAQUD0CqkZA7XEyFARUj4CqEVB7nAwFAdUjoGoE1B4nQ0FA9QioGgG1x8lQEFA9AqpGQO1xMhQEVI+AqhFQe5wMBQHVI6BqBNQeJ0NBQPUIqBoBtcfJUBBQPQKqRkDtcTIUBFSPgKoRUHucDAUB1SOgagTUHidDQUD1CKgaAbXHyVAQUD0CqkZA7XEyFARUj4CqEVB7nAwFAdUjoGoE1B4nQ0FA9QioGgG1x8lQEFA9AqpGQO1xMhQEVI+AqhFQe5wMBQHVI6BqBNQeJ0OxYf0G9Sn0vBd2qc+g561fv3GyVfKNPQHtMbPCWepTAMT2h3DFydkTAe0xBBQgoCiJgAIEFCURUICAoiQCChBQlERAAQKKkggoQEBREgEFCChKIqAAAUVJBBQgoCiJgAIEFCUdPnRYfQqA2PChQ0dOzp4IKACUREABoCQCCgAlEVAAKImAAkBJBBQASiKgAFASAQWAkghoF7sxtC4dX7Vw5syF9x+ffAnoFmWe/0WngYB2r6OzQvPSsctCwyXHJlsCukWZ53/haSCgXevQtSE0L64IMx966aWHZoYVky0BXaLU87/wNBDQ7nT09qtmhNYTaPe08GTy88kwbc+Jl4CuUPL5X3waCGh32p/+JpIu3B/OG05+Dp8XHjzxEtAVSj7/i08DAe1Ox7/1rW+tbz2BvhhuSy/cFr504iWgK5R8/hefBgLavXa1nkCXhjXphTXh8hMvAV2jxPO/+DQQ0O418gQ6L2xKL2wK8068BHSNEs//4tNAQLvXyBNoRtiRXtgRzjjxEtA1Sjz/i08DAe1e459A28P0Ey8BXaPE87/4NBDQ7jXmV5in0gubwnknXgK6Ronnf/FpIKDda+QJdFn4ZnphTbjsxEtA1yjx/C8+DQS0e+0a/RrHnemFO8OyEy8BXaPE87/4NBDQ7rVr9IvEF6YXLgyrTrwEdI0Sz//i00BAu9euMX/KtjX5ua3152oTLgFdo8Tzv/g0ENDuNfIEqt0U5mwcPPbknPDVyZaAblHm+V94Ggho9xp9Ah27NISB/hAuH5xsCegWZZ7/haeBgHav0SdQ7fj9C2bMWPjA0ORLQJco9fwvOg0EFABKIqAAUBIBBYCSCCgAlERAAaAkAgoAJRFQACiJgAJASQQUAEoioHDhj6J2H53qLl/T2M2KnKuO+VGVjoe7Lb33VZ4HJsYDARc6BzSKXlN+l3kC2jxANqBTOmoBBNQ+Hgi4cNGfpf5TFDUvLahNPaArVqx4/oTr2Avozvo5v5KAWsEDAVeyUZlqQCddx15AT3gDqsYDAVey7Vi27PqTtavOmgfIBnRKRy2AgNrHAwFXTmI7CuzK0IdIgvPAxHgg4Eq2HWPC9vzf/8TrXv8Tf/1srXbXH/7Qq7/3lL6R/5sMuwZ+6c2v+bH3X3Vkol2dYPNsOds6enjhB3762783/uiDI7scd6jBKz70Q69+6/tmHWwuH7/2I2993bf/l49cfXz00Mcv+R9vfs2Pf+CG1r9B+UL4xTe+7qc+v2vi4xBQM3gg4MpEAX3ge6Lv/PYoit7wwCej6DuSSz/yQrrOtd/V/Nj+LY9OsKsTbH6igH7l+1tfCPi9QxMcauPPNq/4wVsbyxt+rrXJz6xvHXrve6Lozf+pftX7jjauuvmN6Ro/8B8THoeAmsEDAVcmCugbf/K+40M3fU8UfVv0tnsHh2767ij6RGOVq+vh+fE/+/zvvCGKXr+q865OsPkJAvrAK6LoP3/83/7u3fUD/HXnQ239vij6rg98/i9+JIpe+3h9eXN9t69892c/d9qrouhNTzV39f5XzNpROzyvfunzyTX3vLK+8q/8za+/PnrrKyY4DgG1gwcCrkwU0B9s/Mp9cT0yP7w3ubQgin46+fnsd0avODf5ffn534uinz3WcVcn2PwEAf1AFP1pY39L6v0d7nio90ZRvK3+89DHouiPa7Xh06LorQ8kmySvTE8bTnfV/CbqrCj6+fqPY/81it62oX5hy880v+457jgE1BAeCLgyUUAXNZa31puzuHXptcnPf4qigXTdY3EUXdVxVyfY/AQB/d4o2pFu/9+i6MlOh7q//vpzd+OKIz8avXm4dlf99eemdI3N9deZX0939f7mNemXo74YRa97rnHFrtdHnY9DQA3hgYArEwU0LctQvTlbW5caX9Z8a/Sm1v+BxZui6OMdd3WCzU8Q0LdE0QOZUxt3qM9G0f9pXvGF7//+jbW/jqJPtVb+31H0t+multfGHvGjUfT3zVX+Lup8HAJqCA8EXJkooM1PtaPoVcOtS0mPno2i32mtvCuKfq7jribe/EQB/aMoevMFO0d3N/5QvxhF9409YH15ZevyLVF0arqr9bWxR/yxKLqnecU9UcfjjPuPACEeCLgy4deYUqN/JJReejD79/Nv6biriTc/UUCf+5Fkjz//V5c+nW4x/lD1X76fHXvAH4iip1qXN0TRD6e7an3lKT3id478wp4UudNxxv1HgBAPBFwpGNDbslV7Y8ddlQtobd+/vbn5laSzjnQ81Cuibzs+9oCvj6JDrcsHo+g7sncnPWJ9m9b3QY9HHY8z7j8ChHgg4ErBgD4QRR99YtT6jrsqGdBabegbMz7w3UnafvLFTof6rijaNfaA/3nMK9DNUfQDtQ4BfdPoK9Dnoo7HGfcfAUI8EHClYEC3RNGvTbqr0gFNDD/6D98TRR/rdKgfj6KHmxeHDh8erP1CFN3Wum1lFP1CrUNAf3b0PdD7o47HGfcfAUI8EHClYECHfzB6c+tNxqfmzr23467KBfTief/eXPfhKPrRTof60yg6o3nF/Cj6cu3TUfSZ1sE/k14eF9C/GP0U/l+jjscZ9x8BQjwQcKVgQJOvDk1rXvP7UXR3x12VC+jPRz/c/MT+2Sh6c6dD3RBF37OvsTx8evRtL9buiKJXbU5X2PTKKLqj1iGgK6Lodenv8HveEHU+DgE1hAcCrhQN6KbXRK++IPlY5vgXouhtwx13VS6gfxJFobG/w38URb/V6VDHfyaK3r29fsVQX/1CPaPvqL+EbPxS/8Bbo+hdw7UOAT3+k82/RNrx9uZfIo07DgE1hAcCrhQNaO2ieoZ+6mP/9sl6st64KdeX1gAAAbBJREFUpvOuygX0pvqO3/aZ8PmPvimKXvlwx0M99u1R9IYP/d+/rof0NQ/Vlze+sf4a9PTPfe70+uvP797cdneaR1z16ih67Xv+9kNviN7/6gmOQ0Dt4IGAK4UDWrvotc0vFv1Y5lvtU/8Qqf/bRr7zeesEh3roR5pXvOnmxvITP93a5OfH/VFm69i3vild4xd2TXgcAmoGDwRcKR7Q2tb/89/f+Nqf+q3LB2tZU/4U/huffPdbX/198YevH9nzuEMdmvfe73vVj76//1vN5cGrPvyW173uLX943Zh/D7T92C/0/eIbXvO2vgMTH4eAmsEDgV7lN0N+z7zr8ECgV/nNkN8z7zo8EOhVfjPk98y7Dg8EepXfDPk9867DA4FeVc/QE088sV99GgUdqJ/zqwioFTwQ6FWvaXw1aIX6NApq/qNP6tNAigcCvYqAYsp4IACgJAIKACURUAAoiYACQEkEFABKIqAAUBIBBYCSCCgAlERAAaAkAgoAJRFQACiJgAJASQQUAEoioABQEgEFgJIIKACUREABoCQCCgAlEVAAKImAAkBJBBQASiKgAFDS/wdhhu7VKdg7WgAAAABJRU5ErkJggg==" title alt width="672" /></p> -<p>We see that using the compiled model is by a factor of 13.7 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,iVBORw0KGgoAAAANSUhEUgAABUAAAAPACAMAAADDuCPrAAAC+lBMVEUAAAABAQECAgIDAwMEBAQFBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQVFRUWFhYXFxcYGBgZGRkaGhobGxscHBwdHR0eHh4fHx8gICAhISEiIiIjIyMkJCQlJSUmJiYnJycoKCgpKSkqKiorKyssLCwtLS0uLi4wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6Ojo7Ozs8PDw9PT0+Pj4/Pz9AQEBBQUFCQkJDQ0NERERFRUVGRkZHR0dISEhJSUlKSkpLS0tMTExNTU1PT09QUFBRUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVlaWlpbW1tcXFxdXV1eXl5fX19gYGBhYWFiYmJjY2NkZGRlZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1ubm5vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8fHx9fX1+fn5/f3+AgICBgYGCgoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uMjIyNjY2Ojo6Pj4+QkJCRkZGSkpKTk5OUlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2enp6fn5+goKChoaGioqKjo6OkpKSlpaWmpqanp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKzs7O0tLS1tbW2tra3t7e4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXGxsbHx8fIyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY2NjZ2dna2trb29vc3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr6+vs7Ozt7e3u7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+/v7////Ha1yKAAAACXBIWXMAAB2HAAAdhwGP5fFlAAAgAElEQVR4nO3dCZQc5X3v/cr1Ht/EdhInN7lJfPO+yXW2mzd5AC+AHRITB+I4hjheMXYSL9jGsX3txDbGz2i0IUBIsgBJCAkki1XGICQTIxDggBCLAbMIJLSAJJpFYtMK0ix9zttV1dPLTI+6q1rU7//MfD/nJOqufrqqh67/9/RM94yjMgAgl0j9AAAgVAQUAHIioACQEwEFgJwIKADkREABICcCCgA5EVAAyImAAkBOBBQAciKgAJATAQWAnAgoAOREQAEgJwIKADkRUADIiYACQE4EFAByIqAAkBMBBYCcCCgA5ERAASAnAgoAORFQAMiJgAJATgQUAHIioACQEwEFgJwIKADkREABICcCCgA5EVAAyImAAkBOBBQAciKgAJATAQWAnAgoAOREQA+VPS+Yt3f//p3qxzDO7XxR/QjGud379+/ubGVnY09AD5XnSubtLZe3qx/DOPdiAKfJmFbp4s7OVnY29gT0UAlgMgioHAEVI6BWBTAZBFSOgIoRUKsCmAwCKkdAxQioVQFMBgGVI6BiBNSqACaDgMoRUDECalUAk0FA5QioGAG1KoDJIKByBFSMgFoVwGQQUDkCKkZArQpgMgioHAEVI6BWBTAZBFSOgIoRUKsCmAwCKkdAxQioVQFMBgGVI6BiBNSqACaDgMoRUDECalUAk0FA5QioGAG1KoDJIKByBFSMgFoVwGQQUDkCKkZArQpgMgioHAEVI6BWBTAZBFSOgIoRUKsCmAwCKkdAxQioVQFMBgGVI6BiBNSqACaDgMoRUDECalUAk0FA5QioGAG1KoDJIKByBFSMgFoVwGQQUDkCKkZArQpgMgioHAEVI6BWBTAZBFSOgIoRUKsCmAwCKkdAxQioVQFMBgGVI6BiBNSqACaDgMoRUDECalUAk0FA5QioGAG1KoDJIKByBFSMgFoVwGQQUDkCKkZArQpgMgioHAEVI6BWBTAZBFSOgIoRUKsCmAwCKkdAxQioVQFMBgGVI6BiBNSqACaDgMoRUDECalUAk0FA5QioGAG1KoDJIKByBFSMgFoVwGQQUDkCKkZArQpgMgioHAEVI6BWZZmMldPuyPXsd4mAyhFQMQJqVYbJuP+d7v1b8j3/XSGgcgRUjIBalWEy5jrnrsv3/HeFgMoRUDECalWGyfhGJaCz8z3/XSGgcgRUjIBalWEyPlMJ6DfzPf9dIaByBFSMgFqVMaAn5Xv+u0JA5QioGAG1KmNA/zbf898VAipHQMUIqFUZA3rEtnwnQDcIqBwBFSOgVmUMqFub7wToBgGVI6BiBNSqrAG9Ld8J0A0CKkdAxQioVVkDujzfCdANAipHQMUIqFVZA7ok3wnQDQIqR0DFCKhVWQM6J98J0A0CKkdAxQioVVkDOi3fCdANAipHQMUIqFVZA3pavhOgGwRUjoCKEVCrsgb01HwnQDcIqBwBFSOgVmUN6KfznQDdIKByBFSMgFqVNaD/mO8E6AYBlSOgYgTUqqwBFfwyPAGVI6BiBNSqbAE9wh2V7wToBgGVI6BiBNSqbAH9a+e25jsDukBA5QioGAG1KltAP+Tcw/nOgC4QUDkCKkZArcoW0E85d3e+M6ALBFSOgIoRUKuyBfRLzt2S7wzoAgGVI6BiBNSqbAH9D8X/LicBlSOgYgTUqmwBnezclfnOgC4QUDkCKkZArcoW0FnOXZTvDOgCAZUjoGIE1KpsAV3g3Ln5zoAuEFA5AipGQK3KFtArFH/PjoDKEVAxAmpVtoD+2LmefGdAFwioHAEVI6BWZQvoT537Rr4zoAsEVI6AihFQq7IF9F7nvpjvDOgCAZUjoGIE1KpsAX3UuZPznQFdIKByBFSMgFqVLaAl5z6c7wzoAgGVI6BiBNSqbAF99l3u7/KdAV0goHIEVIyAWpUxoMe6v8x3BnSBgMoRUDECalXGgH7IHfFEvlMgPwIqR0DFCKhVGQN6knOb8p0C+RFQOQIqRkCtyhjQU5x7MN8pkB8BlSOgYgTUqowB/aZzd+Y7BfIjoHIEVIyAWpUxoBOcuynfKZAfAZUjoGIE1KqMAT3buRX5ToH8CKgcARUjoFZlDOhcwV9UJqByBFSMgFqVMaBLnFuY7xTIj4DKEVAxAmpVxoAuE/xFZQIqR0DFCKhVGQO6SvAXlQmoHAEVI6BWZQzonc75fKdAfgRUjoCKEVCrMgZ0reAvKhNQOQIqRkCtyhjQrc6dku8UyI+AyhFQMQJqVcaAPu/cp/KdAvkRUDkCKkZArcoY0L7D3An5ToH8CKgcARUjoFZlDGj5aPf+fKdAfgRUjoCKEVCrsgb0OHdUvlMgPwIqR0DFCKhVWQP6EXfY1srlbRec/ViuMyEHAipHQMUIqFVZA/qvzj1cuTzduTm5zoQcCKgcARUjoFZlDejXnbu7VNr6V86dke9UyI6AyhFQMQJqVdaA9iR/EPQmR0DHEwIqRkCtyhrQc5y7tlQ6j4COKwRUjIBalTWgC5xbUip9jYCOKwRUjIBalTWgS52bWyqdQEDHFQIqRkCtyhrQlc6dWdp6BAEdVwioGAG1KmtA73Lu9NKdjoCOKwRUjIBalTWgjzr3ldJVBHR8IaBiBNSqrAHd4dynS7MJ6PhCQMUIqFVZA9of/zmmbxHQ8YWAihFQq7IGtPxX7q9KJxPQ8YWAihFQqzIH9MPu8K3vI6DjCwEVI6BWZQ7oF5y7wxHQ8YWAihFQqzIH9HTn5jv3LgI6nhBQMQJqVeaAnuvcqc59lICOJwRUjIBalTmgS5072rlvEdDxhICKEVCrMgd0dfID0AsI6HhCQMUIqFWZA7ol7uc7biCg4wkBFSOgVmUO6MC7KgH9+C0EdDwhoGIE1KrMAS1/uhLQaQR0XCGgYgTUquwBPb8S0NUEdFwhoGIE1KrsAX3mGPeZAQI6rhBQMQJqVfaAlp9f/XKZgI4rBFSMgFqVI6AxAjquEFAxAmoVAUV7BFSMgFpFQNEeARUjoFYRULRHQMUIqFUEFO0RUDECahUBRXsEVIyAWkVA0R4BFSOgVhFQtEdAxQioVQQU7RFQMQJqFQFFewRUjIBaRUDRHgEVI6BWEVC0R0DFCKhVBBTtEVAxAmoVAUV7BFSMgFpFQNEeARUjoFYRULRHQMUIqFUEFO0RUDECahUBRXsEVIyAWkVA0R4BFSOgVhFQtEdAxQioVQQU7RFQMQJqFQFFewRUjIBaRUDRHgEVI6BWEVC0R0DFwg+o9ztabN1748XTJ5176d39LW7b4f0r/JgyHbH1F0BA0QECKjZGA3rfVJ+asXnkjUUGND0WAcUrg4CKjc2AbvR+0YZd+7ffNcNPen7ErQT0ECGgcgRUbGwGdKH/4WByYd9Mv2zErUUGdOfSpe2OSECRGwEVG5MBHZhc23aPP2fEPfgZ6CFCQOUIqNiYDOiL3r9UvfjsVVePuAcBPUQIqBwBFQs4oJuXnDFx7j2D1f703Xnh1Emzf5L8xLOvx68etvi55bMmTj7/hj3x5SRnP/bz01se8JP3N999uOeundF7zg82DrbaUf8tM3rPXr5v8L7F0yed+9MD5aSH2y4/s3f61dtrx6oFtPkozV/AiKN2/iQS0HGLgIqFG9A16dvs16f92T03vTp5fXzbFd5fvu5Aw+K1E9Obp20rV3O21fsXk5suiX9K2nz3Zuur913ZakdXnHPdVT3+giU9l65c4n38arfykHqm/mDJNN/7UHlYQJuP0vwFjEBA0R4BFQs2oM/3+vmPv7TtIp/0Z/AiP/OhvX1PXOYnV76i8q45lc0TF69+On3VWN4+0c/Z/PKeh6f7s1+q5mxwpr89vmnvBL9t+N2bDzTZX1zq23Wj95ta7Oj8feXyXZVjPV5Z+VM/oS8OqL/s5cqLzeV+4rPNAW0+SvMXkDows259f8c+3xzQ2Z3fszuV/7oDRR0LLQ3wBGgNdDoDuzvrWmEB/ZE/t1Krcv+cpD+b/LTk9eTgFf66+N++ny/qTT4GelvyOvRq//14cXn3FH/bUM5Wpd/D3+3PGxxx90bL/dyB9IDLW+zo0cqVPd7fnGxMHor305Mlgxf6a5sD2nyU5i8gdcDXre38v8UXmgN6buf3BFCIvZ0tKyyg0/0Dyb/3J/252a9IN2/251cXHNhw/fmVDl0QP/BZ1cXlG/2lQznbnn4PvzB+Idri7jWz/UPJv09dv6bFjiovQMuD3m8pp/8mAb0tXbLOf785oM1Haf4Cqo+ZgAJjkrGAVlLzXHLhuaQ/S+vhmdaw6tkVPfHrxv6e6uJK1GbXfiI5N07ni75376h3jw1MGLpvebQd1d5HrwZ0Q7rkeT9hoCmgTUcZ9gWMxM9A0R4/AxUL9WegL3ifvkl0IOnPxf7M2VXnVV7U3b9vaN1qP3Ggvrj8pJ9U697q+Hv42/yV5eF3b32g5iuNOxoe0CfTJfvjl7iNAW06yrAvYCQCivYIqFioAa29Fnw+6c+y9KeQVWf4h4cu7osrVn/huN7PqnVvZ0/lpjl+Y3n43Zsc6Km+W9901KYdDQ/oI+mSp31PX1NAm44y7AsYiYCiPQIqFmpAyzOqP0Jcm/TnVn9Zunnnhifjn2teM7Rsh5800PCjy5v8JfXuXexv3+7Pid8har778AOtS/597qoVo+1oWEBXpkvW1Bub/v/mozR/ASMRULRHQMWCDegyf178p+oG5iX9KfmeJ+Ktg4v8muRzRfelH2AaXOYvLMfveM8eevP81nr37vHzb0pfFDbfvdkP/fxkXzfEyWy9o2EBnZS8ZN0zLf7xa2NAm4/S/AWMREDRHgEVCzagL/T6C7e8XFo8Me3PpX7avbsHti/10/ZUvj2uRGnBfVteePKhRd7HH1rfPtHPfezlvfHHN/fVu/dSrz+rJ/3cZ9Pdy0uX7qkf6Olev7jUt29NT/xzgdY7GhZQf9aDu3Y9dLafsnvYB+mbjjLsCxiBgKI9AioWbEDLd6TvZ9+T9mfvvOov+cS/IVTe+4Oht7sn3p0sfqg3vXpm/BKw9pvpl3u/OL3UfHfvGz9Of/+E9Lb/HHVHwwK6qidZ0ru2PCygzUdp/gJGIKBoj4CKhRvQ+FfJe+esHapX35r5UybP+Un1pePgxmULz5k4a/ENQ7/c/uzymRPPWrSq/rvwsYd99TOew+7eHNDy08vmTZqxZNPoOxoW0B2ly6dNaP278E0PsvkLGI6Aoj0CKhZwQF85/Q1vvGc2Sg8zI6Boj4CKEdAWXvAH2i8aDQFFcQioGAFt4eYFXdyZgKI4BFSMgLawtOWfBe0QAUVxCKgYAT3UCCiKQ0DFCOihRkBRHAIqRkCtIqBoj4CKEVCrCCjaI6BiBNQqAor2CKgYAbWKgKI9AipGQK0ioGiPgIoRUKsIKNojoGIE1CoCivYIqBgBtYqAoj0CKkZArSKgaI+AihFQqwgo2iOgYgTUKgKK9gioGAG1ioCiPQIqRkCtIqBoj4CKEVCrCCjaI6BiBNQqAor2CKgYAbWKgKI9AipGQK0ioGiPgIoRUKsIKNojoGIE1KocAR14hoCOMwRUjIBalT2guz7sTiOg4wsBFSOgVmUP6CLn3P0EdFwhoGIE1KrsAf18JaCzCOi4QkDFCKhVmQM6eHQloJ8hoOMKARUjoFZlDuhTlX66d68ioOMJARUjoFZlDuidcUDdxQR0PCGgYgTUqswBvdq5dzp3OgEdTwioGAG1KnNA5zj3WedOIqDjCQEVI6BWZQ7oBOfOde5oAjqeEFAxAmpV5oB+2bmfJj8GJaDjBwEVI6BWZQ7oR91hj7+XgI4vBFSMgFqVOaDHuveWPk5AxxcCKkZArcoa0IHD3QdL3yCg4wsBFSOgVmUN6HPOnVyaQUDHFwIqRkCtyhrQTc6dWrqCgI4vBFSMgFqVNaD3OHda6TYCOr4QUDECalXWgK5yblrpscMJ6LhCQMUIqFVZA3qVc+eXSh8goOMKARUjoFZlDehFzi0plU4loOMKARUjoFZlDegs564uJW/DE9Dxg4CKEVCrsga017kbS6XrCOi4QkDFCKhVWQP6TefuKJUeP4qAjicEVIyAWpU1oJ9z7qHK5dOdm5nrTMiBgMoRUDECalXWgH7Muccrlzee/vWHc50JORBQOQIqRkCtyhrQv3NH5jsF8iOgcgRUjIBalTWg73XH5jsF8iOgcgRUjIBalTGgA4e5E/KdAvkRUDkCKkZArcoY0J3OfTLfKZAfAZUjoGIE1KqMAX3Cuc/nOwXyI6ByBFSMgFqVMaCPOPd/850C+RFQOQIqRkCtyhjQu507Pd8pkB8BlSOgYgTUqowBvanA30AaQkDlCKgYAbUqY0CXOTc73ymQHwGVI6BiBNSqjAG91LkF+U6B/AioHAEVI6BWZQzoBc5dnu8UyI+AyhFQMQJqVcaAnuPctflOgfwIqBwBFSOgVmUM6ETnVuU7BfIjoHIEVIyAWpUxoP/u3Jp8p0B+BFSOgIoRUKsyBvSLzj2Q7xTIj4DKEVAxAmpVxoCe7NyGfKdAfgRUjoCKEVCrMgb0BHf4E/lOgfwIqBwBFSOgVmUM6Pvde/OdAV0goHIEVIyAWpUxoEe64/OdAV0goHIEVIyAWpUtoE8794/5zoAuEFA5AipGQK3KFtCNzn0q3xnQBQIqR0DFCKhV2QJ6n3On5DsDukBA5QioGAG1KltAb3XuG/nOgC4QUDkCKkZArcoW0P907nv5zoAuEFA5AipGQK3KFtArBX9PmYDqEVAxAmpVtoAuFPw9ZQKqR0DFCKhV2QI6W/D3lAmoHgEVI6BWZQvoZOeuyHcGdIGAyhFQMQJqVbaAftu55fnOgC4QUDkCKkZArcoW0FOduznfGdAFAipHQMUIqFXZAnqyc3fnOwO6QEDlCKgYAbUqW0BPdG5tvjOgCwRUjoCKEVCrsgX0WOe25DsDukBA5QioGAG1KltA3+mOzHcCdIOAyhFQMQJqVbaAOvc3+U6AbhBQOQIqRkCtyhrQE/OdAN0goHIEVIyAWpU1oCfnOwG6QUDlCKgYAbUqa0BPzXcCdIOAyhFQMQJqVdaAfivfCdANAipHQMUIqFVZAzol3wnQDQIqR0DFCKhVWQNa/F+zI6B6BFSMgFqVNaCL850A3SCgcgRUjIBalTWgy/KdAN0goHIEVIyAWpU1oLfkOwG6QUDlCKgYAbUqa0AfyHcCdIOAyhFQMQJqVcaAHr413wnQDQIqR0DFCKhVGQP6vlxPf3cIqBwBFSOgVmUM6MfyPf9dIaByBFSMgFqVMaBfy/f8d4WAyhFQMQJqVYbJ+EoloNPzPf9dIaByBFSMgFqVYTK+XwnoNfme/64QUDkCKkZArcowGWsOd+/ZmO/57woBlSOgYgTUqiyTsegrK3I9+10ioHIEVIyAWhXAZBBQOQIqRkCtCmAyCKgcARUjoFYFMBkEVI6AihFQqwKYDAIqR0DFCKhVAUwGAZUjoGIE1KoAJoOAyhFQMQJqVQCTQUDlCKgYAbUqgMkgoHIEVIyAWhXAZBBQOQIqRkCtCmAyCKgcARUjoFYFMBkEVI6AihFQqwKYDAIqR0DFCKhVAUwGAZUjoGIE1KoAJoOAyhFQMQJqVQCTQUDlCKgYAbUqgMkgoHIEVIyAWhXAZBBQOQIqRkCtCmAyCKgcARUjoFYFMBkEVI6AihFQqwKYDAIqR0DFCKhVAUwGAZUjoGIE1KoAJoOAyhFQMQJqVQCTQUDlCKgYAbUqgMkgoHIEVIyAWhXAZBBQOQIqRkCtCmAyCKgcARUjoFYFMBkEVI6AihFQqwKYDAIqR0DFCKhVAUwGAZUjoGIE1KoAJoOAyhFQMQJqVQCTQUDlCKgYAbUqgMkgoHIEVIyAWhXAZBBQOQIqRkCtCmAyCKgcARUjoFYFMBkEVI6AihFQqwKYDAIqR0DFCKhVAUwGAZUjoGIE1KoAJoOAyhFQMQJqVQCTQUDlCKgYAbUqgMkgoHIEVIyAWhXAZBBQOQIqRkCtCmAyCKgcARUjoFYFMBkEVI6AihFQqwKYDAIqR0DFCKhVAUwGAZUjoGIE1KoAJoOAyhFQMQJqVQCTQUDlCKgYAbUqgMkgoHIEVIyAWhXAZBBQOQIqRkCtCmAyCKgcARUjoFYFMBkEVI6AihFQqwKYDAIqR0DFCKhVAUwGAZUjoGIE1KoAJoOAyhFQMQJqlWgy7v3CKQ90upaAyhFQMQJqlWgyTnbuc52uJaByBFSMgFqlmYx7nXOHdfoSlIDKEVAxAmqVZjIWVQLqlnS4mIDKEVAxAmqVZjJOjwPqO1xMQOUIqBgBtUozGSfFAT25w8UEVI6AihFQqzST8ZfuPe9xx3S4mIDKEVAxAmqVZDLWO3dS5UXous5WE1A5AipGQK2STMZ/Offt05y7ubPVBFSOgIoRUKskk7HKud7znbuss9UEVI6AihFQq2QBXebczM5WE1A5AipGQK2SBfQu577T2WoCKkdAxQioVbKAbuv4lzkJqBwBFSOgVskCuv8wd0JnqwmoHAEVI6BWyQJaPtYd3dlqAipHQMUIqFW6gH7CuY0drSagcgRUjIBapQvo15y7s6PVBFSOgIoRUKt0AZ3q3HUdrSagcgRUjIBapQvohZ3+QTsCKkdAxQioVbqALnNuVkerCagcARUjoFbpAnq7cz0drSagcgRUjIBapQvoo86d2tFqAipHQMUIqFW6gL7g3Cc7Wk1A5QioGAG1ShfQwXe54ztaTUDlCKgYAbVKF9DyB9w7O1pNQOUIqBgBtUoY0E87t76T1QRUjoCKEVCrhAH9pnNrOllNQOUIqBgBtUoY0Gkd/ioSAZUjoGIE1CphQBd0+KtIBFSOgIoRUKuEAe30V5EIqBwBFSOgVgkDertzvpPVBFSOgIoRUKuEAd3Q4a8iEVA5AipGQK0SBvRF5z7eyWoCKkdAxQioVcKAlt/t3t/JagIqR0DFCKhVyoD+gztiWwerCagcARUjoFYpA/o55x7sYDUBlSOgYgTUKmVAT3fu5g5WE1A5AipGQK1SBvQ8567sYDUBlSOgYgTUKmVAr3LuvMq1dTMu2Hqw1QRUjoCKEVCrlAFd7dz3SqWt/+jcvIOtJqByBFSMgFqlDOjjzn2uVFrknJt8sNUEVI6AihFQq5QB3X+4+2Cp9EkCah4BFSOgVikDWv6Ae8e2Bw8noOYRUDECapU0oF9y7u6FjoCaR0DFCKhV0oCe7dxVXyag9hFQMQJqlTSg1zh3xlEE1D4CKkZArZIGdL1z73XuaAJqHQEVI6BWSQPaf3T88vPfCah1BFSMgFolDWj8P8zpjryOgFpHQMUIqFXagK6uBPSs/yKg1hFQMQJqlTag5bnHfnMvATWPgIoRUKvEAY0RUPMIqBgBtYqAoj0CKkZArSKgaI+AihFQqwgo2iOgYgTUKgKK9gioGAG1ioCiPQIqRkCtIqBoj4CKEVCrCCjaI6BiBNQqAor2CKgYAbWKgKI9AipGQK0ioGiPgIoRUKsIKNojoGIE1CoCivYIqBgBtYqAoj0CKkZArSKgaI+AihFQqwgo2iOgYgTUKgKK9gioGAG1ioCiPQIqRkCtIqBoj4CKEVCrCCjaI6BiBNQqAor2CKgYAbWKgKI9AipGQK0ioGiPgIoRUKsIKNojoGIE1CoCivYIqBgBtYqAoj0CKkZArSKgaI+AihFQqwgo2iOgYgTUKgKK9gioGAG1ioCiPQIqRkCtIqBoj4CKhR5Q3+Cnh3TPO7zv9CHs6OpAo9yfgKI9Aio21gLacfbaIqAE1D4CKhZ+QB/fPWR/mYB2iYAGhoCKhR/Q5vrsXLr0EO2ZgBJQ+wio2FgL6KFDQAmofQRUjICOhoASUPsIqNhYC+hQ9p668syJ8x4ZurnvzgunTpr9k+eHVqxfMGXqBfcOlMs/9vPT+z3gJ+8vl7cs/f6kKbOvfqK+p1pGW+yp/hCevvLMSXOveTK93riT8t4b5049Y97qA+lNzffevOSMiXPvGSSgyI2Aio3RgD40IXlf/p705t1z07fpJ69PV9w14dKVi7xfWS5v9f7F5H6X+GXlcmW576n834R15VEC2rSn2kNY05ts7Lk3vtq0kxfPStcvHCiPuPea9Mr1BBS5EVCxsRnQvVP8D57qe+qS3uTmwYv8zIf29j1xmZ/8QrJiyubKmpt8b195cKa/Pb7b3gl+W/z/rnlu4MCWuX5uuXVAm/dUewg9Zz6wa+eD0/2kF8vDdnKZn/dE//77e/395eH3fr7Xz3/8pW0X+YYv4cDMuvX9AvcND+jMg60eLJcHinpkaGmAJ0BroNMZ2N1Z0JSfA91ZHgreT/wF8Uu+gQuTOG3y05KXmYNX+OuSFTfE13Ynt61Kv4e/2583WF7nZw3GV7b4nnLrgDbvqfYQJiU13XuW/1F52E6m+cfjK9fNiD+j2nzvH/lz+ypX+uc0BrThq1n7Svznauf+4QGdpXgUwJizt7NlxQd0xuwhcePT4M3yDye3rk3idLNfkS7e7M9PVjwVX0l/9rg9/R5+YfxCdP/O9Ivc0pDO5oA276n2EG5IL9zlZwzfydn+lsHawuZ7T/cPJFfuJ6DAmGc2oC2+he+f4NPvsJ9Pbl5ar9K0ZMXLDXedG6fzRd+bfn0Dz2/+2dVnjhrQ5j3VHsK69MJ233Ng2E7+y/tzVm54Kb296d6VVj6XbH2On4EiN34GKjYmfwb6gvd9ybUDyc0X+zOHXqSeN+JtodXx9/C3+SvjLVsvmVSJ2/SFLQI62GJPtYdQfft9f/Jqtmkngw/Oid9Puvih+HVo070rj/FAw2MciYCiPQIqNiYDeqCn+gr0heTmZf7m4Svqd93ZU6neHL+xcvGRHr/w1g17mtI5tHpPiz3VHkLtFejEgWE7qdj18+XTvb++POze/T3VV6DPE1DkRkDFxmRAyzOqPwNdl9x8q78svXHnhidHfrLzYn/7dn9O/J7Tuf6W5IZtzQFNPoG0ucWeagsqKFwAACAASURBVA9hZXrhTn/e8J089Ux8eeC+5Jv75nvPqP4MdC0BRW4EVGxsBnSZnx9/zzy4MLm55HuST7UPLvJrRgb0Hj//pvSl4SS/LblhRUNAn/M+aeClLfZUewjpu/B7zvQ3Dd/J2T174iu7vd8z/N7L/Hn9lSsD8wgociOgYuEHtP7XmOrvwu+c5Jc81ffM5VPSOF3qp927e2D7Uj9tz8iAvtTrz+pJCjjP/3BfeXD71d777UPrBib7BTvLe388tcWeag+h98wHdu168Gx/9v7hO1niF28f6H9ysZ9dHn7vF3r9hVteLi2eSECRGwEVCz+gdVPLtTw+GL+R4yc+6pMfhu6dV/0NoG3lkQEtX+794uTCoz2+56xJvuenZ/qeG4fW3Rq/Z+4nrE9XN+2p9hDunphufGzETp7ura5PXno23/uO9Mo9BBS5EVCxMRrQ8tNLZ05eVNpVfTu+b838KZPn/CR51TgioA97/1C6ZcuS6b0zlj1TfmzW1NuG1g0+MH/q5IseH1rduKfaQ9jxzBXTemeseHHkTsrbr/r+pDPmrNyZLm2+9+YlZ/TOWTvipxBVBBTtEVCx0AN6cE8mUQ0TAUV7BFRsbAb0zqvSP/axJv3ZY5AIKNojoGJjM6C3+Yvjd+H7Zvtr1Q8lNwKK9gio2NgM6AuT/bXP7X9ioe95sv1iowgo2iOgYmMzoOW1k9L3lVarH0h+BBTtEVCxMRrQ8s7rfzDnspUdPmiTCCjaI6BiYzWg4SOgaI+AihFQqwgo2iOgYgTUKgKK9gioGAG1ioCiPQIqRkCtIqBoj4CKEVCrCCjaI6BiBNQqAor2CKgYAbWKgKI9AipGQK0ioGiPgIoRUKsIKNojoGIE1CoCivYIqBgBtYqAoj0CKkZArSKgaI+AihFQqwgo2iOgYgTUKgKK9gioGAG1ioCiPQIqRkCtIqBoj4CKEVCrCCjaI6BiBNQqAor2CKgYAbWKgKI9AipGQK0ioGiPgIoRUKsIKNojoGIE1CoCivYIqBgBtYqAoj0CKkZArSKgaI+AihFQq8QBvepj0w4QUPMIqBgBtUob0HsPc24OATWPgIoRUKu0AT3NOffelQTUOgIqRkCtkgZ04JhKQN3pBNQ6AipGQK2SBnSTc0c5dwwBtY6AihFQq6QBXeHchHfFL0InH2w1AZUjoGIE1CppQGc4t/RzBNQ+AipGQK2SBvQrzt1xAQG1j4CKEVCrpAH9oDti670E1D4CKkZArVIGtO8I94FS6Z8IqHkEVIyAWqUM6FbnPlsqXUBAzSOgYgTUKmVA73TutFLpsb9zbvbBVhNQOQIqRkCtUgb06rSc935nyqaDrSagcgRUjIBapQzoXOcu72A1AZUjoGIE1CplQL1zqzpYTUDlCKgYAbVKGdBTnLu/g9UEVI6AihFQq5QB/ZA7fGsHqwmoHAEVI6BWKQN6lDu2k9UEVI6AihFQq4QB3eXcRztZTUDlCKgYAbVKGNBNzn2pk9UEVI6AihFQq4QBvdO50ztZTUDlCKgYAbVKGNBrnZvRyWoCKkdAxQioVcKAXuTc4k5WE1A5AipGQK0SBvRM51Z0spqAyhFQMQJqlTCg/+7c6k5WE1A5AipGQK0SBvQzzq3rZDUBlSOgYgTUKmFA/96984lOVhNQOQIqRkCtEgb03e64jlYTUDkCKkZArdIFdKdzH+9oNQGVI6BiBNQqXUA3OXdqR6sJqBwBFSOgVukCeodz3+toNQGVI6BiBNQqXUBXdPiLSARUj4CKEVCrdAFd0OEvIhFQPQIqRkCt0gV0Woe/iERA9QioGAG1ShfQrzu3pqPVBFSOgIoRUKt0AT3JuUc7Wk1A5QioGAG1ShfQ97ujOltNQOUIqBgBtUoW0L7D3T90tpqAyhFQMQJqlSygJef+tbPVBFSOgIoRUKtkAb3Hue90tpqAyhFQMQJqlSygHX+OnoDqEVAxAmqVLKDznFvS2WoCKkdAxQioVZLJuMW500537sbOVhNQOQIqRkCtkkzGWuc+/Rnn1na2moDKEVAxAmqVZjKOdscc497T4WICKkdAxQioVZrJ+Lir+GSHiwmoHAEVI6BWaSbjO3FAT+twMQGVI6BiBNQqzWQsiAO6qMPFBFSOgIoRUKs0k3FnHNB7O1xMQOUIqBgBtUo0GR9x7qRO1xJQOQIqRkCtEk3G7R/96B2driWgcgRUjIBaFcBkEFA5AipGQK0KYDIIqBwBFSOgVgUwGQRUjoCKEVCrApgMAipHQMUIqFUBTAYBlSOgYgTUqgAmg4DKEVAxAmpVAJNBQOUIqBgBtSqAySCgcgRUjIBaFcBkEFA5AipGQK0KYDIIqBwBFSOgVgUwGQRUjoCKEVCrApgMAipHQMUIqFUBTAYBlSOgYgTUqgAmg4DKEVAxAmpVAJNBQOUIqBgBtSqAySCgcgRUjIBaFcBkEFA5AipGQK0KYDIIqBwBFSOgVgUwGQRUjoCKEVCrApgMAipHQMUIqFUBTAYBlSOgYgTUqgAmg4DKEVAxAmpVAJNBQOUIqBgBtSqAySCgcgRUjIBaFcBkEFA5AipGQK0KYDIIqBwBFSOgVgUwGQRUjoCKEVCrApgMAipHQMUIqFUBTAYBlSOgYgTUqgAmg4DKEVAxAmpVAJNBQOUIqBgBtSqAySCgcgRUjIBaFcBkEFA5AipGQK0KYDIIqBwBFSOgVgUwGQRUjoCKEVCrApgMAipHQMUIqFUBTAYBlSOgYgTUqgAmg4DKEVAxAmpVAJNBQOUIqBgBtSqAySCgcgRUjIBaFcBkEFA5AipGQK0KYDIIqBwBFSOgVgUwGQRUjoCKEVCrApgMAipHQMUIqFUBTAYBlSOgYgTUqgAmg4DKEVAxAmpVAJNBQOUIqBgBtSqAySCgcgRUjIBaFcBkEFA5Alqwh2Zc33SdgFoVwGQQUDkCWrCPuHfc23idgFoVwGQQUDkCWqx7nXPnN24goFYFMBkEVI6AFuvWSkC/0biBgFoVwGQQUDkCWqw4oB9t3EBArQpgMgioHAEtVhzQYxo3EFCrApgMAipHQIsVB/Swxxo2EFCrApgMAipHQIsVB9Q1vg1PQK0KYDIIqBwBLVYS0FUNGwioVQFMBgGVI6DFSgL6o4YNBNSqACaDgMoR0GIlAb2oYQMBtSqAySCgcgS0WElAZzVsIKBWBTAZBFSOgBYrCeikhg0E1KoAJoOAyhHQYiUB/Y+GDQTUqgAmg4DKEdBiJQH9UsMGAmpVAJNBQOUIaLGSgJ7csIGAWhXAZBBQOQJarCSgJzZsIKBWBTAZBFSOgBYrCejfNmwgoFYFMBkEVI6AFisJ6NENGwioVQFMBgGVI6DFSgJ62Lb6BgJqVQCTQUDlCGixkoC69fUNBNSqACaDgMoR0GKlAb2vvoGAWhXAZBBQOQJarDSgt9U3EFCrApgMAipHQIuVBnRlfQMBtSqAySCgcgS0WGlAl9U3EFCrApgMAipHQIuVBvTS+gYCalUAk0FA5QhosdKALqhvIKBWBTAZBFSOgBYrDejs+gYCalUAk0FA5QhosSoBfbdz0+obCKhVAUwGAZUjoMWqBPRY5ybUNxBQqwKYDAIqR0CLVQnoPzT9RWUCalUAk0FA5QhosSoB/YRzp9Y3EFCrApgMAipHQItVCehnnfvX+gYCalUAk0FA5QhosSoB/WrlRWh9AwG1KoDJIKByBLRYlYB++zD3ofoGAmpVAJNBQOUIaLEqAf3uUY1/kp6AWhXAZBBQOQJarDig73fvrW8goFYFMBkEVI6AFisO6AnuiPoGAmpVAJNBQOUIaLHigH7SucdqGwioVQFMBgGVI6DFigP6OefW1jYQUKsCmAwCKkdAixUH9GvO/ay2gYBaFcBkEFA5AlqsOKCnNf5vehBQqwKYDAIqR0CLFQd0cuP/pgcBtSqAySCgcgS0WHFAZzp3bW0DAbUqgMkgoHIEtFhxQOc7d3ltAwG1KoDJIKByBLRYcUAvce6i2gYCalUAk0FA5QhoseKAXuPc+bUNBNSqACaDgMoR0GLFAV3p3Nm1DQTUqgAmg4DKEdBixQFd7dzE2gYCalUAk0FA5QhoseKA/ty579Q2EFCrApgMAipHQIsVB3SDc1+tbSCgVgUwGQRUjoAWKw5oybnP1TYQUKsCmAwCKkdAixUHdKdzJ9U2EFCrApgMAipHQIsVB7TfuRNrGwioVQFMBgGVI6DFigNafrc7vraBgFoVwGQQUDkCWqwkoMc2/G96EFCrApgMAipHQIuVBPRDDf+bHgTUqgAmg4DKEdBiJQE9yblN8ZUrznmUgJoVwGQQUDkCWqwkoKc491Dl8lLnegmoWQFMBgGVI6DFSgL6TefurFw+2blvEFCzApgMAipHQIuVBHSCczeVSpuPIKCWBTAZBFSOgBYrCeh055aXSjc6AmpZAJNBQOUIaLGSgM5z7opSaQEBNS2AySCgcgS0WElAL3FuQan0PQJqWgCTQUDlCGixkoBe69zsUumfCahpAUwGAZUjoMVKAnqzc1NLpeMJqGkBTAYBlSOgxUoCeo9zp5W2HE5ATQtgMgioHAEtVhLQR537t9LdjoCaFsBkEFA5AlqsJKDPOPcvpeUE1LYAJoOAyhHQYiUBfcm5f0o+xURADQtgMgioHAEtVhLQ8rvc8aWpBNS2ACaDgMoR0GKlAT3OHVX6KgG1LYDJIKByBLRYaUA/5tzjHyegtgUwGQRUjoAWKw3oKc7d/z4CalsAk0FA5QhosdKAnubcysPc4QTUsgAmg4DKEdBipQGd7ty5zh1HQC0LYDIIqBwBLVYa0EXOfcm5fyWglgUwGQRUjoAWKw3oj507xrkeAmpZAJNBQOUIaLHSgP4sfgPJLSCglgUwGQRUjoAWKw1oKQnoKgJqWQCTQUDlCGix0oD2v6PSzyPXElDLApgMAipHQIuVBrR8YiWgn36EgFoWwGQQUDkCWqxqQL9TCehZBNS0ACaDgMoR0GJVA/qjSkBXE1DTApgMAipHQItVDejeE91nBwioaQFMBgGVI6DFqga0vO++A2UCaloAk0FA5QhosYYCGiOgpgUwGQRUjoAWi4AGI4DJIKByBLRYBDQYAUwGAZUjoMUioMEIYDIIqBwBLRYBDUYAk0FA5QhosQhoMAKYDAIqR0CLRUCDEcBkEFA5AlosAhqMACaDgMoR0GIR0GAEMBkEVI6AFouABiOAySCgcgS0WAQ0GAFMBgGVI6DFIqDBCGAyCKgcAS0WAQ1GAJNBQOUIaLEIaDACmAwCKkdAi0VAgxHAZBBQOQJaLAIajAAmg4DKEdBiEdBgBDAZBFSOgBaLgAYjgMkgoHIEtFgENBgBTAYBlSOgxSKgwQhgMgioHAEtFgENRgCTQUDlCGixggmo9ztabN1748XTJ5176d39LW7b4f0hOvahd9DH1vpLJaDoAAEtVtgBvW+qT83YPPJGmwFNHxUBxSuDgBYr6IBu9H7Rhl37t981w096fsSthzCgoyQtDwKKVxIBLVbQAV3ofziYXNg30y8bcavNgO5curRMQPFKIaDFCjmgA5Nr2+7x54y4h82ApggoXhkEtFghB/RF71+qXnz2qqtH3IOAFo6AyhHQYtkP6OYlZ0yce89gtSp9d144ddLsnyQ/8ezr8auHR2b5rImTz79hT3w5idSP/fz0lgf85P3Nd28Wr16/YMrUC+4dGL6n+cn7VI80HOXaGb3n/GDjYKtD9t8yo/fs5fsG71s8fdK5Pz1QTnq47fIze6dfvb32qGoBbX48zV/qcAFMBgGVI6DFMh/QNenb7NenVdk9N706eX182xXeX77uQMPitRPTm6dtK1cjtdX7F5ObLol/Stp89yaV1XdNuHTlIu9XDt/Tjysb5y/aUlu6vnrbiIXxTq4457qrevwFS3ouXbnE+/h1ceXB90z9wZJpvveh8rCANj+e5i91hAAmg4DKEdBiWQ/o871+/uMvbbvIJ1UZvMjPfGhv3xOX+cmVx1neNaeyeeLi1U+nrwXL2yf6OZtf3vPwdH/2S9VIDc70t8c37Z3gtw2/e5PK6inxZ6Fu8r19w/fU/E3185P9xaW+XTd6v6nFIc/fVy7fVXlUj1dW/tRP6Ivv7C97ufJic7mf+GxzQJsfT/OXmjows259v3mV52FA/RjGuQGegEI93BzQ7/YPdDoDuzsrYJcB/ZE/t9Kgcv+cpCqb/LTk9eTgFf66+N++ny/qTT4GelvyOvRq//14cXn3FH/bUKRWpd/D3+3PGxxx90aV1Tck900O1Lyn5oAu93MH0oe2vMUhH61c2eP9zfV9eT89WTJ4ob+2OaDNj6f5S00d8HVru/svCeCQW9cc0NM7v+fezpZ1GdDp/oHk3/uTqtzsV6SbN/vzqwsObLj+/EpdLogfzqzq4vKN/tKhSG1Pv4dfGL8QbXH3msrqp+J/059ANu+pOaCz/UPJv09dv6bFIfel+9hS35dPGlyxzn+/OaDNj6f5S61+dWfUPTJgXvwKFFKDg+pHML4MfwUaz0BnT8GezgrYXUArr8CeSy48l1Rlaf312LSGVc+u6IlfDfb3VBdXUjW79nPGuXE6X/S9e0e9e6Ky+uXkQnygYXtqCujAhKHbyqMdsra+GtAN6ZLn/YSBpoA2PZ5hX+pIAfxwi5+ByvEz0GIZ/xnoC96nbxIdSKpysT9zdtV5lZdq9+8bWrfaTxyoLy4/6SfVarY6/h7+Nn9lefjdm9U+WBQfaNiemgJav63pSuMhhwf0yXTJ/vjFcGNAmx7PsC91pAAmg4DKEdBiGQ9o7RXe80lVlqU/W6w6wz88dHFf3Kb6y8H1flatZjt7KjfN8RvLw+/erCmgw/bUFNADPdX39ZseX9Mhhwe0+gmop31PX1NAmx7PsC91pAAmg4DKEdBiGQ9oeUb1B4Nrk6rc6i9LN+/c8GT8c81rhpbt8JMGGn4geZO/pF6zi/3t2/058fs+zXdv1hTQYXtq/hnoDL8u+fe5q1aMdshhAV2ZLllTb2z6/5sfT/OXOlIAk0FA5QhosawHdJk/L/5TdQPzkqqUfM8T8dbBRX5N8mmh+9IPMA0u8xeW4/exZw+9JX5rvWb3+Pk3pS/1mu/erDmgzXtqDugP/fzkqDfEyWx9yGEBnZS8ZN0zLf5BbWNAmx9P85c6UgCTQUDlCGixrAf0hV5/4ZaXS4snplW51E+7d/fA9qV+2p7KN72V1Cy4b8sLTz60yPv4o+jbJ/q5j728N/5Q5r56zV7q9Wf1pJ/7bLp7eenShjfCmgPavKdy0+8hPd3rF5f69q3piX+C0PqQwwLqz3pw166HzvZTdg/7IH3T4xn2pY4QwGQQUDkCWizrAS3fkb5LfU9alb3zqr+6E//eT3nvD4bexJ54d7L4od706pnxC7taEi/3fnF6qfnu3jd8nL45oM17Kp9TaWD9l5fun5De9p+jHnJYQFf1JEt615aHBbT58TR/qSMEMBkEVI6AFst8QONfEO+ds3aoSX1r5k+ZPOcn1ZeOgxuXLTxn4qzFNwz9cvuzy2dOPGvRqvrvwsce9tVPbg67+8EC2rSn8ubze8/aVF/79LJ5k2Ys2TT6IYcFdEfp8mkTWv8ufNOX0/ylDhfAZBBQOQJaLPsBfeX0N7yd/go6VH/KKYDJIKByBLRY4zmgL/gD7Rd1j4CiOAS0WOM5oDcvKOQwBBTFIaDFGs8BXdriz4K+AggoikNAizWeA1oQAoriENBiEdBXHAFFcQhosQhoMAKYDAIqR0CLRUCDEcBkEFA5AlosAhqMACaDgMoR0GIR0GAEMBkEVI6AFouABiOAySCgcgS0WAQ0GAFMBgGVI6DFIqDBCGAyCKgcAS0WAQ1GAJNBQOUIaLEIaDACmAwCKkdAi0VAgxHAZBBQOQJaLAIajAAmg4DKEdBiEdBgBDAZBFSOgBaLgAYjgMkgoHIEtFgENBgBTAYBlSOgxSKgwQhgMgioHAEtFgENRgCTQUDlCGixCGgwApgMAipHQIs1FND96wcIqG0BTAYBlSOgxaoG9KWPui8PEFDTApgMAipHQItVDegy59ydBNS0ACaDgMoR0GJVA/rdSkCnE1DTApgMAipHQItVDeiHKwH9DAE1LYDJIKByBLRYaUAH3lkJ6NEPE1DLApgMAipHQIuVBrTkYqsIqGUBTAYBlSOgxUoDek8S0IsIqGUBTAYBlSOgxUoDep1z73VuAgG1LIDJIKByBLRYaUAXO3eKc58loJYFMBkEVI6AFisN6DnOzXbueAJqWQCTQUDlCGix0oB+17nrD3OHE1DLApgMAipHQIuVBvSLzv38r+P3kQioXQFMBgGVI6DFSgP6Mece+xgBtS2AySCgcgS0WGlAj3NHlr5KQG0LYDIIqBwBLVYa0He740pTCKhtAUwGAZUjoMVKAvqScx8uLSCgtgUwGQRUjoAWKwnoM879c+laAmpbAJNBQOUIaLGSgG5w7t9KdxFQ2wKYDAIqR0CLlQT0Hue+U3r8cAJqWgCTQUDlCGixkoDe4tzUUulvCahpAUwGAZUjoMVKArrcue+XSp8hoKYFMBkEVI6AFisJ6KXOLSiVvktATQtgMgioHAEtVhLQec5dXirNJ6CmBTAZBFSOgBYrCeh055aXSisJqGkBTAYBlSOgxUoCOsG5m0qlTUcQUMsCmAwCKkdAi5UE9JvO3VG5/CkCalkAk0FA5QhosZKAnuLcg5XLlzk3gYCaFcBkEFA5AlqsJKAnObcxvrJo4iME1KwAJoOAyhHQYiUBPcEd8cTQBgJqVQCTQUDlCGixkoD+jXtvbQMBtSqAySCgcgS0WElAj3TH1zYQUKsCmAwCKkdAixUHtN+5E2obCKhVAUwGAZUjoMWKA7rTuZNqGwioVQFMBgGVI6DFigNacu5ztQ0E1KoAJoOAyhHQYsUB3eDcV2sbCKhVAUwGAZUjoMWKA/rz+O8pDyGgVgUwGQRUjoAWKw7oaud6axsIqFUBTAYBlSOgxYoDeoNzZ9c2EFCrApgMAipHQIsVB/Qa586vbSCgVgUwGQRUjoAWKw7oJc4trG0goFYFMBkEVI6AFisO6PzkD9JXEVCrApgMAipHQIsVB3SWc8tqGwioVQFMBgGVI6DFigM6xbnraxsIqFUBTAYBlSOgxYoD+l3nbq1tIKBWBTAZBFSOgBYrDujXnLu7toGAWhXAZBBQOQJarDign3dubW0DAbUqgMkgoHIEtFhxQE9ybnNtAwG1KoDJIKByBLRYcUBPcIfX/hc9CKhZAUwGAZUjoMWKA/p+9576BgJqVQCTQUDlCGix4oAe7d5f30BArQpgMgioHAEtViWg3znMfai+gYBaFcBkEFA5AlqsSkC/6tzH6xsIqFUBTAYBlSOgxaoE9LPO/Ut9AwG1KoDJIKByBLRYlYB+wrlT6xsIqFUBTAYBlSOgxaoE9B+c+/f6BgJqVQCTQUDlCGixKgE91rme+gYCalUAk0FA5QhosSoBfbdzZ9Q3EFCrApgMAipHQItVCWjF7PoGAmpVAJNBQOUIaLHSgC6obyCgVgUwGQRUjoAWKw3opfUNBNSqACaDgMoR0GKlAb2mvoGAWhXAZBBQOQJarDSgK+sbCKhVAUwGAZUjoMVKA3pbfQMBtSqAySCgcgS0WGlA761vIKBWBTAZBFSOgBYrDej6+gYCalUAk0FA5QhosZKAHratvoGAWhXAZBBQOQJarCSgRzdsIKBWBTAZBFSOgBYrCWjDH6QnoGYFMBkEVI6AFisJ6IkNGwioVQFMBgGVI6DFSgJ6csMGAmpVAJNBQOUIaLGSgH6xYQMBtSqAySCgcgS0WElA/6NhAwG1KoDJIKByBLRYSUAnNmwgoFYFMBkEVI6AFisJ6KyGDQTUqgAmg4DKEdBiJQFd2LCBgFoVwGQQUDkCWqwkoD9q2EBArQpgMgioHAEtVhLQVQ0bCKhVAUwGAZUjoMVKAnpPwwYCalUAk0FA5QhosZKAPtawgYBaFcBkEFA5AlqsOKB/2biBgFoVwGQQUDkCWqw4oB9p3EBArQpgMgioHAEtVhzQrzduIKBWBTAZBFSOgBbr7kpAZzduIKBWBTAZBFSOgBbsRHfEzxqvE1CrApgMAipHQAt231krmq4TUKsCmAwCKkdAxQioVQFMBgGVI6BiBNSqACaDgMoRUDECalUAk0FA5QioGAG1KoDJIKByBFSMgFoVwGQQUDkCKkZArQpgMgioHAEVI6BWBTAZBFSOgIoRUKsCmAwCKkdAxQioVQFMBgGVI6BiBNSqACaDgMoRUDECalUAk0FA5QioGAG1KoDJIKByBFSMgFoVwGQQUDkCKkZArQpgMgioHAEVI6BWBTAZBFSOgIoRUKsCmAwCKkdAxQioVQFMBgGVI6BiBNSqACaDgMoRUDECalUAk0FA5QioGAG1KoDJIKByBFSMgFoVwGQQUDkCKkZArQpgMgioHAEVI6BWBTAZBFSOgIoRUKsCmAwCKkdAxQioVQFMBgGVI6BiBNSqACaDgMoRUDECalUAk0FA5QioGAG1KoDJIKByBFSMgFoVwGQQUDkCKkZArQpgMgioHAEVI6BWBTAZBFSOgIoRUKsCmAwCKkdAxQioVQFMBgGVI6BiBNSqACaDgMoRUDECalUAk0FA5QioGAG1KoDJ2LRu3Tb1Yxjnnt2hfgTj3JZ16x7vbGVnY09Ax5Efed/haQGMTfd5f+uh3B8BHUcIKMY7AorcCCjGOwKK3AgoxjsCitwIKMY7AorcCCjGOwKK3AgoxjsCitwIKMY7AorcCCjGOwKK3AgoxjsCitwIKMY7Aorc9u/bN6B+DIBS/759fYdyfwQUAHIioACQEwEFgJwIKADkREABICcCCgA5EVAAyImAAkBOBHSsutYPXepfM2/KlHm397e/BowJeU7+fKNAQMeo/Wf46qUDC33iwgPtrgFjQp6TP+coENCxXwvtZQAACTJJREFUad/l3lcvLvdT7nr++bum+OXtrgFjQa6TP+coENAxaP+Nl0z2Q+fQCxP8o/G/j/oJLx78GhC+nCd/3lEgoGPQ7vSbkfTK7X7WYPzv4Cx/58GvAeHLefLnHQUCOgb1P/PMM+uGzqEr/Q3phRv8Dw9+DQhfzpM/7ygQ0DFqx9A5tMA/kF54wF908GvA2JDj5M87CgR0jKqdQ7P8xvTCRj/74NeAsSHHyZ93FAjoGFU7hyb7J9MLT/qpB78GjA05Tv68o0BAx6iR51DJTzr4NWBsyHHy5x0FAjpGNXwXsym9sNHPOvg1YGzIcfLnHQUCOkbVzqGF/sH0wgN+4cGvAWNDjpM/7ygQ0DFqR/2THDelF27ySw9+DRgbcpz8eUeBgI5RO+qfJT4vvXCeX3Pwa8DYkOPkzzsKBHSM2tHw22xb43+3Df3G2qjXgLEhx8mfdxQI6BhVO4fKK/z0DX0HHp3uf9zuGjAm5Dn5c44CAR2j6ufQgQXeT+z1/qK+dteAMSHPyZ9zFAjoGFU/h8r9t8+dPHneHQPtrwFjQa6TP98oEFAAyImAAkBOBBQAciKgAJATAQWAnAgoAOREQAEgJwIKADkRUADIiYAiCB+JhvtEt7t8XbKb5R0ubfinKC0Pd0P61Rf5ODA6nggEoXVAo+h1+XfZSUCrB2gOaFdHzYCA2scTgSCcf3Lqv0VR9dLccvcBXb58+dMHXWMvoNsrj/nVBNQKnggEpTkq3Qa07Rp7AT3oDSgaTwSC0tyOpUuvOlS7aq16gOaAdnXUDAiofTwRCMohbEeGXRl6E0nwODA6nggEpbkdDWF7+uv/+w1v/N9ffqJcvuXD//O1v3Z4T+1/mWHHxHe99XW/f9wlL4+2q4Pcvbmcwzr60rzj//gXf8194s7aLkccqu/iD/zP177tfWfsrV7vv/yf3vaGX/xf/3Rpf/3Q/Re+562v+4Pjrx76S5TP+r948xv+6Ns7Rj8OATWDJwJBGS2gd/xq9Eu/GEXRm+74QhT99/jS7zybrrn8l6tv2//ufaPs6iB3P1hAr/mNoQ8EfGjfKIfa8KfVDb91fXJ9/f8ZusufrBs69M5jouit/62y6X37k03XvTld8Zs/H/U4BNQMnggEZbSAvvkPV/cPrPjVKPqF6O239Q2s+JUo+myy5NJKeP7g5G9/8E1R9MY1rXd1kLsfJKB3vCqK/sc/f+9rR1YO8OXWh9r661H0y8d/+19+J4pe/3Dl+ubKbl995KlfOeo1UfSWTdVdHfeqM54svzS7cunb8ZZbX11Z/Jf/9jdvjN72qlGOQ0Dt4IlAUEYL6G8l33JfUInMb++ML82Noj+O/33il6JXzYy/X376Q1H0pwda7uogdz9IQI+Pok8l+1tc6e9gy0P9dRS5bZV/9306ij5aLg8eFUVvuyO+S/zK9KjBdFfVT6KeEUV/VvnnwP8TRW9fX7mw5U+qH/cccRwCaghPBIIyWkDnJ9e3VpqzaOjS6+N/vxlFE9O1B1wUXdJyVwe5+0EC+mtR9GR6//83ih5tdajbK68/X0g2vPx70VsHy7dUXn9uTFdsrrzO/K90V8dVt6Qfjroyit7wVLJhxxuj1schoIbwRCAoowU0LctApTlbhy4lH9Z8W/SWof+ZxRVR9M8td3WQux8koL8bRXc0PbQRhzo1ir5V3XD6b/zGhvKXo+iUocWfj6KvprtaVm484iei6OvVJV+LWh+HgBrCE4GgjBbQ6rvaUfSawaFLcY+eiKIPDi3eEUX/p+WuRr/7wQL6kSh667nb67sbeai/iKLVjQesXF85dPknUXREuqt15cYj/n4U3VrdcGvU8jgj/iNAiCcCQRn1Y0yp+i8JpZfubP79+d9tuavR736wgD71O/Ee/+xzCx5L7zHyUJVvvp9oPOBvRtGmocvro+i3010NfeQpPeIv1b5hj4vc6jgj/iNAiCcCQckY0Buaq/bmlrvKF9Dyru+9tfqRpLNebnmoV0W/0N94wDdG0b6hy3uj6L83fznpESv3Gfo8aH/U8jgj/iNAiCcCQckY0Dui6BOP1K1ruaucAS2XB342+fhfidP2h8+1OtQvR9GOxgP+j4ZXoJuj6DfLLQL6lvor0KeilscZ8R8BQjwRCErGgG6JomPb7ip3QGOD9/3fX42iT7c61B9E0d3ViwMvvdRX/vMoumHotpVR9OflFgH90/rPQG+PWh5nxH8ECPFEICgZAzr4W9Fbh37IuGnGjNta7ipfQC+Y/Z/VtXdH0e+1OtSnomhqdcOcKPpR+YtR9KWhg38pvTwioP9Sfxf+tKjlcUb8R4AQTwSCkjGg8UeHJlS3nBBFP225q3wB/bPot6vv2D8RRW9tdairo+hXdyXXB4+OfuG58qooes3mdMHGV0fRqnKLgC6Pojek38O/+Kao9XEIqCE8EQhK1oBufF302nPjt2X6T4+itw+23FW+gH4sinyyv5c+EkV/3+pQ/X8SRUeWKhsGeioXKhl9R+UlZPJN/R1vi6J3D5ZbBLT/D6u/ifTkYdXfRBpxHAJqCE8EgpI1oOXzKxn6o09/7wuVZL35gda7yhfQFZUdv/1L/tufeEsUvfruloe6/xej6E0f+M6XKyF93V2V6xveXHkNevRXvnJ05fXnr2we9uVUj7jmtVH0+mO++oE3Rce9dpTjEFA7eCIQlMwBLZ//+uoHi36/6VPt3b+J1PsLtc98Xj/Koe76neqGt1yXXH/kj4fu8mcjfilz6NjXvyVd8ec7Rj0OATWDJwJByR7Q8tZv/X9vfv0f/f1FfeVmXb8L/7MvHPm21/66O/Gq2p5HHGrf7L/+9df83nG9z1Sv911y4u++4Q2/++ErGv4e6PBjP9vzF2963dt79ox+HAJqBk8ExqtwMxTuIx9zeCIwXoWboXAf+ZjDE4HxKtwMhfvIxxyeCIxX4WYo3Ec+5vBEYLyqZOiRRx7ZrX4YGe2pPObXEFAreCIwXr0u+WjQcvXDyKj6R5/UDwMpngiMVwQUXeOJAICcCCgA5ERAASAnAgoAORFQAMiJgAJATgQUAHIioACQEwEFgJwIKADkREABICcCCgA5EVAAyImAAkBOBBQAciKgAJATAQWAnAgoAOREQAEgJwIKADkRUADIiYACQE7/P+GgfxCd7GCfAAAAAElFTkSuQmCC" title alt width="672" /></p> +<p>We see that using the compiled model is by a factor of 13.6 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 13.734549 -## Eigenvalue based 1.272901 +## deSolve, not compiled 13.581348 +## Eigenvalue based 1.256428 ## deSolve, compiled 1.000000</code></pre> </div> <div id="benchmark-for-a-model-that-can-not-be-solved-with-eigenvalues" class="section level2"> @@ -136,19 +136,19 @@ smb.1["median"]/smb.1["deSolve, compiled", "median" smb.2 <- summary(mb.2) print(mb.2)</code></pre> <pre><code>## Unit: seconds -## expr min lq mean median uq -## deSolve, not compiled 20.475764 20.494740 20.507391 20.513716 20.523205 -## deSolve, compiled 1.244022 1.244327 1.261983 1.244631 1.270963 +## expr min lq mean median uq +## deSolve, not compiled 21.324080 21.368031 21.460777 21.411981 21.52913 +## deSolve, compiled 1.376772 1.414208 1.461651 1.451643 1.50409 ## max neval cld -## 20.532695 3 b -## 1.297295 3 a</code></pre> +## 21.646269 3 b +## 1.556538 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,iVBORw0KGgoAAAANSUhEUgAABUAAAAPACAMAAADDuCPrAAAC4lBMVEUAAAABAQECAgIDAwMEBAQFBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQVFRUWFhYXFxcYGBgZGRkaGhobGxscHBwdHR0eHh4fHx8gICAhISEiIiIjIyMkJCQlJSUmJiYnJycoKCgpKSkqKiorKyssLCwtLS0uLi4wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6Ojo7Ozs8PDw9PT0+Pj5AQEBBQUFDQ0NERERFRUVGRkZJSUlKSkpLS0tMTExNTU1PT09RUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVlbW1tcXFxeXl5fX19gYGBhYWFiYmJjY2NlZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1ubm5vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8fHx9fX1+fn5/f3+AgICBgYGCgoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uMjIyNjY2Ojo6Pj4+QkJCRkZGSkpKTk5OUlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2enp6fn5+goKChoaGioqKjo6OkpKSlpaWmpqanp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKzs7O0tLS1tbW2tra3t7e4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXGxsbHx8fIyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY2NjZ2dna2trb29vc3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr6+vs7Ozt7e3u7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+/v7////Td9k+AAAACXBIWXMAAB2HAAAdhwGP5fFlAAAgAElEQVR4nO3debzcdX3v8Z91QbStWEt3pb23rV297dR6W+ttbYu3y+1tbb2t195eW5e61CoqfE/OSQgBkhADSSBhSVgEI9IYtEIA9bIEoYhggEBClCUBk0BYspGzzP93Zn5zlt9n5shnJp+T83uf3+v5R8+cc34zv5kO79djzhazOgCgL9ls3wEAUEVAAaBPBBQA+kRAAaBPBBQA+kRAAaBPBBQA+kRAAaBPBBQA+kRAAaBPBBQA+kRAAaBPBBQA+kRAAaBPBBQA+kRAAaBPBBQA+kRAAaBPBBQA+kRAAaBPBBQA+kRAAaBPBBQA+kRAAaBPBBQA+kRAAaBPBBQA+kRAAaBPBBQA+kRAAaBPBBQA+kRAAaBPBBQA+kRAAaBPBBQA+kRAAaBPBBQA+kRAo+zfB8g6/MKh2b4LJeObPQGN8tROQNZYfWS270LJ+GZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0Si8B/fZ3ZuxpB/pBQC3f7AlolB4Cel7tbZfO3BMP9I6AWr7ZE9AoPQT0PbXam784c8880DMCavlmT0Cj9BDQ/12r1U7mq3iUCAG1fLMnoFF6DGjt3Jl76oFeEVDLN3sCGqXXgL6dl6AoDwJq+WZ/7AOa0p4uHz1wwyWL5597xZ0jXT63J6UZvk89nbH7A+gxoH9Tq62Zuece6BEBtXypKElA716Yckt3dH7yWAY0P9fMB/S6Wu0dj7beuWXlTTP0XwDgRkAtXzDKEdDtKa3d9twLu+9YmuY/3fHZuRjQve+u1S5pXNz+oVrtLffN2H8EgA8BtXzBKEdAL0qfH2tdOHhO2tDx2WMZ0GfXr3+xM8YE9NZa7Y+27dzyP5vfDr1l5v4rAFwIqOULRikCOrpg4mN3pSUd15iL3wPdW/+HWu0fN55cI6AoAwJq+VJRioA+k9Kh9sW9V1/TcY25GdDvvLUVzz/+XwQUs4+AWr5UHMOA7rjsjKFVd421+zP8jTUL5y//Sus7nsMD6VZz8FMblw0tWLFpf/NyK2dfSqvzz9ybFrxQvLr11BeXDi65dPtYtxsa+drSwbM3Hhy7e93i+ed+/Ui91cPHrjxzcPE1uyfONRHQ4lmKD6DjrP6nphXQ+k2/0+jnOx//FAHF7COglq9qxy6gm/Mfs1+X9+f5Vfm7Cx5sfu6qlK7cemTKwfcN5Z9e9Fi9nbNHU3qm9anLm98lLV696MH2da/vdkNXLfny1QPpgssGrrj+spSar3Ybd2lg4aWXLUqDW+omoMWzFB9Ah54DWt+28JT1L9QJKEqAgFq+rB2zgD49mFZ/99BjF6dWf8YuTudsOTD8+GfTgn2NTz63svHhoXW3Ppm/aqzvHkordxzef//idPahds7Gzkm3NT91YF56zF69eKIF6ZKdw8/dkNLDXW5oxcF6/Y7Gub7bOPLrad5wM6Dps4cbLzY3pqG9xYAWz1J8ALnhtZO2H3b7hzyguUZAv+2/KjAjGv+5z/ZdKJdnfF07ZgH9Qjq3Uav6yMpWfx5Oi1p3cOyq9OXm2+FvrR1s/RroLa3XodekzzQPrj9/erplPGc35l/D35nOG+u4+lQb06rR/IQbu9zQQ4139qf01dYHW3clpcWtQ8bWpC8WA1o8S/EB5I6kSff5/3/xXhPQLi+jAcymA77DjllAF6d7W2/vafXnq+na/MM70or2AUe2Xbei0aELmnd8Wfvg+g3pivGc7c6/hr+o+UK0y9UnLE9bWm+fuG5zlxtqvACtj6X0SD1/2wroLfkhW9NnigEtnqX4ANr3mYACc1LJAtpIzVOtC0+1+rN+MjyLphy199qB5uvGkYH2wY2oLZ/4juSqZjqfSYMHpr160+i88evWp7uhiZ+jtwO6LT/k6TRvtBDQwlnMA8iNHZy05wm39+QBHb3jxu+1Anqb/6rAjBirj8z2XSiXXb6wHauA7ksp/yHRkVZ/LklnLm87r/Gi7p6D48fdmoZGJw+u70rzJ7p3a/Nr+FvS5+r26t1PVHxn6g3ZgLb/X/VC8yXu1IAWzmIeQKeef4j0xN/Uam85b5QfIqEE+CGS5QvbsQroxGvBp1v92ZB/F7LtjHT/+MWDzYpNvnB8MC2b6N6zA41PrUzb6/bqBUcG0uS3f7vfkA3oA/khT6aB4UJAC2cxD6BTrwE9+Jet3wP98EcJKGYfAbWmK0zRMfse6NL2txDva/Xn5vTZ/MPPbtvV/L7mv40ftifNH53yrcub0uWT3bsk3bY7LWn+hKh4dXuira23T1197XQ3ZAJ6fX7I5snG5v+3eJbiA+jUa0AX12p//Mk385dIKAUCavm6dswCuiGd1/yn6kbPb/VnZxp4vPnRsbVpc+v3iu7Of4FpbENaU2/+xHv5+A/Pb57s3l1p9U35i8Li1Ys+n1a3bmtTM5ndb8gEdH7rJev+Rc1vv04NaPEsxQfQqceAPvaW2m/fvPPq3yGgKAMCavm6dswCum8wrXnk8M51Q3l/rkiLvvn86O71adH+xpfHjShdePcj+3ZtWZtS80fSu4fSqu8cPtD89c2Dk907NJjOGsh/77Nw9fr69fsnT/TkYFq3c/jg5oHm9wW635AJaDrr2889t+XsdPrz5hfpC2cxD6BDjwEdqNXmNS7e+PZGQO+csf8IAB8Cavm6duz+Eun2/OfZd+X9OXB++498mn8hVD9w6fiPu4fubB28ZTB/98zmS8CJv0y/MqV1+aXi1VOa+uv098zLP/fv096QCeiNA61DBu+rm4AWz1J8AB16C+iWt9R+t/Wv2D1wxrsWzdR/AoAXAbV8WTu2fws/uPK+8XoNb159+oKVX2m/dBzbvuGiJUPL1m0a/+P2vRvPGTpr7Y2TfwvfdH9q/46nuXoxoPUnN5w/f+llD09/Qyage3ZeuWhe97+FL9zJ4gOwegvoKbXawMw990CPCKjlq9qc+N9EGknOv7vqZpoe9qy3gL659tvfmrnnHugRAbV8s58TAd2Xjrz4QdOZnYDWah+fuace6BUBtXyznxMB/eqFR3HlWQrof/3mzD31QK8IqOWb/ZwI6Pqu/yyo0ywF9KyZe+aBnhFQyzf7ORHQozILAf1YrfZPj87cMw/0jIBavtkT0FkI6EOrvvDYzD3xQO8IqOWbPQGN0kNAgbIhoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFMAJq+WZPQKMQUAirQkC3be3laN/sCWgUAgphFQjoN9/2lk09HO6bPQGNQkAhrAIBvaJWO7uHw32zJ6BRCCiEVSCga2u1hT0c7ps9AY1CQCGsAgG9uFab38PhvtkT0CgEFMIqENALa7XBHg73zZ6ARiGgEFaBgK6p1eb1cLhv9gQ0CgGFsAoE9IJaLfVwuG/2BDQKAYWwCgR0Va12Wg+H+2ZPQKMQUAirQEBX1mqn9nC4b/YENAoBhbBqBPRTPRzumz0BjUJAIawCAT2vVvtkD4f7Zk9AoxBQCKtAQJfzCrTMCCiEVSOgp/RwuG/2BDQKAYWwCgR0Wa32iR4O982egEYhoBBWgYCeU6t9vIfDfbMnoFEIKIQRUMs3ewIahYBCWDUC+rEeDvfNnoBGIaAQVoGALqnV/qWHw32zJ6BRCCiEVSCgiwlomRFQCKtAQM8moGVGQCGsAgE9i4CWGQGFsAoE9EwCWmYEFMIqENBFBLTMCCiEVSCgZxDQMiOgEFaBgPIKtNQIKIRVIKALCWiZEVAII6CWb/YENAoBhbAKBHQBAS0zAgphBNTyzZ6ARiGgEFaBgM4noGVGQCGsAgEdIqBlRkAhrAIBHSSgZUZAIawCAZ1HQMuMgEIYAbV8syegUQgohBFQyzd7AhqFgEJYBQKaCGiZEVAIq0BATyOgZUZAIYyAWr7ZE9AoBBTCKhDQUwlomRFQCKtAQD9NQMuMgEJYBQL6KQJaZgQUwioQ0E8S0DIjoBBWgYDyCrTUCCiEVSCgpxDQMiOgEFaBgH6CgJYZAYUwAmr5Zk9AoxBQCKtAQD9OQMuMgEJYBQL6rwS0zAgohFUgoB8joGVGQCGsAgH9FwJaZgQUwgio5Zs9AY1CQCGMgFq+2RPQKAQUwioQ0I8S0DIjoBBWgYB+hICWGQGFsAoE9EMEtMwIKIQRUMs3ewIahYBCWAUC+s8EtMwIKIRVIKAfJKBlRkAhrAIB/QABLTMCCmEE1PLNnoBGIaAQVoGAvo+AlhkBhbAKBPSvCWiZEVAIm/sB/c5vEdAyI6AQNvcD+lCNgJYZAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9Ao5ZLuOEAABlXSURBVBBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2BDQKAYUwAmr5Zk9AoxBQCCOglm/2UQFNaU+Xjx644ZLF88+94s6RLp/bk1LQueN93/vW/aESUCgjoJYvFTMa0LsXptzSHZ2fLGdA83tFQFExBNTyBWMmA7o9pbXbnnth9x1L0/ynOz4bGNBpktYPAopKIqCWLxgzGdCL0ufHWhcOnpM2dHy2nAF9dv36OgFF5RBQyxeMGQzo6IKJj92VlnRco5wBzRFQVAwBtXypmMGAPpPSofbFvVdf03ENAgqUBgG1fKk46oDuuOyMoVV3jbWrMvyNNQvnL/9K6zuewwPpVhuZjcuGFqzYtL95uRWpL6XV+WfuTQteKF69qHn0gxeevvCCb47aW1rd+jnVA1PO8sWlg0su3T7W7ZQjX1s6ePbGg2N3r1s8/9yvH6m3evjYlWcOLr5m98S9mgho8f4UH6pFQCGMgFq+/h1tQDfnP2a/Lq/K86vydxc82PzcVSldufXIlIPvG8o/veixejtSj6b0TOtTlze/S1q8ekHj6DvmXXH92pSut7f0pcYHV699ZOLQB9uf6ziweSNXLfny1QPpgssGrrj+spSar4sbd35g4aWXLUqDW+omoMX7U3yoHQgohBFQyxfAowzo04Np9XcPPXZxalVl7OJ0zpYDw49/Ni3Y1/jkcysbHx5ad+uT+WvB+u6htHLH4f33L05nH2pHauycdFvzUwfmpcfs1QsaR5/e/F2om9LgsL2l4hfVTy9Il+wcfu6GlB7ucsoVB+v1Oxr36ruNI7+e5g03r5w+e7jxYnNjGtpbDGjx/hQfau5ImnTf0f1/EsBMeq4R0FP9hx/wHXaUAf1COrfRoPrIylZVHk6LWq8nx65KX26+Hf7W2sHWr4He0nodek36TPPg+vOnp1vGI3Vj/jX8nem8sY6rT9U4elPruq0TFW+pGNCNadVoftc2djnlQ4139qf01cnbSmlx65CxNemLxYAW70/xoeYIKCCilAFdnO5tvb2nVZWvpmvzD+9IK9oHHNl23YpGXS5o3p1l7YPrN6QrxiO1O/8a/qLmC9EuV5/QOPqJ5tv8O5DFWyoGdHna0nr7xHWbu5zyYH4bj0zeVmo1uGFr+kwxoMX7U3yo7Ud3zqQHRwBZjf+aZ/suzLB9jYB+2n/4874CHl1AG6/AnmpdeKpVlfWTr8cWTTlq77UDzVeDIwPtgxupWj7xfcZVzXQ+kwYPTHv1lsbRh1sXmicyt1QI6Oi88c/VpzvlxPHtgG7LD3k6zRstBLRwf8xD7cT3QCGM74Fajv7Vjzag+1LKf0h0pFWVS9KZy9vOa7xUu+fg+HG3pqHRyYPru9L8iZrd2vwa/pb0ubq9etHELxY1T2RuqRDQyc8V3pl6ShvQXfkhLzRfDE8NaOH+mIfaiYBCGAG1HP2rH21AJ17hPd2qyob8e4ttZ6T7xy8ebLZp8uXgg2nZRM2eHWh8amXaXrdXLyoE1NxSIaBHBto/1y/cv8IpbUDbvwH1ZBoYLgS0cH/MQ+1EQCGMgFrTtqjgKL8HurT9jcH7WlW5OX02//Cz23Y1v6/5b+OH7UnzR6d8Q/KmdPlkzS5Jt+1OS5o/9ylevagQUHNLxe+BLk1bW2+fuvra6U5pAnp9fsjmycbm/7d4f4oPtRMBhTACak0bvYKjDOiGdF7z28+j57eqsjMNPN786NjatLn120J357/ANLYhrak3f469fPxH4jdP1uyutPqm/KVe8epFxYAWb6kY0M+n1a2zbmoms/spTUDnt16y7l/U/Ebt1IAW70/xoXYioBBGQC1fAY8yoPsG05pHDu9cN5RX5Yq06JvPj+5enxbtb3zR20jNhXc/sm/XlrUpNX8VffdQWvWdwweav5R5cLJmhwbTWQP5730Wrl5fv37/5ImKAS3eUr3wd0hPDqZ1O4cPbh5ofgeh+ylNQNNZ337uuS1np9OfN79IX7g/5qF2IKAQRkAtXwGP9i+Rbs9/Sn1XXpUD57f/dKf5dz/1A5eO/xB76M7WwVsG83fPfLw+JYlXprQuv1S8ekpTfp2+GNDiLdWXNBo4+cdL98zLP/fv057SBPTGgdYhg/fVTUCL96f4UDsQUAgjoJYvgBF/Cz+48r7xJg1vXn36gpVfab90HNu+4aIlQ8vWbRr/4/a9G88ZOmvtjZN/C990f2r/5qa5+vcLaOGW6jtWDJ718OSxT244f/7Syx6e/pQmoHt2XrloXve/hS88nOJDtQgohBFQyxG/eqn/N5FGpvw4fQZF/VNOBBTCCKjlm32JA7ovHXnxg44eAQUIaAff7Esc0K9eeExOQ0ABAtrBN/sSB3R9l38WdAYQUICAdvDNvsQBPUYIKEBAO/hmT0AJKEBAO/hmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhTACavlmT0CjEFAII6CWb/YENAoBhbC5H9DvvpmAlhkBhbC5H9Cdf0dAy4yAQlgFAvo+AlpmBBTCKhDQDxDQMiOgEEZALd/sCWgUAgphFQjoBwlomRFQCKtAQP+ZgJYZAYWwCgT0QwS0zAgohFUgoB8moGVGQCGsAgH9CAEtMwIKYRUI6EcJaJkRUAirQED/hYCWGQGFMAJq+WZPQKMQUAgjoJZv9gQ0CgGFsAoE9F8JaJkRUAgjoJZv9gQ0CgGFsAoE9OMEtMwIKIRVIKCfIKBlRkAhjIBavtkT0CgEFMIqENBPEtAyI6AQVoGAfoqAlhkBhbAKBJRXoKVGQCGsAgHlFWipEVAIq0BAP01Ay4yAQlgFAnoqAS0zAgphFQjoaQS0zAgohBFQyzd7AhqFgEJYBQKaCGiZEVAIq0BA5xHQMiOgEFaBgA4S0DIjoBBWgYDyCrTUCCiEVSCgvAItNQIKYRUI6BABLTMCCmEVCOh8AlpmBBTCKhDQBQS0zAgohBFQyzd7AhqFgEJYBQK6kICWGQGFMAJq+WZPQKMQUAirQEAXEdAyI6AQVoGAnkFAy4yAQlgFAsor0FIjoBBWgYCeRUDLjIBCGAG1fLMnoFEIKIRVIKBnE9AyI6AQVoGALiagZUZAIawCAV1CQMuMgEJYBQJ6Tq32sR4O982egEYhoBBWgYAurdU+3sPhvtkT0CgEFMIqENBlBLTMCCiEVSOgn+jhcN/sCWgUAgphFQjo8lrtlB4O982egEYhoBBWjYB+qofDfbMnoFEIKIRVIKArarVP9nC4b/YENAoBhbBqBPTTPRzumz0BjUJAIawCAV1Zq53aw+G+2RPQKAQUwioQ0PNrtdN6ONw3ewIahYBCWAUCekGtNtDD4b7ZE9AoBBTCKhDQNbXaYA+H+2ZPQKMQUAirQEAvJKBlRkAhrAIBvbhWW9DD4b7ZE9AoBBTCKhDQdbXawh4O982egEYhoBBWgYBeWast6eFw3+wJaBQCCmEVCOh9//33b+3hcN/sCWgUAgphFQhoj3yzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQhgBtXyzJ6BRCCiEEVDLN3sCGoWAQtjWrQ/N9l0oGd/sCSiA+vy0dLbvgiQCCoCA9omAAiCgfSKgAAhonwgoAALaJwIKgID2iYACIKB9IqAACGifCCgAAtonAgqAgPaJgAIgoH0ioADqhw4emu27IImAAkCfCCgA9ImAAkCfCCgA9ImAAkCfCCgA9ImAAkCfCCgA9ImAApX2xTR+aWTz+aeffv5tI7N4Z+QQUKDKXjgjtS8duSi1rDkym/dHDAEFKuzglSm1L25Mp9/x9NN3nJ42zuYdEkNAgap64YbLF6TxgO6blx5qvn0ozXtmFu+TGAIKVNXz+dfs+Tu3pWVjzbdjy9I3ZvE+iSGgQFWNfO9739s6HtDPpU35hU3p87N2j+QQUKDK9owH9MJ0b37h3nTxbN0ZPQQUqLKJgC5L2/ML29PyWbs3cggoUGUTAV2QduUXdqWFs3Zv5BBQoMo6A7ozzZ+1eyOHgAJVNuVL+IfzC9vTslm7N3IIKFBlEwG9KH07v3BvumjW7o0cAgpU2Z7JX2O6Kb9wU1o/a/dGDgEFqmzP5C/Sn5dfOC9tnrV7I4eAAlW2Z8qfcj7afPsYf8rZAwIKVNlEQOvXpsXbho88tDh9aTbvjxgCClTZZECPXJjS0GBKFw/P5v0RQ0CBKpsMaH3ktlULFpx/++gs3hs5BBQA+kRAAaBPBBQA+kRAAaBPBBQA+kRAAaBPBBQA+kRAAaBPBBQA+kRAIe+dmfWuo73J41o3szHi3r34mTo/uCl/GMfg9Dg6PEeQ1z2gWXZc/zdJQOHCcwR5K/4u9wNZ1r60qn70Ad24ceOTYXfx+56p84O7Gyd/GQEVwHOEOaPYoqMN6FHfnaM807G7C+gfzxHmjGJy1q+/OuqmZhIBlcZzhDkjMDkEFC48R5gzislpv9d88+RHfvH4V//i+x+v17/2lz/9ih9988DE/2bFnqHfOfG4nz/58sPT3tSh89/xK6/60dq7vlGf9jrDl/zJT7/ipD8840D7/ZEr/+qk41/1s391xcjkjY2s+W8nHvcL77hm/B/b3Jt+84Tjf/mUPRNn6jgPAVXAc4Q5Y7qA3v667IdelWXZa25/b5b9YPPS6/fmx1z5w+0f27/h7mlu6t9+fPwn+39+cJrrbPu19gd+6rrW+w/++vhVfnXr+I09+wdZduIPND70hy+0PvTlE/IjfvJb7TN1noeAKuA5wpwxXUBP+KVbR0avfV2WvSR74y3Do9f+SJb9n9YhVzR69Qt/d8qfvSbLXr25603d/tIs+4n3nPbhtzaOfH/36zz6Y1n2w+845e9fn2WvvL/x/o7Gzb/srR/44O+9PMte+3D7xk5+6Rm76oeWNy6d0vzIzS9rHPz7//zHr85OemnrTJ3nIaASeI4wZ0wX0J9qfcV+QaNNP/Ns89KqLPuV5tvHfyh76TnNL7Of/PMs+7Uj3W7qHVn2t61PrGt0eKzrdd6eZbXHGm8PvjvL/rpeH/u9LDvp9uZVmq9Mf28sv7H2r5SekWVvarw58p+y7I0PNi488qvtX/fsOA8B1cBzhDljuoCubr3/aCNVa8cvvbL59mNZNpQfe6SWZZd3u6kfzbJd+Uf+c5Y91O06tzVef+5rfeDwz2UnjtW/1nj9uT0/Ykfjdeb/y2/s5PZH8t+t+lyWHf9E6wN7Xp0HtOM8BFQDzxHmjOkCmgdptJGqR8cvtX5F9KTsteP/A5TXZtl7ut3UG7Ls9sI5Oq7zgSz7RPsDp/74j2+rvz/L/nH84P+bZR/Kb2xDfeqZ35VlH2kf8uE8oB3nIaAaeI4wZ0wX0PYPw7Ps5WPjl5oZezzL/mz84D1Z9uvdbuqdWXbiubsnP955nd/MslunXrPx/vXjl7+SZb+d39jW+tQz/3yW3dz+wM15QO15Oh4NyonnCHPGtL/GlJv806T80jeKfz//hm439cTrm5960z9c+J38453XaXzx/fjUa/5klj08fvnBLPuZ/MbGf+UpP/MPTXzB3ixy1uU8HY8G5cRzhDmjx4BuKsbwhK439dxpJ7Z/Jemsw12v89LsJSNTr/nqLDs4fvlAlv1g8X7lZ25cZ/z3QUfa/2aIOU/Ho0E58RxhzugxoLdn2bsemLR1mpsa/Y8F7/iRZtp+6alu1/nhLNsz9Zo/MeUV6I4s+8l6l4C+dvIV6BMT/+hS4TwdjwblxHOEOaPHgD6SZX/ku6n62N0ffV2WvbvbdX4hy+5sXxw9dGi4/htZtmn8c9dn2W/UuwT01ya/B3pb4V+tmzhP511AKfEcYc7oMaBjP5WdOP69yYeXLr2l201dsPzf2x+5M8t+rtt1/jbLFrY/sDLLvlD/pyx73/itvC+/3BHQv5/8Kfyn8oDa83Q8GpQTzxHmjB4D2vyNo3ntj/yPLPt6t5t6U/Yz7Z/cP55lJ3a7zjVZ9rrnWu+PvS17yVP1G7Ps5TvyA7a/LMturHcJ6MYsOz7/Gv6Z1+QB7TgPAdXAc4Q5o9eAbj8ue8W5zZ/mjJyaZW8c63ZTf5NlqfWJQ+/Msj/tdp2RX82yt+5sfGB0oHGhkdG3NF5Ctr6ov/2kLPvdsXqXgI78UvsvkXb9VvsvkTrOQ0A18Bxhzug1oPUVjXr98rtPe2+jdCfc2/Wmrm0c8cb3pVPe9dose9mdXa9zz6uy7DV/8sn3N0J63B2N97ed0HgN+rYPfvBtjdefP7LD3K/2mTe/Iste+Qcf+pPXZCe/ovXJzvMQUAk8R5gzeg5ofcUr27+P9POFX4afcq3Bl0z8zud101znjte3P/DaL7fef+BXxq/ypo4/yhw/83WvzY/4jT3tT3aeh4Aq4DnCnNF7QOuPfuK/nPDKX/7Ti4frRZPX+o/3vvWkV/xY7S+unjik4zoHl7/9x17+cycPfq/9/vDlf/GG449/w19eNWJvbPLMewd+8zXHvXFg/8QnO85DQBXwHAGdSlCvEtwFvCieI6BTCepVgruAF8VzBHQqQb1KcBfwoniOgE4lqFcJ7gJeFM8R0KlRrwceeOD52Tr9/sbJX05ABfAcAZ2Oa/1G0cbZOn37H32ardPDjecI6ERA4cJzBAB9IqAA0CcCCgB9IqAA0CcCCgB9IqAA0CcCCgB9IqAA0CcCCgB9IqAA0CcCCgB9IqAA0CcCCgB9IqAA0CcCCgB9IqAA0CcCCgB9IqAA0CcCCgB9IqAA0CcCCgB9IqAA0Kf/Dym7UeNOnmHoAAAAAElFTkSuQmCC" title alt width="672" /></p> -<p>Here we get a performance benefit of a factor of 16.5 using the version of the differential equation model compiled from C code using the inline package!</p> +<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAABUAAAAPACAMAAADDuCPrAAAC+lBMVEUAAAABAQECAgIDAwMEBAQFBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQVFRUWFhYXFxcYGBgZGRkaGhobGxscHBwdHR0eHh4fHx8gICAhISEiIiIjIyMkJCQlJSUmJiYnJycoKCgpKSkqKiorKyssLCwtLS0uLi4wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6Ojo7Ozs8PDw9PT0+Pj4/Pz9AQEBBQUFCQkJDQ0NERERFRUVGRkZHR0dISEhJSUlKSkpLS0tMTExNTU1OTk5PT09QUFBRUVFSUlJTU1NVVVVWVlZXV1dYWFhZWVlaWlpbW1tcXFxdXV1eXl5fX19gYGBhYWFiYmJjY2NkZGRlZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1ubm5vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8fHx9fX1+fn5/f3+AgICBgYGCgoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uMjIyNjY2Ojo6Pj4+QkJCRkZGSkpKTk5OUlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2enp6fn5+goKChoaGioqKjo6OkpKSlpaWmpqanp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKzs7O0tLS1tbW2tra3t7e4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXGxsbHx8fIyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY2NjZ2dna2trb29vc3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr6+vs7Ozt7e3u7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+/v7////rs2anAAAACXBIWXMAAB2HAAAdhwGP5fFlAAAgAElEQVR4nO3de5zddX3n8Z/1gmhbsZbele1uW3t126O2FtS1tHahuC3aspZuVbpVrNZqW1tZ134nMwkQIIQYIOFOMFwDRQiiXIUlEC5FSgyXQIiGXMAECJfcyEzmPB57bnP5nXOmvH8nn2TmPb/X8w8zmTnf8zvj5P16nMwlZFUAQE+yyX4AAOCKgAJAjwgoAPSIgAJAjwgoAPSIgAJAjwgoAPSIgAJAjwgoAPSIgAJAjwgoAPSIgAJAjwgoAPSIgAJAjwgoAPSIgAJAjwgoAPSIgAJAjwgoAPSIgAJAjwgoAPSIgAJAjwgoAPSIgAJAjwgoAPSIgAJAjwgoAPSIgAJAjwgoAPSIgAJAjwgoAPSIgAJAjwgoAPSIgAJAjwgoAPSIgAJAjwgoAPSIgEbZugVw8PLLOyb7ITjQZk9Aozy7ATCwsVrdNdmPwYE2ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIapWBA1954+j997Mg/+qM//ovPnbB4+bq982cA6EBANdrsCWiUIgFdt/Qf3lsZ772fPH3Z+r32RwEYQ0A12uwJaJQCAb3tf4yE85BxHT308+fxTBR7HQHVaLMnoFEKBPTYWi4P/sLF/7Zpd+3cllU3nvE3Ix095Oh/OvmsC6+oWXrzfWv33h8OlBcB1WizJ6BRCgT0E5XKFdtzh4e+e/6n3lNp8+4Pn733/nigrAioRps9AY1SLKDPdN7BjrsXfPoD+YS+a9Xe+/OBkiKgGm32BDTKHge0YfMDN1554YULv/rV02b87fsqlRV7788HSoqAarTZE9AoMQEd728IKOIRUI02ewIahYDCAgHVaLPf9wFNaXOX1267+cI5M0+/5L6hLm/bnNJefkyFrtj9HSCg8EBANVoqpkhAHzghNc1d0/nGfRnQ5rUIKKYtAqrRgjE1Aro6pUWPv/jypnvnppnPdbyVgAJhCKhGC8bUCOj56crhxgvbT0vXdLx1Xwb0hSVLXumKBBTOCKhGC8aUCOjuWaOvuz+d2nGCz4ECYQioRkvFlAjo8yntaL34zFVXd5wgoEAYAqrRUrEPA7pm8YkDC+8fbvVn8J5zT5g5/1uNz3gO9qU722787NJ5A7POvGlr/eVGzr6Rzmm+ZUWa9XL+eLtnr53bf+rXVg93u6Oh2+b2n7J0+/ADF82Zefrtu6qNHq677KT+OVdvGr3WaEDzV8m/Ax1X1T80BBSTh4BqhIVW92VAlze/zH5Dsz8vLWz+dtaq+tsuT+myR3eNu/FDA803z15XbeXsyZSeb7zp4vpnSfPH81a1zt7Y7Y4uP/X6q/rS2Yv7LrlxcUr1Z7u1h9R3wtcWz079K6ttAc1fJf8OdCCgsEBANcJCq/swoM/1p3O+v2PdBanRn+EL0mkrtw2uvzTN2lJ744sLaq8euOjOp5vPGqubBtKCNTu3PjwnnbKjlbPh09Jd9Tdtm5HWtR/PX2hWunDD4Is3p/RElzs6c3u1em/tWt+v3fL2NGOwHtB06c7ak82laeCZfEDzV8m/A02Di8as3in7pBzQdfqdApran+zJfggGnte6ts8C+q/p9FqtqkMLGv15Is1uPMDhy9P19V8H/31Rf+PbQJc1nodenb5av3H1pePTspGc3dL8O/x96YzhjuPjLU0LdzcvuLTLHT1W+83WlL7deGXjoaQ0p3GT4XPTtfmA5q+SfweadqUxD+n/XxwrB/QH+p0CiLNNu9k+C+ictKLx64ON/nw7Xdd89Zp0ZusGux6/4cxah86uP/B5rRtXb06XjORsU/Pv8OfXn4h2OT5qflrZ+PWpG5Z3uaP6vyI3nNLaavPXRkCXNW/yaPpqPqD5q+TfgdZjJqDAtDTFAlpLzbONF55t9GfJWHhmj7vVM9f11Z83DvW1blyL2vzRz0gurKfz+dS/bcLjdbtnjJytTnRHo19HbwX08eZNnkszducCmrtK2zvQNLx9zOanZMfIAV2p3ykgebq2xsl+DAY2Cgut7ruAbkmp+UWiXY3+XJhOmt9yRu1J3YOj/7zwnWlg99iNqxvTzNHu3Vn/O/yydEW1/Xj3C+V/M/6O2gPa+r/q5fpT3PEBzV2l7R3oxBeRYIEvImmEhVb3XUBHnws+1+jPNc3PQracmB4eeXF7vWJjTxxXpXmj3Xuhr/amBWl1tf14zq6+NPbp3+531B7QR5o3eTr1DeYCmrtK2zvQiYDCAgHVCAut7sPPgc5tfQrxoUZ/7kiXNl/9wuMb65/X/PrIzTanmbvHfery1nTxWPcuTHdtSqfWv0KUP95+oUcbvz571XUT3VFbQG9s3mT5WGOb/5u/Sv4d6ERAYYGAaoSFVvdhQK9JZ9T/qbrdZzX6syH1ra+/dnhRWt74vqIHmt/ANHxNOrda/4r3/JEvnt8x1r370zm3Np8U5o/nXZnOadzXTfVkdr+jtoDObDxl3Tq7/unX8QHNXyX/DnQioLBAQDXCQqv7MKBb+tO5a3duuGig2Z9L0uzvvLR705I0e2vtr8e1KJ33wNotG1cuSqn+TeubBtLC7+3cVv/2ze1j3dvRn07ua37fZ+54dcmSrWMXero/XbRhcPvyvvrnBbrfUVtA08nfffHFlaek419q+0b63FXa3oEOBBQWCKhG69q++0mku5tfz76/2Z9tZ7V+yKf+E0LVbV8b+XL3wH2NG6/sb/72pPpTwNGfTL8spYuaL+WPpzT+2+kfnNF82zcnvKO2gN7S17hJ/0PVtoDmr5J/BzoQUFggoBphodV9/bPw/QseGqnX4PJzjp+14Futp47Dq685/9SBeRfdNPLD7c8sPW3g5EW3jP0sfN3DqfU9nm3H8wGtPn3NWTPnLn5i4jtqC+jmDZfNntH9Z+FzDzL/DrQjoLBAQDXCQqvT5L+JNJTEn7vqZoIeFkZAYYGAarTZT4uAbkm7XvlGEyGgKBUCqtFmPy0C+u3z9uDwVAvo0HPr1z/y8MMrln39SAKKeARUo81+WgR0Sdd/FlQ0dQK6e9WVs/7qD99ZGWfl3vvzgZIioBpt9tMioHtkcgK6tv342iv+8QOVdp9Yv/f+fKCkCKhGmz0BnYSAHlOpvOeL1z09cnLXI1d95fCRZh585Mc++bmaL8+Yd/m9e+9PB0qLgGq02RPQKAUCenbzr+mHfuIfv9L3z3/9J+9utfN3j5lz7Xf33h8IoI6AarTZE9AoBQK6Ydlx7237y/q7P3riDWv32p8FYBQB1WizJ6BRigR0w4bvf73vo4e02nnEp0+5dvVe+lMAtCGgGm32BDRKsYA2PHL37bcvX7ku/mMPTIiAarTZE9AoPQQU2PcIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBYIqEabPQGNQkBhgYBqtNkT0CgEFBbKGdDHHyl6Qps9AY1CQGGhlAF94P3v/kbBI9rsCWgUAgoLpQzoeZXKlwoe0WZPQKMQUFgoZUBPrlQ+WfCINnsCGoWAwkIpAzqjUvnzgke02RPQKAQUFkoZ0H+uVI4seESbPQGNQkBhoZQB/btK5bCCR7TZE9AoBBQWShnQYyuV3yt4RJs9AY1CQGGhlAE9plI5uOARbfYENAoBhYVSBvToSuVdBY9osyegUQgoLJQyoEdVKpW1xY5osyegUQgoLJQyoEfWArq62BFt9gQ0CgGFhVIG9IhaQB8udkSbPQGNQkBhoZQB/cNaQFcUO6LNnoBGIaCwUMqAHloL6P3FjmizJ6BRCCgslDKg76sF9O5iR7TZE9AoBBQWShnQg2sBXVbsiDZ7AhqFgMJCKQP6u7WA3lbsiDZ7AhqFgMJCaQN6U7Ej2uwJaBQCCgulDej1xY5osyegUQgoLJQ2oNcWO6LNnoBGIaCwUNqAXl3siDZ7AhqFgMJCaQO6pNgRbfYENAoBhYXSBvTSYke02RPQKAQUFkob0MXFjmizJ6BRCCgslDagFxQ7os2egEYhoLBQ2oCeW+yINnsCGoWAwkJpA7qw2BFt9gQ0CgGFhdIG9PRiR7TZE9AoBBQWShvQ04od0WZPQKMQUFgobUDnFDuizZ6ARiGgsFDagM4udkSbPQGNQkBhobQBnVXsiDZ7AhqFgMJCaQPaX+yINnsCGoWAwkJpA5qKHdFmT0CjEFBYKG1Av1zsiDZ7AhqFgMJCaQP6pWJHtNkT0CgEFBZKG9AvFjuizZ6ARiGgsFDagH6h2BFt9gQ0CgGFhdIG9HPFjmizJ6BRCCgslDagnyl2RJs9AY1CQGGhtAH9VLEj2uwJaBQCCgulDej/LnZEmz0BjUJAYaG0Af14sSPa7AloFAIKC6UN6P8qdkSbPQGNQkBhobQB/WixI9rsCWgUAgoLpQ3oUcWOaLMnoFEIKCyUNqAfKXZEmz0BjUJAYaG0Af2TYke02RPQKAQUFkob0COKHdFmT0CjEFBYKG1ADyt2RJs9AY1CQGGhtAH9YLEj2uwJaBQCCgulDeihxY5osyegUQgoLJQ2oO8rdkSbPQGNQkBhobQBPbjYEW32BDQKAYWF0gb0t4sd0WZPQKMQUFgobUDfWeyINnsCGoWAwkJpA1pZV+iINnsCGoWAwkJ5A/q9Qke02RPQKAQUFsob0NWFjmizJ6BRCCgslDegqwod0WZPQKMQUFgob0BXFjqizZ6ARiGgsFDegD5Y6Ig2ewIahYDCQnkDen+hI9rsCWgUAgoL5Q3oPYWOaLMnoFEIKCyUN6B3FjqizZ6ARiGgsFDegN5e6Ig2ewIahYDCQnkDekuhI9rsCWgUAgoL5Q3oDYWOaLMnoFEIKCyUN6DfKHREmz0BjUJAYaG8AT39m0WOaLMnoFEIKCyUN6CVynkFjmizJ6BRCCgslDmg/1jgiDZ7AhqFgMICAdVosyegUQgoLBBQjTZ7AhqFgMICAdVosyegUQgoLBBQjTZ7AhqFgMICAdVosyegUQgoLBBQjTZ7AhqFgMICAdVosyegUQgoLBBQjTZ7AhqFgMICAdVosyegUQgoLBBQjTZ7AhqFgMICAdVosyegUQgoLJQ0oO8hoFMbAYWFkgb0UAI6tRFQWChpQI/4bQI6pRFQWChpQD/03wnolEZAYaGsAf0kAZ3SCCgslDWgK44moFMZAYWFsga0+ggBncoIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02e/daNtgAABuvSURBVAIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLMnoFEIKCwQUI02ewIahYDCAgHVaLOPCmhKm7u8dtvNF86Zefol9w11edvmlIKuHe8/fGzd31UCCg8EVKOlYq8G9IETUtPcNZ1vnJoBbT4qAoppi4BqtGDszYCuTmnR4y++vOneuWnmcx1vDQzoBEnrBQHFNEdANVow9mZAz09XDjde2H5auqbjrVMzoC8sWVIloJjGCKhGC8ZeDOjuWaOvuz+d2nFiaga0iYBi2iKgGi0VezGgz6e0o/XiM1dd3XGCgAKTgIBqtFTscUDXLD5xYOH9w62qDN5z7gkz53+r8RnPwb50Z3tkls4bmHXmTVvrLzci9Y10TvMtK9Ksl/PH8+q3XnXe8Sec/Z3d7fd0TuPrVI+Mu8q1c/tP/drq4W6XHLptbv8pS7cPP3DRnJmn376r2ujhustO6p9z9abRRzUa0Pzjyb+r7QgoLBBQjda/PQ3o8uaX2W9oVuWlhc3fzlpVf9vlKV326K5xN35ooPnm2euqrUg9mdLzjTddXP8saf54Tu3W98645MZFKd3Yfk/fqL3ynEVrR2+6qvW2jhvW7+TyU6+/qi+dvbjvkhsXp1R/Xlx78H0nfG3x7NS/stoW0Pzjyb+rHQgoLBBQjRbAPQzoc/3pnO/vWHdBalRl+IJ02sptg+svTbO21N744oLaqwcuuvPp5nPB6qaBtGDNzq0Pz0mn7GhFavi0dFf9TdtmpHXtx3Nqtz6+/r1Qt6b+wfZ7yv+l+rlZ6cINgy/enNITXS555vZq9d7ao/p+7Za3pxmD9cPp0p21J5tL08Az+YDmH0/+XW3alcY8tGf/TwLYaw5uBvQr+olt2s32MKD/mk6vNag6tKBRlSfS7MbzyeHL0/X1Xwf/fVF/49tAlzWeh16dvlq/cfWl49OykUjd0vw7/H3pjOGO4+PVbn1T42zjQvl7ygd0aVq4u/nQlna55GO132xN6dtj95XSnMZNhs9N1+YDmn88+Xe1iYACDqZqQOekFY1fH2xU5dvpuuar16QzWzfY9fgNZ9bqcnb94cxr3bh6c7pkJFKbmn+HP7/+RLTL8VG1Wz9V/7X5Gcj8PeUDOj+tbPz61A3Lu1xye/M+1o7dV2o0uObR9NV8QPOPJ/+utt6708asGgIc1P7kT/ZD2OdaAf2/+omXtALuWUBrz8CebbzwbKMqS8aej80ed6tnruurPxsc6mvduJaq+aOfZ1xYT+fzqX/bhMcbarfe2XihfqG2e8oFdPeMkbdVJ7rk6O1bAX28eZPn0ozduYDmHk/bu9qJz4HCAp8D1Qj9q+5pQLek1Pwi0a5GVS5MJ81vOaP2VO3B7SO3uzMN7B67cXVjmjlaszvrf4dflq6oth/PG/3GovqF2u4pF9Cxt+V+M/6S7QHd2LzJy/Unw+MDmns8be9qJwIKCwRUI/SvuqcBHX2G91yjKtc0P7fYcmJ6eOTF7fU2jT0dXJXmjdbshb7amxak1dX243m5gLbdUy6gu/paX9fPPb7cJdsD2voOqKdT32AuoLnH0/audiKgsEBANRO2KGcPPwc6t/WJwYcaVbkjXdp89QuPb6x/XvPrIzfbnGbuHvcJyVvTxWM1uzDdtSmdWv+6T/54Xi6gbfeU/xzo3PRo49dnr7puoku2BfTG5k2WjzW2+b/5x5N/VzsRUFggoJoJo5ezhwG9Jp1R/6fqdp/VqMqG1Le+/trhRWl547uFHmh+A9PwNencav3r2PNHviR+x1jN7k/n3Np8qpc/npcPaP6e8gG9Mp3TuOpN9WR2v2RbQGc2nrJunV3/RO34gOYfT/5d7URAYYGAarQC7mFAt/Snc9fu3HDRQLMql6TZ33lp96YlafbW2l96a6k574G1WzauXJRS/VvRNw2khd/bua3+TZnbx2q2oz+d3Nf8vs/c8eqSJVvHLpQPaP6eqrmfQ3q6P120YXD78r76ZxC6X7ItoOnk77744spT0vEvtX0jfe7xtL2rHQgoLBBQjVbAPf1JpLubX6W+v1mVbWe1fnSn/nM/1W1fG/ki9sB9jRuv7G/+9qT11XFJvCyli5ov5Y+nNO7b6fMBzd9T9dRaA8d+eOnBGc23fXPCS7YF9Ja+xk36H6q2BTT/ePLvagcCCgsEVKMFMOJn4fsXPDTSpMHl5xw/a8G3Wk8dh1dfc/6pA/Muumnkh9ufWXrawMmLbhn7Wfi6h1PrOzfbjv9HAc3dU3XNmf0nPzF226evOWvm3MVPTHzJtoBu3nDZ7BndfxY+9+7k39V2BBQWCKhGiF91Sv83kYbGfTl9L4r6p5wIKCwQUI02+ykc0C1p1yvfaM8RUJQKAdVos5/CAf32efvkMgQUpUJANdrsp3BAl3T5Z0H3AgKKUiGgGm32Uzig+wgBRakQUI02ewJKQFEqBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLBFSjzZ6ARiGgsEBANdrsCWgUAgoLZQ3oHR8koFMZAYWFsgb0ExUCOpURUFgoa0D/gIBOaQQUFkoa0CPeRUCnNAIKCyUN6O9VCOiURkBhoaQB/R0COrURUFgoaUArBHRqI6CwQEA12uwJaBQCCgsEVKPNnoBGIaCwQEA12uwJaBQCCgsEVKPNnoBGIaCwQEA12uwJaBQCCgsEVKPNnoBGIaCwQEA12uwJaBQCCgsEVKPNnoBGIaCwQEA12uwJaBQCCgsEVKPNnoBGIaCwUOaAfrHAEW32BDQKAYWFEgf0nZcXOKLNnoBGIaCwUN6ALr6/yBFt9gQ0CgGFhfIG9PpCR7TZE9AoBBQWyhvQGwod0WZPQKMQUFgob0BvKXREmz0BjUJAYaG8Ab290BFt9gQ0CgGFhfIG9K5CR7TZE9AoBBQWyhvQewsd0WZPQKMQUFgob0C/U+iINnsCGoWAwkJ5A7qi0BFt9gQ0CgGFhfIGdGWhI9rsCWgUAgoL5Q3oqkJHtNkT0CgEFBbKG9DVhY5osyegUQgoLJQ3oN8vdESbPQGNQkBhobwBXVfoiDZ7AhqFgMJCaQP6rmJHtNkT0CgEFBZKG9DfKXZEmz0BjUJAYaG0AT2k2BFt9gQ0CgGFhdIG9P3FjmizJ6BRCCgslDaghxY7os2egEYhoLBQ2oB+sNgRbfYENAoBhYXSBvTwYke02RPQKAQUFkob0A8VO6LNnoBGIaCwUNqA/kmxI9rsCWgUAgoLpQ3oR4od0WZPQKMQUFgobUCPKnZEmz0BjUJAYaG0Af3zYke02RPQKAQUFkob0L8sdkSbPQGNQkBhobQB/USxI9rsCWgUAgoLpQ3oXxc7os2egEYhoLBQ2oAeW+yINnsCGoWAwkJpA/rZYke02RPQKAQUFkob0L8rdkSbPQGNQkBhobQB/ftiR7TZE9AoBBQWShvQLxY7os2egEYhoLBQ2oAeV+yINnsCGoWAwkJpA/rlYke02RPQKAQUFkob0L5iR7TZE9AoBBQWShvQgWJHtNkT0CgEFBZKG9Djix3RZk9AoxBQWChtQE8qdkSbPQGNQkBhobQBPbXYEW32BDQKAYWF0gZ0XrEj2uwJaBQCCgulDegZxY5osyegUQgoLJQ2oGcVO6LNnoBGIaCwUNqAnlfsiDZ7AhqFgMJCaQN6YbEj2uwJaBQCCgulDejiYke02RPQKAQUFkob0MuKHdFmT0CjEFBYKG1AlxQ7os2egEYhoLBQ2oBeXeyINnsCGoWAwkJpA3ptsSPa7AloFAIKC6UN6DeLHdFmT0CjEFBYKG1Aby52RJs9AY1CQGGhtAG9rdgRbfYENAoBhYVSBvSQWkCXFTuizZ6ARiGgsFDKgL6/FtB7ih3RZk9AoxBQWChlQH+/FtD7ix3RZk9AoxBQWChlQP+wFtAVxY5osyegUQgoLJQyoB+qBfSRYke02RPQKAQUFkoZ0A/XArq62BFt9gQ0CgGFhVIG9KhaQJ8sdkSbPQGNQkBhoZQB/YtK5d0Fj2izJ6BRCCgslDKgx1QqhxQ8os2egEYhoLBQyoB+ulI5tOARbfYENAoBhYVSBvTzlcrhBY9osyegUQgoLJQyoF+qVD5c8Ig2ewIahYDCQikD2l+pHF3wiDZ7AhqFgMJCKQM6p1L5VMEj2uwJaBQCCgulDOiiSuXLBY9osyegUQgoLJQyoA8f8f6C/xwoAd3HCCgslDKgPdBmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhQUCqtFmT0CjEFBYIKAabfYENAoBhYX1jz762GQ/Bgfa7AkoUCpDKS2Y7McwfRBQoFQIaCQCCpQKAY1EQIFSIaCRCChQKgQ0EgEFSoWARiKgQKkQ0EgEFCgVAhqJgAKlQkAjEVCgVAhoJAIKlAoBjURAgXLZvn3nZD+E6YOAAkCPCCgA9IiAAkCPCCgA9IiAAkCPCCgA9IiAAkCPCCgA9IiAAmVxbRp5aWj5Wccff9ZdQ5P4YKYHAgqUxMsnptZLu85PDefumszHMx0QUKActl+WUuvFpen4e5977t7j09LJfEDTAQEFSuDlmy+elUYCumVGeqz+62NpxvOT+JimAwIKlMBLzb+zN39zV5o3XP91eF66ZxIf03RAQIESGPrBD37w6EhAr0g3NV+4KV05aY9oeiCgQElsHgnoeWlF84UV6YLJejDTBAEFSmI0oPPS6uYLq9P8SXs00wMBBUpiNKCz0sbmCxvTCZP2aKYHAgqURGdAN6SZk/ZopgcCCpTEuL/CP9F8YXWaN2mPZnogoEBJjAb0/PTd5gsr0vmT9mimBwIKlMTmsW9jurX5wq1pyaQ9mumBgAIlsXnsG+nPaL5wRlo+aY9meiCgQElsHvejnE/Wf13Hj3LuKQIKlMRoQKvXpTmPD+56bE76xmQ+numAgAIlMRbQXeelNNCf0gWDk/l4pgMCCpTEWECrQ3ctnDXrrLt3T+KjmR4IKAD0iIACQI8IKAD0iIACQI8IKAD0iIACQI8IKAD0iIACQI8IKAD0iIDC3lFZu6P39C73a9zN0ohH98pX6nzlTc13Yx9cHnuGjxHsdQ9olu3X+10SUEj4GMHemR9r+qEsa720sLrnAV26dOnTYQ/xP7xS5ys31S7+GgJqgI8Rpo18i/Y0oHv8cPbwSvvuIaB3fIwwbeSTs2TJVVF3tTcRUGt8jDBtBCaHgELCxwjTRj45rd/Vf3n6739p/zf+0mfXV6u3/enPvu7H39U3+h+y2DzwngP3+4XDLt454V3tOOvwX33Dj1eOvqc64ZnBC4/42dcd9Psnbmv9fuiyPzto/zf8pz+7ZGjszobOfd+B+/3i4VeP/Aucz6TfOmD/Xzlu8+iVOq5DQB3wMcK0MVFA735L9iNvyLLsTXcfm2U/XH/prc80b3PZj7a+bP+2Bya4q6//5MhX9o/cPsGZx3+99YqfuaHx+1W/MXLk1x4dubMXPpBlB/5Q7VW//3LjVdcf0LzFT/9760qd1yGgDvgYYdqYKKAH/PKdQ7uve0uWvSp7+7LB3df9WJb9deMml9R69YsfO+6P35Rlb1ze9a7ufnWW/dQx//KFg2u3/Gz3M0/+RJb96OHH/dVbs+z1D9d+v6Z29685+G8/d8hrs+zNT7Tu7LBXn7ixumN+7aXj6q+54zW1G/+3v/vgG7ODXt24Uud1CKgFPkaYNiYK6M80/sZ+dq1NP/dC/aWFWfar9V/X/0j26tPqf81++sgs+/Vd3e7q8Cz7y8YbLqp1eLjrmUOzrLKu9uv2j2fZ/6xWhw/JsoPurh+pPzM9ZLh5Z61vKT0xy95R+2XXf86yt6+qvbD211rf7tlxHQLqgY8Rpo2JAnpO4/dP1lK1aOSl19d//WKWDTRvu6uSZRd3u6sfz7KNzdf8lyx7rNuZu2rPP7c0XrHz57MDh6u31Z5/rm7eYk3teeb/a97ZYa3XNL+36oos2/+pxis2v7EZ0I7rEFAPfIwwbUwU0GaQdtdS9eTIS41vET0oe/PIf5Xyuiw7pttdvS3L7s5do+PM32bZl1qv+MpP/uTj1c9m2adHbvypLPt8886uqY6/8tFZ9vetm3yhGdCO6xBQD3yMMG1MFNDWF8Oz7LXDIy/VM7Y+y/545Mabs+w3ut3VUVl24Ombxl7feea3suzO8Sdrv79x5OVvZdm7m3f2aHX8lX8hy+5oveKOZkDbr9Px3mBq4mOEaWPCb2NqGvvRpOZL9+R/fv5t3e7qqbfW3/SOT573vebrO8/U/vK9fvzJn86yJ0ZeXpVlP9e8s5FveWpe+UdG/8JeL3LW5Tod7w2mJj5GmDYKBvSmfAwP6HpXL/7Lga1vSTp5Z9czr85eNTT+5BuzbPvIy9uy7Ifzj6t55dqZke8HHWr9myFt1+l4bzA18THCtFEwoHdn2dGPjHl0grva/W+zDv+xetp++dluZ340yzaPP/lT456Brsmyn652Ceibx56BPjX6jy7lrtPx3mBq4mOEaaNgQNdm2R9od1UdfuAf3pJlH+925hez7L7Wi7t37Bis/maW3TTythuz7DerXQL662OfA70r96/WjV6n8yFgSuJjhGmjYECHfyY7cORzk0/Mnbus212dPf+brdfcl2U/3+3MX2bZCa1XLMiyf63+TZZ9ZuRePtN8uSOgfzX2VfgvNwPafp2O9wZTEx8jTBsFA1r/jqMZrdd8OMtu73ZX78h+rvWV+/VZdmC3M1dn2VtebPx++L3Zq56t3pJlr13TvMHq12TZLdUuAV2aZfs3/w7//JuaAe24DgH1wMcI00bRgK7eL3vd6fWv5gx9JcvePtztrj6aZanxhh1HZdmHup0Z+rUsO3hD7RW7+2ov1DL627WnkI2/1N99UJb97nC1S0CHfrn1k0gb39n6SaSO6xBQD3yMMG0UDWj1zFq9fuXj/3JsrXQHrOh6V9fVbvH2z6Tjjn5zlr3mvq5nHnxDlr3piP/z2VpI97u39vvHD6g9B33v5z733trzzx9b0/a4Wlde/rose/0HPn/Em7LDXtd4Y+d1CKgFPkaYNgoHtHrm61vfj/QLuW+GH3eq/1Wj3/N5wwRn7n1r6xVvvr7x+0d+deTIOzp+KHPkyje8uXmL39zcemPndQioAz5GmDaKB7T65Jf+6wGv/5UPXTBYzRs79W/HHnzQ636i8pGrRm/ScWb7/EN/4rU/f1j/D1q/H7z4I2/bf/+3/enlQ+13NnblZ/p+6037vb1v6+gbO65DQB3wMQI6TYF6TYGHgFfExwjoNAXqNQUeAl4RHyOg0xSo1xR4CHhFfIyATlOgXlPgIeAV8TECOtXq9cgjj7w0WZffWrv4awmoAT5GQKf9Gt9RtHSyLt/6R58m6/KQ8TECOhFQSPgYAUCPCCgA9IiAAkCPCCgA9IiAAkCPCCgA9IiAAkCPCCgA9IiAAkCPCCgA9IiAAkCPCCgA9IiAAkCPCCgA9IiAAkCPCCgA9IiAAkCPCCgA9IiAAkCPCCgA9IiAAkCPCCgA9Oj/A1FX2XCkTa+RAAAAAElFTkSuQmCC" 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>This vignette was built with mkin 0.9.41.9000 on</p> <pre><code>## R version 3.2.2 (2015-08-14) ## Platform: x86_64-pc-linux-gnu (64-bit) |