Appendix-E-Ensemble-Kalman-Filters.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "[Table of Contents](http://nbviewer.ipython.org/github/rlabbe/Kalman-and-Bayesian-Filters-in-Python/blob/master/table_of_contents.ipynb)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Ensemble Kalman Filters"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>\n",
       "@import url('http://fonts.googleapis.com/css?family=Source+Code+Pro');\n",
       "@import url('http://fonts.googleapis.com/css?family=Lora');\n",
       "\n",
       "//@import url('http://fonts.googleapis.com/css?family=Open+Sans');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Vollkorn');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Karla');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Poppins');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Arimo');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Roboto');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Lato');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Domine');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Chivo');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Cardo');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Arvo');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Crimson+Text');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Ubuntu');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Fontin');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Raleway');\n",
       "//@import url('http://fonts.googleapis.com/css?family=Merriweather');\n",
       "\n",
       "\n",
       ".CodeMirror pre {\n",
       "    font-family: 'Source Code Pro', Consolas, monocco, monospace;\n",
       "}\n",
       "    div.cell{\n",
       "        width: 850px;\n",
       "        margin-left: 0% !important;\n",
       "        margin-right: auto;\n",
       "    }\n",
       "    div.text_cell_render{\n",
       "        font-family: 'Lora';\n",
       "        //font-family: 'Open Sans';\n",
       "        //font-family: 'Karla',verdana,arial,sans-serif;\n",
       "        //font-family: 'Roboto',verdana,arial,sans-serif;\n",
       "        //font-family: 'Lato',verdana,arial,sans-serif;\n",
       "        //font-family: 'Domine',verdana,arial,sans-serif;\n",
       "        //font-family: 'Chivo',verdana,arial,sans-serif;\n",
       "        //font-family: 'Cardo',verdana,arial,sans-serif;\n",
       "        //font-family: 'Arvo',verdana,arial,sans-serif;\n",
       "        //font-family: 'Poppins',verdana,arial,sans-serif;   \n",
       "        //font-family: 'Ubuntu',verdana,arial,sans-serif;\n",
       "        //font-family: 'Fontin',verdana,arial,sans-serif;\n",
       "        //font-family: 'Raleway',verdana,arial,sans-serif;\n",
       "        //font-family: 'Merriweather',verdana,arial,sans-serif;\n",
       "        //font-family: 'Crimson Text', verdana,arial,sans-serif;\n",
       "        //font-family: verdana,arial,sans-serif;\n",
       "        //font-family: arial,sans-serif;\n",
       "        line-height: 125%;\n",
       "        font-size: 130%;\n",
       "        text-align: justify;\n",
       "        text-justify:inter-word;\n",
       "    }\n",
       "    div.text_cell code {\n",
       "        background: transparent;\n",
       "        color: #000000;\n",
       "        font-weight: 400;\n",
       "        font-size: 12pt;\n",
       "        //font-style: bold;\n",
       "        font-family:  'Source Code Pro', Consolas, monocco, monospace;\n",
       "   }\n",
       "    h1 {\n",
       "        font-family: 'Open sans',verdana,arial,sans-serif;\n",
       "\t}\n",
       "\t\n",
       "    div.input_area {\n",
       "        background: #F6F6F9;\n",
       "        border: 1px solid #586e75;      \n",
       "    }\n",
       "\n",
       "    .text_cell_render h1 {\n",
       "        font-weight: 200;\n",
       "        font-size: 30pt;\n",
       "        line-height: 100%;\n",
       "        color:#c76c0c;\n",
       "        margin-bottom: 0.5em;\n",
       "        margin-top: 1em;\n",
       "        display: block;\n",
       "        white-space: wrap;\n",
       "        text-align: left;\n",
       "    } \n",
       "    h2 {\n",
       "        font-family: 'Open sans',verdana,arial,sans-serif;\n",
       "        text-align: left;\n",
       "    }\n",
       "    .text_cell_render h2 {\n",
       "        font-weight: 200;\n",
       "        font-size: 16pt;\n",
       "        font-style: italic;\n",
       "        line-height: 100%;\n",
       "        color:#c76c0c;\n",
       "        margin-bottom: 0.5em;\n",
       "        margin-top: 1.5em;\n",
       "        display: block;\n",
       "        white-space: wrap;\n",
       "        text-align: left;\n",
       "    } \n",
       "    h3 {\n",
       "        font-family: 'Open sans',verdana,arial,sans-serif;\n",
       "    }\n",
       "    .text_cell_render h3 {\n",
       "        font-weight: 200;\n",
       "        font-size: 14pt;\n",
       "        line-height: 100%;\n",
       "        color:#d77c0c;\n",
       "        margin-bottom: 0.5em;\n",
       "        margin-top: 2em;\n",
       "        display: block;\n",
       "        white-space: wrap;\n",
       "        text-align: left;\n",
       "    }\n",
       "    h4 {\n",
       "        font-family: 'Open sans',verdana,arial,sans-serif;\n",
       "    }\n",
       "    .text_cell_render h4 {\n",
       "        font-weight: 100;\n",
       "        font-size: 14pt;\n",
       "        color:#d77c0c;\n",
       "        margin-bottom: 0.5em;\n",
       "        margin-top: 0.5em;\n",
       "        display: block;\n",
       "        white-space: nowrap;\n",
       "    }\n",
       "    h5 {\n",
       "        font-family: 'Open sans',verdana,arial,sans-serif;\n",
       "    }\n",
       "\n",
       "    .text_cell_render h5 {\n",
       "        font-weight: 200;\n",
       "        font-style: normal;\n",
       "        color: #1d3b84;\n",
       "        font-size: 16pt;\n",
       "        margin-bottom: 0em;\n",
       "        margin-top: 0.5em;\n",
       "        display: block;\n",
       "        white-space: nowrap;\n",
       "    }\n",
       "    div.output_subarea.output_text.output_pyout {\n",
       "        overflow-x: auto;\n",
       "        overflow-y: scroll;\n",
       "        max-height: 50000px;\n",
       "    }\n",
       "    div.output_subarea.output_stream.output_stdout.output_text {\n",
       "        overflow-x: auto;\n",
       "        overflow-y: scroll;\n",
       "        max-height: 50000px;\n",
       "    }\n",
       "    div.output_wrapper{\n",
       "        margin-top:0.2em;\n",
       "        margin-bottom:0.2em;\n",
       "}\n",
       "\n",
       "    code{\n",
       "        font-size: 6pt;\n",
       "\n",
       "    }\n",
       "    .rendered_html code{\n",
       "    background-color: transparent;\n",
       "    }\n",
       "    ul{\n",
       "        margin: 2em;\n",
       "    }\n",
       "    ul li{\n",
       "        padding-left: 0.5em; \n",
       "        margin-bottom: 0.5em; \n",
       "        margin-top: 0.5em; \n",
       "    }\n",
       "    ul li li{\n",
       "        padding-left: 0.2em; \n",
       "        margin-bottom: 0.2em; \n",
       "        margin-top: 0.2em; \n",
       "    }\n",
       "    ol{\n",
       "        margin: 2em;\n",
       "    }\n",
       "    ol li{\n",
       "        padding-left: 0.5em; \n",
       "        margin-bottom: 0.5em; \n",
       "        margin-top: 0.5em; \n",
       "    }\n",
       "    ul li{\n",
       "        padding-left: 0.5em; \n",
       "        margin-bottom: 0.5em; \n",
       "        margin-top: 0.2em; \n",
       "    }\n",
       "    a:link{\n",
       "       font-weight: bold;\n",
       "       color:#447adb;\n",
       "    }\n",
       "    a:visited{\n",
       "       font-weight: bold;\n",
       "       color: #1d3b84;\n",
       "    }\n",
       "    a:hover{\n",
       "       font-weight: bold;\n",
       "       color: #1d3b84;\n",
       "    }\n",
       "    a:focus{\n",
       "       font-weight: bold;\n",
       "       color:#447adb;\n",
       "    }\n",
       "    a:active{\n",
       "       font-weight: bold;\n",
       "       color:#447adb;\n",
       "    }\n",
       "    .rendered_html :link {\n",
       "       text-decoration: underline; \n",
       "    }\n",
       "    .rendered_html :hover {\n",
       "       text-decoration: none; \n",
       "    }\n",
       "    .rendered_html :visited {\n",
       "      text-decoration: none;\n",
       "    }\n",
       "    .rendered_html :focus {\n",
       "      text-decoration: none;\n",
       "    }\n",
       "    .rendered_html :active {\n",
       "      text-decoration: none;\n",
       "    }\n",
       "    .warning{\n",
       "        color: rgb( 240, 20, 20 )\n",
       "    } \n",
       "    hr {\n",
       "      color: #f3f3f3;\n",
       "      background-color: #f3f3f3;\n",
       "      height: 1px;\n",
       "    }\n",
       "    blockquote{\n",
       "      display:block;\n",
       "      background: #fcfcfc;\n",
       "      border-left: 5px solid #c76c0c;\n",
       "      font-family: 'Open sans',verdana,arial,sans-serif;\n",
       "      width:680px;\n",
       "      padding: 10px 10px 10px 10px;\n",
       "      text-align:justify;\n",
       "      text-justify:inter-word;\n",
       "      }\n",
       "      blockquote p {\n",
       "        margin-bottom: 0;\n",
       "        line-height: 125%;\n",
       "        font-size: 100%;\n",
       "      }\n",
       "</style>\n",
       "<script>\n",
       "    MathJax.Hub.Config({\n",
       "                        TeX: {\n",
       "                           extensions: [\"AMSmath.js\"],\n",
       "                           equationNumbers: { autoNumber: \"AMS\", useLabelIds: true}\n",
       "                           },\n",
       "                tex2jax: {\n",
       "                    inlineMath: [ ['$','$'], [\"\\\\(\",\"\\\\)\"] ],\n",
       "                    displayMath: [ ['$$','$$'], [\"\\\\[\",\"\\\\]\"] ]\n",
       "                },\n",
       "                displayAlign: 'center', // Change this to 'center' to center equations.\n",
       "                \"HTML-CSS\": {\n",
       "                    scale:95,\n",
       "                        availableFonts: [],\n",
       "                        preferredFont:null,\n",
       "                        webFont: \"TeX\",\n",
       "                    styles: {'.MathJax_Display': {\"margin\": 4}}\n",
       "                }\n",
       "        });\n",
       "</script>\n"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#format the book\n",
    "%matplotlib inline\n",
    "from __future__ import division, print_function\n",
    "from book_format import load_style\n",
    "load_style()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "> I am not well versed with Ensemble filters. I have implemented one for this book, and made it work, but I have not used one in real life. Different sources use slightly different forms of these equations. If I implement the equations given in the sources the filter does not work. It is possible that I am doing something wrong. However, in various places on the web I have seen comments by people stating that they do the kinds of things I have done in my filter to make it work. In short, I do not understand this topic well, but choose to present my lack of knowledge rather than to gloss over the subject. I hope to master this topic in the future and to author a more definitive chapter. At the end of the chapter I document my current confusion and questions. In any case if I got confused by the sources perhaps you also will, so documenting my confusion can help you avoid the same.\n",
    "\n",
    "\n",
    "The ensemble Kalman filter (EnKF) is very similar to the unscented Kalman filter (UKF) of the last chapter. If you recall, the UKF uses a set of deterministically chosen weighted sigma points passed through nonlinear state and measurement functions. After the sigma points are passed through the function, we find the mean and covariance of the points and use this as the filter's new mean and covariance. It is only an approximation of the true value, and thus suboptimal, but in practice the filter is highly accurate. It has the advantage of often producing more accurate estimates than the EKF does, and also does not require you to analytically derive the linearization of the state and measurement equations. \n",
    "\n",
    "The ensemble Kalman filter works in a similar way, except it uses a *Monte Carlo* method to choose a large numbers of sigma points. It came about from the geophysical sciences as an answer for the very large states and systems needed to model things such as the ocean and atmosphere. There is an interesting article on it's development in weather modeling in *SIAM News* [1]. The filter starts by randomly generating a large number of points distributed about the filter's initial state. This distribution is proportional to the filter's covariance $\\mathbf{P}$. In other words 68% of the points will be within one standard deviation of the mean, 95% percent within two standard deviations, and so on. Let's look at this in two dimensions. We will use `numpy.random.multivariate_normal()` function to randomly create points from a multivariate normal distribution drawn from the mean (5, 3) with the covariance\n",
    "\n",
    "$$\\begin{bmatrix}\n",
    "32 & 15 \\\\ 15 & 40\n",
    "\\end{bmatrix}$$\n",
    "\n",
    "I've drawn the covariance ellipse representing two standard deviations to illustrate how the points are distributed."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAuoAAADaCAYAAADwgmciAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xt8nNV97/vPWs8zo9HdF1l32ZZsycIXfJGFLAROuBhC\nSAlJm9KEDZScsF+7J+yG5rTZyd49JN07KU12smnTJmn3aZuSJmxIUiChAQIBDEYY2ZYtY0uWdbEl\nW5YlXyTrPtLMs9b5Y2lkGxtzieWL+L1fL7+wRzOjNY+H8ff56ff8lrLWWoQQQgghhBCXFH2xFyCE\nEEIIIYQ4kwR1IYQQQgghLkES1IUQQgghhLgESVAXQgghhBDiEiRBXQghhBBCiEuQBHUhhBBCCCEu\nQRLUhRBCCCGEuASd96D+/e9/n5UrV5KZmUlmZiZXX301zzzzzGn3+drXvkZBQQEpKSlcd911NDU1\nne9lCCGEEEIIcVk770G9qKiIb33rW+zYsYP6+nquv/56br/9dnbv3g3AN7/5TR5++GG+973vsW3b\nNrKzs9mwYQMjIyPneylCCCGEEEJcttSF2Jl07ty5/NVf/RX33Xcf+fn5/PEf/zFf/vKXAYhGo2Rn\nZ/Od73yH++67b7qXIoQQQgghxGVhWnvUjTE89thjjIyMUFNTw/79++np6WHDhg1T94lEIqxfv57X\nX399OpcihBBCCCHEZcWfjifdvXs31dXVRKNR0tPTefLJJ1m6dCmbN29GKUVOTs5p98/JyaG7u3s6\nliKEEEIIIcRlaVqCenl5OTt37mRgYICf//zn3H333bzyyivv67kGBgbO8+qEEEIIIYS4cDIzM9/X\n46al9cX3fUpKSli9ejXf+MY3WLVqFQ8//DC5ublYa+nt7T3t/r29veTm5k7HUoQQQgghhLgsXZA5\n6sYYxsfHKS4uJjc3lxdeeGHqa9FolE2bNlFTU3MhliKEEEIIIcRl4by3vnzlK1/h1ltvpaioiKGh\nIX7yk5/wyiuvTM1Sf+CBB3jooYdYsmQJpaWlfP3rXyc9PZ1Pf/rT7/jc7/fHBkIAbNu2DYC1a9de\n5JUI8duR97KYSeT9LGaKt76Xz0f79nkP6j09Pdx111309PSQmZnJlVdeyXPPPceNN94IwJe+9CWi\n0Sj3338//f39VFVV8fzzz5Oamnq+lyKEEEIIIcRl67wH9R/+8IfveJ8HH3yQBx988Hx/ayGEEEII\nIWaMC9KjLoQQQgghhHhvJKgLIYQQQghxCZKgLoQQQgghxCVIgroQQgghhBCXIAnqQgghhBBCXIIk\nqAshhBBCCHEJkqAuhBBCCCHEJUiCuhBCCCGEEJcgCepCCCGEEEJcgiSoCyGEEEIIcQmSoC6EEEII\nIcQlSIK6EEIIIYQQlyAJ6kIIIYQQQlyCJKgLIYS4YOJxQzxuLvYyhBDisiBBXQghxAURjxtqa6G2\nFgnrQgjxLkhQF0IIIYQQ4hLkX+wFCCGE+GDwfU1NjZn6vRBCiHOToC6EEOKCkYAuhBDvnnxiCiGE\nEEIIcQmSoC6EEEIIIcQlSIK6EEKI8+ZyGb94uaxTCPHBJkFdCCHEeXG5jF+8XNYphBAS1IUQQggh\nhLgEydQXIYSYIRLV4Ys1WeVyGb94uaxTCCEkqAshxAyQaOcAqKkxFzyAXuyThPfqclmnEOKD7bx/\nUj300ENcddVVZGZmkp2dzW233UZjY+MZ9/va175GQUEBKSkpXHfddTQ1NZ3vpQghhLgApOdbCCGm\nx3kP6q+++ir3338/mzdv5uWXX8b3fW688UZOnDgxdZ9vfvObPPzww3zve99j27ZtZGdns2HDBkZG\nRs73coQQ4gPBtXNATY1Ui4UQYqY4760vzz777Gl//td//VcyMzOpra3l1ltvBeBv/uZv+MpXvsLt\nt98OwCOPPEJ2djaPPvoo99133/lekhBCfCBIb7oQQsws0/6JOjg4iDGG2bNnA7B//356enrYsGHD\n1H0ikQjr16/n9ddfn+7lCCGEmAa+ryWkCyHEeTbtF5N+4QtfYM2aNVRXVwPQ09ODUoqcnJzT7peT\nk0N3d/c5n2vbtm3Ttk7xwSHvIzFTvNv3srUuQCs1s/vHPyivc6aSz2YxUyTey6Wlpb/1c01rUP/i\nF7/I66+/Tm1tLUqp6fxWQgghzsJaTWNjBgDLlg1eFiH2/QTuy/F1CiHEO5m2oP4nf/In/PSnP2Xj\nxo0sWLBg6vbc3FystfT29lJYWDh1e29vL7m5ued8zrVr107XcsUHQOIMV95H4nL3Xt7L8bhhdNT9\nfvXqS6+H/K1jHU8fM/nu13upv07x9uSzWcwUb30vDwwM/NbPOS2fZF/4whd4/PHHefnll88o+xcX\nF5Obm8sLL7wwdVs0GmXTpk3U1NRMx3KEEOID6/1Og4nHzbSPWjwfYx0T65SpN0KImei8V9Q///nP\n8+Mf/5hf/OIXZGZm0tvbC0BaWhqpqakAPPDAAzz00EMsWbKE0tJSvv71r5Oens6nP/3p870cIYT4\nwHi7TYfea3A9H5snvd1a3mljpPcyQeZib/IkhBDT7bwH9R/84AcopbjhhhtOu/2rX/0qDz74IABf\n+tKXiEaj3H///fT391NVVcXzzz8/FeSFEEK8N6eG1qqqAN9XFyS4ni14v12APtvtZwvlEriFEMI5\n70HdmHf348sHH3xwKrgLIYQ4XyybN1u0VmeE5FOdKwy/XYCORgMAIhFvsuXEUlfnBgW834r2bxPK\nZX67EGKmm/bxjEIIIaZfIrTG41BXd2a7iatkW4yxaK3fMVi/9WvRaMAjj1gA7rwzTn29+7oxBq3P\nbLV5u0r5+Q7WEtCFEDOZBHUhhJgh3KZDTEuVOR5PhPxTR+0qqqvB9999X7wEayGEePckqAshxAxz\nttDswrua/PXeJ8DU12uWLjVUVEBamn/KyYA3dZ/3+rznw8X6vkIIcSFIUBdCiA+A9xrMz/4YTSRy\n5vO51hrXFnMhp6/I1BchxEwnQV0IIcSUs4Xfs/WWnxrm43FLc7ML6lVVAO5C099m8sy5KuXTPd/9\n3a5DCCGmmwR1IYSYwX6boBmPW+BkWD/1OU8P84qyMksQ2MmLWaG52VJWZrn22jOr7++0nnNVyt+6\nc2lin7zpCNJSsRdCXGwS1IUQ4jJj7blDLiQq3e89aJ6cHvPuRy/6vqa62t1/61aL1vZt13a+Z72f\nq+IuwVoIcbmToC6EEJeReNzQ2JgBwOrV56o2v//2EPecice7Krnvn/71s02WSYx+rKqyVFXxtkHc\nGPO2JwGJkP12lfJ3GvF4PqvgMqddCHGxySePEELMYFVVlpqad241OduGSFVVbiRjXZ0669d9X2Ot\nZWBwnNdfN4ClosIQiXj4vuKtXPCF6uqT02feuo7aWqaC9rlGPJ7r9Rhj3vXme+/knb6XEEJMJ6mo\nCyHEZcT3NcuWDU79/tQ2j1PbVlyvuAvEieo0uDAcjbq+8sQOo29XgfZ9hdaKw/39/OyVblq6jrOn\n8xjNB45x5MQIw2MTjERjGGPxtCY1nMyc9BTKFmSQm5JPeUE+d38in4J5GVPfG9zOptNZqT591rsQ\nQly+JKgLIcRlIhF0lTJTfz5byHYB3Z6xa2g8bti0KWDjRsjPt9xzD1OVb2MM0ShEIjAcneCFbe38\npn4fv6nfz77D/edcl681cWMYjI4wGB2h4+hRoB2A//p/YP2VC/j09cspilxBWiT5bdtyzl+riQR1\nIcTMIEFdCCEuA6eG8pSUxIhEy9uH0rffNfRUrsUloPZ1w599s5Xdx3axpb2NiVgwdZ9ZaRGuXlbE\n8uJ5lBVmMXEii6z0dK7/UJjMtAhYiE7E6e0boW9ojH2H+3jqxW6aDx2i6VA3r77ZyatvdhLyfs0n\nq9ZSsuRq9u1NxlX731rFfz/Tac78qcL7fS4hhLiUSFAXQohL1NtNL7FW09SUxuiooqrq9HnlpwdV\n77TH+b7m2muhsvJk68uB3gH+99P1/MMvGzg2NASAUq4K/pGrFvHhVcWMHcnF066/3E2DcReOpiYp\nmBzwEgn7dLVnApl8/Ooccr2lgGVwdIxNzXv5zZu72NHRyeOvv8Gv6htYl381n6q+iqoq/7QLVd/P\nMTrb3HchhJgJJKgLIcQl6GwBNBHAd+w42Tpytskq55pb7vuatDRNZ88Jvvz/vchPNzZijEvbiwvm\n8NmPrOaej6wkPyv95DqOJZ7LhfRYLCAU0pOB3UxNejl1Ta4vXgEp3LJhNX/JaupbDvPl//0bNu7s\n4Df7X6I7tpcbb/1dSvJnn6/DJoQQM4oEdSGEuEwkwrZShmXLBlm9+tyh/Gxhf2h0nL969DW+89PN\njMcCQr7m9z+8jM99dDXe0AKUUmTPOv25EicIIyOGxkaLUoqlS09fWzxuicVcWI9GFZHImZXtqqUF\nPP8//wMPP9LG//zlr2jqOsS6z/8jT/z3O7hmxfwzXu+7mYcurS5CiJlMgroQQlyC3imAKmVOC+Tx\nuGXzZheU37obKLhWlR8+u4P/+o8v0dM3DMBnbljBQ/fdwPyczNNC/ambHVVUxKcuOP35z6Gvz1JZ\naamo8IhE3FSZxP1bWgxKWTo7NSUlhurqk205iTX6vuKLf1jKPZ/8j9z90FM8v62dm7/0Y57/1n+g\n5pSw/l7moUtAF0LMVBLUhRDiEvVuAmgi0BpjaGkBdxGpnbqItKbG8PL2Dqr+7xdoaO8BoOqKAr79\nn25i3dLC075Hot89cYFqLBbn3/4NolHNpz5l8H0Ih6G93f335AmBwZiA2bMDwmFLX5+bOrN5s1tP\nZaWhvl7R3GwpK7Nce60mZ04av/qrz3Dft5/mX55r4KNffpQXvn0XVxbnTr72S2dyi+x0KoS4WCSo\nCyHEZSpRpXbBWlFeDqCmQm77oT7+9Acv8FRtMwDzMtL5zudv4I4PL2fzZkVt7ckdTF31Wk1tjlRV\nFfDKK4p9+6CkxLBlC2RlWQoKoL3dm/z+lsQOprGYZft2RSymuOEGQzyuqKtzF6Zq7Sr6xkBwcpgM\nvqf5xz/9HYZHx/n5q3u4+U8f5c8/8jnSwunceaelqkq97e6mF8r53OlUCCHeKwnqQghxkb2fiq21\neipcu0q4ngzO7nmert3LZ77+FMPRKGEvRHXhOv7LZ9ax4frIWb7/yQtBo1FLJOIq2p4H1dWWVasM\nP/iBxlrFAw/ANde4E4HEvPaKCovWkJ4OsRgkJWna2mB83LJ4sQU8KivdfRKtMuBCr+dpHvnyJ9h3\ncIzt+zv4lzd+xReu+33eeEMRCunTNmsSQogPGgnqQghxEZ2Pim2igl5Xp4jF47zUvpGH/s9rAKzI\nX8wfrPkIV61M55prvKl+8YqKxKZJiro6hTEBExPwk58oystd+N6/32KtJSlJEQpZPA927fKoqUm0\npViamy3RqCUctpSUKIqKFBUVGq3dn9etg6Skkychxpip/vfE601JDvHEX36c1f/xH9jd3c5I5g5C\nocqzHiv3ei9cVVsuVhVCXEwS1IUQ4hKXCKinUspMzTVP9JT3Dw/z3378C5p69qG14it/sJ7rFl2N\nNYp169zc9HjcsHFjQHs7+L6lsFARDrvn9KbGrlsmJgxdXQprYckSxbp1llBIASfHPFZVBcRiltZW\nd7Kgtfu670MopPA8Q0ODBhQVFYaKCvfs9fVnvsaCrAy+94WP8pmvP8H/+y8v0fqvy5ibkXLWCTZV\nVcEFbYmRgC6EuFgkqAshxEX0ThVbF1Dt1LxyUCQn+yjlHpOoToezuvnjH/2Mg0cGyUxJ5t51H8fr\nKiFWZElOVtTXu+8Tj7tgffy4QSkYH4dly1zArq52myFt3w47d2quv94CmoqKxK6famqd8bjbNGnd\nOkUoBMYoVq40hMP2lNeRqLwb/vVfXXvN0qUKY+zkrqneKa8R8v1l3Limgd9s38d3n6jjG5+74Yzj\ncbaKvBBCzFQS1IUQ4iI7V9iMx117iTGwdKkLu7t3Z6CUZcUKi7Xw7/U7+NtfP8dELOCKgnz+4vc/\nyYnedLq6QKmTU1zc91JTYTkW05NVcjX1vbZuhfZ2RXm5OqUC7taX6BePRgPq6gDcugDWrLE89pi7\n3z332Mn2GD15f0tTE1jr1u9OOJgM+6e/9gfv+dBkUN/Cn3yqmqzMlKlj5E401FRQF0KImU6CuhBC\nXKISU13Kylx7S1WVuxBz0yZ3Qej4RJwfbn6WHz2/E4D/9Dtr+eTKDYR9n1uuT1Tp1eT8clfBHh6O\nMzZm2L/fY+lSF8a3boUgcCMU29sVixZZKioSLSonq/nRqOGNNyxg6ewE34eCAoXnnRqcLfH4meMV\ny8stWitiMQBDXV0i/JtTQrgFCrhq8SK2tLXz97/Yxpc/c81pzxOJeNIzLoT4wJCgLoQQl4hTL5ZM\ntLw0N8OiRW7SSiTiEY0GzJ8fJR4EfPobG3l2SyvJST7/8MWPcddNK095Dv8tzxGwZk3Ak08quros\nc+eCq3i7KncibJeXQ1WVmyDz0kvua1/8osXzDFu3BtTVKQoLLWNjLqBXVkJDg6KhQfMHfxCwdaud\nDP5xPA8mJhTt7W4Ge0kJeJ5+y2s+OeIxMUXmtoo1bGlr57GXG7l2wTW4kwWL1lraXYQQHyjT8mm3\nadMmPv7xj1NYWIjWmh/96Edn3OdrX/saBQUFpKSkcN1119HU1DQdSxFCiMtCNBqwaZPr1X7rxaOe\npyYr464/Ox4EfP+VTTy7pZWM5GQ2PvyH3HXTSoDTLr5MCIcD6uosTz9ticUMycmwYoUbu1hfr9Fa\nsWaNJRZzX0/sIJqba8nJMUxMGF57zfDyywqtLQsXKkpLFaWlbhRjonXG8xTt7ZqODsOmTZbaWkNH\nR5y2NktfnwEs1dWKa6/VVFXBqlUBmzdbamst0WhivYo/umMR6ZEIjR1H2H/kyHQfeiGEuGRNS0V9\neHiYFStWcM8993D33Xef8fVvfvObPPzwwzzyyCOUlZXxF3/xF2zYsIGWlhZSU1OnY0lCCHHJcgEc\nWlqgrMzgJqe4qSqrVlmSkvRUlT0eBHx/4ybqDxwkJRzhO3ffyZrSvNOC+cnpKC6MDw9Dd7eb6lJU\nBF1d0NHh+tQhABQNDYr9+wM6O6G11fCZzyhKSixdXZYf/hBGRy3Z2Zryctcr72orrq2mqirR6qIp\nL3d96wcOQBC458/Pd1NgOjs9Kivd69u82dLSApGIYXRUYUxiR1UFhPjQ0iv49+07aD3RxH/4xIdJ\nnAxINV0I8UEyLUH9lltu4ZZbbgHgnnvuOePrf/M3f8NXvvIVbr/9dgAeeeQRsrOzefTRR7nvvvum\nY0lCCHGJU5NtJydbX1yYdbfX1Bgslr976YnJkB7muW9+hurleQBs2uSCenW1C+BBYNm8Gfbtc8+d\nmQnz5lnCYVDKopSlpcVj4UJ3kegbb1jCYcPYmKK727J9u6G21icSsfT3Q2am4o47LA0N8Nd/bRkb\ns5SUWPbu1SxYoKmpUaSl+dTUGPr7obPT0NkJPT2K7m7N4sUwMWF4/HHFkiWJHUoVCxe6XU0TxyAR\nxP+v28v49+07eL3p4Gk/JTjbBahCCDFTXfAe9f3799PT08OGDRumbotEIqxfv57XX39dgroQ4gPn\n5IjGt58NHh0PuPdbT/JUbTPJoRB/uuF6qpcX4PuaaDSgpcWF3cpKJueZJzY0clVwN+9cMTFhGB2F\nnBxLEAR0dHisW2fxPMv4uCYvLyAeN7z5pkdxsaGoyGKtJhz2gIAdOyy9vVBYaBkehpERy9Gjhq4u\nzZ13BsTjlr/7O8XAgGL1akss5kYwZmaC51l27VIEgeWuu/TkDqeKV14xtLa6nUurqgIA1pTmA1DX\n1M34eBzP07/1xlBCCHG5ueBBvaenB6UUOTk5p92ek5NDd3f3OR+7bdu26Vya+ICQ95G4XKSlaVav\n1sQDww33v8mWjk5SIz5/97kqls0P0dCwHYAgCDMykkUsBo2Nxzh0KI2DByPk5o6zatUIWk9gjM/Y\nWDZ9fUn4vmVoaIJo1M0/f/HFKM3NYQoKYmRmhtm2LZW8vDHmzh1h69YQCxeOU1Iywj/902wOHUph\n7dooq1cPsnt3mHnzNGNjPseOWXbu7AM8xsdziESgvLwXpcZZscIHNLW1swiHkwiHh2hqOoFSBms1\n9fU5HDsWJiVlmK4uw/79SeTlTTAvLZWjwyM88evXWJyXwf79GQCkpAxOPdZaPTVTHpi6PfF7cfmQ\nz2YxUyTey6Wlpb/1c8nUFyGEOM/OV1BUymBswF/8dCdbOg4RCfn87efWsXxBJomKubWa5uZ04nFN\nXt4ooInFfIzRdHUlEwSWFSsmABgaCjE46JOWZsjIMExMeBw5EmJ42Cc3d5zkZEVXV4i8vDhJSXEO\nHEhieDiJ9vZkUlM11ioKC8eIxeCf/imbFSvGKCiYYGBAUVAQZ+PGLEZGQoTDrrWloyNCcfE4LS3J\nWKvIyoqSlRVlyZJRGhvTUcqydOkwRUWjLFgQxfehoyMJpTQ9PRGK52VxdHiElu4TlOansWzZ4NRx\nsVaza9csurpSWLRojIkJM/V8TU1pACxbNnja34EEeCHE5eaCB/Xc3FystfT29lJYWDh1e29vL7m5\nued87Nq1a6d7eWIGS5zhyvtITKdTt7pPbBAE7+4iyFPHMyZ86e9f4PmGQ6Qlh3nmr+7k2ivnn/Ze\n7uub4Je/tCiluOWW2fz854o9e2D16oCJCY9YbC5lZUUA/OpXASMjUFxsmZhIIjXVbVQ0OqpZvDiV\n+nrXHlNUZFm1KkJjo2VoyNLVpTlwIIvKSrj9ds2TT1rmzLEMDobp73d96OPjrvUmFoO0NDh2TBGN\nJrN4cQE7drjXc+ed7gLT7dstExOasjJYs0azZg1Eo26Oe24urFjhxkhmHMwDOvFT55zx/200GrB9\nu+HYMUtRkaWoSBMKeVx5pRsdCbB69elTcE79e5HWmUuHfDaLmeKt7+WBgYHf+jkveFAvLi4mNzeX\nF154gYrJbe+i0SibNm3iO9/5zoVejhBCTJt43L7r7e5PD5Luvo+9tJv/+fjreFrzP37/DqqXFhKP\nn2ztGB6O8/TTMDoKpaWuxx0Uc+caFi4Erd00mccfV5SUKPLyNFVVbjOjJ55QHD3q+tmLi11Y3r8/\nTlqaJRSCSMR9LSnJBW/fV5M7ilrS0y2FhZZIBBobFdZasrPdLqj5+ZasLMP27YqmJk1eHixcaPE8\nAMWjj1oOH4b16wMqK72pC2fr6zXGGKqrFb7vs2BBwNw9bgrYtl1DZ1xE6vtu99RIxDJ/vqa6Wk9t\n6iQbIgkhZoppCeojIyO0tbVhrdvR7sCBA+zcuZM5c+ZQVFTEAw88wEMPPcSSJUsoLS3l61//Ounp\n6Xz605+ejuUIIcQFc/LCUEiMFDzV2armp3M7e+7u6OGz3/oFAJ+/+SZWFy+cCv4HDmQSj8PwsCUa\nVaxebRgagiee8PjEJ+Ls3g1tbR6xmKGnB0Ihw/i4prxcs2aNoqEBUlMNkYirpltrmZiwFBdbensV\ns2drtNZceaXhd383jtYQiWh274bHH4d9+xSep7n5ZsOJEzA0ZMnOVqSlKfbt02zZ4irs8+aB1m4s\npNaasbE4fX3uz6GQor5eUVUVTO1iqnUibGvWrTM8WZsOwNHBwbc5zlBV5WbMnx7izzy2p/69SIAX\nQlwupiWob9u2jeuuuw6l3IfvV7/6Vb761a9yzz338M///M986UtfIhqNcv/999Pf309VVRXPP/+8\nzFAXQswIpwbBU8Ph2armpz6mqsptAPT8y1G+8KPHGRuP84cfWcW3/2Tt5EzyM33qUxbP83jsMYUx\nAc88Y+nu1hQVucCfnm7R2nDsmGHePI+6OqZaTBYtUnR3G3bssLz8sqtSe55i9WpFLGbYv9/w0kua\n4WHNbbe5Sv7x44rUVEtVlWHNGg8wbNyo2LZN8+EPG+bPdzPa8/Lgrrs0KSkedXVubT//ucIYM3lM\nfNraDMZAZaWioiIxi939uxGJaBYVuX8TVCh62rE69WTHfw//iklAF0JcbqYlqH/oQx/CmHNfrPPg\ngw/y4IMPTse3F0KIS8Y7tbtEowbfV0Qi3lRl+aEnf0FHzwkqy/P52/98C0HAVAtNVZUlJWUAYyLU\n1xtaWjzuvBM++cmApibDv/2bJiPDcvPNsHs3HDxoaWrShMOWvLw4g4Me4+OWykqAEEEQ0NMDxlgW\nLLAsXaq46ip49FFFVxdo7cYlvvEG5OQE1NRYurs9jh5VPP10wMKFAYWFHseOWQ4cUMyfb1m4UAOa\nlBRNJOJRVRUwOmrZsweMgfp6TXa2YWxM09bmNnryfUtJiUHrEDU1bnfUJUvccXLjJk8es3ca0/jO\nP7UQQojLg0x9EUKIC6iqylXG43F45RVDba1lwQLFnXdCJOKxrfd13mhtY3Z6hMf+39+jfqtPLBYQ\nCrmdP31fYa1m+/ZZdHZaFiyIs3mzYtMm1x+emRmQnq54802P9nbDggUBsZjmyBHN3LmQlOSmvXR2\nKm67LWDpUjdTfft2SxDAwICbda5UnOPHISdHs379BC0tPrW1Hu3tlsrKgMZGzZ49ivR0j+XLDWvW\nKPr6FG1tirIySFTGExs37dunWLbMEIvBpk0exlhKSxWeZ+nstOzaZWhshGXLAlatUjz6KLzR7I6Z\ntXbqueJxy9laihLeTZAXQojLhQR1IYS4ABIB0hg7WSG27NtnOHHCTVkB2Np8iP/2Ty8C8E9/ehvz\nszPZ1xTQ2grl5ZbKSgNorPXZsSOVuXPdBJfWVk1ysmV83DI0pElKUigVcPAgzJrlc/31cfbuDdi+\nPUR+vmXNGrfR0aOPGjxPkZxsaGvTxGLgeQH5+a5tJiVFoTU0NoYxxjJvnuH4cc2JE4oFCwIOHPAI\nAsjN1axbp3nsMejuthQXG4JA8dprlnXr3Gvt6bEkJVk++lHNunWn95RXVhomJhSHDyuKik4eM63P\nPH6gqKqyZ/SlCyHETCRBXQghLjgXzMvLNeXlhjVrNAMjUe74i58Tixt+t6qSOXYJ0ajB8yy+7+aS\nb92qJi+KRPoLAAAgAElEQVTsNBQWjpGVFSEU0pSVKVautDz1lE9Wlp0cd+imuoTDhuZmn6NHLUVF\ncXp6fMDgeYbWVlfRnjMHsrMDjh3TtLVpMjMtnqdZvjxgdFTR1OQuFL3qqoCBAcWRIwqlfK66ynLr\nrTB3bhiARYvitLdDc7MFAo4e9fA8t57WVsOWLYqhIbj3XncB6MnKt8e992rGxw0NDZqGBsVnPhMw\n+HgADeB7pwfyc4V0uWhUCDGTyKeYEEJcAO5iUdcbbgxTF1Fq7dPQoPi9L73I/p4TLM7J4cPzr8eY\ngPp6S0eHJRy2eJ6duqBU6zgf+UgPd9/tEQp5hEKK5GSPsjKPvDxNcrIiHPZYuVKxaBGcOAFHjrhR\ni5WVExw86JOe7lplrNXMmWMpKoKUFIjFFF1d7oRgYsInL08xb55hYgKOHfMoKTGkpEAQKPLyFGlp\nPr6vJnvoNWvXuqp+erqlpCQ+teYFCxSZme5YxOOWeNxNBTPG9ekDpKZ6JNpakpI0IxPjAMzNSJ6a\n8vLWGeiuHeb0a6LcRabyz5sQ4vInFXUhhJhmid7qxAWhYCdHEbp+69bDh6ndtwNPaf7rJ29j4Tyf\n5csDfvlLRWcnFBYa9u3zWLwYVq0KaGlxYd1dfGppboZo1FJZaTEGDhxQ5OcbPv7xgJdeUhQVxYlG\nFXv2+FhrWb8+4OhRw9CQZtEig1I+sZhi1aqAoaGApCSP+fM1nZ2u/eWKK+DAAbexEcDAAGRmGvbu\ntezcCZ/7XIAxPpGIx803W+rrLe3t0N3tNjLav18DihtvdK9369ZEW4shCCyPP66Jx+Gee6Cmxn0P\n3/eYnTMKwLxZqZO3nR6+pR9dCDHTSVAXQoj36N1OFTk1oLvedABFdTVTm/OsWRPjy//l11jg9qvW\n8qmPZhMElsce03R1GQoKXLhNSTG0tEBjoyY9PZMlS4bYvNkQi1mKi6G2Fo4eNcyda+nsVLS0WDZt\ncpsh5eV5zJtnGBmBINA0NRl27/bxfcsVV4AxmtJSsFbR0+Mxd66bg37kCBw5oqmoMBQWBrz8sqal\nxY1ETEqy9PXBwIDhxz+GykrLqlUBW7cq2trcOMiODo9QyFJSEjAy4sK6++U2Z2puVlgLqamWeFwR\nj5/ee9511O3qNz8nc3r+IoUQ4hInQV0IId6Dd1vFPf3iUbeZj7sI0gV0gGg04IG/bOSNPV3MTk3h\nng9dQ309xGIunOfnu2p2R4fi+HE72XNuSU11F5S+8YY7Ybj3XsuOHdDY6C5U9X3IyTG0tGiMsaSm\nWm6/3cfzXL/5c8+5dS1eDEVFmp4e8DzL5s2uGl9WFmfLFp/xcUt5eYzNm0MEgYdShkgEkpMVBQWa\n0tIYv/61T1uborIyzvbtHq2tLnwbo7niioCJCU1/vyIpybB3r6agIODAAUV5uaGkRBMO68lRkW5s\nY+K4AnT0nABgYe6ssx5j6UcXQsx0EtSFEGKaVVRYIhE1FdATFfnhsQn+baeb8vI/PnsdyaEITU2W\ncFixaFHAFVdY2tpCLFwIY2PuAtBZsxJBXnPsmJ2ajFJQ4P6bmqpZvBhWrLC8+mowOVJRsXevpqrK\nY+nSOM8/b8nJsdx+u6WjA/r6LNGoZXwcysri9PYqQiHD2Jiivd0jPz/g+HFNXp5hzx6f5GTFjTfC\n3r1JFBUZjLG0tiqUUixdqliyBPbuddXzXbsUx465tWVkuDaY48chHFaUlEB1tZshf2qfeeKnELvb\nXEV9wTkq6hLQhRAzmQR1IYR4D95NFTcROquqXOCsr1fU1Ly1FSbgV3teY2BsmIqyfO772Go2b1aU\nlgZAwMaNimefVSxfHvDJT2oaGjThMPT3Q1oagCI93V0oum2bZWxMkZXlcdttlt27Fb/8pU9xsaWg\nwHLwoObgQcPYGExMuAq7UvDEEx6NjVBUFJCdbSgqUvT3K3bv9li4MAAU+/dr5s+PkZ8P4XCIBQss\nAwOWPXsUWmuuvdawc6fiwAHN+vUBS5bAP/yDJhwOmJhQpKUZrrnGsGWLT2qqoqzM9aBPTLggn9jk\n6dTjGo9DEAQc7u8HYEHO2SvqQggx00lQF0KI9+iddhtNtMZUVcHJjX9OXkwai8X5xfP9/H19HQDf\n+8ItaO3mg8fjivp6TShkGBhwLTMnKRYvtvT3K2prMyktNbS3a7Zu1cyZE5CUZGlq8mhrg6GhODt2\nKObOtSjlquO7dgUcOaJITVUsXGg4ftwF9oEBeO01TXl5jDlzPJSypKd7VFTECQLF1q0+OTmWRYss\nc+YYursVb75pKS6O09sLg4OK3NyAQ4fU5AmBndyYCLq7FYODHrGYZflyWLNGEwSWXbtOHptTj2s8\nbqirU3QeO8ZQNEp+Vjp5c9POw9+aEEJcfiSoCyHEu/B+tqX3fVdJd48HNz9dUVmp+H9++CJxE3DX\nhiupKMtn48Zg6uvt7YpFixTr1xuWL7c0NHgYYygsZHKzJDh8OInkZM26dXD4sGHHDnfBalJSnOJi\nS3c3dHRo5s9XLFpk6O/XgGXxYrfL6d69HsnJAWvXBsTjcOiQZtOmELNnK0pKFMnJ0NzsA5asLHcC\nMT7u2m583zIyAo2NYK1m4UK3QVJ9PaSlKTIyIC1Ns2JFwMGDlr173fz0+npLd7flxAlISlKsWQPg\nnfXY1e/rAOCG1cUo9fY7kZ6vvyshhLgUSVAXQoh3cPLCUEN1tSUSOXu4hLO3xiSqxO7xsPHNTt7s\nbiU9Ocy3/tMGolHDpk1u+smVV8KiRW4jokOHFK2t4HmGxYsVSrnnXbhwjNbWZECjtSUvDw4fdjuU\nHjlimDMnYHzcIxyG4WFDf78lLc0yPKzo7tbMmmVQynD4sCYlxTI87DEyYsnIMIyNeeTkWIaHYWgI\nFi0yZGRYeno0R45ojh3TjI8rRkYMvg/Llxuuu07xzDMaz7PMng0ZGa5q3tGhGR2FcBjS0iwjI4rh\nYUtXl8L3LaWlcd4a1BPH71vP7AfghjXF7+vvCmRkoxDi8iefYEII8S4Y48Yj1tWdrNiebbMdePsN\nd9zsdMVf/mQTAH9+13py56Th+4qSEsvEhJu+snKlobXVfa/jx6GkxLJ0acD27YZt2wzNzWnMmgVL\nlwbs3w8HD2r+8A8VVVWG7m7YtMmjqChg1qyAcNhOtp8Y8vMNAwOKsTE3xzwtDQoL3VjG8nJLQYEi\nPx/S0gx79lj6+gx79yreeMMjEglYtCgOKEIhyM8PyMqyHD8OL7+smD8fVq3SrF3r8alPKbKyLD09\nlv374ehRt/HSF75gyc31WLjQhfZXXtH098eIRoPTjtNYNM7LDZ0A1CxbcMZxfLvjLoQQM41U1IUQ\n4h34vqa6+uQcdHhvldtTq+yNHUfY9OYB0lPC/NHH1wIQiXjcfrvhqacCtm9388UHBxU5OVBRoais\n1NTWBjQ1ucp0Xl6ErVtTcNdaGmIxhecpSkstdXWa0VHFsmVxRkcVmZkBfX0+R44ojIlTXR3j4EGf\nw4c15eWWvDyfEycMIyOK48c18+dbDh6E8XGw1jI+rkhNNTQ0eGRlWZYvj3H8uKKtzV1YOjTk1vsH\nf2C5/XZ3DILA0tioiUYNy5ZZPE8xMuLT2Ai5uZbcXI0xoLXlZz9zPz24555gavrLvzzRzdjEBLkZ\nczjYls7CPHPaTyfOddxlZKMQYiaRoC6EEO9CJOKdFgDfa0U3ERq/99RWAO6+aSXJ4dDU83ieoqND\n099vWboUqqpcL3tSUiJsakpK4qSkGIaGPHwfhoctnudaWl580TIwoFAKMjJcD3l6umVw0MMYRXq6\nq3739ITwfU1+vgvNsVhAa6vH0JCirMxw7Jiir0+TnBxw+LC7sLSsLM7wsEd3t0dfn2LevGBqV9XC\nQkhKAs/TbN1qqK1VZGcbQiHXv/6xj7nXtmULbN4M8+cr7rjDY906mJgw/OxnZ4bpho4OAKrKFk71\n5L+fYy2EEJc7CepCCPEunRoA364X/a33O1Xf4Bg//s2bAFxVsJZNm8xkEFVUVEB2tiY723LHHScD\nel0dgMHzAtLS3M6iQ0Mxfud3Rhkens2OHZpYTHHokJsGk51tWbIkRnt7iKIi9/y9vSc3XurrU9xy\nS5xZsyz//u8hUlIgPR0WLbIsW+bGNWptSU2FEycs0Sj09fncdZfhkUfcZkazZilGRmDdOqis9Bkb\ncxfCPv205vhxS16e5dprNaDZvRtaWhQlJQGRiNtQaWLCkJHhk5bmc+edcYCpano8bmk56vrT15QU\nY4w7ITjXcRdCiJlKgroQQrxPpwbFeNywaZMLkNXVFt9XZwTJHz63g7HxODesLqZo7pzJWy3GuBaO\ne+89GUhfe81tJHT8uOuN933L6Khi/nzNqlX9vPTSbOJxjwULLLGYC/Bz5ypWrLDs3h0mHleTbTKW\nj34UXnsNamt95s8PeOEFD2uhsDDO+LirwiclaVat0lRUKILAtdDMnRvj6FG3y2hXlyYnB44dUxw6\n5BGPw4EDEI/HePVVxfi44uabDeGwIhpVjI0prrlGUV9vSEuLY61HZaXlwAHD00+7STFVVcHkbqSW\nioo49fWKnbvHeaOpC60Va4qL0VpPjno82ebyTvPr3+2JkxBCXOokqAshxHkQj1taWsCNWLRo7Z3W\nQx0Ehu//YhsAN5RfBVjWrIGtW6G93V3cWVPjgnptraGuzjI2ZohELNYq4nEXbn//9xXNzTGOHAmR\nnAwZGW4To6QkRW6ue8zBgz5z5xp27YLf/EZTWgp5eZoNGwyHDyu6u91FqvPmGQ4OdZEeLCQtzbJl\ni8H3PSorFVpDbm6YG28M+Md/hCefhCuuMMyfr9HavZ7eXk1bmwvss2cbliyxJCX5/OIXcOiQZdWq\nOC+8AD09ipqagI99TNPaqunpUSxcaBkfdycpzc2aWAxCIWg8vI+4MVSU5fGRG5NPmz9/rmsB3tq7\n7o4j7/g4IYS4lMknlxBCnAe+rygvd79O3eQoUdV9bksb+7r7WZCTycqiYlpaoKEhUUFXwMlNgkCR\nl2eIxVyrytVXw/r1iuJid3+lDBUVQyxdqjh0yP3KyjIMDRl++lMfYwzLlsUZG9MEAfT0wNCQYWDA\nMjoKGzZMsGhRQH+/x479h4jFFJFInF274rS1BdTXG3p6DJ5nSU521fPRUbeWwkLLwoWKzs4QsRiA\nG/0YiVhefNEjFnOPS0qy7Nnjxjz6vmVwEJqaFKWlcPXVbmLNo49atm61pKQEeJ7ry9877E5m7r5p\n5eT0nPfeoy6EEDOFVNSFEOI8ONk7feZOpDU1hr99cgsAt1y5lkOTc8TBbVi0Zo1h+3ZFXZ2mpsZV\n1let8njqKRfcr77azSX/X/8L3njDkJc3G6UCbr7ZEotZmpoUw8OaiQlLOAyZmYrU1DBXX23p7w+I\nRjXhMOzapfB9VzXvtwfYdriLHf2bmTPHsr92AVn+fG65JUZbW4jmZsvQkKWy0udDH1L098c4dEjz\n2muaggJLfj7k5Sm6uy3j4xCPK7Q2eJ4iM9P12oPPJz5hGRoK6O31aG5WFBYa0tIsTz1lOXrUXfQ6\nMqLIzYW9XUd5ZWcnaclh/vAjq95yXM/dwuL7mqqqYPL33tRxf6fHCSHEpUyCuhBCnCenB8KTU2H2\nHjjOr7e2E/Z9bl65gt4uTXGxYdUqQ329pqPDta8UFwdEo4pIRJOZGeKOOxLBUzEyEmAtDAxYkpOT\nGR+HH/1IMTGhKCmBkhKPSMQyb16MUMjNOs/NVeTn+5SUGPbs0fT0GDo7PU6cgFWrCimePZ+kN+G2\nK9fR3u4zPg7z5mkgoKlJ0dWlGB6OMzGh2L49RHZ2nLExQ1OT5tZbA8LhEACtrYZwGA4edP3u114L\nEKKiwl0IWl/vE4nEaW01bNmiiEQCjh5VrFwJ8+ZBZ6c7kfj+U4lq+pVkpCadcVzP1XOe2FQKTra6\nSEAXQlzuJKgLIcQ0OLUS/PDPWgBYf0U5N34oFYCtWzXbt4PWhn37FDk5bqThT37ixiRee+3JSSiu\n19rjP//nGL/5Deza5ZOSAv39MGsWlJbCunWu5aa93aO2FmIxKCkJGBjQRKOalSsNCxe61pdYDLq6\nNMPDHkXpRbS0hCgrM8ybZ+js1Bw96jZfysoK+NGPPMbGYNasgGPHwBjIzjZ0dbme+Wuu0USjbvb6\ngQOKEydcdbyyMmDrVoXWiooKV8UPhSy+H9DVpcnPt8yfr6mu9ojHLS9vivKj591EnM/fftUZx1N2\nHBVCfBBJUBdCiGmSCJMv7+wA4KrFi6iv11RUGILAbQRUUaGIxVwbzKljCIGpcYWO6xfPzbW0tk4Q\nCsFtt0WYmFDs3+/heYbKSk1xsaajw3DsmCU52TAx4arcfX2Wnh6PzMyAuXPj7N0bor9fsbRsAcuW\nWU6c0Bw5YmltVZNrt2RkBOzf72OtJRRSBIHiyisNRUWaffs0QWDYudNjzhzFddcZ6urcjqpdXdDX\nZ2hv10xMKGIxxYc/rKmpcWMo6+vdK6qqUlM96C+8uZtoLMaHVy6krHAu8fh7C+MytlEIMRNJUBdC\niGk0EQt4dWcnAGuKizHGUF/vpqLMn2+JRHw+/GGmAnlVlXvcyEjA1q2u5xuguVlhjKK6WpOScoR9\n+zJ59VWN71s8L87rr3toDdXVHsuXW375S8OuXSGWLzdkZcGRI5qxMRgc1IyNhSgsDCY3MvKJxRTG\nQH+/wvfB8xKr91m/PqC11W2mNDCgGB/XXHGFYs8eaG/XlJVZQiH304GGBkhPD2hvV4TDkJ9viMc9\nPE8Rj7uRlWlp/lSgBjeZxRh45k23EdQN5Yn58vqMynlV1dnHXk6tVgK6EGKGkaAuhBDnSaICfmqY\nrNvTxUg0xhULsrj5ulSCwNLQoNi3T2GMYc2agNRUb7K/2lJRYamvh4aGOLt2KXJyLDfdpNHazTd3\nVXefY8eSGRszzJ5t6O31SE11s8bHxxW7d0NPj2ZoyG1aNDKiSU5W3Huv4dlnobcXlNIcOqQIhy0d\nHRprASwlJQFdXQpjNMuXGxoaNN3dltRU1y//5psexrjdSBcsgOpqzfi44Y03DEEA0ahiaEiRkwO/\n93uacDixcZM74aiqCqaOT6LnfNu+dtq7+5idkkF1WelZj6tre1FUVZ2cqf7W3WElqAshZpqLGtS/\n//3v8+1vf5vDhw+zbNky/vqv/5prrrnmYi5JCCHel8SGRy0tUF5+so/6xe1ul81rly/k8cdd1fwT\nnzDEYpZXX1U8+ijcdZcLnB0d7vZQyDIxYVDKBfRVqxSrVrmAv3lzwK9/nUt3d4TcXBeYQyE3PaWz\nM05jo6a/31JcHJ9sVwGtAw4f1mzdCocOwfHjitWr4zQ2hohEFB/7mCElRVFX59pjSkoCUlIMLS16\ncsa7ZXRUk5lpSE9X9PR4LFpkWbfObUb0k59Yjh5VrF9vWLnSY+dOCAIIhzX19RpjAiYm3AjGujoX\nphNV9aoqy7eecReR/u66NdRcnSjn26npLSdZNm+2aK2oqgqmTm7crqtnVuCFEOJyd9E+0R5//HEe\neOAB/vzP/5yGhgauvvpqbrnlFrq6ui7WkoQQ4rz7Tf0+95v+BRw/bujqsmzfrqis1OTluSpzEFii\n0YDWVsPmzZbWVhgc9CgtheuvVyQluZDd3GyBAK0VSUkGrdVktdxtHDQ4aOjuNvT1GerqPLq7NSMj\nEI0aurpg82YYH3cbDQ0PWzIzLQsWQFqaa385ckSTkgKxmKajw21ClJwMBQWG9HQLeKxdCzk5MDjo\n1l5f7yr0QQCJf1JWroRQyKO+XhGLBbS1WY4eNXieIRo1GON+8lBbC//880M8s6WVpJDHf/+jNfi+\noq7O/Tp1yktNjWsL0lqCuBDig+OiVdQffvhhPvvZz/LZz34WgO9+97s899xz/OAHP+Ab3/jGxVqW\nEEK8L76vufZaqK4+2foyPDZB3Z5DaKVYXriAK0otQaAJhTxSU10lvb4etm41bN9uOXIEFixwFeO8\nPM2iRYrKSsXoaEBzs6GnBzxPUVAwTkXFMBkZ2Rw+HLBzp2J0VOF5AampAX19HsPDiuRkVwmfPduQ\nmWnJybGkpBgGBjyamkIsW2ZYtkxjjGXvXsW8eQGVlQHPPBPi+HEoKYlRXq45dEgzOKiwVjMx4dpm\nPM9V/cHn+usNS5YYdu+Ghx5SZGbCunUBoZDHmjXuxKOlJdG643Y3BYUxlu8++xwAf3rH1eRlpZ/R\nznLq8fX9U2eje6fMrVdT9xFCiJnkogT1WCxGfX09f/Znf3ba7TfddBOvv/76xViSEEL8VhIBMxI5\n2a7x6s5O4oGhckkBV5RG8DzFNde4CzZ9XxOJJO5pyc1VaG0pLoacHMWKFe6izq1bDXv3xunudicA\n7e2aoaFkYjHLvfdqtm+39PUFdHQoxsY8kpMNaWmKDRssIyNu5nlWlqaiAnbv9ujrc60tx4/DiRNu\nE6TMTBgZMXR1KcJhjw99aJydOz2amjxyc2F0VE9dlOp5HmVlAbW1isceU3zqUzF+9jNFU5PG8wKi\nUUtOjpuRvn2728hp0SJLYaFhYEDR2qopKHCz4dtGGtjbfZiCrHS+8hnX9vhO01tOvU2CuRBiprso\nQf3YsWMEQUBOTs5pt+fk5PDiiy++7eO2bds23UsTHwDyPhLnm7WaxsYMAJYuHUYpg1KGR59tAmBZ\nQTJdXZ1Yq0hNHULr+NRjU1Jc2MzI0CxdGubVV2fT0BBm48YA3zekpir27EllYkKzenUU349z9GgS\nkYiiubmdbdsymT1bk5YW5vBhyzXXjDB/viE/f5R9+9KYPz9MTs4EkcgIAwOZgCY/P874eBodHT7W\nxli5coRZsxSHD6exY4fHvn2aefMCZs8OJuetxxgZiRIKTXDokEdpaYzR0XQmJiyvvTZCW1sqWVkB\nsZhPEHjMnj3I88/Hqa/PYM6cGJmZhokJTVbWOIsWBWRkDLJp8xj/5e83AvBHNy1iT+ObWOuOhVJn\nr6qLmU0+m8VMkXgvl5aeeXH8eyVTX4QQ4rdkrZ4ao9jYmI5SlmXLBuk4MgzA0sJZXHFFlP3702hr\ny6C09MTUY5UyWKsnw32U7Owow8NJjI2FSEsbJxSyZGe7Cnh6+gStrUkcPeoRjSZz6FCYjAzYt8/D\nWpgzJ6CuLoW8vBhbtmTQ3Z3EokUxWluTCYVSWLlylOPHQ7S1pdDZGSIaVSxdOkFl5QnAtcmMjKQy\nPq6JxQzNzSFmzTJUVUV5880IIyMRysrGiMUMc+fGOHEiRGdniKuvHiIWU+zbl8asWXHAEI9DTs44\nBQUTHD6cRE9PiJycKGVlg2gd5x9faKV/ZIJVxbO5aVX+aSc7y5YNSlgXQgguUlDPysrC8zx6e3tP\nu723t5fc3Ny3fdzatWune2liBkuc4cr7SJxPidGBhYWGykqor3dV4dWr4UT0DQBmRVbS3DyPHTsg\nI0Nx661FNDS4FpmKCsPWre4iyaoqS1FRjO9+122GdO+9EX75SxealywJUCodYxQTE24UYkZGEqOj\nbj6678PcuZCVpQFDNApz51qs1Rw5EmLevBCQTCymmDMn4MQJRUYGXHVVEkNDiwmFFOXlcTIzAzIy\nLP39PtGoIiPDY/bsNKx1ffUrVkT4/9u796iqrnvR49+51gI2bxRQBARBJfhGEbdIMKb3WGua1LRp\n0hirbdqOmj5MmuRk3N6epk1OTk7bk47Te5umqacn5zRpY549TZM2idokJipgEUKiKL7wwUPwCQq4\ngb3XvH9M2ICiYsT37zOGI3uvNddec6+sscdvTX7zN2fMiOf3v9c0N2vGjo3gU5+y2LTJpqkpwIwZ\nGseJwbIsbrvNpOscP+7nL3+B1tZYJk8eyfa6Q7xc9BaWpfjvH9xBzpgk/H6XtjaC107SWq4d8tss\nrhYn38vNzc3n/ZmXJFAPCQkhNzeX1atXc9tttwW3r169mttvv/1SdEkIIc6LZZmc84IC8962FXsb\nzY/08NhYjh40ueDDhxMcfe9e/Gj7dlNOcd062LvX4tgxkw8eGmqRmqqprFS0tCiOHLHxeFyyswMc\nOACgiI9XJCZqIiIgOdlm8mRTK33YMJfhwwM0NDiMGKEZMcKlvl7R0uISF6eYOlWRlAQhIQ47d4Lr\n+jl8WBMZaR4mGhtthg83izKVlYUwZYpLUpLJYQdNdrbCslza2hTl5YqZMzWuawG6qySkycUvLtZU\nVVmEhmqys0Frzdf/dRX+gMs3b84lZ4wZnJGVRYUQ4lSXLPXlgQceYMmSJeTl5VFQUMDTTz/N/v37\nWbp06aXqkhBCnLPTBZiHm9toOdFBTGQYc28M7ypfaCZRejymYonPZ8obOo6mrk5jWdDeDlOnKsaM\n8fOnP1mEhChmzNDs328zdCiAxbFjATo7bZSyiIx0OXIETpyAzk6zKujhw5oDBxSHDzvcfLNLIGDR\n2Gjx8cdQV6dIS9NMnqwBM1o/ZoyiogIOHoT9+80IfWKiwu+3iI/XpKUpRo2yGT/e5ZVXXKqqLG6/\n3SUkRLFmjUJruP56RV6epqxMddU5N6PpYKrYZGaC16t48+/b2VhdTbTHw6NfnRO8Xr1LMQohhDAu\nWaB+xx13cOTIER5//HH279/PxIkTeeuttxg5cuSl6pIQQnwi/QWX9YePA5CaEM3f/25G0AsKVJ+2\nZWUWnZ0B0tPNSHpHB4wZowkLC/Dee2b10lGjYNkyiy1bzHHjx7s880w7o0Yp8vJCKS5WNDSA60J0\ntBlZtyzweEwpyK1bHUJCzOTQ664LcPiwTXu76U9JCfh8MGuWxuOxycgIEBbm4vc73HEHbNoEYPH5\nz2s2boSVK03N9MTEAOXlLnV1NsnJplINQGkpXQs+mRVUe5esBMWJjg4efHoVAP/yjTkkxUcBvVce\nRRYtEkKIXi7pZNJ77rmHe+6551J2QQghLoiDTSbhOj4momuhIrNgj+OYwNTvN8GrbcO0aYq8PBPo\ngtaNQsYAACAASURBVGLPHpcDBxQxMZpp0yA83PxUb9mi6exUJCZ24DjQ0mJx4oTLsGEBYmPhwAGH\n+HiNbftxXcXhwxaNjZCS4qekxIx0Z2QESEoKYeJEKC3VtLWB36/IyjL1zsPDLUaOVAQCLhs2QGur\nSdFZuxaUghtvdKmuVrz+us348Zpx4xRg8u0tyyIrK0Bnp1mBtCfodlm3TvP4/7zJnoYmcsYk8e1b\n805bM10IIYQhVV+EEOICaG71ARAX5WH0aN21VeHzBdiwwYxoT5rk55VXFFVVMH68IiREMWlSgNJS\nh5gYP0OGuJw44bB+vVkoKCIC9u5VHD4cRnS0y969Jn2mpcUmPl4xa5aio0OzalUI4eEuI0a4DBli\noTUcOwauq4iMVHz2syaP3HVdwKxGOmUKrF+vOHBAMXJkgD/+UVFbCzk5JphuazOTT8eNM4sXxcdr\nMjMV3auROo6ioAD8fiv4/Xp7u+Jj/rZpM5GeEF58+DbQ9BpF78ntl9F0IYToIYG6EEIMMr/f7VqF\n00yqDAkxkyyLi7u3aaqqwOczixyZ/HUT3Nq2YtYszbvv2tTXQ0aGpro6QH29IjVVkZ2taWrSXbXU\nzcTQzk6L1FSFZcHBgxqPxyxSNHQojB0LFRVW16qkAZqaHFau1JSXK/x+k68eCFiEhlokJmqGDw/Q\n0qKor1dkZQW46SZFaKhFXZ3Jdf/4Y8XBg5rkZI3Xa1FeDqDx+/vm30NP0L2j7jBPrnwLgKfuu4nr\n0hJOGU3vL0CXvHUhxLVOAnUhhBiAMwWNvfeZfGsdTHcxukeYNZZlkZvr4rpme1aW+W9Ojqa42OX5\n52H0aEhNVSQnQ1wc7NljKqhkZIDPZxEXZ6q57NtnU1jokpZm0lreeQfCwixiYwNERmpSUjSrVztE\nR2sKCgK89ZZDRITG7zcPCLGxmowMF9u2cBxFVpYLKDo6LJKS/Bw8CP/5n5CQYCq63HknfPSRorNT\n4brmocKyLFzXpMmAxusN9Fmdtc3XyZf++VXafJ0snjuZr3wmJ3it+puE2zuAl7x1IcS1TgJ1IYQ4\nizNNdjx5n99vgvTaOhOoK6W60joU3QG74zgUFpq2GzaYHO4NGzTFxdDS4mJZinHjzDGdnRbx8WaS\n6aRJ8MorsH17JD6fGSWfN88iLs5hwwYYPTpASoqmuNji0CHF7t1+Ghs1x44pZsywSUhQDBum+Nzn\nTA764cNQUeEQHQ1ZWR1s3WoBigkT/CQkuOzcaRMSAlq7xMfbRETYFBYq8vJcHEfhOIrcXPPXg4oK\n3VXmUZGfH8BxFLat+PoTr7Op+gBjU4fy1Pdu6nNdTw6+u6+l65qa9N1pNb3393ecEEJcrSRQF0KI\n8+D3a1zXjJQDXSPTmkPavG/v9AcDS7/fxedzcRxTZcVxwOsNUFzssnOnJjXV5fhxAEV1tUmNuesu\nyMuD8nKLjz7SZGdrmpvbse1Ohg+PYfNmm+uvN7nkPp+io0ORnGz6VFtrk5Dgkp4OEyc6bN9u4fdD\naKhZTTU2NkBysmbXLnjmGWhtdUlOhjVr4PBhi7lz/YSHW+zbZzF6tKa93cVxbKKiHPx+l7VrXaqr\nzfc5flwRFmZSfYqLTUnGdXvX8+K7m4kKD+W1x+4kOiKs13XrP+h2XZft2+lT4rHnLxWmjYywCyGu\nFRKoCyHEWZwuTcNMDO0dVNr4/S75+QrPsBh4CRqOtAAmMF2zJsDatZCWplm0iGCKSFWVoqkpwJgx\nGsuySU2FhgZFRISmtFR1pZiYdqNGQUrKCbS2OX48luPHNSEhAerrXY4e1VRVKeLjXRISNK2tikOH\nFDU1is2bNaNH93wny9JorUhJMTXcjx41AXxcnHn4OHhQUVzsEBMDUVGaDz6Ad97RzJ4d4Prru9N6\nNPv2mXKQyckwZowiN9eUnfz7jp08vOJdAP7zwc8xflRi8NynC7odxyI/39STN+UdZfRcCHFtk0Bd\nCCEGoL80jeJi3VU33KSB9ASgiozrogGoO3T8LJ+ryMjQbNumaG01CwOFhkJSkkap7omoivHjO2lv\nN5VhAHbvDiEtzYxA19drNm7UxMe7pKVZHDlisWcPREa6hIZqTpwwk0xnz4Zt22wqKhQ+n8vBg5ro\naJeWFovkZE1ODlhWGDfd5GfjRrdr8SRFZmYAjwdqaswCRmbE3CIvT9HRYb5vQYHC4zHXKCbpAD/5\n+Z/QWnPThOtZUJA94Ovc32TU7teycqkQ4lojgboQQpyH0aM1Xm9Peka3YUMisSzFgaOtdPoDhDg2\nc+bAzJk9+d1+v9s1iqyoq7NwXU1aGtTVQVUVJCXBV7+qCQRcfvUrRVhYgHHjLHw+hzFj/AwZoqiv\n12jtcvCgzaFDFo884mffvhBKSjRRUZqtW22OH4edO122bjVlHmfPhiFDNEeOaPbvB601tbWKiRPN\niHZZmcO+fS5ZWZrJk10++shh6FC4/XYIC7OD5RcdRxEa2jNp1udzefXNQzzw3B842upj/oyxvPBP\ns/tMLjXHnTnoPl0gLgG6EOJaI4G6EEJ8QpZlqq04Tk/d8J68aoeRiTHsbWxmR+0Rxo9KxHEsPB6T\nWtKd+pGb66eiwiIz02XKFCgvh5YWTXu7KcPYPaGysdFUa4mJCfD3vw9h2LAAN94YoKnJlHAcNcpM\nGt2zJ5S8PKsrb96lutqkwNg2HDli0lgmTlRUV9uEhweoqbFobYXkZJf9+/vWP9+1S7N/v6KtzfzV\nIDLS6vMd/X5TU92yAmzfrqhvOsz/XfM8x3yt5GaOYn7a5ykrUxQWnppTLkG3EEKcnfxSCiHEJ9ZT\nyaU77aX3Yj952SkAbNhae1Ibk7Liui6lpS47d7pYFpSWKnbtssjIMDnfSpmf6LAwi+xsTVqaZuzY\nALatqa+3UcplyxbF++87LF0KCxZYWJaDx2ORn2/h84UQH6/wel1mzlRMnqxIT1eEh9ssWmRx992m\nCkx4uCm12N6uuhYuUtx5J3g8ZqLp2LHdVVjo8x0dRzF6dAClNJt2Hub/rVnBMV8r0zJG8diXbifU\nCbk4/xuEEOIqJSPqQghxDnpXK+mdvnHyAj4A3nEpvPr+Fkq21HL3/KlA92qgppJLIKBZsQICAVNX\nvbra5KsXFNhdKSU6OFofEwM1NbBtWwhZWcfo7NRMmRLHBx+AbUNoqGLaNIVtm9VGTXlEyMlR1NSE\nUFWlGD3aT12dRVmZJjcXoqMdZs3ys2uXmRCakaGCVVZs26W+XjFkiIvfb1NWZuH1mmoy5joo/H5N\nQ4OmsvoIL1a9QLOvhTk5o/jzY3cS4QnB7++p2iKEEOLcSaAuhBAMrEb3mUoEOo6F1xvoem1ysmeO\nTwVg/eaaYDvLMhMyS0sBNA0N4Lqm7OKRIzahoSaAz8tTlJYqNmzQTJrk0tZmJptqDenpbXR0uAwZ\nksiDDwZYvz7Ak08qoqJcoqMVY8dq8vJg717Tr/R02LlTU10NEREua9YotmzRZGbagMWoUS7bt1t0\ndPRNfdFac+KEqTjjui5+P10VWaC42JRrXLOhiZe2vUhLZws3TEnnL/+6kMjw0K7r8An/ZwghhAAk\n9UUIIYIB+Pr19DsyPtDP2LBBsWGDCn5G3nXJxEaGUbnnIFX7DnW1VHR2unR0aECRmqq7qqVokpJc\n9u93ee01PyUlLpalefddzauvalJS4OabFXfdZXH8eID2djOyHRZmUVNj09KiUSqA1i47dmgqKqCj\nw6xkmpurcRzN5s02Pp9Le7up275rl5/qaggJUShlctJbWwPB79TRofD5ICTElEwsKTHHddvf1MT/\nVK+gpfM4BRNH8tef3BUM0v1+9xNfSyGEEIYE6kIIcQ68Xk1BwcAmQ4aFOnyhcBwAL7yzCcexyMkJ\nUFNjUkby8hR33GHh9bqsXRtCe7smOjrAunWKqirNpEma4cM1TU1w6JAZzbZtM+qttSmxCGai55Qp\nioQERWSkJjLSD8DixYpFi2DLFpvwcJg4UZOW5pCX5xIebnLMs7M1YDN6NBw9Cq++aurDO44iPd1i\n9OiePPwdO0w5xtxcCE9o5PsvPM/RtuPkjx/J2//25T5B+ukefCSAF0KIgZM/TAohrnkDqdHdu0Z6\nQcHAP2Ph/5rIf79dwYq/beZ/f6mQigr48ENFXJxJcamosAkJMXXTwSYuThEZafLOIyJs7rwzwCuv\nwK5dpspKSAhYlsOhQ+G8/LKfjAyHvDzo7LRoaICtW12amhStrQFyc81PfFUVgMWdd5pAf8UKsz09\n3Uw6dRyFz2dRU2Mqy4CpZ/7FL7qUlJg67rm54Lpm3ysfbOZbv/gr7X4/E1JT+etPFhIVHtpvAO73\nm5KQsrqoEEKcOwnUhRCCMwfo5/MZN07NIDUxhp31R3jid1XMmZDNhAmmpGL36DjYLF7sdpU8tEhL\nc7Esk+e+aZON67rMmqUJDzcpMnV14WzZEkV8vCYlJYDjOMyZY3HsmGbnTjPqHRGhWL9eY9uQlQWg\nCAszE0C1thk71mXmTIKTVaOiHBYtMmkvHo9ZYbWiwmLPHk1WFng8FjO8fh78zSqWv1EGwGdyJrPi\nkfkMifF0BeEmLaagwDzM+P06WCGm+yFGCCHEwEmgLoQQp9F3BJjgSPq5jAQ7tsVDX5rFfb96m+fe\nX8tDX80mN9dMFO2upNJdd71bXp4fU/HFwnVdsrMhL8+cs6QEbNtl1KhOhgwJxbZNuUSvVxMaajFz\nphmpz89XvPKKCZIXLTIBuQmaFV/4Qifl5VBerrAsF/NXAhePx8bnCwRTX0CRnQ25uZqaA83c+dgf\n+XtVHaEhNt+dN4/PTs0hOtLuulaaqioTqHu9PQ8Afa6FrC4qhBDnRAJ1IcRVZSDVWz6pgX7myX1Y\neksuT7xURPWBA7xeUsWts8YBuqtN3/KFPl+A5583rxctcrEshetqysoU4AIuSUknSEjoJCVlOCEh\nNqBZt85l1y5T+zwlRREaapGd3V3v3ATSrqsJBDR/+5tm40ZFRoZLeLgZxfd6wefz88c/mr7feqvC\n6zX9emrFHh77459obmsjfXgsLz78RdoPJZ90bRRZWT0TTXsWdHLxeKzgd5QAXQghBk4CdSHEVWOw\nc6DPNgLc30NB3xQQ04ewUIeHF89m6b//hWW/fIuwE2kMjYogN9ftM5LemwmqASwCAbCsAO++a7ZP\nnWrR1NRJWFh3RRfFunXQ0KBxHGhrA9s2JRodR1FcrNm+HbKzYcoUza9/rQgENBkZioiInpHv9es1\nH36oCQnR+P0WmZmaNbvX8y9/eB9Xaz49fTQrfvgF4mMjTvnujmNRWNj7G5gRdtflpO1CCCEGSgJ1\nIcQ15VxH3E9u1ztn/eSHAr/fxedzuyZvdqeAmNff+Ow0Xl5TyTvlu/m3P7/BV6bdgesqCgv7PlB4\nPDaLFvkpK1NUVJgqMaWlms5Ol0AAWlqgtjaUoUPbsSwoK1N4vXD99WaRokDATEStrjZVYrze3r03\nI+0jRgQYORJuuMHuKg1phIYqxoxROA7UNR3i179fTfnuPSjgnxYV8ujdc7Dt04+M996Wm+vHdXUw\n114IIcS5k0BdCHHVGMgI+PmMuPc+3us1ddBP3ue6dKWAqD552pal+N33b2XS156mZMdOspPXMW7c\nDcFje/e5u646uJSWatauVaSmKkaPDhAertm2zaK5OYyRIzVVVeaceXmQl6e68sxdtm83+7xeKCy0\nyM/XXdVdNLW1KriwUu9rUFgIY8ef4F/+8D6/fbMc19XEhIfzg88v4MGvjg0G6QO5TmVlFuB25eBL\nsC6EEJ+EBOpCiKvKxcqBdhzVZ3Jpd7BtWVZXSsupuedx4RE8+Nlb+PHLr/K79z9g4ngPBXj7fXhw\nXY0J1hVpaZCerikttbBtSE3twLI6yM1VuC50dGheeknT0aFYvNjkhGdlucF+mv4FgnXXk5OhpkZR\nWgpz5phtHZ0Blv9lI48++z7Nre1YSnFr3nSefHA2SfGRn+i6WpZ13quTXsg5B0IIcbmTQF0Icc04\n36ojZzq+e58pSWj2dQfePl+AZ581kzkLx2fxj7fczL+9/gb/+PQqtNZMH+5Fqb5VUizLVHzJyzO5\n5uXlFoGACd4jIwNkZLTi8ZiR8rVrA+zZo4mI6Kkkk5+vgkG6Ob9LXZ3m+uvhtttMxRfbVpw4EeDJ\nFTt5etXfqD1yBIC5uZnc5Z1LxrBhxEWa73QuAfdgVXeRuutCiGudBOpCiGvKYEwwPfO+/uuFmxFy\ns6qn1zuFpNQOHvj1Sh76zWoW/UMDv7r3JhzHE/yc3Fw/ZWWaF19UjB6tCQmxKSjopLJSsXdvKG1t\nUfh8kJPjsm6dRXR0gOhozc6d0NnpdrXve/6WFpfNmzVhYTbTp2tW/KWef/z9Gsp37wHgupHx/Pu3\n5zHfO4ZAQOPzuV0VaDRf+UoAj2fgKSwSVAshxPmTQF0IIc5T7/SM/kaTHUcxfnzPiHlxsWbasDx+\n+IVwnnj9Lzz/t02UbKnliXvmcuv12QQCmpISza5d4DhmYaTcXE1pafcqpp34/WbE3bYVaWkm1zwp\nqTv3HFzXBUxKjuMo7rpL8fLLmq07OnjvhUo2P1XBnoMHAYiL8vDDL8/m2wumE+LYKNU9Gt+dfnNp\nSN11IcS1btAD9d/+9re88MILfPjhhzQ3N7Nnzx7S0tL6tGlqamLZsmW88cYbAHzuc5/jySefJDY2\ndrC7I4QQpzUY+c/9pWf093ndueSlpZrt200O+acmTiBj2DB+8fZrbN5zgC/86GVmTRjJw1+eze6d\n6aDh85/XRERYXRNTA2Rluezb53LoUBiZmS6OY7NokepKTzEPAxs2mHP6fC5lZQp/wKUtbDevbNrM\n+1uq6Ah0AhAbEc68yTl8eXY+MRERlBRrLKvne3g8Nl/5Ss9qpRfi+p2NBOhCiGvZoAfqbW1tzJs3\nj1tvvZX777+/3zYLFy6ktraWVatMfubXv/51lixZwp///OfB7o4QQvTrYuc/W5apc97WpsjO1uTn\nm/NZG4bxiy9/g8qj5Tz+/AcUVdYw//88z/DYWOZOnoB6cxxZKUnMyDO10LtXDPX7Lfbt6ynBaKqs\nmGo0rgv7Dh3k5SerWbd5DzsP7aWtoyPYl5xR6Xxu+jTuvSuLjytstm9XxGT3P3J+unQXyR8XQogL\nb9AD9fvuuw+AsrKyfvdXVVWxcuVKioqKmDFjBgDLly+nsLCQHTt2MHbs2MHukhBCXDADT89Q+HxW\nMEj3eOzgiHSIY/OdW/O4e34Ov3i1mGf++iE1B5v5w9oi/rC2CE9ICNOykgjrTCImLJLsMc0c90ez\nd3s8W492srGhg81VPmoPH+HQHw6xveYwbe2dfc4+Pj2B22+YwKK5k8hIGgKY1Uoty5ST9Hp7ykle\nbkG3VH4RQlyrLnqOenFxMdHR0cycOTO4raCggMjISIqKiiRQF0JcFIOZ/zyQ471euoJhq1fuet8+\nxDhh/NOi2cxOK+SjvXtY/fEmPty9l/1NTRRV1gA1APx509n7lBAdzdy8dP4hN5NPTc1g7/YYADKS\nespJbtigcF2X/Hx1ThNF++v7hSIj90KIa9lFD9QbGhpITEw8ZfuwYcNoaGg447EbN268UN0S1xC5\nj8TFpLVFZaUJkidMOIZS/VeF6d1+794Y4m3Fj76UiWWl0dzmo3JfM0UfddLS3k5IRAvNbR1YShER\nZhMe6hAZ5jAyIYL0YdH4Do8gMiys63wBDtZWs3t3DForwsOPY1l+tLbYt8/MC9q0qfmM/dLaBMcD\n6ftA2p0LrS127zbXLyLi7NdPXLnkt1lcLbrv5cEYfB5QoP7www/z+OOPn3a/Uor33nuP2bNnn3eH\nhBDiWtM7wFXKZfz4Fioro9m6NYIJE44RFxnKrOzhxAb6BvxaW2htsWVLVN/t6e1AezCo7f7Mbdui\n2bIligkTjgHg9w+sbwN50DjXB5KBUsoN9leCdCHEtWZAgfr999/P4sWLz9jm5Moup5OUlMTBrpJg\nvR04cICkpKQzHjt9+vQBnUOI/nQ/4cp9JC6EM+VRT516+n09qR0ar5dg7vqJE93H9hzX/TkVFeVo\nbdHWlgNo0tM1lmX1advfeXp/JkBb26nnOPkYv99MgD1Tu+62Z/s8Ifojv83ianHyvdzc3Hzenzmg\nQH3o0KEMHTr0vE8GkJ+fT0tLCyUlJcE89aKiItra2pg1a9agnEMIIS6ms+VRd78/fTCvqaoy1VoK\nC93T5n/3fm1G4TWgyM8n+Pp0HMfC6w10vbaDfe2/P72/k8Lr1cFVTs/0+VLzXAghBteg56g3NjbS\n0NDAtm3b0FpTWVnJ0aNHSUtLY8iQIWRnZzNv3jyWLl3K8uXL0Vpzzz33cMstt8hEUiHEVet0wXz3\nKqSuq7GsngmdZwp2u9Nd0tM1+fmmZOPZJlx2Tx7t3WagAfXZgvSB9FkIIcS5G/Rf1d/85jdMnTqV\nxYsXo5Ti5ptvZtq0acHFjQBeeOEFpkyZwmc+8xnmz5/P1KlTee655wa7K0IIcdF4vZqCgnMPVv1+\nl9JS6OiA3NwzVzUxqSg9edqW1b0Q0uAzI+R8ou8khBBicCit9aVbH3oAeuf3yMql4nxIHqS4EPqO\nlJ85qO0v9cXnC/Dssy719TBnDhQW2mdJRYGIiAoApk7NGUBazann7yYBuLgcyG+zuFqcKUf9k8aw\nF708oxBCXKv6C4wdR5GdrXCcM+eYn0ypvqPvAwnYu+unS11yIYS4MkigLoQQ5+F8J1F2p5j0N2Gz\nd9DdezLo5s39lymUIFwIIa4uEqgLIcR56j2a7fefvULKyaPeJhA/tU3voBsITgaNiLA+cU1xqc4i\nhBBXDgnUhRDiPHUH6Bs2QFWVJitLU1h4trKHgzvq3R38FxSY92crpSiEEOLyJ4G6EEKch96Bt+u6\nnEue+Zl0j3z7/SZ3vfdIeEVF39H0c5nQ2vuY7vP09/5cjxdCCDH4JFAXQohBocjPt4J1zU8XwJ5r\n6snJtc97O7mKy0D1l1ZzLqP8kgsvhBAXhwTqQghxHvoG3vZZWvccc760tvqMop+c8iIj3kIIceWT\nQF0IIc7ThQqGzzb6blJtAPquMnq2Ee/+PvdcRvllQqoQQlwcEqgLIcRl7EyBsGV98nz4/oL38zle\nCCHE4JNAXQghrlj9B+oy4i2EEFcHCdSFEOIKo3V3uop5f7pFkoQQQlzZ5JdcCCGuIH6/S2VlDJWV\nMcCpQfr69aaCyyetCCOEEOLyIYG6EEIIIYQQlyFJfRFCiCuI41hMmHAs+PrkfZKbLoQQVw8J1IUQ\n4gqj1OnTWiRAF0KIq4f8ogshxEXm97uSQy6EEOKsJFAXQoiL6GJP+JSHAiGEuHJJoC6EEFcpqQIj\nhBBXNslRF0KIi0gmfAohhBgoCdSFEOIiu1gBujwUCCHElU0CdSGEuIpJgC6EEFcu+QUXQgghhBDi\nMiSBuhBCCCGEEJchCdSFEEIIIYS4DA1qoH706FHuvfdexo0bR0REBGlpaXz729/myJEjfdo1NTWx\nePFi4uLiiIuLY8mSJTQ3Nw9mV4QQQgghhLiiDWqgXl9fT319PT//+c/ZvHkzzz//PB988AF33XVX\nn3YLFy6koqKCVatWsXLlSsrLy1myZMlgdkUIIYQQQogr2qBWfZkwYQKvvvpq8H1mZiZPPPEEt9xy\nCy0tLURFRVFVVcXKlSspKipixowZACxfvpzCwkJ27NjB2LFjB7NLQgghhBBCXJEueI56c3MzYWFh\nREREAFBcXEx0dDQzZ84MtikoKCAyMpKioqIL3R0hhBBCCCGuCBc0UG9qauJHP/oR3/zmN7Esc6qG\nhgYSExNPaTts2DAaGhouZHeEEEIIIYS4Ygwo9eXhhx/m8ccfP+1+pRTvvfces2fPDm5rbW3llltu\nYeTIkfzsZz87/56CTDgV56U7rUruI3Glk3tZXE3kfhZXiwtxLw8oUL///vtZvHjxGdukpaUFX7e2\ntjJ//nxs2+aNN94gNDQ0uC8pKYmDBw+ecvyBAwdISkoaaL+FEEIIIYS4qg0oUB86dChDhw4d0Ae2\ntLQwf/58lFK8+eabwdz0bvn5+bS0tFBSUhLMUy8qKqKtrY1Zs2adY/eFEEIIIYS4OimttR6sD2tp\naWHu3Lm0tLTw2muvERUVFdw3dOhQQkJCALjpppuoq6tj+fLlaK1ZunQpmZmZvPbaa4PVFSGEEEII\nIa5ogxqov//++3zqU5/qs01rfUoOe3NzM8uWLeP1118HYMGCBTz55JPExMQMVleEEEIIIYS4og1q\noC6EEEIIIYQYHBe8jroQl4M5c+ZgWVbwn23bp6yY29TUxOLFi4mLiyMuLo4lS5ZIFQJxWfr1r39N\nZmYm4eHhTJ8+nXXr1l3qLglxRo8++mif32DLskhOTu7T5pFHHiElJYWIiAhuvPFGtmzZcol6K0SP\ntWvXsmDBAlJTU7Esi+eee+6UNme7dzs6Oli2bBmJiYlERUWxYMEC6urqBnR+CdTFNUEpxde+9jUa\nGxtpaGhg//79LF++vE+bhQsXUlFRwapVq1i5ciXl5eUsWbLkEvVYiP699NJLfO973+OHP/whFRUV\nzJo1i/nz51NbW3upuybEGWVnZwd/gxsaGti0aVNw389+9jN+8Ytf8NRTT7Fx40aGDRvG3LlzaW1t\nvYQ9FsLMv5w0aRK//OUvTymQAgO7d++77z7+9Kc/8dJLL7Fu3TqOHTvGzTffzICSWrQQ14A5c+bo\nZcuWnXb/1q1btVJKFxcXB7etW7dOK6X09u3bL0YXhRgQr9erly5d2mfb2LFj9Q9+8INL1CMhzu6R\nRx7RkyZNOu3+ESNG6J/85CfB9ydOnNDR0dH6P/7jPy5G94QYkKioKP3ss8/22Xa2e7e5uVmHhobq\nF154IdimpqZGW5alV61addZzyoi6uGa8+OKLJCYmMnHiRB566CFaWlqC+4qLi4mOjg6WDAUopugV\nRwAABExJREFUKCggMjKSoqKiS9FdIU7R2dlJWVkZc+fO7bP905/+tNyn4rJXXV1NSkoKmZmZLFy4\nkN27dwOwe/duGhoa+tzXHo+H2bNny30tLmsDuXc3btyI3+/v0yY1NZVx48YN6P4eUB11Ia50ixYt\nIj09neTkZCorK/n+97/Ppk2bePvttwFoaGggMTHxlOOGDRtGQ0PDxe6uEP06dOgQgUCA4cOH99k+\nfPhw3nnnnUvUKyHObubMmfzud78jOzubAwcO8Nhjj1FQUEBlZSUNDQ0opfq9r+vr6y9Rj4U4u4Hc\nu42Njdi2TXx8/CltBhJfSKAurlgPP/wwjz/++Gn39y4L+o1vfCO4fcKECWRmZjJjxgwqKirIycm5\nGN0VQohr1rx58/q8nzlzJhkZGTz77LN4vd5L1CshLn8SqIsr1v3338/ixYvP2CYtLa3f7bm5udi2\nzY4dO8jJySEpKYmDBw+e0u7AgQMkJSUNSn+FOF8JCQnYtk1jY2Of7Y2NjXKfiitKREQEEyZMYMeO\nHSxYsACtNY2NjaSmpgbbyH0tLndJSUlnvXeTkpIIBAIcPny4z6h6Y2NjcH2hM5EcdXHFGjp0KFlZ\nWWf85/F4+j32448/JhAIMGLECADy8/NpaWmhpKQk2KaoqIi2tjZmzZp1Ub6PEGcTEhJCbm4uq1ev\n7rN99erVFBQUXKJeCXHufD4fVVVVJCcnk5GRQVJSUp/72ufzsXbtWrmvxWVtIPdubm4ujuP0aVNb\nW8vWrVsHdH/bjzzyyCOD3nMhLiPV1dX86le/Iioqis7OToqKili6dCnp6ek89thjKKVISEhgw4YN\nrFixgqlTp1JTU8M999zDzJkz+c53vnOpv4IQQTExMfz4xz9mxIgRRERE8Nhjj7F27Vr+67/+i9jY\n2EvdPSH69dBDD+HxeNBas23bNr773e+ya9culi9fTmxsLIFAgJ/+9Kdcd911BAIBHnjgARobG1m+\nfDmhoaGXuvviGtba2srWrVtpaGjgmWeeYfLkycTGxtLZ2TmgezcsLIz9+/fz1FNPMXnyZJqbm/nW\nt77FkCFD+OlPf4pS6swdGIxyNUJczmpqavQNN9ygExIStMfj0WPHjtX333+/Pnr0aJ92TU1NevHi\nxTo2NlbHxsbqJUuW6Obm5kvUayFO7+mnn9YZGRna4/Ho6dOn63Xr1l3qLglxRnfeeadOSUnRYWFh\nOjU1VX/xi1/UW7du7dPm0Ucf1cnJyTo8PFzPmTNHV1ZWXqLeCtFjzZo1WimlLcvq8+/uu+8Otjnb\nvdvR0aHvvfdenZCQoCMjI/WCBQt0bW3tgM6vtB5ItXUhhBBCCCHExSQ56kIIIYQQQlyGJFAXQggh\nhBDiMiSBuhBCCCGEEJchCdSFEEIIIYS4DEmgLoQQQgghxGVIAnUhhBBCCCEuQxKoCyGEEEIIcRmS\nQF0IIYQQQojLkATqQgghhBBCXIb+P5FCo25Z78OAAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x73e9cf8>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "from numpy.random import multivariate_normal\n",
    "from filterpy.stats import (covariance_ellipse, \n",
    "                            plot_covariance_ellipse)\n",
    "\n",
    "mean = (5, 3)\n",
    "P = np.array([[32, 15],\n",
    "              [15., 40.]])\n",
    "\n",
    "x,y = multivariate_normal(mean=mean, cov=P, size=2500).T\n",
    "plt.scatter(x, y, alpha=0.3, marker='.')\n",
    "plt.axis('equal')\n",
    "\n",
    "plot_covariance_ellipse(mean=mean, cov=P,\n",
    "                        variance=2.**2,\n",
    "                        facecolor='none')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## The Algorithm"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "As I already stated, when the filter is initialized a large number of sigma points are drawn from the initial state ($\\mathbf{x}$) and covariance ($\\mathbf{P}$). From there the algorithm proceeds very similarly to the UKF. During the prediction step the sigma points are passed through the state transition function, and then perturbed by adding a bit of noise to account for the process noise. During the update step the sigma points are translated into measurement space by passing them through the measurement function, they are perturbed by a small amount to account for the measurement noise. The Kalman gain is computed from the \n",
    "\n",
    "We already mentioned the main difference between the UKF and EnKF - the UKF choses the sigma points deterministically. There is another difference, implied by the algorithm above. With the UKF we generate new sigma points during each predict step, and after passing the points through the nonlinear function we reconstitute them into a mean and covariance by using the *unscented transform*. The EnKF keeps propagating the originally created sigma points; we only need to compute a mean and covariance as outputs for the filter! \n",
    "\n",
    "Let's look at the equations for the filter. As usual, I will leave out the typical subscripts and superscripts; I am expressing an algorithm, not mathematical functions. Here $N$ is the number of sigma points, $\\chi$ is the set of sigma points."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Initialize Step"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "$$\\boldsymbol\\chi \\sim \\mathcal{N}(\\mathbf{x}_0, \\mathbf{P}_0)$$\n",
    "\n",
    "This says to select the sigma points from the filter's initial mean and covariance. In code this might look like\n",
    "\n",
    "```python\n",
    "N = 1000\n",
    "sigmas = multivariate_normal(mean=x, cov=P, size=N)\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Predict Step"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "$$\n",
    "\\begin{aligned}\n",
    "\\boldsymbol\\chi &= f(\\boldsymbol\\chi, \\mathbf{u}) + v_Q \\\\\n",
    "\\mathbf{x} &= \\frac{1}{N} \\sum_1^N \\boldsymbol\\chi\n",
    "\\end{aligned}\n",
    "$$\n",
    "\n",
    "That is short and sweet, but perhaps not entirely clear. The first line passes all of the sigma points through a use supplied state transition function and then adds some noise distributed according to the $\\mathbf{Q}$ matrix. In Python we might write\n",
    "\n",
    "```python\n",
    "for i, s in enumerate(sigmas):\n",
    "    sigmas[i] = fx(x=s, dt=0.1, u=0.)\n",
    "\n",
    "sigmas += multivariate_normal(x, Q, N)\n",
    "```\n",
    "\n",
    "The second line computes the mean from the sigmas. In Python we will take advantage of `numpy.mean` to do this very concisely and quickly.\n",
    "\n",
    "```python\n",
    "x = np.mean(sigmas, axis=0)\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can now optionally compute the covariance of the mean. The algorithm does not need to compute this value, but it is often useful for analysis. The equation is\n",
    "\n",
    "$$\\mathbf{P} = \\frac{1}{N-1}\\sum_1^N[\\boldsymbol\\chi-\\mathbf{x}^-][\\boldsymbol\\chi-\\mathbf{x}^-]^\\mathsf{T}$$\n",
    "\n",
    "$\\boldsymbol\\chi-\\mathbf{x}^-$ is a one dimensional vector, so we will use `numpy.outer` to compute the $[\\boldsymbol\\chi-\\mathbf{x}^-][\\boldsymbol\\chi-\\mathbf{x}^-]^\\mathsf{T}$ term. In Python we might write\n",
    "\n",
    "```python\n",
    "    P = 0\n",
    "    for s in sigmas:\n",
    "        P += outer(s-x, s-x)\n",
    "    P = P / (N-1)\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Update Step"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In the update step we pass the sigma points through the measurement function, compute the mean and covariance of the sigma points, compute the Kalman gain from the covariance, and then update the Kalman state by scaling the residual by the Kalman gain. The equations are\n",
    "\n",
    "$$\n",
    "\\begin{aligned}\n",
    "\\boldsymbol\\chi_h &= h(\\boldsymbol\\chi, u)\\\\\n",
    "\\mathbf{z}_{mean} &= \\frac{1}{N}\\sum_1^N \\boldsymbol\\chi_h \\\\ \\\\\n",
    "\\mathbf{P}_{zz} &= \\frac{1}{N-1}\\sum_1^N [\\boldsymbol\\chi_h - \\mathbf{z}_{mean}][\\boldsymbol\\chi_h - \\mathbf{z}_{mean}]^\\mathsf{T} + \\mathbf{R} \\\\\n",
    "\\mathbf{P}_{xz} &= \\frac{1}{N-1}\\sum_1^N [\\boldsymbol\\chi - \\mathbf{x}^-][\\boldsymbol\\chi_h - \\mathbf{z}_{mean}]^\\mathsf{T} \\\\\n",
    "\\\\\n",
    "\\mathbf{K} &= \\mathbf{P}_{xz} \\mathbf{P}_{zz}^{-1}\\\\ \n",
    "\\boldsymbol\\chi & = \\boldsymbol\\chi + \\mathbf{K}[\\mathbf{z} -\\boldsymbol\\chi_h + \\mathbf{v}_R] \\\\ \\\\\n",
    "\\mathbf{x} &= \\frac{1}{N} \\sum_1^N \\boldsymbol\\chi \\\\\n",
    "\\mathbf{P} &= \\mathbf{P} - \\mathbf{KP}_{zz}\\mathbf{K}^\\mathsf{T}\n",
    "\\end{aligned}\n",
    "$$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This is very similar to the linear KF and the UKF. Let's just go line by line.\n",
    "\n",
    "The first line,\n",
    "\n",
    "$$\\boldsymbol\\chi_h = h(\\boldsymbol\\chi, u),$$\n",
    "\n",
    "just passes the sigma points through the measurement function $h$. We name the resulting points $\\chi_h$ to distinguish them from the sigma points. In Python we could write this as\n",
    "\n",
    "```python\n",
    "sigmas_h = h(sigmas, u)\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The next line computes the mean of the measurement sigmas.\n",
    "\n",
    "$$\\mathbf{z}_{mean} = \\frac{1}{N}\\sum_1^N \\boldsymbol\\chi_h$$\n",
    "\n",
    "In Python we write\n",
    "\n",
    "```python\n",
    "z_mean = np.mean(sigmas_h, axis=0)\n",
    "```\n",
    "    \n",
    "Now that we have the mean of the measurement sigmas we can compute the covariance for every measurement sigma point, and the *cross variance* for the measurement sigma points vs the sigma points. That is expressed by these two equations\n",
    "\n",
    "$$\n",
    "\\begin{aligned}\n",
    "\\mathbf{P}_{zz} &= \\frac{1}{N-1}\\sum_1^N [\\boldsymbol\\chi_h - \\mathbf{z}_{mean}][\\boldsymbol\\chi_h - \\mathbf{z}_{mean}]^\\mathsf{T} + \\mathbf{R} \\\\\n",
    "\\mathbf{P}_{xz} &= \\frac{1}{N-1}\\sum_1^N [\\boldsymbol\\chi - \\mathbf{x}^-][\\boldsymbol\\chi_h - \\mathbf{z}_{mean}]^\\mathsf{T}\n",
    "\\end{aligned}$$\n",
    "\n",
    "We can express this in Python with\n",
    "\n",
    "```python\n",
    "P_zz = 0\n",
    "for sigma in sigmas_h:\n",
    "    s = sigma - z_mean\n",
    "    P_zz += outer(s, s)\n",
    "P_zz = P_zz / (N-1) + R\n",
    "\n",
    "P_xz = 0\n",
    "for i in range(N):\n",
    "    P_xz += outer(self.sigmas[i] - self.x, sigmas_h[i] - z_mean)\n",
    "P_xz /= N-1\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Computation of the Kalman gain is straightforward $\\mathbf{K} = \\mathbf{P}_{xz} \\mathbf{P}_{zz}^{-1}$.\n",
    "\n",
    "In Python this is\n",
    "\n",
    "```python\n",
    "K = np.dot(P_xz, inv(P_zz))```\n",
    "\n",
    "Next, we update the sigma points with\n",
    "\n",
    "$$\\boldsymbol\\chi  = \\boldsymbol\\chi + \\mathbf{K}[\\mathbf{z} -\\boldsymbol\\chi_h + \\mathbf{v}_R]$$ \n",
    "\n",
    "Here $\\mathbf{v}_R$ is the perturbation that we add to the sigmas. In Python we can implement this with\n",
    "\n",
    "```python\n",
    "v_r = multivariate_normal([0]*dim_z, R, N)\n",
    "for i in range(N):\n",
    "    sigmas[i] += dot(K, z + v_r[i] - sigmas_h[i])\n",
    "```\n",
    "\n",
    "\n",
    "Our final step is recompute the filter's mean and covariance.\n",
    "\n",
    "```python\n",
    "    x = np.mean(sigmas, axis=0)\n",
    "    P = self.P - dot3(K, P_zz, K.T)\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Implementation and Example"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "I have implemented an EnKF in the `FilterPy` library. It is in many ways a toy. Filtering with a large number of sigma points gives us very slow performance. Furthermore, there are many minor variations on the algorithm in the literature. I wrote this mostly because I was interested in learning a bit about the filter. I have not used it for a real world problem, and I can give no advice on using the filter for the large problems for which it is suited. Therefore I will refine my comments to implementing a very simple filter. I will use it to track an object in one dimension, and compare the output to a linear Kalman filter. This is a filter we have designed many times already in this book, so I will design it with little comment. Our state vector will be\n",
    "\n",
    "$$\\mathbf{x} = \\begin{bmatrix}x\\\\ \\dot{x}\\end{bmatrix}$$\n",
    "\n",
    "The state transition function is\n",
    "\n",
    "$$\\mathbf{F} = \\begin{bmatrix}1&1\\\\0&1\\end{bmatrix}$$\n",
    "\n",
    "and the measurement function is\n",
    "\n",
    "$$\\mathbf{H} = \\begin{bmatrix}1&0\\end{bmatrix}$$\n",
    "\n",
    "The EnKF is designed for nonlinear problems, so instead of using matrices to implement the state transition and measurement functions you will need to supply Python functions. For this problem they can be written as:\n",
    "\n",
    "```python\n",
    "def hx(x):\n",
    "    return np.array([x[0]])\n",
    "\n",
    "def fx(x, dt):\n",
    "    return np.dot(F, x)\n",
    "```\n",
    "\n",
    "One final thing: the EnKF code, like the UKF code, uses a single dimension for $\\mathbf{x}$, not a two dimensional column matrix as used by the linear kalman filter code.\n",
    "\n",
    "Without further ado, here is the code."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false,
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAuAAAADXCAYAAABbPuQ1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xd4VNX28PHvzGTSe29AAumhBwhFuoBIFVDwRfGiiIqK\nXPGHqEhRFBUUwYuiclVEERApEREBAQWRFkgCJIQSCKQnk2QmUzPlvH8MjHLpEAjI/jxPHsiZc/ZZ\nZ3Z4WDmzz1oySZIkBEEQBEEQBEG4JeT1HYAgCIIgCIIg3E1EAi4IgiAIgiAIt5BIwAVBEARBEATh\nFhIJuCAIgiAIgiDcQiIBFwRBEARBEIRbSCTggiAIgiAIgnALiQRcEARBEARBEG6ha07At2/fzqBB\ng4iMjEQul/P11187XrNYLLz88su0aNECT09PwsPDGTlyJGfOnDlvjNraWp5//nmCgoLw9PRk0KBB\nFBYW3vjVCIIgCIIgCMJt7poTcK1WS7NmzZg/fz7u7u7nvabX68nIyOD111/nwIEDpKWlcebMGfr2\n7YvNZnPs98ILL7B69WqWL1/Ojh070Gg09O/fH9ETSBAEQRAEQfink91IJ0wvLy8WLFjAqFGjLrlP\nTk4OycnJHDx4kOTkZDQaDUFBQSxevJgRI0YAUFBQQKNGjdiwYQO9evU673i1Wn294QmCIAiCIAhC\nvfPx8Tnv+5u+BlytViOTyfDz8wMgPT0di8VyXqIdGRlJYmIiO3fuvNnhCIIgCIIgCEK9uqkJuNls\nZuLEiQwcOJDw8HAASkpKUCgUBAQEnLdvSEgIJSUlNzMcQRAEQRAEQah3TjdrYKvVysiRI9FoNKxb\nt65Oxjx27FidjCMIgiAIgiAIN1NsbOwlX7spd8CtVisjRozg0KFDbNmyxbH8BCA0NBSr1YpKpTrv\nmNLSUkJDQ29GOIIgCIIgCIJw26jzO+AWi4Xhw4eTnZ3Nb7/9RlBQ0Hmvp6Sk4OTkxKZNm857CDMn\nJ4dOnTpdduw2bdrUdbhXZd++ffV6fqF+iHm/+4g5vzuJeb87iXm/O93Keb9cIZFrTsB1Oh3Hjx9H\nkiRsNhunT58mMzMTf39/wsPDGTZsGOnp6fz4449IkkRpaSlgf/rT1dUVb29vnnjiCSZNmkRQUBD+\n/v5MnDiRli1b0rNnz+u/SkEQBEEQBKHeWSwW0tPT0ev1tGzZ8ryVEILdNS9B2bdvH61atSIlJQWj\n0ci0adNo3bo106ZNo6CggLS0NIqKikhJSSE8PNzxtWLFCscY8+bN44EHHmDEiBF07twZb29v0tLS\nkMlkdXpxgiAIgiAIwq1TWFhImzZtaN++PT169CA6OppffvmlvsO67VzzHfCuXbue11Tnf13utXOU\nSiXz5s1j3rx513p6QRAEQRAE4Tb1yCOPkJmZSUhICEFBQRw6dIhhw4Zx9OhRwsLC6ju828ZNrwMu\nCIIgCIIg/PPp9Xq2bduGTCbj0KFDZGVl0bFjR7RaLTt27Kjv8G4rIgEXBEEQBEEQbphcLkcmkyFJ\nEhqNBrPZjE6nA0ChUNRzdLcXkYALgiAIgiAIN8zV1ZXBgwcDkJSURFhYmGM5Srdu3eo3uNvMTWvE\nIwiCIAiCINxdvvjiC4xGIz///DMmk4moqChWrlyJv79/fYd2WxEJuCAIgiAIglAnfH19Wb9+PcXF\nxej1eho1aoSTk0g3/5d4RwRBEARBEIQ6JSqeXJ5YAy4IgiAIgiAIt5BIwAVBEARBEAThFhIJuCAI\ngiAIgiDcQiIBFwRBEARBEG6J6upqcnJy0Gq11z3G37uuFxYWsn79eiyWWmpr1axfv5Jhw/pjMBTX\nRbg3jUjABUEQBEEQhJtKkiSmTp1KYGAgSUlJBAUFsXDhwgv2q6ioYPPmzUiShNVqJTv7IK+99jIm\nkwqDoZhff/2BLl3aU1NzGLU6g8LCHezevZrq6s1oNNtp1gxmz34Ui6WyHq7y6okqKIIgCIIgCMJN\ntWjRIt58803AXiFFp9Px9ttv4OEhZ8mS71i16lOsVgOnTx9n3brVtGolQ5JqUSiqaNHCg5qanQA0\nbapg1arXMJnyAGjc2J3nnx8IWAFwc3Opl+u7Vtd8B3z79u0MGjSIyMhI5HI5X3/99QX7TJ8+nYiI\nCNzd3enevTvZ2dnnvV5bW8vzzz9PUFAQnp6eDBo0iMLCwuu/CkEQBEEQBOG2Y7PZyMrK4qef1tG8\neWM+/fQN2rdPZubMYcTG1pCb+z3vvfcoRmMuZvNpGjZ0ZsqU4UiSHrAQEOBFjx6tHOMpFArk8sun\nrxUV5WRmZpKTk4MkSTf5Cq/PNSfgWq2WZs2aMX/+fNzd3S94/d1332Xu3LksWLCAffv2ERwcTK9e\nvdDpdI59XnjhBVavXs3y5cvZsWMHGo2G/v3737ZvkiAIgiAIgnB1zGYzNpsNs1mHVpvPmDGP0KCB\nCi+vU/j5VfHZZ+NRKEwYDFrARmRkUJ2de+XK73niiSeYNm0aSUlJjBw5ErPZXGfj15VrTsD79u3L\nzJkzGTJkCDKZ7ILX582bxyuvvMLgwYNJSkpi8eLF1NTUsHTpUgA0Gg1ffPEFc+bMoUePHrRs2ZIl\nS5aQlZXF5s2bb/yKBEEQBEEQhFvOvm7bwr33duWPP5aiVm+jtvYQ69e/RfPmcdhsNj76aD7/93//\nx7fffgtAu3Zt6+z8+/bt5euvlyBJEBcXh4uLC9999x1vv/12nZ2jrtTpQ5gnT56kpKSEXr16Oba5\nurrSpUsXdu60r93Zt28fFovlvH0iIyNJTEx07CMIgiAIgiDcGTZv3szGjevQanOort7K4sUvkJjo\nA/xVrWTQoIH06nUvFouV3NxcJAlGjBhOx46d6iyOvXv3AfDgg8N47733HDd/169fX2fnqCt1+hBm\nSUkJMpmMkJCQ87aHhIRQVFQEQGlpKQqFgoCAgAv2KSkpqctwBEEQBEEQhDpgtVoda6+3bt3KiRM5\nPPLIA9hsevT6HECHyWRfSuzpeeESZblcwfPPj2fo0KGUl1cQHh5OUFDdLT0B+/pwAJ1OD4BarT5v\n++3kjqqCsm/fvrv6/EL9EPN+9xFzfncS8353upvm3Ww28+mnn/L7778D0L17d5588kmcnC6fDu7Z\n8wfr16fx5pvjsFprqKzMw2arJjNzGQDBwQDOHD9+7Kri8PBwR62uRq2uvpHLuUBUVCNkMlj/88/s\nzyhl9Q/pAHTs2PGCeb4V8x4bG3vJ1+o0AQ8NDUWSJEpLS4mMjHRsLy0tJTQ01LGP1WpFpVKddxe8\ntLSULl261GU4giAIgiAIwlkzZszgl19+cXx/8uRJysrKmDZt2nn7GY1GNm7cyJAhPbDZykhONhMf\n3xOTKR+AyEhfIiN9b2nsl2KzSejNJlRaDSW1zkS2nEp+Xgyl6gPI5QcYOXIkDz30UH2HeYE6TcCj\no6MJDQ1l06ZNpKSkAPZJ3L59O++//z4AKSkpODk5sWnTJkaMGAFAQUEBOTk5dOp0+XVAbdq0qctw\nr9q535Lq6/xC/RDzfvcRc353EvN+d6qreS8tLWXHjh0oFAq6du2Kn5/fZfevra3lu+++4/jx4zRq\n1IiRI0fi5uZ2QzFcjby8PH755RecnZ1Zs2YNVquVwYMHs27dOhYuXEh2djYdO3ZEqbSh15/hu+9y\nCAlphVLpBsTf9PiuhiRJaPQ6ilTlbPptP1mH5ajVyZQXplByKh5VSUOQ7Mtk2tzTFJVqKr6+5/+i\ncCv/vZ9bAnMx15yA63Q6jh8/jiRJ2Gw2Tp8+TWZmJv7+/jRo0IAJEyYwa9Ys4uPjiY2NZebMmXh5\nefHwww8D4O3tzRNPPMGkSZMICgrC39+fiRMn0rJlS3r27Hn9VykIgiAIgnALbd68maFDh6LRaAAI\nDg7mp59+umRyZzQa6d27N9u3b3ds+/jjj9m2bRve3t43NdbS0lIA4uPj6dOnDy++OJ7u3dug16so\nL09n6dKvCAmpIDzcG7AxY8aomxoPgMVqwWAyYag1Of6stVgwmqxUVSuoqlJSXa1EXe1M/mlXTub5\nUFYUT1VpT2qqxl8wnlxhJijiOKFRubTsUIyvb6uLnPX2cM0J+L59++jevbujBOG0adOYNm0ajz32\nGF988QWTJk3CaDTy3HPPUVVVRWpqKhs3bsTDw8Mxxrx581AqlYwYMQKDwcC9997LkiVLLlrWUBAE\nQRAE4XajVqt58MEH0Wg0pKamYjAYyMrKYtiwYRw/fvyi66o/+OADtm/fTlhYGI899hjLli3jwIED\nvPHGG8yZM6fOY7TZbMhkMqxWMx9//CGtWjUkOFjPnDn/D7X6JGr1HtzdnfH3NzB79uhzR9V5HABW\nm41jBUWk/WJl794gatSuGHRemAyhmAweGPWe1J770+hxxfHkCjO+QcX4h5whuEEuPr67adNaT5/u\nqYQHBBHoW3fVVW6Ga07Au3btis12+cmZOnUqU6dOveTrSqWSefPmMW/evGs9vSAIgiAIQr3LzMyk\nurraUUbZYrEQGRlJfn4+J0+evOgDeOnp9ocC33rrLUaPHk1qaioPPPAABw4cqJOYzjU0tFpr+fe/\nnycpqREjRnTDalUzcWIfMjP9+eyzhaxZcwIAZ2cZzz777EUbK16N4uIicnJycHZ2oVWrlnh4eF4Q\nz7EzKlasduLXXxtwPKsbRp3PlQeW2XBx0+HipsPVXYuLmxYPn0r8QgoICS/GyekQkuUgwwf3o1lc\nAmtWriVtbRoBbR4iPrLRdV3LrXZHVUERBEEQBEG4HZxbt61SqaisrESn01FTU3Pea//rXNm9tLQ0\nunXrxurVqwEIDAy8YN9zZf8utTqgoqICjUZDdHQ0NpuVjz6aQ2lpAS+/PBKLpYoXX+yJh4crVmsl\nAJGRQURG3k+TJo3Zu3cPICM1NZW4uLjruv6tW7fy0UfzsVisAAQE+PPGG2/QoEFD1FojH38JP65t\nQl52H6xmF8dxPoHFNEpMx8uvAld3PV5eJnx8avHxNePra8HX14KPtwUXZwVKhRPOTkqUTk54uroR\n6h+Av1ci7733I3/8kcmyilIyE5McVV0SEhKv61rqg0jABUEQBEEQrlGrVq1o2bIlGRkZhISEOFYH\n9OrVi4iIiIse89xzz/H111+zZs0a1qxZA4CPjw8vvTSBtWuXsXr1aubPnwpIrFy5gc2b/2DRonnI\nZEp+/XU7mZmH+Pe/x2KzGVi//keOHz/Bv/89AqtVx8CBUbi4xGGxVADg5XXxu9oJCQkkJCTc0LWX\nlZU6ku+UlNaoVCpOncrn/fmf4OwzlTXfdaKq7K/3IKThUaKb7aFZm0y6dXClaVRjAr19UDr5X9fy\n47Fjx3L6dD5nzhRQXGzvITN48CBHAZA7gUjABUEQBEEQrpGTkxPr1q1j5MiR/PbbbwAMGDCAr776\n6pJJZXx8PNu2bWPu3LcoKDhKZWUNCxdOJjq6iqAgOcnJg6mtzQNg4MB4BgyIQ6/PBKBRIz0+PsHo\n9fblKvffHwvEYrXaa2l7e1953XRdyc09isVipUWL5kydOo2Dxwp4ZbqetLWjqKkKtscTUELLbmkk\nt9lH26a+NIuOoWHQPY5mPjfCz8+PDz74gD179pz9FKAxSUlJNzzurSQScEEQBEEQhOsQERHBtm3b\n0Gq1KBSKSy49kSSJkpJ8OnXqwt693/LRR08C0t/2sOLpeeGxf0/kw8ICCAsLuGCfW8lssVBZo6HK\nbMQvuj0nqzrw4BgLe7c+jL7GXn7RL+QMbe79ga69crgnOYn4BvfjdBM6Ubq4uNK5853bP0Yk4IIg\nCIIgCDfAw8MDm83C//3fRF566QV8fb2QJAu9ew/iq6/m4u9vQ6msZOvW95Gkuu3+eCnZ2Yc5cSIP\nb29vUlPb4ep6bbXGLVYrJVUqjp+pJPOwlZwjrpw+GU55QRLlhf3Rqc//ZSAo8gRt7v2evveX0Ck5\nmQZB/W64up0kSfz662Z27tyJzSbRvn0qffrc94+omicScEEQBEEQhGskSRJPPvk4jz/+AImJ/lit\nNcTEyDEYdqFQ2JPd2bMfw82tCIvFnm55eLjekti++OK/rFmz1vF9REQEb7/9FpXVvqxK82TrthDU\nVe64e5hwdzfh7lmLu0ctHh61mC02Tp1yo7jQl+qKzhck2ucoXQwEhp8kKPIkDeP3EB+TzugHhhLi\nn3zN8ZpMRlasWMGxY8fx9PRkwID+JCYm8e2337BixfeO/fbv309hYSFPPDHm2t+U24xIwAVBEARB\nEICysjIkSSI42L6Oed++fbi6utK0aVMAtmzZglJpo3XraMzmEh59tDURESas1ioAhg49f0lEkybh\nt/YCsJc6XLNmLU5OTnTs1oOMw0HknO5Kp+6dKC+Muebx5Aoz3v5l+ASWEBB2Ct/AHJKSq0ht606Q\nrze+7h6EBybj5nJ9nSXNZjOvv/46R47kOrbt3LmTiRMn8v339uR77NgncXJyYuHChaxdm8bgwYMJ\nCLiwcgzIee21L3nqqSdITr74g7C3C5GAC4IgCIJwV5AkCYPBgJubGzKZjE2bfiY7ez/DhvVDkows\nWrQQLy9XHnmkL2AjM3MzXl7uRESoACgq+oOgIHeMRiMAyclR9Xcxf2OxWsgv0bAnXcYvm4JQW5ej\nqYhn7/vRGPV/ddh0UpqIjMsiKmkffsGF1JrcqDW6n/2y/12yyfEJKME7sAR//yIMNZlEhQbx8JAH\n2bFlC98s+QZbYEvua/tGncS+c+cfHDmSi7+/P2PHPkl6ejqbNm3m22+/QZIgIiKc/v0HALBhwwby\n8k5SVlbmSMCXL/+N+PimdOzYDSenAB55xJvGjZvh4uJ7udPWO5GAC4IgCILwj3T8+HHKy8tp27YV\nFksN8+d/RHl5Ka+9NhqbLZOoKB2NGlnR6/cDMHas/Q621WpPuAcPtpe1s9ns9b379r361ubZ2dms\nX78enU5LfHw8Q4cORal0vqb4LRaJ0nIoKJYoLpUoKZVRXOpEWZmSigoXqirdUFe7U13ph7o8FEm6\n8GFHb/9SGiXtIyopnUZxOTQM86ZBUAj+Xj5YrTasNiMWmw6rzeaoPR7mH0BkYFM+X7iDLXv30WP0\naMIDAmkQGQlAbW3tNV3H5ZSWlgFwzz330LFjJxITk9i0aTNVVdXI5XIKC4tYtmwZTk4K8vJOYjDI\nKC620KJFNEplILGxzkRENMTdPQqAzp0711lsN5NIwAVBEARBuGOd6/5of2BvA7/9tpVXX30Wq1VP\nTs5vnDhxnNjYKsDKv/6Vglwux2Ipx2rVERUVRExMZJ3HlJ6+jzffnOmoDZ6evp9Dhw4xY8YbKM5W\nBLHZoKjIhaPHncjOlXPipMvZtdfeVFd5oKvxwGhwB+nqyvbJ5Fb8Q/IJDM8nIPwUAWH5+PjlgvUY\nOpWKe9t1YWj/4ddUkaRZs+Zs2bKVb75ZQnZ2NocOHQSgadNm1/iO/KW4WEVGxgn69u0IgELhhU7n\nxp49+4mJ2c7WrTvR6dxo2jSWDh06smjRN3z++RrCw/1xcYmldesuyGRxeHklI5PJ6N079LpjqU8i\nARcEQRAE4Y5hs9mQyWTYbDZ27NjMrFnvsHz5HCwWFRERKvr2jUKvtyeKHTpE0qFDJGDv1lgXNaiv\nxn//+19sNhu9evWiefPmLFq0iKysg2zafIBj+Z35aV04WelNMNe6XH4gmQ1Xdw2uHjW4utfg6lGD\nu5caN89q3LzUuHuq7X96VdOksZbIEG9C/AKQ1ZpZvWwFJ3YexdvbmyFDhjJw4MBrrh7So0cPcnOP\nsGHDL+zevRuANm3a8NBDD13X+6JQ+GE0uqBSVeHn1wuQaNs2CIXiT3btKmL37q+prTVjtbqwYMF8\nOnToiF6fzMqVPxAc3JC+ffsyYsSI6zr37UYmnfvV8TalVqsdf/fx8amXGPbt2wfYf+iEu4eY97uP\nmPO7k5j3a2M0GnF2dr4pyazBYGDnzp1069YNSTJTWlrMlCnT+OSTD5AkK6dP5zNkyKPs3PkdFksF\nJpOW2loznp4X7/p4KRUV5fz003qsVis9e/akUaNGdXodQ4cOwWy28O3Sbygsg3fmlJJ5oB1FJ1Ox\nWZWO/dy9qvAOKMU7oBSfgFK8/e1/9/CuxNWjBhd3Ha7OCjzd3PB0dcPTzR0vd3e83DzwdnfHy90D\n77PbnBQX3lOVJKlOSvYdP36MwsJCgoKCSExMRCa7+rk3m20MHTqDtWt/IDAw6qLxGI1GPvjgA/bt\n24e/vz/jxo2jdevWNxz3xdzKf++Xy2Hr/A64zWZj2rRpfPvttxQXFxMWFsbIkSOZMWPGef9Yp0+f\nzueff05VVRWpqaksWLDgjutiJAiCIAh3i6ysLB599FGysrJwc3Pj+eef5+2333YsqbiU8vJyAgPt\nD8xZrVa+//57x11Mg8HA5MmTeO+9adhsOiorS1iwYBYtW1qQJAs2m5H7749Bq/0TAF9fKz/99CZm\ncyEAzs5KnJ2VFz/xJRw+fIg333wTvd4AwI8/pjF+/At07979msY5R5IkyqqrKFdXoVLryD7ijsxv\nJurSeDp2S6CsoIkj6ZbJrETEHCSmxU6aNN+Fh08Vzk5K/L28OJN3ivLCIowqHdrTesxGEw8/NJwh\ngwdfV1z289VNveyYmFhiYmKvev/qai1yuQtBQYl4e0eycGH0JZNvAFdXV1599dU6ifVOUecJ+Dvv\nvMMnn3zC119/TdOmTcnKyuKxxx7D1dWV1157DYB3332XuXPnsnjxYuLi4pgxYwa9evXi6NGjeHjc\nulaqgiAIgiBcWVlZGffeey/l5eWAPXF+7733cHZ25s0333TsJ0kSa9asoX//fkiSAb1eRevWndi/\nfxVKpRyr1cwPP3xBnz4RgBWbzUJUlAyNZidyuRxXV/jssxeQJPtDfm5uznTv/teDjwqF4qIdI6+W\n1Wpl9uw56PUGYmJi8PBwJzMzi48++oiWLVvg5+d/VePUWsycLCnmWOFpduyScWBnO4ryUqkojMZi\nPn9ZiUxmJTIuk9gWf9C4+W6iIuVEh4XTOLQtYf4BeLi6sX9/OusWfomfny8vv/wyJ07k8fnnn/P9\n8mU8MGjQJRPXAwcOsHTpUlQqFQ0aRPLEE2No2LDhdb8/dUEmUzJv3kY6dbqX4cMTkMlkNG/evF5j\nuh3VeQL+559/MmDAAO6//34AGjZsyIABAxxrhwDmzZvHK6+8wuCzv9UtXryY4OBgli5dypNPPlnX\nIQmCIAiCcJVMJhMrV64kPz+fmJgYhgwZwtq1aykvLyc1NZUtW7bw66+/MnDgQD777DNiYmLo168X\n7u42rFYNK1cuonlzKz4+roDEgQOfARVYLPbxFy4cj83210fzo0b1rpO4bTYrcvnl78ZXVqqorKzE\n09ODp59+CrlczqeffkZubi4nT568ZAJuMJkoV1dTpCojr7iIg7kmcvZ14sjeZ1AVR523r3dACcGR\nJwhucJygBicIDjtKsK+FLm1SiY3oi5f7hTcaq6rsdcSTk5uSlJRMfHw8n3/+OTqdHovFjJOTkuLi\nYsxmMxER4Tg5KcnKymTGjOnYbPaVxBUVFUye/DIffvghwcEh1/EOXj+tVk96+gl69x6Cq2sj5s3r\nfcvW29+p6jwBv+eee/jkk0/Izc0lPj6e7OxstmzZ4rj7ffLkSUpKSujVq5fjGFdXV7p06cLOnTtF\nAi4IgiAI9USr1dKrVy927drl2HbvvffSp08fAGJjY/n5558JCQkgNTWWwEAXTpzYSWGhRESEvWPi\nvHlPnz3y1jxitmvXLhYtWkRZWRnBwcGMGTOG9u3bX3RfNzf7WnGdTsepU/l4enpSXFwEgMLZmdKq\nSmoMeqq1NVRoqilX278qq2xUl4dTURjN0f1jOHOsuaM6iau7hoSUP+nY9QTNmlfSMNwJXw8v/Ly8\n8PVohdKp7RWvoWHDRo5r+eabJeTn5wMQFhaKTqdn1qy3yck5AkBgYCCvvvoKq1atxmaT6NOnN4MG\nDWLhwoVkZR1kw4YNjBr12EXPY18TLkcmc+O//11L796tiIy8WEOby7PZbPz5ZzadOjVHJnPGbPZh\n48bNPPBAnEi8r1KdJ+Avv/wyNTU1JCUloVAosFqtvPbaazz11FMAlJSUIJPJCAk5/7ezkJAQioqK\n6jocQRAEQRCu0syZM9m1axcREREMHz6ctLS1ZGTso2FDfzp3jqG8fBdffJGF1VqBXF5EdHRrxo/v\nX2/xHjp0iHfemeW4C1xWVsY778xi5sy3HN0rz7HZbJhsVjrf14vjp/P5ad+fuHi4E9GhM2G2IBat\nL8OoM2LQeqHTBFBd1ozq8giqyiIuaMcuV5iJabqH/oPzeHCwgcbhgcjlDYAG13UdcXFxDBw4kLS0\nNEfrdWdnZ55/fjxz5swhJ+cIbm5uuLm5UlFRwRtvvImfn73RTNeu3YiMbEC7du3IyjroePBPp7M3\nC/LwcEcu9+TZZ+cwbNgQBgwYilLphb//SYKCuuHh4YLVWsPYsROZMuVRgoM9ARkLF67hX//qj5ub\nBzKZE+PGvcPcuVPx8gpAJnNlwYJ5dO/+FF5efvj7w6efdrnotQkXV+cJ+LJly1iyZAnLli0jKSmJ\njIwMxo8fT3R0NKNHj76hsc89uVpf6vv8Qv0Q8373EXN+dxLzDr///jtBQX706NGWkSPbEBtbw+ef\nf0Vt7UHuvz+JH39M41zttODgYPr168fRo0fZvv13MjMzsdkkEhMTuffenigUCtLT09m8+Ve02hrC\nwsIZOnTIBTfgbsTKld8jyWS073oPKR3asy89nfwzp/lhy0YOlxZgtlowW6zUmAyodQZUZWGUnupP\nSX4cJafiUatCMZuuvKZcrjDjG1hIeAMVjaL2UFX2JmFB8OB9L4ER8vLUVxzjSrp27Yq/vz95eXm4\nuDiTkpKCJNnIyspCLpfz8suT8PDw4IMPPqC4uAQ3N1/MZhmffPIpLVu2ZtWqP4AQAgObceYMzJ69\nlKSkRPr1G4TN5sSYMRPw8vIiK+soAC1btqS4WEVxsf38nTr1paYmCpMJJElOcbETp0754+7ujiRJ\npKb25sS74VWZAAAgAElEQVQJK66uGkDD22+/x9GjeTd83fXlVvx7j4299IOrdZ6AT5o0iUmTJvHg\ngw8CkJyczKlTp5g1axajR48mNDQUSZIoLS0lMvKv4velpaWEht6ZxdQFQRAE4U6Vm5tLbu5Bhgzp\nSseOvthsetzdT6NSHSUraw8uLmaCgrzp1q0r8fHxFBScwdXVjfj4eJydlaSlpbFt22+O8U6fPk15\neTnJyUl8++1Sx/bjx4/zn/8s4KWXXsLHx/tioVyWxWZFbzJRYzRQrddRpddBgxBaN+mPRS5j98mj\n2Hx88JV3oVTtT95Gf3Qaf7TVAZQXNKY0Pw6TwfOCceUKM24eNbh51uDhpcXLS4ePn44GkVU0jqqh\nsnwTGenf0bNHV/r378/p06f58MNiamsDLhLl9ZPJoHnzZjRv/leTG51OD4Cnpw9paTl4evri5BSH\nzeaJQpGIwaBkzx49W7duR6vVERsbS5s2D2K1KnnxRfvSX6sVQLpiKefU1NS/7Q+jRz8O4GgmdKd0\nmLxT1HkCrtfrL1j/I5fLHRMYHR1NaGgomzZtIiXF3uLVaDSyfft23n///cuOXV81WkWN2LuTmPe7\nj5jzu9OdNO9paWm89957VFRU0Lp1a+bMmUN4ePg1jVFQUMAPP6xk7NgRWCxVODkpkcmciYy0MHBg\nF3bt2kBW1n6ysuzt2ZVKJx5++GGioxtfUIpOq63ht9/syfeECS/g6urG++/PYf/+/ahUFQAMHz6c\n3r17MWfO++Tk5HDy5EmGDRt20dh0RgOl1VWUV1dSrdWi1mtR63SodVr0JqNjP0kCbXUgZWe6UHam\nCWVnYlAVRaGr8b1s50hP30oax50iuXkpTZrkkRinJqVlI7zd3VE6/W9KpAT82bMnmsx0G7//vh2L\nxUp29mEAmjdvcU2l+a5EkiT0ehMeHq6AgnXrDnD4cAEJCUNZv/539u37FU9PT/Lz8wkICGD9+q9Q\nKBQsXryYwsJCEhISeOSRR1Aqr60s492mvuqA/686T8AHDBjAO++8Q1RUFMnJyezfv5+5c+fyr3/9\ny7HPhAkTmDVrFvHx8cTGxjJz5ky8vLx4+OGH6zocQRAEQfhHWLVqFUOHDnV8n5uby+7du0lPT8fX\n1/eix5SXlzNmzBh27PidVq1iGDiwB337dqG29hg63V4AIiPdiIy0PyiYlJTM669PZfHirygtLSMs\nLIwnnxxDdHTji45fVVWFJEFQUCA9evQ8G2c0R48eo6amBoCOHTsSFBRM8+bNyMnJQVOjQaVRo9Zp\nqdbVUF5dTWl1JWXVVY4kW5LApPdEp/FDpwlHp/ZHr/FDp/GnuiKMstMxGLQXu2Yb7p4qQsJMBAbp\nCArSExyiJyFBT/d7jDSOtjn2PH7cBLgS4H35u/Ft27Zl8OBBrFmzlt9//x2AJk0aM2bME5c97vJk\n5OQUsHPnEZ5+eiQymSvr12/nhx9+YcmSz3Fy8qZ//3bcf78Ns9lMZubAs7/UqAgJCeGHH37A399e\nsWX8+PE3EIdQX+o8Af/Pf/7D66+/zrPPPktZmf0f71NPPcXrr7/u2GfSpEkYjUaee+45RyOejRs3\nihrggiAIgnAJ5+ptT5w4keHDhzN69GgOHz7Mt99+y7PPPnvevjqdDrlczrhxj6HVHsfX14ROt49l\ny/ahVh/kmWeeueR5UlJSHJ9QX0lwcDDOzs6Ul1fw6acLcXNz4+jRYyicFDROiIeCMyxcsZSAsBAq\nNRpaDOxDgbPEgh9XYql1RqsOQF0RglrVEnVFKOqKUDSqUNSqECy1rpc9t6eXgabNKmmTUoOvTyY+\nPpkkJbnStm2L8zo1Wq0W5HLFVTelMZlMaDQa/P39UCickMlkPP74E9xzzz1nSxX60apV6yveaZYk\nCZVKQ2CgLzKZkuzsMt5550u+//5rnJy8iIxUkZQUgbd3G2QyGcOHN2XEiHGO4wMD/7r+PXv2cOjQ\nIUwmE02bNsXd/do6f16NLVu28OGHH1JZWUlKSgpvvPFGvXUgvxuIVvRX4U76eFKoO2Le7z5izu9O\nd8q8h4WFUVJSwrFjx4iJiWHy5Mm8++67TJkyxZGcS5KE1VpL7973MnJkFxYvno2XlwezZr1DeXkZ\nr78+FblcxooVK3B2drnCGa/Oxk2bWLxsKZ5BAbj7+eLm44Wrl32ttU7tR/GpeCqKotBVB6JV+6NT\nB6BV+2PSe112XKWzFl8/DY0aQXhYLRHhFsLCTEREGGjRopqGDQ1cLqcuKCjgww/ncuzYMZydnenX\nrx+PPjrK0bnz+PFjAI5lJJIk8c0337Bq1Q9YrTY8PT0YP3487dt3uKr3wWis5c8/D9OjRyoKhR/5\n+VWMGTOZP/7YhkLhSW1tLVVVVde8ZOhWWLduHYMGDXIsFwZo1aoVf/zxB25u19/46Hb0j21FLwiC\nIAhC3WvatCklJSWMHj2azp078+mnn+Lj40ObNs2ZPPlFGjYMYuTI+7BYSvnmm3+Tnr4Xi8VMgwYN\nCQ8PJyQkGBcXZ0ymWgwGw3Un4FarlUJVBfllxZwqLeaMqpT4Hvdgs8pRFTfiaEYCxacSKDkVj0Z1\n6eIKCoUFX38NoREaoqN0xMbUkpu9kjNnNuDqegInJ3vy8uCDw3j00VHXFKNWq+X116egUlUCYDLV\nsmrVamQyOY89dvEa2evW/cj339tLAHp5eVJTo+Xdd99lzpw5NGkSc9H3YcWK3xg58gGUykBkMmeW\nLVvCkCGvIZfLadFCxt69fRz728sI3p7J7KuvvorNZuOZZ56hf//+jBs3jgMHDrBixYpLvl/CjREJ\nuCAIgiDc5iRJYtasWZw8eQK9voQVKz4nIsKJBx5IpWNHVxISUvDwcMVsLgDA2VlJkyZNUCjkZGVl\nMWfOHCorVZhMtQQHB+PtffWfKEuSRFl1FSdLisgrKSS/rASzxYLFrKT0dCxFJzpSlJdE8cnEC0r6\nuboZiU0soFlzFdGNaolqZKNxIxvh4bX4+9fy95oNeXknmDDhbQIDXZg0aRLV1Wrmz5/PqlWrePDB\nB8nNzWX9+p8xGAwkJiYwbNgwlErni8a8d+9eVKpKoqOjmDXrHY4cyWH69Bls2PAzo0aNuuhylK1b\ntwLwwgsv0KNHDz78cC5bt25jx44djgR89eod9O7dFT+/hjg5+ZGbuxGlshVeXva7+T/8sOaq39fb\nyZkzZwCYMmUK4eHhDBkyhLlz53L69Ol6juyfSyTggiAIwh3l2LFjvP/++xQWFpKUlMTkyZPx8/Or\n77Dq3LkVolarheXLv+SXXzbw88/vsHFjGmq1miZNmtC9ezdkMggIuPBBwuDgEJ566mkWLvzE8fCg\nh4c7L7300hXXQ+uMBo4VnjmbdBehMxqQJKgojObEwWEUHm9KaX4cVsv5CXBYRDUpbSq4p0MN7dur\nSUioQXGJ7vAFBWc4ffo0fn5+JCQkUFlpv1sdExNDmzZtkSSJr776Eo2mhl9/3cJnn33qqEGekZFB\nTk4O06dPv2j7eaPRAEBERCTu7u6OZSZGoxF7h84Lr99isdff8/f3QyaT4e/vT20t6HQWnJyCUCqD\nyc//A0jG29veufI///n4su/jnSIuLo49e/bwzDPP0KdPH5YsWeLYLtwcIgEXBEEQ7hg5OTl06NDB\nsbZy3bp1pKWlsXv3bryvUM3idlZVVcVvv/3GgAH9sFr1ZGbu5/nnJ/LLL//FYimjZ89gevf+19mH\n9YZf9bj33XcfcXGxHDx4CKVSSWpqOwICLt563Fhr4siZfA7l53GypAhJkpBsMkpPx3I8swMnsjpc\nsKQkOqacTh2q6dJZQ4cOKkJDTVcV17fffsvy5csd37dq1ZInnhgDwOHDh/n668Wo1Wo0mho8PNxZ\nteoHJAkGDOhPQkIiCxcuJCMjk/T0/bRte2Gr94SERGQy+OOPHcyaZaWgwH6HNzEx8bwHNP+uZcuW\nHDlyihkzZpOS0oyjR8soKvLFySkFH59UZDIZM2e+c1XXd6eZN28ePXv2JC0tjbS0NMD+s3OpcpHC\njRMJuCAIgnDHmDJlCmq1ml69evH4448zffp0jhw5woIFC3jllVfqO7zz2Gw2R1+M2tpa9uzZQ6dO\nnbDZLBQWnmTMmKdZvfq/2GwGVKpCNm1ayj33uABWGjSwsGLFq44lJU5Ol7iNfBUaN25C48ZNLtgu\nSRJ6k5G84iIO5+dxvLgAm82GzSqnKC+ZE1ntOZHVHp36r4Tdz1/HfX2LuP++KlJTK/HzM19zPHv2\n7GH58uXI5XJatmxBbu5RDhzIYNOmjTzyyEi++eZbVq78AeBsJZdxzJ49B4DRo0fj5KTkwIH9bN78\nK2VlpUiS5Lijr9HosFisREdHM2bMGBYtWsSff/4JQEhICH36DCU9/QTt2qXg4tKQL79Mw9f3DJMn\nT+Cppz7gl1/KyMnJ5eTJk5SWVjB9+nRefXXKNV/jnaZ9+/bs3buXhQsXUllZSZs2bXj22WcdD6wK\ndU8k4IIgCMId48SJEwDMmDGDDh06UFZWxgsvvODYfr3sJeNUBAQEnG2IoufNN9/gjTdeA+zfP/30\neL7+ehFgL/P31FPPs2TJIsCGRqNhyJCRbNy4CkmyUl2tolmzruTn78Jms6DRVPPmm5P57ruZ2Gx6\nnJ3NvPrqYAyGQwD4+sKMGY8CFsCecF9L0q3RaFixYjlnzpwhICCAIUOGOrpNG0wmTpYW/a32thb1\n2S/L2baHZpMLp3PbkncwlVOH22DU//VpQmiYloEDSunfv4R27SovuaTkamVlZQIwdOgQHn10FBkZ\nB5g6dRqZmZnMn/8RjRs3Zv/+Azg5Kbjnns7ExcWxZMkSSkpKmTFjLu7uQZw6VYi3dwxHjkhs376M\n+fPfQy53Ji1tBbm5R5k27UUeeiiJ06c9yMg4yNixI2nXrgs7d2aiVpvx9b0HqzWdhx4aS2pqqqOk\n4ObNO0hPT6e4uJjExERiYi58+PKfKikpifnz59d3GHcNUYbwKtwpJaqEuiXm/e7zT5jzwsJCxo0b\nx59//om3tzfPPvssEyZMuOoayLe7Bx54gDVr1tC3b1+efPJJpkyZQnZ2NjNnzuS111677LFGoxFX\nV9ezpfqsvPTSRN59dwbZ2RmYzVoGDhxNTs4GZLJazOYaPv10JU89NQCZTIbZbGHLlv306dMOsK8X\n3rJlP717239WrFYrhw+fonlz+51mSZKw2Wy35A6iVlvDxIkTKS4ucWzzCw1m2KhHKNfVUKgq5+//\n1VvMyrM1t8OoLg+jKC+JM7ktsZj/qooS3VhD/35l9O9fTKtW1Zct93dlLiiVschkSmQyGSUlJVRU\nlOPj40tkZAOqqyspLCzEzc2N6OgmWCwWnJ2VgAyNxoBKpSE4OIATJ05hMpmxWm1IkoS/vz+JiYnI\n5fJL/nzbbDZsNhtOF3S5xNEo6NwDlMLdoa7m3dnZ+YLO7//rcjmsSMCvwj/hP2Xh2ol5v/vc6XOu\n1WpJSUnh6NGj521/9913mTRpUj1FVbcOHjxIx44d0Wq1jm0xMTHs3bv3gm6QaWlp9O7dG7ncjMVS\nQ1RUcw4c+BEPDxlWq45Fi37g0Ud7cfr0ybPj1F1b8aslSRI1NRrc3NwuWdHjSr5b9h1pv2wgonEU\n8c2bUVRdiezs3fNakyuFx5MpONqCiuJGqMvDqKkOvGi79uCQLJ4co6BfvzJiY7UXvH59XHB1bYmH\nR9A/5pdAQZAkCaPRiIuLy2WTcFEHXBAE4S7w448/cvToUWJjY1m3bh2//fYbY8eOZc6cObdtAl5R\nUcGECRPYvn07np6ePP300zz33HOXTNaaNWvGnj17eO+99zh06BAdO3bk9ddfx8PDmYcffpAZMyYR\nERGI1apj9eoviY3VExTkCUgcPLgIhaICi32VB0880ffWXehFZGdn8/77cygvr0Aul9O/fz9Gj378\nsnfNrTYbBRVlFJSXUVKlorRKRYVFR+K9XQAo0mgoL47ldG5LzuS2pCgvEZv1/I6NCoUNV9c8nJ1z\nadFCRvPmRjIyJmE2H6Nr1xnExraqs2tUKmNF8i3848hkMlxdXTGZTLi6Xr5j66WIBFwQBOEfQqVS\nAfYHquLi4ggODmbs2LFUVlae90Dg7cJkMtGrVy8yMjIc28aPH49Op2Py5MkX7H/o0CGCgoJo3DiE\n//znNV58cQqjRrUD9qFWWxk9uj0uLqcwGOxLMWbPfvzskfYPem+nB8rKykqZMWMGBoPB0RwnLe1H\n3NzcGDnyEcd+kiRRoakmr7iQvJIi8ktLqLXYH3y0WeVUl4dTUdic8sJoKgqjKStoglH31/ptudxG\nq1YqevRQkZJSRePGOkJDq3j4YXt1ixUr1iCTKZg9O4Tt24+hUlXc0HVZLFZ++mkXgwZ1Q6FwB9xF\n8i38I93oz7VIwAVBEP4hUlJSAFi2bBleXl4cOHAAgNatW992yTfAhg0byMjIID4+juXLP+ftt99j\n9+69/PTTMv7972dYvfongoKC6dKlE2aziq++mkfXrkl06GB/MO6tt0aeHcmekLZufefULN61axcG\ng4GWLVswfuJE0jMO8O2yZew5kk1k1n50RgNao4HCinJqDHqAsyUBYzhxsD0Fx5qhKm6E1XxhN0sX\n1wIC/DcxYUI4gwdLF6lU4kxQUCDl5RXMnfshUVFR7N69B4CIiIhrug5Jkpgz53smTnwGD48Q5HJP\nfv31Ox55ZDrOzi5n624LgvC/RAIuCILwD9GhQwdeeeUVZs2axccf2xuEBAYG8t///reeI7uQJEns\n2rWT0FAfRoxoRkSEmunTH+Kxx37Cai2jsnILHh5nUCiqUKtNgI3JkwfWd9g3xGK1UFZdRUlVJfk6\nDQk978EpwJ+P19lL7sXcY3/A87eDBxzH2KxyCk8040RWe/IOpp5XEhAgPKKG5s1qiI9XUV39E3L5\nn0RF1fLAA4MJC6u9ZCzjxj3LW2/NZNu2bY5tvXrdS0JC4hWv49NP1zF06CAaNIhDofAlPDwPZ+ck\nPDzsa/CXLv2rvre4+y0IFycScEEQhH+Qt99+m549ezqqoDz44IOEhYXVd1iAvVLIrl27aNu2GSZT\nAX36hLFxo5otW9bg5mYmP/80Tk7QoEEkSqUTHToknT3SVq9xXw1Jkqi1mDGYTNQY9FRrtVTraqjW\n1lClrXGU/vt73QOPgAA0lcGUF0Zj0PpgMnhQa3THZPBwfJWeivufkoA6+vcrpU+fUlq1qsbHx/K3\nKJLPfl1ZSkoKs2fP5tdff8VoNJKUlEzPnj0vmjCXlFQCTjRokIBSGUJkpApPz+Z4eNjLHL7wwoTr\necsE4a52UxLwkpISJk+ezPr166mpqaFJkyZ88skndO7c2bHP9OnT+fzzz6mqqiI1NZUFCxaQlJR0\nmVEFQRCEq9GzZ0969uxZ32Fgs9mYM2cO48Y9gUJhwmzWMHXqRL788v9wd3emadMmPPLIML7/fiVr\n19q777m7uzF+/Ph6i9litaLWae1Js7YGtV6H2WLGYrX+9WWzYrFYMJnNGGpNGGtrMdSauFxRMatV\ngaYijPLCxpQXNKGsoDHlBY0x6a9cCi0mpoZ+/Uro16+Yli3VN1gS8C9NmsTQpMnF61zbm9socHIK\nZNWqHTRunExycmtkMhmPPfavuglAEO5idZ6Aq9VqOnXqRJcuXfj5558JDAwkLy+P4OBgxz7vvvsu\nc+fOZfHixcTFxTFjxgx69erF0aNH8fDwqOuQBEEQbqklS5bw1ltvUVZWRosWLfj4449JTLzyR/t3\nOkmSePnlSTzzzKMEBblgtWqoqcmlrGwz3t7uACxffn6t7kcfHUVycjLZ2dm4urrRpUtngoNDbnqs\nOqOBsuoq9mVlkpF9CJyVuP3pge1/1srbbHJkMttlk16bVY5B54NBa//SVvujqQyhpjIYjSoETWUw\n2uoAJOnCh0D9/PRERRUQFqajYUNPfH2t+PiY8fEx4+1toXFjXR2WBLwaMrZty+XHH3exaNGXODl5\nMHVq6i08/51j8eLFjB49+qKvyWQyfv75Z3r37n3N4+3atYt27do5tut0Ou677z727NnDd999x5Ah\nQ5gxYwYzZsy46Hlnz57Niy++eO0XJNxSdZ6Av/vuu4SHh/Pll186tjVq1Oi8febNm8crr7zC4MGD\nAfsPXXBwMEuXLuXJJ5+s65AEQRBumaVLlzJq1CjH99u2baNr165kZGQQHh5ej5HVDavVel5jk1mz\n3qJ791SaNo3EYqmieXM3zOaDGI32ZRMvvDD4imO2bp1C69YpdR6r2WKhSquhqqaGksoKThWcQVtr\nQm+x37kGsFqcqDG3orownJrqQLRVQdRUBVFTHYC2Kgh9jR8ymRUn51qclCaUziYUylqclLVYap0x\naH3OWyJyaTacXU7j4Z5BcPBRJk5sQ/v2ZkJDjXV2R/ty1GotPj6ewNnKKhVqgoJ8Ha8tXfo7EyY8\nj7NzKH36dKVHj9EolZ43P7A7nEwmY8aMGTRu3PiC11q0aHFd4/2dXq+nb9++7Nmzh2XLlvHAAw+c\nt++CBQvw9j7/5+/cw9jC7a3OE/C1a9fSt29fRowYwdatWwkPD2fMmDE8++yzAJw8eZKSkhJ69erl\nOMbV1ZUuXbqwc+dOkYALgnBHe//99wH7MrvRo0fz8MMPs3PnTpYsWcLLL79cz9FdWU5ODl5eXo5q\nGMuWfUdUVENat26Kzabn8cef5f77uzBwYGesVj1Nmzrj5VWIyWSvdnHffW1vabySJKE16ClXV1Ou\nrqZCXU25ppqqGo2jegjYE+2qsggqSxKoLGlIZUlDVCUNUFeEIdkuX55QkhSYTW6YTW4YLvK6TCbh\n719LYGAtgYEmgoNNNGyop1EjPbt3f01u7lratPFmwIDeLFu2jFOn8tHrHyIs7JGLjHbjJEkiLW0n\nAwd2Ran0xWp1o337e8nPzwQkrFYL3bq14+TJvYCEUilHLj+Bh0cCcrkclwsLqwiX0bt37/PuWNeV\nc8n37t27L0i+zxkyZMh5KwyEO0edJ+B5eXl8/PHH/Pvf/+aVV14hIyPD0VRh3LhxlJSUIJPJCAk5\n/yPGkJAQioqKLjv2uS519aW+zy/UDzHvd58bmfOSEnsN6ri4OMrKyoiPj2fnzp1kZ2ffNj9L9vW9\nMuRyOatXf4+fnwc9eqQgSQa+/no1jRuH0qFDEpJkprT0AJLki0xmT9aef74TcrlETs7vAAQHO6HV\nVnL8eOUti71Kr+N0ZTllGjVqgw6z1Xr2NTDqvNBUhlBVmkBlSSSVpQ2oKo1EXRF60SUgMpmV4NBS\nGjaqxGbNRlW+m9RUf9q29ePQodXodEfw9fWjS5f78fdvgNGoxGh0wmhU4uZmxtdXj4+PEYXi4uu/\n9+79ATe3fLp0eZrg4GDatm3HqVP5HDt2jOPHj13yOtVqNUePHkWSJGJj4/Dz873kvjKZgo8//pnH\nHx+Ot3cQ4MmPPx4mJqYvHh5KwML69T+TnV3sOOann9aTnV3u+L5Pn/vYv3//Fd79a9eoUaPrblTy\nTyGXy3n66afp06cPU6ZM4dixY8TExPD+++/Tp0+fix5jMBi4//772bVr1yWTb6H+1dTUcOjQoUu+\nHht76e66dZ6A22w22rVrx1tvvQXYP4I5evQoCxYsYNy4cXV9OkEQhNtKQkICRUVFTJkyhXbt2pGW\nZn+4sL7WgOt0OvR6PSEhIchktSxZ8hUGQw1PPjkIm01DQoKEm5sBozEHgOHD7V0QLZYqANq3P/cf\niD3BvJp64haLBZOpFnd3d8fyCovFyv796ahUKvz8/GjdOgVnZ+XlBzp35nNJt6qc05XlaE1GKksi\nyT/S9ewa67PrrSuDMZvcLjqGTGYlJKyMBg0rqKrciEzaS6NGZWi120lKiiI2NpbNm3/FRWmgffsH\n2bDhZ2pq7GuvCwvzOXbsMOPHjyci4tLLiMrLy8nIyMBsNhMTE0tcXCz+/v6cOpXPpk2b0en07Nix\nA4CAgIBLjnPs2DG+/PJLjEb7MhlnZ2cee2wUcXHxODkpUSg8mDx5EWPHjiI2NglwIyFBj+3/s3fe\n4U1W7R//PEmapnvvCXSwZ4Gyl1RAFFEcbEVUfOWnCAiigKAIThBQgRdRQWXIK4iyBESmbJBZCt0t\n3W2aJml2nt8fgUKlTAs4ns915Wp7ctbzJE2+5z73uW97NDab41Dn5Mlv/uH6pXCAdwqNRlOVBOtK\nrnyNf/vtN3766SdeeOEFPDw8mDdvHgMGDCA7OxsfH59q7fR6PX369GHfvn03FN+lpaXVEkzJZLKr\n+pP4a1LrAjwkJOSqL5oGDRowb948AIKDgxFFkcLCQsLDw6vqFBYWEhwcfN2+ExISanu6N8Ulq9W9\nGl/i3iC97v8+auM1X7ZsGZ07dyY9PZ309HQAhg0bxsSJEykoKOCNN97g9OnThIaGMmnSpFrfui4q\nKiI7O5tmzRpitWr4+uuvycrKYOzYAYiigZdf7oFK5XTxS9uF6OjaO/Bos1lZsmQJGzduwm63Exoa\nwvjx44mMjGTKlKkkJydX1T1y5AjvvDMTFxcXLpQWk1GQh8Fkwi7asdvt2EURm93xe25JEWqdFoPO\ng3PHunP2UFeKsmtOuqNyqcTLu5CQkCLaJyoJCytg1aoJBAYW8N//fozNZuf551/HYDAwZswk5szZ\ny9mzKZw9mwJAu3aJWK1WtFodDRrUZ+jQYaxbt44DBw6wf/8BJkyYUOO4p06dZPbs2ZhMjtjb27b9\nwuDBg3juuecZP34cqamppKamAo4d36eeegoPj8u+u2lpeXh6uuLl5cJbb71NWZmVRo0aEhDgTUpK\nHq+/vpqPP/6QRx55EoXCjUWLuhEeHo7zRX+RFi3+mgcl/+mJeERRrNGKLQgCBoMBpVIJwNmzZ0lO\nTqZOnToAdO3alWbNmrFixYpqxklRFBkxYgR5eXk3FN+iKNKoUfWwk/7+/hQVFdXGpUncBB4eHtf9\nvsmZlt8AACAASURBVNBoNNd8rtYFeIcOHUhJSalWlpKSUnUQs06dOgQHB7N169aqgwJGo5Hdu3dX\n+U5KSEj8tTCZTMyYMYOdO3fi6urKyJEjGTBgwL2e1l+SsLAwjh49yqpVqygsLKR58+b07duXsrIy\n2rVrR1ZWVlXdjRs3smPHDtq1a3fb45WXl3P8+HHat2+DxaLm+PHtbN/+K9HRjwB2HnmkAdAAUXT4\nQ7u53Tl3gGXLvmb9+g0AODkpyMvL5803p5GU1JPk5GR8fX3p0aM7u/buocxs5LPvV4CLikqTEZPB\nFbPRtaovQbjk0iFSmN2Qs4e6kXmmFXabw2qucjXQrUcGLZqZiK1nITKykjVrZnHkyJYqq3t6uhfP\nPjuLM2f0pKeXXAyfJ2K12qhfP562bRP5+OM5rFq1isrKSlq1akXPnj2ZPXs2AD163Efjxo2xWi0c\nOHDgmmnaRVFkzpw5mExmWrduTUhIMOvXr+fbb5fTrl17Zs+ezZo1ayguLiYsLJzHH3+M3347h49P\nEJ06dUUu92b9+t20bt2KOnXCyMryxmz25umnP6Rr1/to0qQpBQUaoqObVB2MrFevXq2/fn8FNu5f\nweYDq25c8Tbo1fYJ+iQOrNU+BUFg/vz51K9f/6rnLolvgO7du1eJb4AmTZrg6elZtUi/kqKiIlQq\nFZGRkTcce/Xq1Xh7X3ZPunJMib82tS7AX3nlFTp06MDMmTN54oknOHr0KPPnz+fdd9+tqjNmzBhm\nzZpFfHw8sbGxzJgxAw8PDwYOrN1/DAkJiT+P3W6nf//+bNq0qars559/ZsGCBYwaNeoezuyvi5eX\nF88991y1skWLFpGVlUWzZs2YPXs2ixcvZuXKlUyfPp3NmzffdN8Wi4WDBw/Stm1rLJYycnJO8O23\ni2jUaARgp0ULL1q0eJh7kbxm27ZtALz99ts0atSQN954g+Tks/x+6iTeocF07pWEzVVFZKe2iEBR\nuS9pe9qRdiKRvLSGNfpoX4kgs5PYIZOnhhbRp08xLi6Xr3H37t0cPboFd3dXHnroIY4dO8bZsyl8\n/vlipkyZwvvvf1BlgW/atAnjxo1HEARCQkLp27cvADExDneb8PAIAL7//ntsNitbt269WB5OTej1\neoqLS1AoFLzxxuvIZHLy8ws4dOgQmZkZdO7clfr1OxEeDo8//hhyuRcBAX54efni4dHw4j17H4C0\ntDTOnMnE1dWVqKgYsrNzqs5HeXjcOGa4xN0nISHhhjtZERERV5X5+PigVqurlQmCwMKFCxk/fjy9\nevVi165d13Vf69ixo3QI829KrQvwhIQEfvjhByZNmsSMGTOIjIzknXfeqfZFPWHCBIxGI6NHj65K\nxLNlyxYpBriExF+QnTt3smnTJry9vVmyZAlnzpxhypQpTJo0ieeff17yLb1JcnJyABg6dCjdu3fH\n3d2dlStXVpVfj+3bt9OpUyfsdj1a7QWmTXuVb7+djExmIyQEZs58ir9Ctkiz2YTcSYHdWcGBlGQ8\nY+vQNCYSpasLPkCuXoMmy4W0E/1IO96OgqzLVkO5woqPrwa4+H4ShYte5wL+/gaefLyAAQPyCQ42\n1Th2ZmYmAL1792bgwEF06dKFUaNeIDMzEz8/f9599110Oi2CIODufn0h269fP377bS+ZmVksWLAQ\nAD8/XwYNGlxjfRcXF0RRjlZrZ+/efYSF1eXw4RJycrwJDOyMp2cn4uMjMBqNuLs7XGeSknrX2Ffd\nunXp2bMnW7duJS7usptNly5d/hWx5P+pXOmnfSU1JW+Kj49ny5YtdO3alaSkJHbv3k10dPQdnqHE\n3eaOZMLs3bs3vXvX/OFyialTpzJ16tQ7MbyEhEQtckkg9ujRg0ceeYT+/fszffp0ysvL0Wq1V8Wg\nlaiZ+Ph4AObNm4erqysrVqwAHIc2t27dSocOHXB2dsZuN/PYY48xb94s/PxcsdmMfPzxDOrUGYOH\nhwwQWbFiEmC7dxfDxVjS6lJKKioo1WspKCulyQM9wUnBT0f2A1BpDaAoJ4ai7HqOnzkx6Ct8q/pw\ncjLTqWsejz5cxv33F/4hrfqtcWkb/vDhI3Ts2JHdu3dXKxcEoZrP9fVwcXFh6tS3WLlyDaJoICgo\nmLp1m7J58ymGDeuHIDixZ88J1q3bzpw5byGTOdOmzXCWL/+el1/+moqKSioqKmjZsiVt2/bGycnp\npmNCC4LAd999x6hRo1i7di2iKPLggw+yePHif8Vit0/iwFp3E/k70qxZM9avX09SUhI9e/Zk9+7d\nNzwnJ/H34o4IcAkJiX8Ol4Tjhg0bmD17NmfOnMFqtRIcHCxtid8kaWlpDBkyhI0b15OTk8rEiWOJ\njQ2kX7+WzJw5kpkzPyE0VENQkAdgZdSoLshk5zAYHP6cX3zxysWerp3q/E4giiJ6o5EyrYYybUXV\no0hdRnF5GcIVVj2LyZmi3KYUZsVRkBVHUXYMWvXVW+NOSg3h4XsYNtTI00+74OZWOwuJbt26sm7d\nD2RmZjJmzCtV5Y8//vgN2wqCjOxsDT/99BMTJvwHudyFCxdOkZfnzaeffoMgKMnIyEKvP4G3d0cA\nOnVqQNOmvXB3d/jpzpmzmIiIBnzxxRfI5c48+uijfPTRRzg53Vyklyvx9vZm5cqV2O12RFG8pvVU\n4p9Nhw4d+P777+nXrx9JSUns3LlTinDyD0IS4BISEtelTZs2jBw5ks8//5xx48YBju3UBQsW/Css\ncreD2WzGYrGgUBix2bRMnPh/vPzyk3z11cts2LCWkyfPEh4exEMPPYifn5WPPhp5qSUACQlXH+i6\nE9jtdgxmEzqDAZ3RgLZSX01ol2krMFstV7Uz6D1QFzamrCCCopwYCrNiKS2IvCqhjbPKRJMmahJa\n6WjRQkPz5uVER1fiiGSopDat+O7uHrz77rssWbKE9PR0PD09efTRR0lMdBxwtdlsFBWVExLihyAo\nycys4JVXPmTz5jUolf74+Ojp0MGMl1dzAHr0iKdHj0er+o+Nja0W09fLywsvL6+qv2UyGePHj2f8\n+PG1dk03E/JR4t4iiiKbN2/m/PmrY7q3bduWmJiYW+7vSnr16sU333zDwIED6dWrF7/88gvu7lKG\n0n8CkgCXkJC4LoIgsGjRItq1a8eOHTtwdXXlqaeeIjEx8V5PrdaxWCyUl5djtVqrUq3fKna7nfHj\nRxMfH8ATTzju0cKFL1Q9379/f+5mTg2TxUyJRkNJxaVMkWo0ej16owG9yVijDyqA2ehCRVloVXxt\ndWE4ZYXhlBVEYNBdnRRGEKzExhWRmKgnoZWGVq3UxMTouJvGW5nMmYkTJyEITuh0dsaM+YhOnZ5A\nJnOmrEzH4MFvcOjQHhQKD5o0sbNyZQdcXUOx2fLx9vaWQo5K3DKXUtHXxPz584mJiUEQhBqNFTWV\n11Tvscceo6Kigueee45+/fpVOxAv8fdFEK/16fsX4coYildaG+4mUjzofyfS6/7vYu7cuUycOBGT\nyYSHhwfz589n+PDhN9W2srKSw4cPk5AQg8mUgdFYhJPTvbFvmCxmzmRnkpydQaG6rFo6dnCkZNdr\nfKnUelGp9aFS601lhTeVOm/0Gl+06gAqygIxVV7bvUgu1+HllUaPHj7Ur69h8+bxCMJuliz5mODg\nkDt9iYBjoXPgQDKJiU1QKLyxWl1o2fJBkpN/x8XFB1F0HF5NSkq6YV/S//qdw2g0/uszYUr8c7nR\n+/t6GlaygEtISPzr+e677xgzZgwA7u7uaLVannrqKcLDw+nRo0eNbXQ6HXPnzmXr1s14egoUFxex\nZs1MnJyc7rr4ttvtpBfkcSL9PGdzs7DabJiNKsoKIlAXOSzX6sJw1EVhVJQEY7ffeH7OzlZ8fIow\nGk+gVGbh4pJCs2YWxo7tyowZT6HX61Cp6pGSYkahyCE4OIiAgLsTDk0QnHF2DuG//11IYuLTeHsH\nIggC2dm51SyINyO+b4f09HT27duHq6srPXr0kA4iS0hI3DKSAJeQkKhGVlYW58+fJyoqqprP6z+Z\nZcuWAfD8888zcuRIVqxYwezZs/niiy+qBLjVaq3aMjaZ1ISFRdO6tT9OTuWUlmoQBJHZsz9iwoSJ\nd9w3XhRFNHodReVqMovyOZWRhs5oAKCiNJCj2/tz5kAPbNark3IIgp2AQC1mcyoy4QLt20cSG+tC\nQICJwEAjkZEGwsMrKSg4zOTJkwFo2LAB6ekZFBcb2bcvm0mTJjJr1ixSU9MA8PX1ZdKkSXf4sKCM\nyZOX0aNHEv37D0MuV/DTT9W34u/GmYRvv/2WESNGYDY7/PUjIyPZsmVL1WFlCQkJiZtBEuASEhKA\nQ9RNnTqVGTNmVJWNHDmShQsX/i2jMGRlZbF3715UKhU9evSotv0niiIffPAeDzyQRN26oURHu+Ln\n50pQkAWZ7CwJCe6Ehnrg61tJRcVpZDIFHTr0Z8GC6cTH+6DRFNKggQ6DQcfEia+hUjkzc+ZM9u79\njQsXLlwzYcutYrfbKdfrKNNWOHy4y9UUlasp1pRfdTiyrCCcI788QsqRLlWHIes3KKN+fCVhYYVs\n2fI+fn6ZrFjxOm5uAmPHvkJaWjpjxkynefMWlJaWUF5eTnBwMG5u7mzadAyABx/sy7PPPseJE8eZ\nPHkKR48eYeTIkSxatIizZ88ik8lo2LARrq6uV83/VsnPL8VgMFG3bjiCoGD9+kMYDFaGDBmEUhnM\npEmNCQkJQaG49cgitUFmZmaV+E5KSiInJ4fk5GQGDRrEkSNH7smcJCQk/p5IAlxCQgKA1atXM2PG\nDARBIDExkcOHD/P555/ToEEDxo4dCziE67p160hMTCQgwB+73cJrr73G008PIj6+KQqF6i8RGWX1\n6tUMHToUk8mRtCUyMpLhw4fSuXNr2raNw2otIyzMgN1+Er3+Ak2auHD0aCU7d35NdvZe0tLSiIqy\n0qFDBGazI1X0tm3vIggComigvLwcmQwCAgJo164dgiAQGRlJamoa5eXqmxbgoihSaTKi0eupqNRR\nUalHrdVSejHsX7lOi/06x3REEYpy6nFk2wDSTrYFUYZcbufRx3N46aVU4uN1ANjtNnJzf6WoqJg3\n33wVNzc30tLScXZWEhUVxfz586uyPTo5KRg1ahQKhUPEV1RUIIoiFRVa4HJCEU9PL9q0aXtLr4vd\nbkej0ePj4wnI2b79BFlZpYwa9TRyuQvHjv1MWZmGli2HI5MpiYvzZOfOnaxcuYsuXbrc8x2Z/fv3\nYzab6dGjBz///DM6nQ5vb2+OHj2KTqeTolNISEjcNJIAl5CQABxxvgFmzJjB66+/zrJlyxg+fDgb\nNmxg7NixWK1mDIZUjh3bRFhYOQpFEGCnc+dgnJ3TqKgox2YLZcGC75kyZeo9C6F24cIFhg0bhlwu\n5777OhIS4kxRUTpHjiznoYeiMJkc87r//ssH7h5+uD8ZGZns2rWLlJQUAB54oA89etxXVefKhUVQ\nUBAqlYqiomLmzZuHi4uK1NQ05HIZYWE1i29RFMkvKyU1L4esogI0eofgttquHYpPFMFsdEVX7odW\nHYCmNJiKkqArfgZhtbgAIJdbGDIkm9Gj04iOrn7wUiaTM2HCRKZPn14VLk2pdGLcuHH8/PMWtm7d\nikwmIyQkhAsXLvDJJ5/wf//3EgqFnJ07d/H777+j0VQAjoyMN0tlpYn8/FLq1YtCofBmzZq9HDhw\nkvnz5yKXu9K6dX3i4yvw8HBkeBw27PmqtkeOHKFfv34UFxcD4OTkxFdffcWgQYNuevza5pKVPycn\nB41GQ2pqKjabDblcjlJ5tbuPhISExLWQBLiEhARw2bJZUVFR7afJZGLKlAm88kof7HYt//d/fS+2\ncGQt7NChMQCiaKai4hTu7hqMxguoVGF3VYQnJydz5swpPDzsdOwYga+vGyNH9qdhw0gGDRqIxWLF\nz69mNwm5XM64ceNo1aoVanUZrVq1Iioq+ppjOTs7M3r0i8yePYdffvmlqnzkyGerJcowmEyk518g\nNT+X1Lxc9Bf9tI2VblRW+GCsDMOo93A8Kj0w6D2orPBBp/FDV+6HXuOHxeRy3et2dlbj7/8FwcEf\nkZT0CNHR99dYLy4ujgULPuPYsWPYbDYaNWpMUFBQVUbO116bSGJiOxYs+IxNmzaTk5PDhAkTmTdv\nLhpNBTKZQN++D/LII49ccy5FReX8+utxBg16CJnMg/T0DJYu3c3nn49AJlPw9NNtGTHi8kKmTp06\nNfYjiiJPPvkkxcXFtG7dmqCgINavX8+IESPo3Llzrbn43Crdu3enTp06nDt3Dh8fn6oQjsOHD5cE\nuISExC0hCXAJCQnAkTHwyy+/5L333mP58uXk5OQQGhrIqFF9UCjM2O3aG/bh7e3O8OHdqaz8HZMp\nh++/P0Ljxq3vSMzwsrIydu3ahUpl58KFZCor1djtFho29EWvP4+7exBxcSEUFORjsTgWC0ql8zX7\nEwSBiIhwIiLCryu+L9G5cxdCQkLZv38fdrudhIQEGjVqjN1uJzUvl0PnzpBekFcl0sxGF1J/78HZ\nw125kNrkpq9TpbLg4VGMp2cxdesaadvWg6NHvyYnZwv/+U9XBgzoyoYN51i0KJ8jR46QlFSzAAeH\n20iXLl2rlV3ycLlk4b+0aBJFkcTERBISWlFSUoqnpyeCIOe3387QsWNzZDJXcnIqeOGF6Wzb9j2C\n4IK7uxadLg9v704AdO3agq5dry3Yr0VRURGpqam4urqye/dunJ2d6dSpE3v27OHIkSP3TIC7u7uz\nZcsWhgwZwoEDB1AoFDz11FPMmzfvnsxHQkLi74skwCUkJABHxrXPPvuMsWPHkpeXR/PmdRk1Kon7\n7298W37dNlspnp5qnJxy0evTUSi8Wbt2C926dSMoKOiq+larlQ8++ICxY19EEERAzqRJb/L+++8j\nk8kQRZGff/6Znj27YTaXUVR0mvfeG4eTUw7WiwcS69SJpnfvNwkNDSEvL59BgwZVCcyuXbvg7Hxt\nAX47XJkd0WAysS/5JIfPJaPWORYrdpuM7JSWnD3UhfRTbbFZHOM7Ka2Ehenx87Xi52fBx8eCj48Z\nHx8LAQEmQkMNhIYa0Wp/54MPXqvyZdfpQKnsTf36atTqo9jtDQHQah2+3jLZrR+WbdcukczMTN57\n7z1CQ8PIzs5GECAxsS16vYW33/6a2bOnI5e7odGYWLFiJX37vopcrsDDw8rq1a1wdQ0DIDIyiIkT\nJ/65mwq4ubkBjvjqJ0+eJDw8nMzMTAA8PK4dn/xuEBMTw/79+6msrESpVN52wiYJCYl/N3f8k2PW\nrFm88cYbjB49upqVYNq0aSxevBi1Wk3btm359NNPadiw4Z2ejoSExHUYNWoUAwc+yMKFs2nY0I/2\n7Zv+qf66d28BgMFwGoDjxzfTsqUvHh5xyOVeJCU9yKpVC3B1FbDbdZSXn6Gk5BecnRVYrSI+PhWU\nl/+CICixWAQ+/fRdWrSwolDA4cObEIR0nJ3d6N37fg4cOEhGRiarVn3H9OlvMXv2RyQnO6J0dO3a\nhRdeeOF6U70tHH7dJRw5f5aTmWlV/twGnQfHdz3Iqd+SqmWNbN++hMcfz+XBB/Px9LTesP8XX5yN\nyWQiISGBevXqsXbtGjZt2sTgwYPYt28/33zzDRs3bqSsrAyALl06X9WH0WhGqVRctGzLmDNnNc8/\n/yhubg7XlsxMGT17PsiBA/tQq42kpLgwb97btGkzFLncg44dK3B3b4xMJsPNDVat+l9V305OToSF\nhf2ZW1gj7u7uDBs2jGXLltG6deuq8saNG9OhQ4daH+92qI2oLxISEv9e7qgA379/P4sXL6ZZs2bV\nyt977z3mzJnD0qVLiYuLY/r06fTs2ZNz585VWT4kJCTuHufOneOTT2bz9tsjsVoLGDmy2x0ZZ9y4\nRwETlZUnAZg+/XFEMRWz2fFR9Oqrj1fVVSgERo16CFE0IYom5HL46qvx1eYM8OSTT9KvXz/atGnD\nlClTOX/+PEFBQbz33vuYTCYUCjly+c191FlsVioMBny1FXi6uqKooZ3OYCA9P5e0/AukF1xAbzRe\nfq7cj2O/9uP0viQsZkd2tJgYHY8/nsOAAReIiDBcc+z9+/dx9OhRFAoFnTp1Ij4+npycHAAmTXoN\nJyclRUWF/PrrDjw9PXnmmREsXbqMsrIynJwcrhCJie2qXF4UCn8UCi9at+7Ntm0/EB4ehUymwsvr\nNP7+PXBxcQjw6OhknnnmGUpKSigsLESr1dK5c+cqy+5TTz11U/eutlm4cCEqlYply5ZhsVjo3r07\nX375Za3vYkhISEjcC+6YANdoNAwZMoQvv/ySadOmVXtu7ty5TJo0iYcffhiApUuXEhgYyPLly3n2\n2Wfv1JQkJCSuoLi4GB8fHyyWUjw9C+jYMRSLJf+uzqFx45oP4d0M7u4OV4QTJ07QpUuXqnTiHh6X\nQ8HdjFgr12k5dyGbc7nZZBbmYxdFfj7tiIHt6qzC09UNT1c33FQq8stKKFCXXd1HcYgjBvehbths\njo/V++4r5KWXUklMLONGHjzffvsNq1Z9V/X3xo0bGTt2LF5eXmg0Gn76aT3169fn5EnHwsXHx4fE\nxHb06tWL0tIy/Px8cXZWAQKvvbaU3r0f4tFHH0QQBNLTs6q5EE2YMKHa2P/3f/8HOEI1RkZG3vB+\n3S1cXFxYtGgRCxYsQBTFv2UsegkJCYlrcccE+HPPPcfjjz9+VciqjIwMCgoK6NmzZ1WZSqWic+fO\n/Pbbb5IAl7jn6PV6fvllG5GRXths/5wdmS1bttCyZUs8PZ2wWnX06dOHxYsnEx6uQqG47C7yd+H+\n+5PYvHkThw4dYtiwYVXlDz3U77rtDCYTBepS0gvyOJ+bTZFGfc26lSYjlSYjBerSq57Tqv3JOtOS\nrOS2ZJxpgWiXIQgiDz98gZdeSqVJk4qbuo78/DxWrfoOmUzgiSeeQK1Ws3nzz3z22WcMHjyYxYs/\n56uvvqqqHxNTj1atHCEUnZ1VBAcHk59fRnR0HZydo3nrrRYEBwdXie6/Qlz2P8O9CmcpISEhcSe5\nIwJ88eLFpKenV4W3upKCggIEQbjqEFZQUBB5eXnX7feSheteca/Hl7jzCALo9Tl89dVcXn31QZyc\nAnjrrXW0bJlIcHDwvZ7eLXH48GEiIkIJC3NHFLWsXPkFFRUtCA/3BuzMnz8cozGH1NR7PdPLJCcn\nc+KEw8rbpEnjG54LGTVqFGvWrKWoqBAvL2969+6Fl5cnqamOWNd6k5EyvQ51pQ61Xo+6Ukel2XRV\nP1aLE9lnW3D+9w4UZjbA3UuNl38OXoF5eAfk4ROYh6dvIcUX6pJ5JoHcs60pzL1sLVYobPR+4CRD\nhhwiMrIc4Kbv6+nTDv/4unXr0bp1a0QRdu/eg16vJygomEceGcDWrXvx8VFRp05dWrVqz4wZ3zJi\nxMMIghNHj6ayZs0u3nqrCZAJQH7+3d3J+KcgfcbXPlFRUahUqns9jTvC0qVLefrpp9m/fz9t2rSp\nKtfr9fTq1YuDBw+yYsUKTp48yfTp069qLwjCxYPnY+/mtCVqEa1Wy6lTp675/PWSh9W6AD937hxv\nvPEGe/fulSwXEneN48ePo9WW0bVra0TRxpo1GzGbrQwc+BggsG3bDiwWO0lJfQAoKSlBoVBUpSdf\ntWoV7ds3ICIClMpyJkx4CACrtQSd7hxKZSMEIaTKv7Y2MBqNlJaW1tohtoKCAmw2G1FRAYCGzMwD\nKJV++PmFAvDcc5d2o+y1Ml5ts3PnTtat+7Hq7wMHDvDQQw/StWvXa7aJjIxkzJiXq/42WS1klxVT\noCmnQKNGZzJes63NqiA7pRmpxzqScbotJsPlQ3Wa0kAupMdfd74uLmZat86mffsMOnRIx99ffxNX\neTVeXl5YLJCdncWpU2coLtaSmiqjQYN4/P1b0KBBfXbsqGDixLcQRTnl5RWEhjoDDbHbRZo1q0ez\nZtcOPSghIXHn+OMOU2VlJb179+bgwYOsXLmS/v37c/LkSQRB4NNPP8XT07Na/VatWt3N6Ur8hah1\nAb5v3z5KS0urWa5sNhu7du1i4cKFnDp1ClEUKSwsrBbLtbCw8IYWxoSEhOs+f6e4ZBW5V+NLVEcU\nRX75ZQv79u1m7NinsFq1aDRmNBoz4eGO8G8DBzbDbrcTHOxwL2jXzgebTSQ6ugSZzI1t27Yhlyt5\n4YWRADRv7k1UlJ3IyCAgAKDKijphwlAAnJxApwvAzc0Lb29vbhWz2Ux6ejpxcbFYLBp++WUj69f/\nyIcfvo5C4cHu3cc5fz6H0aNfAqC8vByz2UxgYCAAu3btQqerICmpG3a7mUWLFlNRoeGVV57Bbrfw\n228HcXKyXIwSITJ+/K3HX76p67BaUGu1lOu0lOkqUGu1qHValE4KWsXUp27IrS8odDot69evB2DA\ngEcRBBmrV69m/foNPP7443h6etXYzmqzkl1cREZBHhkFF8grLbm6jllJRVkgmtJgNCXBlBcHoS4M\noCi7IWbT5S/DJk00dOx4gvbtM3B3jyEtzZ3UVLeqn9nZrkRFVdKzZxH33VdIu3ZlODtfWsyE3sK1\nGvjkk5+YPHk0MpkLwcEG0tKWUKdOIO+/v5XSUi0VFQaGD3+SNm0cr+EDDwyu1kf37t1vejyJGyN9\nxt85jMZrL4L/aVwS3wcOHKgS31fyyCOPVH2eS/wz8PDwuO7nhkajueZztS7A+/fvXy1sFDhO0cfF\nxfHGG28QFxdHcHAwW7durVr5GY1Gdu/ezUcffVTb05H4h1BeXs6mTZt4+OGumM0FBAcX06lTIAZD\nMgDNmoUDlxd0gYE+1drHx0cAjmyNNpuZZ55xhGszGBxbR/fff+NwexZLHitXfo3J5MmkSdNv6Fsr\niiJpaWnUrVsXi6WClJRjvPbadL7++nVE0URcnMiAAQ05fnwz8fHxBAaWIZPZUat3Ipd7s3z51ciH\nNwAAIABJREFUj5SUlPLyy4Ow201UVh7FZDKg0QiASK9e4chkkRgMZwAYNOjS/13tWemvvJaTGans\nOvU7Zdpr+zYnZ2cSHRRC9+YJhPvf/BeNw3pvJyIinGHDhgNw8OABsrKyKSgoqBLgdrudfHUpGfl5\nZBTmkVNceFUqd1GEwqxY0k60J/1kW8qLry2Ovb3P4uHxNRMnevLEEw2rFl0xMT507Fjd71sUueFh\nymtRXq7D3d0NF5dwAgKCCAnJxtOzNTKZDE9PKCgo5p133uHnn38mKCiSIUOG8Nxzz93eYBISEncd\ng8FAnz592L9/f43iW0Lij9S6APf09LzKb9PNzQ1fX18aNGgAwJgxY5g1axbx8fHExsYyY8YMPDw8\nGDhwYG1PR+IOIIoidruN6dOnMmjQEOLjHa/r77//TmxsLO7u7oiiyIoVK+jV637c3Z2w2Qy0aNGB\ntWu/IioqHBDo0KEvq1Z9Tnh4OIKg4IUXxvPOOzMIDg5HJpOxY8cOEhMTkclMVFamsmfPWu67zw1B\nEAgO9iI4uGar6J1k+PAeiCKo1btRKHx56qnXmDTpdRIS2gEO1xY/Pz9sNhPl5Tn07t2T9u29ycw8\ni7+/Ly++OAJRNPHbb3uZPXsOZrMZcCSQmTZtGpGRMdjtWux2LU884TgUabUWAdC6dczFWTgEtq9v\n9a3MO0VphYYNB/eSWZhPpdaL0vwm6Mr90ZX7oSv3R1vuj07tj7OrjsQ+y4EzfPHzT8SHR9KtWSsC\nvX1vOIafnx8Aubm5bNqyGbPdjtZuJTC2LmeLCzhZlEdFpZ5ijRrjxXt2JXa7jPyM+qSfaEfGyQ5o\nyi6PqVDYiYioJDJSz9mz36NSpfHqq31JSLDz+ecvc/r0aYKDx91wjrcrvmUyT/7znw95660ZtG3b\nAkEQGDdufLU6SqWS6dOn1+gnKiEh8ddGr9fTp08f9u3bd13xXVpaWi2aj0wmw8fHp8a6Ev987koK\nrz9aCidMmIDRaGT06NFViXi2bNkixQD/i/PDDz/g6elEQkI0Fkse9es7IQjHKSu7gEzmwocfvsfY\nsc8SF1cPu93IyZPbaNLERkiIDyCyZcu7uLiUUVnpcAtZufI1vLzyqax0CMw+fWKAY6jVychkKubN\ne4c6dcbj7m5FqRSZPn3ovbv4KxAEsNs1mM0apkzpj79/Pmr1DhQKP7p1e4zVq98nIEDBuXNnCA7O\n5fz5bADy8wuYOXMmY8eOY968uVgsVmJjYyktLSEjI5OPP/6Y6dPfusdXdxmrzcqe0yfYe/o4ZrOM\nQ1uf5Mi2R7HbnK7ZZs38mcS32kGHh5aSQjYpudk0ia5HdFAI7i6uuLu44OHiiquzCplMhtFsIq+s\nlPzSYjoOeAiN0cCh4gsA1G3vsOj/nlnzaUaTwZWclGbkn29HxukWaMovZ0gMCTHQt28+ffvm06aN\nGoXCsWiZNm0JR48e5ccfF7Nrlxc5ObmoVCoaN25cW7cNgB9+2Idc7sPAgU+hVPqxefOv0pkYCYl/\nIKIoMmLECPLy8q4rvkVRpFGjRtXK/P39KSoquhvTlPgLclcE+Pbt268qmzp1KlOnTr0bw0v8CfR6\nPUqlgNmcjyhmIIoCZrPD7zUpyeH3JIpGbDYjc+c6tsxNpnQAxo279EHkED+urtVjMvv5Vbfgdu/e\n8uJvFux2C4sXvwxYqtU5fvw4X331ZdWZgREjnql18XQrREQ43Czsdh1ms44dOz68+IydTZs2Y7PZ\n6dGjB6NGPc8XX3zJpk2b+PHHH7FYrLRo0Zzp099CrVYzfPhwjh8/gc1m+0vEO84oyGPjwd8o1WrI\nz4hn+6oXKStwRP1ISCglKspAaKgBhSKb7dv/i92eilr9IHl5r5FypCvpp9rSptdKmnVez8nMNE5m\nplXrXxAEXJTOVF55SFIu4OxWc3ZBUQSbRUl5SQgFqW25cK4tacnR2GyX71VUlL5KdLdsWU5NenfM\nmDG8/fbbnD9/Hq1Wi5ubK6+++mqVBf52KS2tIDU1n44du+PkFExCQjjOzipUqoA/1a+ExL+NOxk1\nsxbP0FejqKgIlUp13Tj6giCwevXqaueHlErlnZmQxN+CuyLAJf6eHDiwn9dfH8eqVa8Ddjp1irlh\nm+thNptwclLedlzilJSzTJv2JjabYwGQmprG1KlT+eCD96lX78/N7U5QWemIitGoUUOcnR1W1k2b\nNmGxOFwodDo9VqsFjcYRtk4mkyGT1e63jyjaKSkpQRAE/Pz8b3jvKyr1bDt2iFOZaZiNKvZtGMmJ\nPX1AlBEdrWHu3FO0b192sW+RUaNG4eqaT6NGjQgM3MmmTcvIyppDWVk/9q57mjP776Nz/yWExZxC\nrricel0UxWri22pWUpRbj/yM+hRm1cdQEYTV7ILFpMJkdMZoUGK3V1fUcrmdxMRS7ruviB49imjU\nqOKGX97e3t588MH7ZGRkYjIZiY6OxtX11nfe7HY7OTlFREVF4OTkh1qt4LffTtGnj8PFpHnzmz+U\nKSEh8fdFEAQWLlzI+PHj6dWrF7t27apyt/0jHTt2lA5hSlQhCXCJqxBFEbNZQ2ysiSVLXuLPhq07\nfvw4n376CQUFhbi7uzFo0CD69n3wlvtZv349Npudbt26MmzYcL788gt27drNhg0beOmll2/cwV2m\nfv367N9/gC+//JKUlBT2798PQIsWLSkrK+P8+fMMHToUg8GRnrxLly4IQu25KRQUFDBr1kwyMjIB\niI+P47XXXsPPz/+qulablX3Jp9hz+jgWq5XMMy3ZsXoUWnUgMrmd0S+e49VXz6NSXX4v6PV68vPz\nUSgUzJjxNnK5Ao3mTVSqh+nefQnLlz9ORkYE6xZOQxDsePpo8PIrxt27EFfvAty8ytAUh1KS04j8\n7Ehs1ut/HDk72/DzM9O5czH33VdE167FeHlZr9umJmQyOfXq1buqfNeuE7i7exMYGIdKFcWAAc/w\n8cfj8PdXAVYyM/OJjg5BJnNHq1Xw/PPT2LdvPwqFksREgcTE3rc8FwkJiercKSv1nSQ+Pp4tW7bQ\ntWtXkpKS2L17N9HR0fd6WhJ/cSQB/ic4c+YUb7/9NsuXr/zbZ5u7xE8//URRUSoPP1wfUbTi7l6z\nS8DNkpWVxVtvTcdicQglnU7Pf/+7GDc3d7p163ZLfWm1OgBat26Dn58frVolsGvXbrRa7Z+a452i\nX79+nDp1isOHj/Dzz1sAaNq0KUOGDKZjx468//77FBYWAtCpU0dGjXq+1sa22ay8/fZb5OTk4uLi\ngijaSUk5xzvvzOSjjz6sEvqiKHIuN5stRw+g1mkpyqnL/k2DyDrjcC9q3LiMefNO1pjVUaVyRi6X\nYbVa+f3344SEhJCVlQVAly65vPDCTj77rC7ffBNFbq4LmjIfNGU+QNxVfQmCSKNGGhIS1CQkqImL\n0+HubsXNzYq7uxW5vJLffz9IRUUF9erVJT6+fq3dKxBwcgomJeUIERFe+Pr64O4exzvvfEy9ek2R\nyQSsVh3PPNORbds24uMTga+vwOHDx2pxDhISEn9nmjVrxvr160lKSqJnz57s3r37b5e8TeLuIgnw\n28RiMRAYqGbIkHaYzWUolb5/exFutRoJCjLj7i5DFG/dslgTv/66HYvFSocO7Rk//lV+/PFHvvzy\nSzZv3nTLAjwmph5Hjx5lyZLPOXfuHDt37gD4S7qfAMjlCiZPnszhw0coKMgnMDCI1q1bI5fLiY2N\n5b//XURZWRnOzirc3d1rdezMzExycnLx9vZm4cKF2GxWnn32WdLS0zmfkYG7lxc6QyUHz50hPf8C\nZQXh7N/0AmnH2wPg6mrm1VdTGTUqveoA4x9RKJx44IEH+PHHn6pF76hbtw5NmjRGobDzyiupvPJK\nKhaLQH6+ipwcV3JyXMjJceXCBRUREQYSEtS0bFmOh0fN7zm1uozXX59MTk5uVdmAAY9WhSu8XY4e\nPc/u3elMnvw2SqU3kyYlVMuEeGWCDLncm+PHr53tTEJCQqJDhw58//339OvXj6SkJHbu3ClFOZG4\nJpIAv0UKCwvZvPlHHnwwFtDStm09dLpDuLg0Jj/fVOPW9l8ZURSZO3cuAwbcj6trAXXrKoG6tda/\n0ehI+x0dHY1cLqdu3boXy289OcNjjz3G8eMnSElJ4YcffgAc/tWPPFK78VZTU8/zySefkpmZiYeH\nB4MHD6ZXr1631ZdMJq+WovhKBEFWoztIbWCxWnHz9cG/bhRr9u1ErdMS36sbcicnVu7bUVVPUxLE\ngc0vc+5IZ0RRjlJp4ZlnsnjppTT8/a8O9/dHnn56BM7OKrZt24rZbKZp06a88MJ/UCiqR0pxchKJ\njDQQGWm45WtZuHAROTm5BAcHUb9+ffbs2cP//vc9zZo1o1mz5tdsV16uY9++M/Tp0xW53JWUlDze\neGMO69YtRRCcadiwES4uF3B2lr4gJSQkbo8/Zkfu1asX33zzDQMHDqR3795s27btHs1M4q+OJMBv\nAVEUqajIIzPzIHZ78BXlFs6d28JLLy1i+/a9KBRX39bKykqUSgVmcwmnTuXQrFlznJ2dr6p3N9i3\nbx9BQQGEhLhisZRit2dRXn4AlerPRYKoiUaNGrJx40a+++47srOzOXny1MXyRjdoeTXOzipmzpzJ\nnj17KCwsJCQkmI4dOyKX197buLCwkMmTJ1NZ6RCKGo2Gzz77DGdnJd261X72QatVQCYTa4zYcbOI\noojBZEKt15JdVEBGQR7ZRYU06NkZ0S5w5DhoSuPQV/ig1/ii1/igr/BFX+FDSW5d7HYFMrmNwYNT\nmTA+g5CQm18cyeVyhg4dytChdy5E5NmzjmRLU6ZMISIiEnd3D9avX09ycnI1AV5ZaeJ//9vFiBED\nkct9UCqtHD26nyFDHK9b8+Zmvv66PW5uIQC4ukJ4eO0tNiUkJP591LTz/dhjj1FRUcFzzz1Hv379\nrmmEkfh3I4h/XL79xbgyjaeX191PvAKONMV2u51GjQIvZh201VhPFEWUyig8PBohkymqylatWsr6\n9T8wd+5/EEUT48cv5sUXx5GQ0Pmuuq3Y7TaMxlzmz/+ERo2CSUy82he3thFFkUWLFrFx48aqsgYN\nGjBt2jRcXFxqdazk5Gx+/VWHu7sTPXrUITRUccshrVauXMHy5SuIi+tBw4YjKS3dwO7d3xAbG1tr\nmVqtVoFffw3g2+VBbN0ahs0m4Otfga+/Gt+Acnz8yy4+1AQGqfEPLkflLCKTyZDLZIiiiM5gQGuo\nrHrY7XZEEXTl/hRmx1CUHUNhdhxFOfUwG68d5UMQbASFriSu3vt88cV4vLy8r1n3XjFq1PPk5eUz\nevSLdO7cmenTp3P69BlGjHiG7Gw7L730DCpVAILgweTJ7/HBBx/WuAi+GaSU5P9OpNf9zmE0GlGp\nVPd6GhISd4Qbvb+vp2ElAX4TLFmyhP/9bxlffz3+xpWBykoXZs5cyUcfTcJszsNoLEEmE/4Q31mO\nShWDXB6OUul8R2M/2+12fvhhDd2718FqzbultkajAaVSiUz25+Z39uxZcnKy8fX1o3nzZn/aam23\nw/nz7hw96s2xYz5s3w7Z2WGI4mXXB2dnKyEhJkJDDYSGGgkNNRAefvkREVGJu7sNUYTsbBf27/fj\nq69KOHUqCqPRsTgRBDve3tuIjf2B//3vAVxcbi8ijCjCiRNefPddKKu/D0Fdeguh7wQ7bp5lePoW\n4elXhJunGotJhbHSDZPBDZPBHZPBDaPOE2Pl1dkxfXzLiYouoV4dgYgIG9u2LcFsPs/w4V3o1SuK\n+fMnkpJyjhdeGEXv3n1u6/ruJD/+uI7PP1+CIAgUFzsTFhaIm5s/CxYs45tv1vGf/7xYa58NkhD7\ndyK97ncOSYBL/JP5MwJcckG5CZo1iycgoOtN15fJNLRu7Y9efwqZTIaTU0232YbRmMLnn3+G0ejO\nm2++fces4RUVZaxb9zUJCUOvSoZzLdLSUpk9ezY5Obk4OSno2/dBhg8fdttCvH79+tSv/+cjV2g0\nCubOjWHp0mgqKv6YkdGGX2AmdtGGThOAyeRJZqaCzMxri10vLzPOznaKiqr/AzkpDfiGZFOcWxe1\nOomDB5No0MBEv375DBx4gbZty27Kuq7RKFi5MoKlS8M5f/6yddk7MJf6CTuIb7ULV49ydBo/tOoA\ntGUBjp8XHxWlQWjL/dFrHI/8jIbXHc/ZRUd8gwLatzPSrk0lLVuWExxsqlbn1KlP0Osr6devD0FB\ncqKiokhJOUdlZeWNL6gG1q3bT2RkFK1a1cdu19y4wS2wc+dxmjbtxqhRUSxduo7CwnOoVF7MmbOE\niIhGTJp0665MEhISEhIS9xpJgN8EcrmcoKCbP6ilUinp16/DTdV9+umumExQWZmFi0skOp0OT8+r\nrZi3i9msxW4/zUcfjbzpNmp1GVOnTkWr1SGTybBYrKxduxalUsngwYNvax5Wq8CBA754elqIi9Ph\n7HxrlmSTScaXX0Yze3YsarUje1hYWCUNGxciuP6GKvA4IXWzUTpf9l82G1XoNH7oyv0wVARhN4Ri\n1oWj1wRSVuxNfp47Go2jL08vE5GxKXiFHSGs3mn8wzKQy20Y9B6cP9qR5EPdKMqOY/nyaJYvjyYo\nWMv9SaX07FlEx44luLtXd0s6ftyLL76I5vs1oZiMjn8zlZuGuJa7qd96B4ERqchkAjGh4UQFNkch\nlyETZMjlZuRCHnJZIQiO3QuzWaSwyI28XE/y8zzRqF3x8hQJ9Bfx9xXZtGEpmRmHUSjKcXIqwM3N\nhRHD36VOnTo13ssGDRpw+PAR3nxzKjExsezZs+di+fXFvQMZ33yzCxcXL4YOHYRc7klIiIyAgDC8\nvZtiMOQxZcoknn++11WZTm+G8+dzsdmgceOWODn5U1ycSXh4DMOGDWbYsCn3NFOoVqvFxcXltt1b\nJCQkJCQkLiG5oNwEx47tw2w+cMcjnJhMSjp0GMmJE0fx8Lhx1sLrodVqefrpocye/QwuLnbWrl3L\nzz9vwWQy0rhxY5577nk8PT05ePAgS5d+RXFxMaGhoTz77LNkZWWzcOFCGjRowIwZMzh27BgzZszA\nx8ebpUuX3fJczpzxYMyYZhw75ljEyOV2YmN1NGigpWHDCho2rKBePT2hoYarXDzsdli7NoyZM+uT\nne2ISd42sYhHh2/AqNpJ8cUskreDq7ML3so66HVKym3HEWSX/xUEQGUDFHIMokNclxWEk3yoGymH\nu6LXXD6w6qS00aF9GffdV4ibm5Vly6KqrhUgIu53GnfYTJ3Gh5DLbXi6utGiXhzN68Xh5fbnwg/+\n8ssvzJ07F29vb5544gkOHjzAsWO/07RpU2bMmFFjm5KSYiZPnkxeXn5V2aBBg3jyySevOU5JiYaQ\nkHhcXOLIzi5GpVIRHh5+VT273c6SJUt48sleWK3Z2Gxl153/uXM5FBfr6Nq1G3K5N2vWbEehcGfg\nwEH3JKxnTa4IR44cYciQIZw9exalUsno0aN57733JCH+D0JyQblzSC4oEv9kJBeUfwjOzmb27ZuH\n2XwQjSaAjAw9P/64k6lT37wpMSKKIhaLBScnJxQKHUOGtMPFRWTFihWsXLmqqt6uXbvJzc1l+PDh\nzJz5Dna7Q3imp2cwdepU+vRx+AGHhITg5OREZGQEACaT6epBr4PZLPDxx7F8/HEsFouMoCAjnp4W\n0tLcOXvWk7NnPVm7NqxaG19f80WfbQNhYUaOHPHmxAmH60ZcnIZ+Q37E6r2OHJMV/jAdk76Ssqxc\nzAYjzu6u1I2LxcPXB7VWi9lquWp+lSYDlaYzAFxKQCkIAk2i69G5SQvKLibJ8QsO5mRGKic90vAN\n/pr2D3xDUW5dss60IjO5FYXZsezYEcCOHQFVfTu76GjQ5hcad/gZn0CH331cWCQtY+OJCQlH9mfC\nnlzBhQuO2Nj335/EAw88QIsWzRk16gVyc3Ov2cbfP4CPP/6YQ4cOo9NpiYmJJTY29pr1i4rMDBz4\nLkePHkculxMTc+2FsEwm49lnnwXAbg9jw4ZV/PLLZqZPfxmQ89NPv7Jlyx4WLHgPQVBgNMopLy/C\n0zMBQRAYNuza87gX5Ofnc//991NaWopcLsdsNjN79mxcXFyuucCRkJCQkJC4EbUuwGfNmsXatWtJ\nSUnB2dmZxMREZs2adVXYuWnTprF48WLUajVt27bl008/pWHDm9kC/2fj8Be3Y7UWolSqadRISUXF\n78jlrmzYsJNDh04xa9ZMZDIVp0+foaCggG7d2mOzGfjkk0/JycnhrbdexGLJp2PHhoiiw/oNMG7c\nOKKjo5g2bRrp6RksX74Cu12kd+/eDB48iAULFrJ3716Ki4sBRxIdrVZLZmYGAE2aNLnp6zh2zIuX\nX25OcrLDDeGppzJ48eUDqNyMWMxKMtJ8OJ/iw7kUH84m+5CT40ZBgQtlZUrKypScOnVZ5AUFGeg3\neDOu0Ssw2s1whZHcSaGgQUQ0oR7erF72Dfnnz+HsrKJr3wcYPGgwMpm8KkxfvrqU3JIicosLyS0p\nwmSpLsobR9ejc5Pm+Hs6BP8lAe7j7kHnJi3o1Lg5eWUlnEhP5ZQqjaDINNr0+g6DzpOssy3IPN0K\nY6UncS13EdtiD05KMzJBoEmdWDo0alrVb23i7+8Q/Xv27KFRo0bs3u1wJwkIuH58cZXKhU6dOl3z\neYvFisUi4OvbkPj4Ohw50u+WXT9kMhn33/8oUVGN8PJqCsCAAfH07/8sbm4Ov/zu3aNuqc+7zbp1\n6ygtLaV9+/Zs376dbdu20bdvXz7//HNJgEtISEhI3Da1LsB37drF6NGjSUhIQBRFpkyZwn333Udy\ncjLe3g4B8t577zFnzhyWLl1KXFwc06dPp2fPnpw7d67qi1kCgoJ8CArywWLJxWKBtm29aNCgBRrN\nDkBBXl4qubl5aDQWwM7Qoc2RyVpisVyo6sNsNmMymZHJBDp0aI9C4URsbCylpQfQ6x2p3Vu2bIGn\npxfNmjVj7969yOVynnnmGZYsWcL27WA09sTDIw6z+RFefdWT0lKHSBZF8Pc3ExBgIiDAVPX7gQO+\nLFxYF7tdoE4dHf8Zv55Kl00s2/mHA3pKcGsCrZpAa0HAy9UThT0UuyEcsy4YvSaASms5LhGrEGV6\nbFcI7yBvXxIbNKZBRDRKJ8dhzDYfNMNmsyGTyartGAiCgKtKRb2QMOqFOCzuoihSrFGTW1KE0Wwm\nNiyCAK/r+/kLgkCYXwBhfgH/3959x1Vdfw8cf33u5bLhiiBTBAQHzhS/qLk1c+TMlZZmpUmpOFJL\ny1Fuzfy5V1lqOXJlZpqm4sIcgQsXOUHFzZR57/39QdxEUFGBq3iej4ePup/7vp977n1f5dwP533e\nvF49kMirURw5H8k/V6IoX2MX5WvsMo5Vq9RU9/Ontn9litnaPfnk51GjRo3YtOk3oqKiGTVqNJCZ\n+D5trX6WadM24u1dhb592z5TKYi5uTlVqlQx3ra0tOT27dsMGjSIw4cP4+joyCeffPLUGx0VtKyF\nqaVLl8bCwgJ/f/9sx4UQQoinUeA14ElJSWi1WjZs2MAbb7wBgLu7O8HBwXz22WdAZg2Ns7Mz06ZN\nM/76OsvLVANeUD76KIgrV65StWoVPDxKsmXLFvR6PY0bN2bHjh04Ohandu1X2bVrFwkJCXzwwftU\nr96VgQN92b0794V8j6NSGWjdaT9+dRaQlo+dMZyLOdCgcnXKe3oVeI3wP/9EAuDn9+iyiKSUZI5f\nPMfxC/8Qfy+JKj5lqO1fCVsr6wKNL0tCQjzLli3jwoULaLVa2rVrT6VKlZ7qXIpijqWlPxqNKxqN\nJt/f46SkJAIDAzl58mS24+vXr6ddu3b5+lxP48Fa4MOHDxMYGIjBYKBly5acOHGCy5cv07p1a379\n9VdThirykdSAFxypARdF2XNdAx4fH49er8fBIfPq4oULF4iJiaFp06bGMZaWltSvX5/Q0NAcCbh4\ndgMGDGTMmNEcPXqMo0ePAfDOO2/Tpk0brlyJ5syZs/z2228AVK0ayMWLgwgOLs+9e2bY2KZTs/ZF\nLG3jMLe+i9r8NmhuodfcwGAAfWoJ9KlOpN9zJCWpOMkJ9mQYMihfZwmOJc+Qdt9Va3MzDQ62duj0\nOnR6/X1/dKSkPXrbc2etAw2qVKO8p7dJFuc9io2lFbXKV6JW+adLenOj1+vQ6/U5tnTPjZ2dPR9/\n3PeZns9gMNCt2yTmzfsWP79Sz3SuR1m5ciUnT57Ez8+PRYsW8csvvzBjxgxGjRr1XCTgD6pRowbT\np09n8ODBxs2kKlWqxIIFC0wcmRBCiBdZgSfgAwYMoHr16tSuXRuAmJgYFEXBxcUl2zgXFxeuXn2y\nTWJE3pQvX55Zs2axd+8+YxeUSpUy67knTJjI/v2h3Lhxk7i4WvzwQzdOn878lhZY7ygVX/s/bLV3\nn+n57aysqVm+EtX9ymFpbp7rmLSMdO4mJHA3MZ67iQncSYjnTkI8apWKV0qXxb/U0yXe6enpLF26\nhN27d6PT6QgICODDDz/E5t/uIzdv3uDu3Vjc3d2wLcBSkbxKTU1h3rz57N69C51OR+XKlRk4cKCx\n1js9PZ3ExAS0Wu0zb470HzOsrcvxzTcL8PX1z6dz5i5rcWjbtm1p2LAhvr6+zJgxgytXrjzmkaYz\nYMAAmjVrRlhYGMWKFaNx48ZyRU8IIcQzKdAEfPDgwYSGhrJv3758uWqZ9WvCwpa19iyrJOFFVbny\nf1do738tdnY+LF3agw0bMmt13d3v0qTztxTz3vfI8xn0BhTVw+e1mLUN/q4lKeVYArVKRfTlS4+N\n0QwoobGkRHFLKO6ceTA9g3Pn/nnsY3OzYsVKDh06ZLy9c2cIly9fpk+fPvz882rjZ8rMzIw332xP\nrVq1cpyjMOd92bJlhIcfMd4+duw4I0aMYNCgQYSEhLBt259kZGRgY2NNp06dstVXPyi1Kp0mAAAg\nAElEQVQ1NZ2zZ6/yyivlURRzYmPjWbs2hPffbw7oiY9PYsOGMN5/P5iMjMx2gX///XeBvj4Li8yN\noL799lsyMjI4cOAAAJ6enib7+52b3GIpWzZzd9QTJ04UdjiikDxPn8GiwsvLS76wiiIrISHhkT8T\nHtVhrMAS8EGDBvHzzz8TEhKCl9d/nQ5cXV0xGAxcv349Wx/h69ev4+rqWlDhiIc4etSdMWNaEhNj\nj5mZjg6dQ/EIXEA6icYxTrb2OFjbYG1uwaZffiX25i1cHZ0wU6u5fPUKtsXs6fL226Qb9CSlpWIw\nGPB0cMRV6/BMX7xiY2NJTU3F0dHxqXouJyXd49ChQ6hUCkFBQVhbWzNnzlzOnTvP2rXrOHz4MCqV\nihIlSnD9+nVWr16Ns7MLpUs/Xd37s0pJSSE8/AiKktmxxt7enokTv+batev8/vvv7Nq1GwBLSwuS\nku6xdOkygoP7U6pUZslIRoaOP/44Qtu2r6FWFyM9HX78cSPVqnVHURQ0mkScnO6gVtcEQFFisba+\nQ3q6Jk+7euaHxo0bU69ePfbs2cOMGTMAsLGxYejQoYUTgBBCCPEcKJAEfMCAAaxevZqQkJAc2b+P\njw+urq5s27aNgIAAIDPx2LNnD9OmTXvkeU21QCY8fD/w+MV4heX8+XMsWvQt0dFRODk58c477xAQ\n8GTvTUaGwjfflGHatLLo9QrVqt1l4IjfOBqzgXRdhnHca9UCqe1fCUVRiIq6zHfnLlC8uAP/N/0b\nFEVFv359iYqKxtvJGX//ChgMBnS6jDzVLj/MvXv3mDZtmvHKtYNDMYYNG0bFik9WYx0Tk7nRjJ2d\nHc2aNUNRVPz666+cPRvJxYsXARgxYjiBgTWZN28umzdv4erVK7z++utA3hdh5pfY2DiuXLGldGl3\nmjTph4WFDeHhX1O2rDfXrinY23tjbe3O118PYtWqNezadYBBg37iyJHt2NhoURRLfvwxksqVuxmv\nNO/a1T7bczRo0CDb7caNGxfKa7vfjh07+Pbbb41dUHr37o2fn1+hx5EbWYz3cpJ5LzgpKSmPHyTE\nC8rOzu6R/27cvwjzQfmegPft25cff/yRDRs2oNVquf5vL2VbW1tji8GBAwcyceJEypUrR5kyZRg3\nbhx2dnZ07do1v8Mpcq5du8rw4SNITk4GIC4unrFjx/Lll19RtWrVPJ0jKsqKoKDqHDxYHEUxMGDA\nWWq1XEXoqf/KDyw0Gt6s04gyHp7GY5aWVgAkJiYSExODRmPO3bt3/73PknXr1rJ69RqSkpLw9PRk\nwIBgypYt99A44uPjmDNnDuHhRzAzU9OoUSN69nyP+fPncejQITQaM2xtbbl7N5Zx48Yxd+5cHByK\n5/m9KlGiBFqtlri4OL788itsbW05ezYSMzMzzM0zvyBkNQHK+m9hL/CcMGEVLVq0oG7dxtjb26NW\njyIk5CwNG7anRIkS3L2bwLFjGdy7Z8fJk1EEBb2Bu3srdLrTbNu2ms6dO1O8eEVjj+758xcWavxP\nw8zMjKCgIFOHIYQQzywpKYkpU6Zw+PBhDh06xK1bt5g0aRLDhg0zdWjiOZfvCfi8efNQFIUmTZpk\nOz569GhGjRoFwLBhw0hJSaFfv37GjXi2bt0qPcDzYPPmLSQnJ1O9ejU++uhj1q1by+bNW/jll/V5\nSsDXr3dnyJAqxMdrcHVNZsjITei0Wwk9dcM4pridPW81aIqTNvvGMU5OTgQEVOfvv8MICvrIeLx8\n+XJERETwww9LjMeioqIYOXIUM2fOzLHgFkCn0/Hll18RGflfffXGjb9x714ye/dm1p7PmDEDd3d3\nRowYwcmTpwgPD6dx4yY5zvUwarUZn3wymPHjxxMWFgZk9sju27cvN25cZ8WKlUyaNAk3Nzeio6+g\nKBAYGJjn8z8tg8GASmWNpaUv77zjhpeXD1ZWmVvbr1u3jubNm3PkSGYduIWFBUuXLuXw4cMcP36c\nZcuWERkZyfbt29Hr9XTu3PmJN8gRQgiRP27dusXYsWPx9PSkevXqbNu2zdQhiRdEvifger3+8YOA\nUaNGGRNykXfx8fFA5q9KXVxcePXVV9m8eQsJCQmPfJxeD8OHV2Lx4sz65oDap3j1za+5lH4bbv03\nrrSrOx3qNsbq3xKG+ymKwtChQ5k7dy6hoaEYDAYCAgIIDg5mxIjPAfj4449p0qQxY8eO48iRI+zZ\ns5uOHTvlONeZM2eIjIxEq9UyZcpk7t69y2efDWf79u2o1ZnbtFtYWKAoKsz/7Zyi1z95y/pXXqnG\nrFmzOXjwADqdnldeeQUfHx90ugxu3LjJ9u3biY6+gkaTeVXW379gd2M9dSqayZPX8ssvv2NmZk71\n6tnrzV955RVOnjzJtm3bSEtLo27duvj4+NCqVStOnDjBpk2b2Lp1KwDDhw+nffv2uT2NEEKIQuDu\n7s7Vq1dxdXXl0qVL+PiYZg2RePEUeBtC8Z+bN29y8eJF7OzsKFu2zFO1kfP19WXHjh2sXLmS27fv\nEBoaCmTu1Pcwej0MHVqFpUu9MNOkU7ftd1SqswXuq7ZQKQo1y1eiySs1UKlUDz2XtbUNQ4YMRa/X\nYTBgvPqakpJZElOmjB8ajTne3l4cOXKE5OTc6/+SkpKAzH+83NzccXZ2wdrainv3kqlatSphYeEM\nGDAAB4fiREVF/bujYuUneq+yuLq60qZN22zH1GozBgwYwNtvdzO2IcxqTZjfDAYDimKGubkX//tf\nQ2bObIaZWe7tGAGKFy9Oly5dsh2ztLTk119/5eDBg1y9ehV/f3/jroxCCCFMQ6PRSAMJ8VQkAS8k\nmzf/zoIFC42/IahYsSIjR36BtfWTld20bNmCv/8+TFhYOGvXrgUyW7h179491/EGA3z2WSWWLvVC\nrUmlVe/xeJY9Zrxfa2NLdb9yvOJbFrsn2LnxwS8PFSpU5MaNEMaPn4C/f3n279//7/HcryiXLl0a\nMzMzTp06xeTJk4iNjePevWScnJwYOHAQEyaM5/TpMyQmJmFlZcWQIUNwds5ZyvKg3buPYWtrRfXq\neVs46eRUwthj+0HHjl3Ay+u/5zx+/DylSjmj1eZM1NPSUjl27BhJSfcoU8YPd3cP432DBy+ia9cP\naNGiAoqiULr00/UbV6lUubZJFEIIIcSLRRLwQnD27Fnmz5+PwQAVK1bg8uXLREREsHDhIgYOHPhE\n51KrzRg5chR79+4lOjoaJydHGjRoYFwgeT+DAYYP9+f7731Qm6XR6oOJeJY9hqIolHH3JKBMeXzd\nPB55xTuvevXqxeXLlzh//gJ79uwFoEOHN42dbh7k6OhI3759mTVrFvv2ZV7Ft7GxZujQIRQrVozJ\nkydz8eJFkpOT8fLyysPVaQWNxgOtVk18fALW1lXJyIhn3br1ODtbERj48MWguTEzc2b//ts4OtbE\nzq4uYGD16l/p2fNt3Nz8yci4y08/fccbb9QkOTmRUaNGEhWVucmMSqXw4Yd9aNmyNZaWZRg/fjYe\nHiWfux08hRBCFJ74+HiGDRvGqlWrsnXHUBSFGTNm0K9fPxNGJwqbJOCF4OjRIxgM8NprTQgOHsC5\nc+cYNGgQR46EP9X51Go1DRo0IC1NQaMx5NrD2WCAT4eX4fvv/FCp02n5/iRKlT9CWY9StPhfbbT5\nXG5hb2/P1KlTCQsLIy4uDm9vH+PGJQ/TpEkT/Pz8OH78GGZmGv73vxo4OjoBoCgqfHweXlaTJTY2\nic8++54lS5ZhbV2CJk2qZbvf3f0qTk4O2NiUIiPjLj/88D2Bgb74+rrlONfZs1GcPHmTt98OwtLS\nmb59BwNgYeEAwJgx4/nggw/4448/UKlU+PuXo02bTixcOISoqGicnZ3x9CzJwYNhjBz5IxUqdCUg\nwA9bW0m8hRAiL8aMGVOg/zWVhIQE6tevT+PGjdm1axfh4eH06dOH/fv34+TkhKen5+NPIooUScAL\ngVqd+TbHxsZiMOi5e/fOv8efrnuFTgcTJ5Zn7lxfihdPo169W9Ste4v69W/h6ZmMwQBDP/NiyeLy\n/ybfk/GuEEZt/8qPrfF+FhqNOTVrPlmJhJeXV7aNmp6EmZkTHh6BfPhhcaysnHK9wpzV0zuTC1pt\nRRwdq2NpaUF6+nV27txFYGB5bGxcsbd3xNLyPFZWuXdtad26NQcPHjQeO3LkGMHBkzl+/DL29qWo\nUaMpQUHtmD37F6ZOXc7evWHUqFH/qV6bEEK8jB5MlPP7tqkMGzaM6tWr88033wBQpUoVZs2aRVRU\nFNWqVXvMo0VRJAl4IahduxYrVizn8OG/ef+jIDIMBoqVdKNG7dqEnztLhi4DnU6PjaUlpd08sMml\nnCTLnTsaPuxTjV0hmUni9euWrFlTkjVrMncV9SgZh7PbHcIP+aBSZdCi51R8K4fRKrAer/g++or0\ni2Lu3E24upanV683UKlUDyTZj9ajRw/j/+t0pfnpp6+pX78rWq03Dg4qKlV6NdfHhYWFcfDgQZyc\nnAgLCyM2NpaAgADWrVuHg4MDkZF3eecdfywsAti6dQL37iVjZ/d0td5CCCGKjjt37rB48WJOnz6d\n7Xh6ejoZGRkPeZQo6iQBf0p6vY5Lly6RlpaGt7c3FhaWDx3r5ubOkE8/ZX3oLqwdHYzHb6Jj4197\ncoz3cCxBGQ9Pynh44urgiKIoJKemsiUkheGDW3D7hiNWtrE0f/drrO1iiTpblejIylz5pxJXorVc\nidaiqHQ0e/drKtY4Sud6zfFyyVly8eJRMDf3oWvXQTg7uz7zlXy1Ws3atb/kaWxW+8eSJUtSsmRJ\n3Nzc0Gq13Lp1iw4dOvDtt98yaNAQhg79jIyMDEqUKEHr1q2fKT4hhBAvvn379uHh4ZGtRWFiYiKn\nT5+W3VdfYpKAP4U7d24zduxYzp07D4C9vR3Dhw9/6FbpF2KusvdSJFbFHbgZ7UNaijU22tvYaO+g\nMU/LMf7K7ZtcuX2TkGNh2FlZ42Brx59byrN95cdkpFvg7BlJy/cnY+eQ2cC7uGs0VetvQq9XceuK\nD1fOVaCExwVeCYima8M2ONjZF9ybkUd6vZ5bt+Jxds7cyfL8+SskJiZTpUpmnXdycioajRlmZjnL\nctLS0hk2bDH/93+zsbPzwt6+8GuqX3nlFWxsbDhy5AjNmzcnISGBW7du4eHhwYwZM7C2tmbevHmk\np6dToUIFfvrpJ5ydnQs9TiGEEM+X1NRU3N3dsx1btmwZzZo1o1SpUiaKSpiaJOBP4euvv+bcufPY\n2dlha2vDtWsxjB8/nvnzF2Bv/1+yq9fr2X3iCLuPZy62PLjlLQ7+8Va2c1nb3KOYYwLFHeMxt7kF\n5texso3FyjYWa7s4rO1iOXWwMUd3Z15N9Q/8k4adFlDCwZKyJStioTEnNT2N1PT0zP96pJBW7QAe\nTiVoULk1luYP7zddmPbu/Yd16w7z/feLAT0xMVu4du0qderUA3T89NNCzp//h6++6oden0h4+FEU\nRUeVKuXQaivQsWMvtFoPk3UScXR0ZNmyZXTr1s24EY6DgwOrVq3C2tqaGTNmMG3aNJKTpfRECCFe\nJnPmzCE2Npa7d+8CsGPHDtLT0wEIDg6mXr16fPrpp+j1elQqFZcvX2b+/Pn89ttvpgxbmJhiMBie\nfHvBQnR/qx6tVmuSGMLD95OWdgBfX1/u3bvHW2+9hUqlYv6ihSTr0pkzYyaRp88w4rPPjIsQE5Lv\nsX5fCBevXwMgbEc79v3aE5VKT/Xqsdy4YUlMjAVpaXlbiKlSp9Pq7RW82/MS5T29KKEtliMZTUtL\nZfny5Rw7dgyNxpwmTZrQtGlTk7e/02hKYWPjj6KoH7nw9P6P4tq1a9DrU2jfvj1mZjYmeQ2HDx8G\nyPYrwvPnzxMSEoJaraZZs2ayAUMRk9uci6JP5r3gpKSkYGn58BLNosDHx4fLly/net+FCxcoVaoU\n69evZ9OmTbi7u3Pjxg1GjBghV7+LgMd9vh+Vw8oV8CekUmUmgi7lfFmwZTPpaWbYVypLQKWy7Lhw\nmqO3Y9Da2BB96wZJKZm7QB7b04J9v/ZEUQzMmnWUzp0z+0UbDDBu3CJCQi6RluZBeroL6enOODsH\noNWW5/DfN9DpXTEzS8XLZwj2+svUKjcZc/Oc28QbDHomTJhAWNh/rQ1PnTpFfHw8HTt2fOrXq9fr\nURTlqRLgX34JRa93onfvVnl6/P1jctu+/nlQunTpR+46KoQQ4uVy4cKFx45p37497du3L4RoxIui\nYPrRFWEZBgMBrV7H3LElP06cwbcjl7Dv1x6k3LMhXacj5u5tzkRfNibfJw80ZtfaPgBMnXrcmHwD\n3Lx5nUOHNqLVRjBzZk3mzfPGy2saGk0XnJ3bUMG/KYMGDOTHpevwLhnBuXPn2b0756JNgFOnThMW\nFo6NjQ2jR4+iT5/M51y5ciU6ne6pX+8HH8wgJOQa5uY+aDQezJr1OxER18j87qbin3+iSUxMzvE4\nlcqOmjXb0qBB3pJvIYQQQoiXhUkT8Llz51K6dGmsrKyoUaMGe/fuNWU4j/XP1WgW/L6e0yfeYPX/\nTSb2pge6dAvCdrzJ0nHzCdvRlox0jXH85RON2Lkqc2ersWMjePfdS9nOFxeX2VnD1dWVatWqU61a\ndWNJw82bNwHo2rUbNWvWol69egDcuHE919hu374NgL+/PwEBNWjZsiUWFuakpaWRlJT4VK/XwsKP\nNWs207Hjh9jbV8Levhr16nWgTJlWODg0RattwKJF+4mKssbSshzm5l6MHfszFy6kY29fi8qVa1Gu\n3JPtQCmEEEIIUdSZrARl1apVDBw4kPnz51OnTh3mzJlDixYtOHXqFCVLljRVWDmcPHmSP7dv455l\nBNdT9Oxa25uTf2X2na7b5ADvvhXFDz/UZ98+Z/b9+h7nDnXivT6hWFunMe+Hluj1Kj777DRBQedz\nnNvd3Q1zcw3R0VeYMWMGANHRVzA31+Du7kFkZCTffruIGjVqsH37jn8f457jPAClSmXuohUeHsbK\nlSuIiYkhNTUNrdb+iRcF7tp1lOrVX8fRsXy2q9eKovDaa68Zb6vVtnz33dJsj23b9j18fV/BzKxo\n1/wJIYQQQjwtk10Bnz59Ou+//z7vv/8+5cqVY+bMmbi5uTFv3jxThZTDokWLqFKlCkeu7+L89WKs\nmTGJk3+9jpkmjc+/CmH9yhu0a2fB+vUHWLnyLypWjCPmmi0Tx7zOyGGtyMhQ0b9/JIMHR+Z6fhsb\nW/r27YdKpbB9+3a2b9+OSqWib99+9O7dG3Nzc/bv/4tZs2YTGxtL5cqVjFfCH+Tl5c2bb7ZHp9Oz\nfPkKduzYiVqton///ijKk0yzmjNnkrl9W/VUpSONGzemePHiT/w4UfSEh4fToUMHAgMD6d69O5cu\nXXr8g4QQQoiXgEm6oKSnp2Ntbc3KlSvp0KGD8Xi/fv2IiIhg586dxmOm6oISGRmJv78/Op2O19p8\nxZ5tg0hNtsXe4SorVpwgMCA1x2P0eli71oMJE8oRHW3De+/9w+TJp3hcHhsZGcmhQ5lbnP/vf4GU\nKVMGgEuXLvH775uIj0/A19eXtm3botFoHnoeg8HA/v2hHD16DHNzDQ0bNsTX1y/Pr1lRNFhbV8XS\n0lXqtpHOCM8iPDycOnXqkJz83/oAV1dXwsLCcHN7fjeFkjl/Ocm8F5yXoQuKeHm9cF1Qbt26hU6n\nw8XFJdtxFxcXtm/fboqQcjhw4AA6nY4WLVrQsE5n/vzVFlfPP/F07UjVSt8COTuRKIqBjIyplCr1\nC8WL+3HhQiQ7dwbRuHGTRz5XmTJljEn3/by8vPjoo4/zHLOiKLz6ah1efbVOnh+T5f/+bz116rSi\neXNJvsWzGzduHMnJybRr144BAwbw2WefceDAAWbNmsWECRNMHZ4QQghhUi9UG8KsqxSF4fr1zMWO\nERERfPbZJYYP38rOnQMwM1O4dOlirtug79+/n9Wr1wDg6nqJhIRUZsyYQVpaOn5+voUW+5NQFDXm\n5p5UqFAfjcaev//+29QhPXcK83NXVJw5cwaAZs2aYWtrS6NGjThw4ABHjx59Id7PFyFGkf9k3vOf\nl5eXXAEXRVZCQgInTpx46P25XVzNYpIacCcnJ9RqtTHJzXL9+vXnZmOTwMBASpYsyeXLl+nVqzf7\n9g1BUQzUqlU71+QbICwsDIAOHTrw5ZdjqFevLgYDHDkSnut4U7p1K5ERI1ah0QSg05WibNlKFCtW\nzNRhiSLC29sbgNmzZ7NmzRp+/PFHIPOHsRBCCPGyM8kVcI1GQ0BAANu2bctWA75t2zY6dXr4BiyF\nXZ+3Z88eevbsSUzMRTQaDS1aNKNnz56YmeVeh23+77bvfn6++PmVwdMzc5crW1tb/Pwe/i2osOh0\nOhRFhYWFJ5aWPjg5vU716rVMHdZzS+pCn97cuXOpVasWERERREREAFClShWmTJmCra2tiaN7OJnz\nl5PMe8FJ+XdPDCGKIjs7u0f+u3F/DfiDTFaCMnjwYHr06MH//vc/6tSpw7x587h27ZpxA5nngbe3\nNyEhIRw+vJeMjMOPLSMJCAjg1KnTTJ8+nd9++41Tp04BUL16QGGE+1gjR/5Is2Zv8uabr6AoCrVq\nSfItCkapUqX4+++/mT17NlFRUZQvX57g4ODnOvkWQgghCovJEvDOnTtz584dxo8fz7Vr16hUqRKb\nN2/G09PTVCE9lFqtJi+bSXbo0IHLl6PYvXs3EREnAejSpQuvvvpqAUf4eBpNKUaPno6Li5ssshSF\nws3NjfHjx5s6DCGEEOK5Y9JFmEFBQQQFBZkyhHylVpvxySef0KlTJ27duoWHh4fJa9p/++0vmjfv\ngqNjZbRaSbyFEEIIIUzNpFvRF0WKouDl5UVAQECBJ98qlT2KYs6pU5dITU3PJRYrzp1LJyHBSq56\nCyGEEPloyZIlqFQqVCoV+/bty3WMn58fKpWKxo0bF3J04n7Jycl8+eWX7N6929ShGEkC/gLat+8E\nq1YdpVixuhQr1oQpUzaRmFgCjcYVRTHnxIkLmJk5YW9fi7FjpzyyDY4QQgghnp6VlRXLly/PcfzA\ngQOcP38eKysrE0Ql7nfv3j2+/PJLQkJCTB2KkSTgLxhFscTPryFVqjRApVKjVpuxceNmypWrh719\nDWxtGzBu3HpSU73RaGTBmxBCCFGQWrZsyerVq9E9sFhs+fLl+Pv74+v7fO4Dklf37t0zdQjPzASb\nvj+WJODPidjYRMLDIx/6IYmJuU1SkgX29rWpWLE29erVyzFGURQsLCz588+duLg8v9t9CyGEEEWB\noih07dqVO3fu8McffxiP6/V6Vq1aRbdu3XL9uT5r1iyqVKmClZUVLi4u9OrVi9u3b2cbs3HjRtq0\naYOnpyeWlpZ4e3szbNgwUlNTs427ceMGvXr1olSpUlhaWuLq6krLli2NndgAVCoVX331VY44vL29\nef/99423s8pqdu7cSXBwMK6urtjZ2Rnvj4+PZ/DgwXh5eWFhYYGvry/jxo1Dr9cbx1y6dAmVSsWU\nKVOYN28evr6+2NjY0LRpU6KiogCYMGECpUqVwtramrZt23Lnzp0csW3dupWGDRtiZ2eHnZ0dLVq0\n4OjRo9nG9OzZEysrK65evUq7du2ws7PD2dmZoUOHGt/3S5cu4ezsjKIojBkzxlg2lPW6k5KSGDJk\nCKVLl8bS0hJnZ2caNWrE3r17c8SUn16onTCLLgWttiLTp3/PqlWz0esTSUiIYfPmXbRvXxdQ+OGH\nUGrWfI1OneSqthBCCPG8KFmyJHXr1mX58uW0bNkSyNzX5ObNm3Tr1o2VK1dmGx8UFMT3339Pz549\n6d+/P1FRUcycOZNDhw5x6NAh454i33//PZaWlgwYMACtVstff/3F9OnTiY6Ozlby0qFDByIiIujf\nvz/e3t7cvHmTXbt2cfbsWfz9/R8Z+8PWh/Xv35/ixYvzxRdfGHtZp6Sk0LBhQy5fvsxHH32El5cX\nBw8eZMyYMVy+fJmFCxdmO8eKFStIS0ujf//+3L17l8mTJ9OxY0eaN2/On3/+yaeffsq5c+eYMWMG\ngwcP5ocffjA+dvny5XTv3p3XX3+dSZMmkZqaysKFC6lfvz6HDh2ibNmyxvgNBgPNmzenZs2aTJs2\njT///JNvvvkGPz8/+vTpQ4kSJZg/fz5BQUG8+eabvPnmmwDG30wEBQWxZs0a+vXrR4UKFbh7965x\n5+a6des+8v17FpKAm1B6egYxMYmUL98US0sXfv99K5D5q5KkpGucObMNa+sAFEVhypSWD92BUwgh\nhHhZHDlyhL59+xIREYGrqyujR4+ma9euJo2pW7dufPLJJyQnJ2NlZcVPP/1EzZo18fHxyTYuNDSU\nhQsXsmzZMt5++23j8ebNm1O3bl2WLl1Kr169gMwk1NLS0jimd+/e+Pn5MXLkSKZOnYqHhwdxcXHs\n27ePr7/+msGDBxvHDhs27Jlej729PSEhIdnyjm+++YazZ88SHh5uXFvWq1cvvL29GTlyJEOHDs22\n5uzKlSucO3fOeAU9IyODiRMnkpycTHh4OGq1GsjcBX3lypUsWLAACwsL7t27R//+/Xn//fdZtGiR\n8XwffPABZcuW5auvvjLurgyQnp5Oly5d+PzzzwH48MMPCQgI4LvvvqNPnz5YW1vToUMHgoKCqFKl\nCt26dcv2Wjdt2kTv3r2ZOnXqM71nT0oyOhPav/8y8+fvwsrKNdu3UEVRcHd35+uvp2Nt7Y6VlZsk\n30IIIV56Fy9epGHDhoSGhhIXF8eZM2fo1q0b69atM2lcnTp1Ii0tjV9++YWUlBQ2bNjAO++8k2Pc\nzz//jJ2dHa+//jq3b982/ilbtiwuLi7s3LnTODYr+TYYDMTHx3P79m3q1KmDXq8nLCwMyFwAam5u\nTkhICHfv3s2319O7d+8cecfq1aupW7cuxYsXzxZ7kyZNMBgMORY4duzYMVv5SvMZUNIAABFnSURB\nVM2aNQHo3r27MfnOOp6enm4sT9m6dSuxsbF07do12/Okp6dTr169bO9RlqwvLVnq1avH+fPn8/Ra\ntVotBw4c4OrVq3kan1/kCrhJKFhY+NKuXUvat5fEWgghhMiLxYsXExcXR5MmTVi6dCmzZ89m4sSJ\nTJ8+3VhaYAoODg40a9aMn376CbVaTXJyMl26dMkxLjIykoSEBFxcXHLcpygKN27cMN6OiIhg6NCh\n7Nq1i+Tk5GzjsspCzM3NmTx5MkOHDsXFxYWaNWvSsmVLunfvTsmSJZ/qtSiKQunSpXMcP3v2LMeO\nHaNEiRKPjR3IsbGiVqsFyBFX1vGsLxCRkZnr4V577bVcn+f+5B1Ao9HkeD8dHBzy/IVk6tSp9OzZ\nk1KlSlGtWjWaN29O9+7djWUuBUUS8EL27bdbsLX1pm/fN6Q3txBCCPEEspKqRo0a4e7uTqtWrZg4\ncWKui/gKW7du3ejRowdxcXG89tprODo65hij1+txcnJi1apVuS7OdHBwADIXO2YtQJw4cSK+vr5Y\nWVlx5coV3n333WyLHgcMGEC7du3YsGED27ZtY9y4cUyYMIFNmzZRv379R8b8YOeWLLm1TtTr9TRu\n3JgRI0bkGvuDSfuDifLjjmedU6/XoygKS5Yswd3d/ZHxA89cIdCxY0fq16/Phg0b2Lp1K7NmzWLK\nlCksWbKEt95665nO/SiSgBciRdHQpcvH2NiUkORbCCGEeEI1atQAYMqUKdy4cYPNmzdnO25Kbdu2\nxcLCgtDQUJYsWZLrGF9fX/78809q1qyJtbX1Q8+1c+dO7ty5w/r167MtBPzzzz9zHe/l5UVwcDDB\nwcFcvXqVqlWrMn78eGMC7uDgQGxsbLbHpKenc+3atTy/Pl9fXxISEmjUqFGeH/M0fH19MRgMODk5\n5dsGRo/LuZydnenduze9e/cmPj6emjVrMnr06AJNwKX+oZAYDCqsravh7V0BZ2dnU4cjhBBCvHC6\nd+9Oly5diI+PZ+bMmURGRuLv71/oC+hyY2Vlxfz58xk9ejTt27fPdUyXLl3Q6XS5tgTU6/XGJFmt\nVmMwGLJd6TYYDEybNi1bMpmcnExKSkq287i7u+Ps7Jwt4fb19c2xC+SCBQseegX8YbEfOnTI+KXn\nfomJiaSlpeX5XI/SrFkzihUrxoQJE0hPz7nL961bt574nFlfdh4sS9Hr9cTHx2c7Zm9vj4+PT44v\nLPlNroAXgqNHzzFnzjbWrv1NrnwLIYQQT0mlUrFixQq6du1KREQEbm5udOrUCVtb07TofbAU4/7O\nJrmpV68effv2ZerUqRw9epRmzZphYWFBZGQka9euZezYsfTo0YM6derg6OhIjx496N+/PxqNhjVr\n1pCUlJTtfGfPnqVx48Z06tSJihUrYmFhwaZNmzh9+jTTpk0zjuvVqxdBQUF07NiRpk2bcvToUbZu\n3ZprPffD9iMZOnQoGzdupF27dvTo0YOAgACSk5M5fvw4a9eu5fjx45QqVSqvb91Dn9POzo758+fz\nzjvvUK1aNbp27YqLiwuXL19my5YtVKpUicWLFz/R+S0tLalYsSIrV66kTJkyODo64uPjQ7ly5fDw\n8KBDhw5UrVoVe3t79u7dyx9//EH//v2f6rXkVb4m4Hfv3mX06NFs27aNS5cu4eTkRKtWrRg3bhzF\nixc3jouNjaV///5s3LgRgDZt2jBr1ixjIX5RExj4Bt7erSX5FkIIIZ6Roii0bduWtm3bmjqUPP9c\nv3/crFmzCAgIYP78+XzxxReYmZlRqlQp3nrrLWPJhYODA7///juffPIJY8aMwdbWlo4dOxpb6WXx\n9PTknXfeYfv27axYsQJFUShbtiyLFy/m3XffNY7r3bs3Fy9e5LvvvuOPP/6gfv36bNu2jSZNmuR4\nDQ97TZaWluzatYuJEyfy888/8+OPP2JnZ0eZMmUYNWoUrq6u2c6R23kedu4Hj3fu3BkPDw8mTJjA\nN998Q0pKCu7u7tSpU4c+ffo81TkXL15McHAwQ4YMITU1lXfffZeFCxfSr18/tm3bxsaNG0lPT8fH\nx4dp06YRHByc63nzi2LIx/05IyIiGD16NO+99x7+/v5cuXKFjz76iJIlS7JlyxbjuBYtWhAdHc13\n332HwWDggw8+wNfXlw0bNuQ4Z9ZKX8BkCXp4+H7S0g488Xayer0eK6sy2Nr6S/L9Ajp8+DDwfNQW\nisIhc/5yknkvOCkpKdl6WQtRlDzu8/2oHDZfr4BXrFiRNWvWGG+XLl2aqVOn0rp1axITE7G1teX0\n6dP88ccfhIaGEhgYCGTWIdWrV4/IyMhsTdxfZBkZOtq0Gc2vv/6BnZ0k30IIIYQQIlOBL8KMi4vD\nwsLCWAC/f/9+7OzsqFWrlnFMnTp1sLGxITQ0tKDDKTRWVh6sWrUu269khBBCCCGEKNAEPDY2llGj\nRvHhhx8a+zTGxMTkWvTv7OxMTExMQYZT4NLS0lm3bh+Wlv7Y2VXHx+fJSlaEEEIIIUTRl6cSlJEj\nRzJ+/PiH3q8oCjt37szW8D0pKYnWrVvj6enJ5MmTnz1S/qvTK2xZPeP/+SfyEaMUVCpXduw4R4UK\nMVhaFmz7GlF4TPW5E6Yjc/5yknnPf15eXlIDLoqshIQETpw48dD7H1VWnacEfNCgQXTv3v2RY+5v\nPZOUlESLFi1Qq9Vs3LgRc3Nz432urq7cvHkzx+Nv3LjxQpZr3LmTQHy8Qvny9dHp7Bgy5FNThySE\nEEIIIZ5jeUrAixcvnq2N4KMkJibSokULFEXh999/z7HTU+3atUlMTOSvv/4y1oGHhoZy7949Xn31\n1Uee21Qr1MPD9wPg5/fgNxk1W7ee49q1NN56q2GhxyUKlnRGePnInL+cZN4LzoObxAhRlNjZ2T3y\n3437u6A8KF+7oCQmJtK0aVMSExP55ZdfSEhIICEhAchM4jUaDeXLl6dZs2b06dOHBQsWYDAYCAoK\nonXr1i9MB5SzZ6MoX74atrYV6Nq1hbQYFEIIIYQQeZavizD//vtvDh48yMmTJylbtizu7u64ubnh\n7u7O/v37jeNWrFhB1apVad68OS1atKBatWosXbo0P0MpQJZMmPALN28Wx9xcK8m3EEII8Qj5uN2I\nEM+NZ/1c5+sV8AYNGqDT6R47TqvVvkAJd+abfOtWIv7+Plhb+/Hbb01NHZIQQgjx3DM3NzduViIX\nrERRYTAYSElJwcLC4qnPka8JeFF16lQkCxduoWXLz+UfECGEECKPVCoVFhYWpKammjqUHLJKZO3s\n7EwciShM+TXvFhYWxhbbT0MS8DwoV64iU6d+I8m3EEII8YRUKtVz2Yowq32cLL59uTwv817gO2EW\nFc/yLUcIIYQQQogsklUKIYQQQghRiCQBF0IIIYQQohBJAi6EEEIIIUQhkgRcCCGEEEKIQiQJuBBC\nCCGEEIVIMTznW1TFxcWZOgQhhBBCCCGemlarzXZbroALIYQQQghRiCQBF0IIIYQQohA99yUoQggh\nhBBCFCVyBVwIIYQQQohCJAm4EEIIIYQQhUgS8MeYO3cupUuXxsrKiho1arB3715ThyTy0cSJEwkM\nDESr1eLs7EybNm2IiIjIMW7MmDF4eHhgbW1No0aNOHnypAmiFQVh4sSJqFQqgoODsx2XOS96YmJi\n6NmzJ87OzlhZWVGpUiX27NmTbYzMe9Gi1+sZOXKk8ed46dKlGTlyJHq9Pts4mfcX2549e2jbti0l\nS5ZEpVKxdOnSHGMeN8dpaWn079+fEiVKYGtrS9u2bbly5UqBxSwJ+COsWrWKgQMH8sUXX3DkyBFe\nffVVWrRoQXR0tKlDE/lk9+7d9OvXj/3797Nz507MzMx47bXXiI2NNY6ZPHky06dPZ86cORw+fBhn\nZ2eaNm1KUlKSCSMX+eGvv/5i0aJFVK1aNdtxmfOiJy4ujjp16qAoCps3b+b06dPMmjULZ2dn4xiZ\n96Jn0qRJzJs3j9mzZ3PmzBlmzpzJ3LlzmThxonGMzPuLLzExkcqVKzNz5kysra1z3J+XOR4wYADr\n169n1apV7N27l/j4eFq1akWBLZU0iIeqWbOmoU+fPtmOlSlTxjBixAgTRSQKWmJiokGtVht+++03\n4zE3NzfDxIkTjbeTk5MNdnZ2hoULF5oiRJFPYmNjDb6+voaQkBBDw4YNDf379zfeJ3Ne9AwfPtxQ\nt27dR46ReS96WrVqZejZs2e2Y++++66hdevWxtsy70WLra2tYcmSJdmOPW6O4+LiDObm5oYVK1YY\nx0RFRRlUKpVh69atBRKnXAF/iPT0dP7++2+aNm2a7fjrr79OaGioiaISBS0+Ph69Xo+DgwMAFy5c\nICYmJtvnwNLSkvr168vn4AX34Ycf0rlzZxo0aJDtuMx50bRhwwZq1qzJW2+9hYuLC9WqVWPOnDnG\n+2Xei6a6deuyc+dOzpw5A8DJkyfZsWMHb7zxBiDz/jLIyxwfPnyYjIyMbGNKliyJv79/gX0OzArk\nrEXArVu30Ol0uLi4ZDvu4uLC9u3bTRSVKGgDBgygevXq1K5dG8isGVUUJdfPwdWrV00RosgHixYt\n4vz586xYsSLHfTLnRdP58+eZO3cugwYNYvjw4Rw5coR+/fqhKAoff/yxzHsR9emnn5KQkECFChVQ\nq9XodDo+//xz+vTpA8jf95dBXub4+vXrqNVqHB0dc4yJiYkpkLgkARfiX4MHDyY0NJR9+/ahKIqp\nwxEF5OzZs3z++efs27cPlUp+Cfiy0Ov1BAYGMn78eACqVq3K2bNnmTNnDh9//LGJoxMFZeXKlSxb\ntoyVK1dSoUIFjhw5QnBwMD4+Prz33numDk+8xOSnz0M4OTmhVqu5fv16tuPXr1/H1dXVRFGJgjJo\n0CBWrVrFzp078fLyMh53dXXFYDDI56AI2b9/P7dv36ZChQpoNBo0Gg27du1izpw5mJub4+joKHNe\nBLm5ueHv75/tmL+/P5cvXwbk73pRNWzYMIYOHUqnTp2oWLEib7/9NoMHDzYuwpR5L/ryMseurq7o\ndDpu37790DH5TRLwh9BoNAQEBLBt27Zsx7dt20adOnVMFJUoCAMGDDAm32XKlMl2n4+PD66urtk+\nBykpKezZs0c+By+o9u3bc/z4cY4ePWr8U6NGDbp27crRo0cpW7aszHkRVKdOHWMdcJYzZ84Yv3DL\n3/Wi6d69ezl+06VSqYxtCGXei768zHFAQABmZmbZxkRHR3Pq1KkC+xyox4wZM6ZAzlwE2NvbM3r0\naNzc3LC2tmbs2LHs2bOHxYsXo9VqTR2eyAd9+/Zl6dKlrFmzhpIlS5KUlERSUhKKomBubg6ATqdj\n0qRJlCtXDp1Ox+DBg7l+/ToLFiwwjhEvDgsLC0qUKJHtz/Lly/H29qZHjx6AzHlR5OXlxVdffYVa\nrcbd3Z3t27fzxRdfMGLECGrUqAHIvBdFp06dYtmyZZQrVw5zc3N27tzJ559/TteuXY0L7mTeX3xJ\nSUmcOnWKmJgYvvvuO6pUqYJWqyU9PR2tVvvYObawsODatWvMmTOHKlWqEBcXx0cffYSDgwOTJk0q\nmLLUAumtUoTMmzfP4OPjY7C0tDTUqFHDsHfvXlOHJPKRoigGlUqV48+XX36ZbdyXX35pcHd3N1hZ\nWRkaNmxoiIiIMFHEoiA0atQoWxtCg0HmvCj6/fffDVWrVjVYWVkZypUrZ5g9e3aOMTLvRUtiYqJh\n0KBBBm9vb4O1tbXB19fX8MUXXxhSU1OzjZN5f7GFhITk+vP8vffeM4553BynpaUZgoODDU5OTgYb\nGxtD27ZtDdHR0QUWs2IwFFSHcSGEEEIIIcSDpAZcCCGEEEKIQiQJuBBCCCGEEIVIEnAhhBBCCCEK\nkSTgQgghhBBCFCJJwIUQQgghhChEkoALIYQQQghRiCQBF0IIIYQQohBJAi6EEEIIIUQhkgRcCCGE\nEEKIQvT/tE6qzk8hwWgAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x83306a0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "from numpy.random import randn\n",
    "from filterpy.kalman import EnsembleKalmanFilter as EnKF\n",
    "from filterpy.kalman import KalmanFilter\n",
    "from filterpy.common import Q_discrete_white_noise\n",
    "import book_plots as bp\n",
    "\n",
    "np.random.seed(1234)\n",
    "\n",
    "def hx(x):\n",
    "    return np.array([x[0]])\n",
    "\n",
    "def fx(x, dt):\n",
    "    return np.dot(F, x)\n",
    "    \n",
    "F = np.array([[1., 1.],[0., 1.]])\n",
    "\n",
    "x = np.array([0., 1.])\n",
    "P = np.eye(2) * 100.\n",
    "enf = EnKF(x=x, P=P, dim_z=1, dt=1., N=20, hx=hx, fx=fx)\n",
    "\n",
    "std_noise = 10.\n",
    "enf.R *= std_noise**2\n",
    "enf.Q = Q_discrete_white_noise(2, 1., .001)\n",
    "\n",
    "kf = KalmanFilter(dim_x=2, dim_z=1)\n",
    "kf.x = np.array([x]).T\n",
    "kf.F = F.copy()\n",
    "kf.P = P.copy()\n",
    "kf.R = enf.R.copy()\n",
    "kf.Q = enf.Q.copy()\n",
    "kf.H = np.array([[1., 0.]])\n",
    "\n",
    "measurements = []\n",
    "results = []\n",
    "ps = []\n",
    "kf_results = []\n",
    "\n",
    "zs = []\n",
    "for t in range (0,100):\n",
    "    # create measurement = t plus white noise\n",
    "    z = t + randn()*std_noise\n",
    "    zs.append(z)\n",
    "\n",
    "    enf.predict()\n",
    "    enf.update(np.asarray([z]))\n",
    "    \n",
    "    kf.predict()\n",
    "    kf.update(np.asarray([[z]]))\n",
    "\n",
    "    # save data\n",
    "    results.append (enf.x[0])\n",
    "    kf_results.append (kf.x[0,0])\n",
    "    measurements.append(z)\n",
    "    ps.append(3*(enf.P[0,0]**.5))\n",
    "\n",
    "results = np.asarray(results)\n",
    "ps = np.asarray(ps)\n",
    "\n",
    "plt.plot(results, label='EnKF')\n",
    "plt.plot(kf_results, label='KF', c='b', lw=2)\n",
    "bp.plot_measurements(measurements)\n",
    "plt.plot (results - ps, c='k',linestyle=':', lw=1, label='1$\\sigma$')\n",
    "plt.plot(results + ps, c='k', linestyle=':', lw=1)\n",
    "plt.fill_between(range(100), results - ps, results + ps, facecolor='y', alpha=.3)\n",
    "plt.legend(loc='best');"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "It can be a bit difficult to see, but the KF and EnKF start off slightly different, but soon converge to producing nearly the same values. The EnKF is a suboptimal filter, so it will not produce the optimal solution that the KF produces. However, I deliberately chose $N$ to be quite small (20) to guarantee that the EnKF output is quite suboptimal. If I chose a more reasonable number such as 2000 you would be unable to see the difference between the two filter outputs on this graph."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Outstanding Questions"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "All of this should be considered as *my* questions, not lingering questions in the literature. However, I am copying equations directly from well known sources in the field, and they do not address the discrepancies.\n",
    "\n",
    "First, in Brown [2] we have all sums multiplied by $\\frac{1}{N}$, as in \n",
    "\n",
    "$$ \\hat{x} = \\frac{1}{N}\\sum_{i=1}^N\\chi_k^{(i)}$$\n",
    "\n",
    "The same equation in Crassidis [3] reads (I'll use the same notation as in Brown, although Crassidis' is different)\n",
    "\n",
    "$$ \\hat{x} = \\frac{1}{N-1}\\sum_{i=1}^N\\chi_k^{(i)}$$\n",
    "\n",
    "The same is true in both sources for the sums in the computation for the covariances. Crassidis, in the context of talking about the filter's covariance, states that $N-1$ is used to ensure an unbiased estimate. Given the following standard equations for the mean and standard deviation (p.2 of Crassidis), this makes sense for the covariance.\n",
    "\n",
    "$$\n",
    "\\begin{aligned}\n",
    "\\mu &= \\frac{1}{N}\\sum_{i=1}^N[\\tilde{z}(t_i) - \\hat{z}(t_i)] \\\\\n",
    " \\sigma^2 &= \\frac{1}{N-1}\\sum_{i=1}^N\\{[\\tilde{z}(t_i) - \\hat{z}(t_i)] - \\mu\\}^2\n",
    "\\end{aligned}\n",
    "$$\n",
    "\n",
    "However, I see no justification or reason to use $N-1$ to compute the mean. If I use $N-1$ in the filter for the mean the filter does not converge and the state essentially follows the measurements without any filtering. However, I do see a reason to use it for the covariance as in Crassidis, in contrast to Brown. Again, I support my decision empirically - $N-1$ works in the implementation of the filter, $N$ does not.\n",
    "\n",
    "My second question relates to the use of the $\\mathbf{R}$ matrix. In Brown $\\mathbf{R}$ is added to $\\mathbf{P}_{zz}$ whereas it isn't in Crassidis and other sources. I have read on the web notes by other implementers that adding R helps the filter, and it certainly seems reasonable and necessary to me, so this is what I do. \n",
    "\n",
    "My third question relates to the computation of the covariance $\\mathbf{P}$. Again, we have different equations in Crassidis and Brown. I have chosen the implementation given in Brown as it seems to give me the  behavior that I expect (convergence of $\\mathbf{P}$ over time) and it closely compares to the form in the linear KF. In contrast I find the equations in Crassidis do not seem to converge much.\n",
    "\n",
    "My fourth question relates to the state estimate update. In Brown we have \n",
    "\n",
    "$$\\boldsymbol\\chi  = \\boldsymbol\\chi + \\mathbf{K}[\\mathbf{z} -\\mathbf{z}_{mean} + \\mathbf{v}_R]$$ \n",
    "\n",
    "whereas in Crassidis we have\n",
    "\n",
    "$$\\boldsymbol\\chi  = \\boldsymbol\\chi + \\mathbf{K}[\\mathbf{z} -\\boldsymbol\\chi_h + \\mathbf{v}_R]$$ \n",
    "\n",
    "To me the Crassidis equation seems logical, and it produces a filter that performs like the linear KF for linear problems, so that is the formulation that I have chosen. \n",
    "\n",
    "I am not comfortable saying either book is wrong; it is quite possible that I missed some point that makes each set of equations work. I can say that when I implemented them as written I did not get a filter that worked. I define \"work\" as performs essentially the same as the linear KF for linear problems. Between reading implementation notes on the web and reasoning about various issues I have chosen the implementation in this chapter, which does in fact seem to work correctly. I have yet to explore the significant amount of original literature that will likely definitively explain the discrepancies. I would like to leave this here in some form even if I do find an explanation that reconciles the various differences, as if I got confused by these books than probably others will as well."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## References"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "- [1] Mackenzie, Dana. *Ensemble Kalman Filters Bring Weather Models Up to Date* Siam News,  Volume 36, Number 8, October 2003. http://www.siam.org/pdf/news/362.pdf\n",
    "\n",
    "- [2]  Brown, Robert Grover, and Patrick Y.C. Hwang. *Introduction to Random Signals and Applied Kalman Filtering, With MATLAB® excercises and solutions.* Wiley, 2012.\n",
    "\n",
    "- [3] Crassidis, John L., and John L. Junkins. *Optimal estimation of dynamic systems*. CRC press, 2011."
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.5.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}