docs/distributions.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Introduction\n",
    "Note: this tutorial needs to be updated - some parts may not reflect the usage of Distributios in the master branch.",
    "\n",
    "We use instances of the `Distribution` class to represent information about the set of coefficients required by each module. This information includes:\n",
    "\n",
    "* Boundaries for the coefficient. For example, the Normal distribution can take on any value. The Beta distribution is bounded by 0 and 1, and the Gamma distribution must be a positive value.\n",
    "\n",
    "* The expected value of the coefficient. This can be used as the starting point for the fitter.\n",
    "\n",
    "* The range of probable values. For example, a Normal distribution centered at 10 with a standard deviation of 1 suggests that the coefficient represented by that distribution is very unlikely to have a value of 0 or 100.\n",
    "\n",
    "A coefficient can be multidimensional. For example, the FIR filter coefficient is a two dimensional array with the shape determined by the number of inputs and the number of taps of the FIR filter. As shown below, the distributions also support specifying multidimensional information."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "\n",
    "import numpy as np\n",
    "import pylab as pl\n",
    "\n",
    "np.set_printoptions(precision=2)\n",
    "\n",
    "from nems.distributions.api import Normal, HalfNormal, Gamma, Beta, Uniform\n",
    "from nems.model import Model\n",
    "from nems.modules.api import WeightChannelsGaussian, FIR\n",
    "\n",
    "# Create some fake data that will be used to generate the priors\n",
    "signals = {\n",
    "    'pred': np.zeros((18, 10000))\n",
    "}"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Specifying distributions\n",
    "\n",
    "Distributions can represent a single coefficient"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x7fd659042eb8>"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEKCAYAAAD9xUlFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzt3Xd8V/XZ//HXlU0mI2EkYW/CJgqi\niFYFRAUXCoqjRamt43bUW2t/1aqtWkf11lsraL2tWnFSRUVBBUHLkLAJM4QVyGIFSAhZ1++Pb0hT\nDBmQk/Md1/Px4GHGSfKWR8I753POuT6iqhhjjDEAQW4HMMYY4z2sFIwxxlSxUjDGGFPFSsEYY0wV\nKwVjjDFVrBSMMcZUsVIwxhhTxUrBGGNMFSsFY4wxVULcDtBQ8fHx2qlTJ7djGGOMT1m+fPleVU2o\n6zifK4VOnTqRlpbmdgxjjPEpIrKjPsc5tnwkIm+ISJ6IrKvlmPNEZJWIpIvIAqeyGGOMqR8nrym8\nCYw52TtFpDnwCjBOVVOACQ5mMcYYUw+OlYKqLgT213LIdcBMVd1ZeXyeU1mMMcbUj5t3H/UAWojI\ndyKyXERuPNmBIjJVRNJEJC0/P78JIxpjTGBxsxRCgCHAJcBo4Pci0qOmA1V1uqqmqmpqQkKdF8+N\nMcacIjfvPsoC9qpqIVAoIguBAcBmFzMZY0xAc/NM4VNghIiEiEgkMBTY4GIeY4wJeI6dKYjIDOA8\nIF5EsoBHgFAAVX1VVTeIyFfAGqACeF1VT3r7qml6x8rKSd9ziB37Csk9dIzSsgqCgoT46DCSmkfS\nLymOuMhQt2MaYxqRY6WgqpPqccwzwDNOZTANV3C0lNlrs/l8zR6WbTtASXlFrcd3TYjioj5tGduv\nLf2S4hCRJkpqjHGCzz3RbJyx98gxpi3YyrtLd1JYUk7n+ChuGt6RIR1b0K11DG3jIggPCaK8Qtl7\n5Bjb9xaxOusgSzL38fr3mby6YCspibH84uzOXDYgkbAQG6tljC8SVXU7Q4OkpqaqjbloPGXlFby1\neAfPf72ZotJyLu3fjl+c3Zn+yfX/rf9gUQmfr8nmzUXbycg7QsdWkdw/uieX9GtnZw7GeAkRWa6q\nqXUeZ6UQuHbtL+LOGStZtesgI7rH88hlKXRrHX3Kn09Vmbcxj6e/2sSm3MOc2aklT13Vjy4Jp/45\njTGNw0rB1GpOeg6/+XA1AE9c0Y9L+zfeb/XlFcqHabt4YvYGissquPeiHkwd0YWgIDtrMMYt9S0F\nW/gNQG/8sI3b3llOl/goZt81gssGJDbqMk9wkDDxzA58c+9Izu+ZwFNfbuSm//uRfUeONdrXMMY4\nw0ohgKgqT87ewGOfr2dUnza8/8uzaN8y0rGv1zo2glcnD+GJK/qxdNt+xr74PSt2HnDs6xljTp+V\nQoBQVZ76aiPTFmYyeVgHXrl+CBGhwY5/XRHhuqEd+OevhxMeEsyk6Uv4Yk2241/XGHNqrBQCxP98\nu4VpCzyF8Pj4vgQ38fp+SmIcn9x+Nn2T4rj93RW8umArvnY9y5hAYKUQAD5YtosXvtnChCHJPDau\nr2u3ibaMCuMftwzl0v7teOrLjTz55UYrBmO8jD285ucWb93HQ/9cy4ju8Tx5ZT/X7wCKCA3mxYmD\naBkVxvSFmRwtKefRcSmu5zLGeFgp+LFd+4u47Z3ldIqP4n+vG0xIsHecGAYFCY+OSyEiNJjpCzMp\nKavwisIyxlgp+K2SsgrunLGSigrl9RtTiWvmXYPrRITfXtyL8JAgXpqXQbOwYB65rI89AW2My6wU\n/NQzczayatdBXrl+MJ3io9yOUyMR4d6LenC0pJzXf9hGTEQI943q6XYsYwKalYIfmr8xj9e+38YN\nwzoytl87t+PUSkT43SW9OXKsjJfmZRATEcLUc7u6HcuYgGWl4GcKikp54OM19Gobw+8u6e12nHoR\nEf50RT+OHCvjidkbaRMbwfiBSW7HMiYgWSn4mce/WM++whLeuPmMJnk4rbEEBwnPXTOAvMPHuP/D\nNSQ2b8YZnVq6HcuYgOMdt6OYRjF/Ux4fLc/itpFd6JsU53acBgsPCWb6DUNIbtGMqW+lsW1voduR\njAk4jpWCiLwhInkiUusWmyJyhoiUi8jVTmUJBIXHyvjdzLV0bx3NXRd0dzvOKWseGcYbN58BwC/e\nXMaBwhKXExkTWJw8U3gTGFPbASISDPwZmONgjoDw8vwM9hQU8+SV/QgP8Z1lo5p0io/itRtT2X3g\nKHe9t5LyCnvq2Zim4lgpqOpCYH8dh90JfAzkOZUjEGzbW8hr32dy5eAkUv1kHT61U0sevzyF77fs\n5Zk5m9yOY0zAcO2agogkAVcAr7qVwR+oKo9+lk54SDAPXtzL7TiN6tozOnD90A68umArn6/Z43Yc\nYwKCmxeaXwAeUNXyug4UkakikiYiafn5+U0QzXfM35THd5vyufvC7rSOiXA7TqN75LIUhnRswf0f\nrmFjziG34xjj99wshVTgPRHZDlwNvCIil9d0oKpOV9VUVU1NSEhoyoxerbxCeerLjXSJj+Km4Z3c\njuOIsJAg/nr9YGIiQvjl28s5VFzqdiRj/JprpaCqnVW1k6p2Aj4Cfq2qn7iVxxfNXJHF5twj/GZ0\nT0K9ZNidE1rHRvDK9YPJOnCU385ca+O2jXGQk7ekzgAWAz1FJEtEpojIbSJym1NfM5AUl5bz/Neb\nGZAcx8V927odx3GpnVpy36gefLEmmxk/7nI7jjF+y7EnmlV1UgOOvdmpHP7qnSU72FNQzLMTBgTM\nZNHbzu3K4q37ePSzdAZ1aE7vdrFuRzLG7/jvmoMfO1xcysvzMxjRPZ7h3eLdjtNkgoKE568dSGyz\nUG5/dwWFx8rcjmSM37FS8EFvLd7BgaJS7h8deGOm46PD+Z+JA9m2t5Dff1rrw/LGmFNgpeBjikrK\neP37TM7vmUD/5OZux3HF8K7x3PWz7sxcsZtPVu52O44xfsVKwcf8Y8lODhSVcsfPfHe+UWO464Lu\npHZswe8/Xcfug0fdjmOM37BS8CHFpeVM/z6Ts7u1YkjHFm7HcVVwkPCXawZSUaHc98EqKmw+kjGN\nwkrBh7y/bBf5h49xZ4CfJRzXoVUkj1yWwpLM/bz+Q6bbcYzxC1YKPqK0vIJpC7ZyRqcWDO3sH0Pv\nGsOE1GRG9WnDs3M2syHbxmAYc7qsFHzE7LXZ7Cko5raRXQPmuYT6EBGevLIfsc1Cufu9VRSX1jlK\nyxhTCysFH6Cq/O2HbXRJiOL8nq3djuN1WkWH88zV/dmUe5hnbcy2MafFSsEHLNt+gDVZBUw5pzNB\nQXaWUJPze7Vm8rAOvP7DNpZk7nM7jjE+y0rBB7z+fSYtIkO5clCy21G82kNje9OhZSQPfLyGohJ7\n2tmYU2Gl4OW27y3k6w25TB7WkWZhvr3NptMiw0J4+ur+7NhXxNNf2TKSMafCSsHL/d+/thEaFMQN\nZ3V0O4pPGNalFTed1ZE3F21nqS0jGdNgVgpe7FBxKR8uz+KyAYl+uauaU/57TC/at2zGf3+8hqMl\ndjeSMQ1hpeDF/rliN0Ul5dw03M4SGiIqPISnrxrAjn1FPGN3IxnTIFYKXkpVeWfJDvonxwXs4LvT\ncVbXVtx4Vkf+b9E2lm3f73YcY3yGlYKX+nHbfrbkHWHyUDtLOFUPjOlFUvNm3P/haltGMqaenNyO\n8w0RyRORGofei8j1IrKm8s8iERngVBZf9PaSHcRGhHDZgES3o/gszzJSf7bbMpIx9ebkmcKbwJha\n3r8NGKmq/YHHgekOZvEp+YePMSc9h6uHtLfbUE/T8G7xTB7Wgf9btI2VOw+4HccYr+dYKajqQuCk\ni7mqukhVj/+ULgHsyaxKH6TtorRcuX5YB7ej+IUHxvSiTUwED368lpKyCrfjGOPVvOWawhTgS7dD\neIPyCuXdpTsZ3rUVXROi3Y7jF2IiQnn88r5syj3MtAVb3Y5jjFdzvRRE5Hw8pfBALcdMFZE0EUnL\nz89vunAu+FfGXnYfPMqkM+0soTFd1KcNl/Rvx0vzMsjIO+J2HGO8lqulICL9gdeB8ap60sdPVXW6\nqqaqampCQkLTBXTBB2m7aB4ZyqiUNm5H8Tt/uCyFZmHBPDRzre3UZsxJuFYKItIBmAncoKqb3crh\nTQ4WlTA3PZfLByYRHmIXmBtbQkw4v7ukNz9u38+MZTvdjmOMV3LyltQZwGKgp4hkicgUEblNRG6r\nPORhoBXwioisEpE0p7L4ilmr91BSXsGEVLvm7pQJQ5IZ3rUVT83eSE5BsdtxjPE6oupbp9Gpqama\nluaf/XHpS9+jCl/cNcLtKH5tx75CRj2/kHN7JDD9hiG2k50JCCKyXFVT6zrO9QvNxiN9TwHrdh9i\nwhA7S3Bax1ZR3HNRD75en8tX63LcjmOMV7FS8BIfpmURFhzE+IFJbkcJCLec05mUxFgenpVOQVGp\n23GM8RpWCl7gWFk5n6zazUUpbWgRFeZ2nIAQEhzEn6/qz/7CEp76aoPbcYzxGlYKXuDbDXkcLCq1\npaMm1jcpjlvO6cyMH3exeKttyGMMWCl4hX+u3E1CTDgjuvv3Mxje6O4Le9ChZSQP/XMtxaU2SdUY\nKwWXHSwq4btNeYwbkEhwkN0F09SahQXz5JX92La3kJfmbXE7jjGus1Jw2ey1OZSWK5fbBWbXnN0t\nnqsGJzNtQSYbsg+5HccYV1kpuOyTlbvpmhBF36RYt6MEtP93SW/imoXy4My1lNsIDBPArBRclHWg\niB+37+fygUn2AJXLWkSF8fBlfVi96yBvLd7udhxjXGOl4KJZq/cA2LMJXmLcgERG9kjgmTmb2H3w\nqNtxjHGFlYJLVJVPVu5mSMcWdGgV6XYcA4gIf7y8L6rw+0/W4WsjYIxpDFYKLtmQfZjNuUe4fKDt\nwexN2reM5L5RPZi3MY/P12S7HceYJmel4JJPV+0mJEi4pL+Vgrf5+dmd6Z8cx6OfpXOwqMTtOMY0\nKSsFF1RUKJ+u2sPIHgm0tLEWXic4SHjqyv4cKCrlydkb3Y5jTJOyUnDB0m37yTlUzPhBdoHZW/VJ\njOXWEV14P81GYJjAYqXggi/W7qFZaDAX9m7tdhRTi/+6oLuNwDABx0qhiZWVV/DVuhx+1rs1kWEh\nbscxtWgWFswTV9gIDBNYnNyO8w0RyRORdSd5v4jIiyKSISJrRGSwU1m8yY/b9rP3SAmX9mvndhRT\nD+d0txEYJrA4eabwJjCmlvdfDHSv/DMV+KuDWbzG52uziQwL5ryetnTkK2wEhgkkjpWCqi4E9tdy\nyHjgLfVYAjQXEb/+9fn40tEFvdvQLCzY7TimnmwEhgkkbl5TSAJ2VXs9q/Jtfmvptv3sLyzhkn5t\n3Y5iGshGYJhA4WYp1DQBrsZzcxGZKiJpIpKWn5/vcCznfL7Glo58lY3AMIHCzVLIAtpXez0Z2FPT\ngao6XVVTVTU1IcE3dyfzLB1lc2HvNkSE2tKRL7IRGCYQuFkKs4AbK+9CGgYUqKrf/qQtztzHgaJS\nLunv15dN/J6NwDD+zslbUmcAi4GeIpIlIlNE5DYRua3ykNlAJpABvAb82qks3uCLNdlEhQUzsodv\nnukYj+ojMJ6YvcHtOMY0OseenlLVSXW8X4Hbnfr63qS0vII56Tlc1MeWjvzB8REYry7YyuWDkhje\nNd7tSMY0GnuiuQks3upZOhprD6z5jbsv7E7HVpE8NNNGYBj/YqXQBL5Yk010eAjn2tKR34gI9YzA\n2L6vyEZgGL9ipeCwsvIK5q7P4cLerW3pyM+c3S2eq4fYCAzjX6wUHPbj9v0cKCplTF97YM0f/W6s\njcAw/sVKwWFz03MJDwmypSM/ZSMwjL+xUnCQqjI3PYdzeyTYmGw/ZiMwjD+xUnDQut2H2FNQzKg+\nbdyOYhxkIzCMP7FScNCc9ByCg4QLe1sp+DsbgWH8hZWCg+ak53Bmp5a0iApzO4ppAjYCw/gDKwWH\nZOYfYUveEUal2FlCoDg+AuNgUSmPfbbe7TjGnBIrBYfMXZ8LwKgUuxU1kPRJjOXX53dj5srdfFP5\nPWCML7FScMic9Bz6JcWR1LyZ21FME7vj/G70ahvDQ/9cS0FRqdtxjGkQKwUH5B4qZuXOg4y2paOA\nFBYSxLMTBrC/sIRHP0t3O44xDVJrKYjI3Gov/9b5OP7ha1s6Cnh9k+JsGcn4pLrOFKo/hjvBySD+\nZE56Dp3jo+jeOtrtKMZFtoxkfFFdpWBP4TRQwdFSFm/dx6iUNojUtA21CRS2jGR8UV2zF7qIyCxA\nqr1cRVXHOZbMR83fmEdZhTLalo4M/15GevHbLYzt144L7el24+XqKoXx1V5+1skg/mLu+hxax4Qz\nMLm521GMl7jj/G7MTc/hoX+u5YxOLYmLDHU7kjEnVevykaouOP4HWA+sP+FttRKRMSKySUQyROTB\nGt7fQUTmi8hKEVkjImNP/X/FfcWl5Xy3KZ+L+rQhKMiWjoyHLSMZX1LX3UciIo+IyF5gI7BZRPJF\n5OG6PrGIBAMvAxcDfYBJItLnhMP+H/CBqg4CJgKvnMr/hLf4YcteikrKbenI/ET1u5Hmpue4HceY\nk6rrQvPdwDnAGaraSlVbAEOBs0Xknjo+9kwgQ1UzVbUEeI//XI4Cz4Xs2MqX44A9DUrvZeak5xAT\nEcKwLq3cjmK80B3nd6NPu1h+O3Mte48cczuOMTWqqxRuBCap6rbjb1DVTGBy5ftqkwTsqvZ6VuXb\nqvsDMFlEsoDZwJ01fSIRmSoiaSKSlp+fX8eXdUdZeQXfbMjlZ71aExZizwSanwoLCeKFiQM5fKyM\nBz9eayO2jVeq61+vUFXde+IbVTUfqOtqWU2L6if+FEwC3lTVZGAs8LaI/CSTqk5X1VRVTU1I8M4d\nzNJ2HOBAUaktHZla9WgTwwNjevHNhlzeX7ar7g8wponVVQq1zf+tazZwFtC+2uvJ/HR5aArwAYCq\nLgYigPg6Pq9XmpOeQ1hIECNt201Th58P78Twrq147PP17NhX6HYcY/5DXaUwQEQOicjhyj+Hjr8O\n9KvjY5cB3UWks4iE4bmQPOuEY3YCFwCISG88peCd60O18Gy7mcuIbvFEhdu2m6Z2QUHCsxMGEBwk\n3PP+KsrKK9yOZEyVum5JDVbVWFWNqfwTW+31WpePVLUMuAOYA2zAc5dRuog8JiLHH3q7D7hVRFYD\nM4Cb1QcXWtP3HGL3waO2dGTqLbF5M/54eV9W7DzItIWZbscxpkqtv9aKSARwG9ANWAO8UfmPfb2o\n6mw8F5Crv+3hai+vB85uSGBvNDc9hyCBC3q3djuK8SHjBiTy9fpcnv96M+d2T6BfcpzbkYypc/no\n70AqsBbPheDnHE/kg+ak53JGp5a0ig53O4rxISLCHy/vS3x0OHe/v5Li0nK3IxlTZyn0UdXJqjoN\nuBoY0QSZfMr2vYVsyj1sY7LNKWkeGcYzE/qzNb+QJ2dvcDuOMXWWQtW834YsGwWSues9T6eOskFn\n5hSN6J7Az8/uxN8X77C9F4zr6nv30fE7jvpXuxvpUFME9HZz0nNJSYylfctIt6MYH/bgxb3o0y6W\n+z9aTU5BsdtxTACr791Hx+84Cqn2cmxtHxsI8g4Xs2LnAbvryJy28JBgXrpuEMfKKrj7/ZWUV/jc\nTXjGT9g8htPw9fpcVGGU7cVsGkHXhGgeHZfCksz9vDI/w+04JkBZKZyGuem5dGwVSc82MW5HMX7i\n6iHJjB+YyAvfbiFt+36345gAZKVwig4Vl7Jo615Gp7S1bTdNozl+m2pyi2b813urbG9n0+SsFE7R\n/I15lJYro23pyDSymIhQXpw4iNxDxTw4c41NUzVNykrhFM1NzyU+OpxB7Vu4HcX4oQHtm3P/6J58\nuS6Hd5budDuOCSBWCqfAs+1mnm27aRx164gunNczgcc/W8+arINuxzEBwkrhFCzaupfCknJbOjKO\nCgoSnr9mIAkx4fzqnRUcLKprWr0xp89K4RTMWZdLdHgIZ3W1bTeNs1pEhfHy9YPJO1zMvR+spsKe\nXzAOs1JooPIK5ZsNuZzfqzXhIcFuxzEBYGD75vz+0j7M25jHXxdsdTuO8XNWCg2Utn0/+wpLGGNP\nMZsmdMOwjowbkMhzczexKOMnO+Qa02isFBpoTnouYSFBnNfTtt00TUdEePLKfnRJiOau91bafCTj\nGCuFBlBV5qTn2LabxhVR4SG8OnkwRSXl3PHuCkrKbBtP0/gcLQURGSMim0QkQ0QePMkx14jIehFJ\nF5F3ncxzuqq23exrS0fGHd1ax/DUVf1J23GAxz5PdzuO8UOO/borIsHAy8BFQBawTERmVW7BefyY\n7sBvgbNV9YCIePV+lnMqt928sLfdimrcM25AIum7C5i2MJOUxDgmndnB7UjGjzh5pnAmkKGqmapa\nArwHjD/hmFuBl1X1AICq5jmY57R9tS6HMzu3pGVUmNtRTID77zG9OLdHAg9/us4G55lG5WQpJAG7\nqr2eVfm26noAPUTkXyKyRETG1PSJRGSqiKSJSFp+fr5DcWuXmX+ELXlHbO8E4xWCg4SXJg4iqXkz\nbntnBdkFR92OZPyEk6VQ0/yHE5+8CQG6A+cBk4DXRaT5Tz5IdbqqpqpqakKCO3f9zEn3bJNopWC8\nRVxkKNNvTOVoSRm/fHs5xaXlbkcyfsDJUsgC2ld7PRnYU8Mxn6pqqapuAzbhKQmv81V6Dv2T40hs\n3sztKMZU6dEmhuevHciarAIemrnWJqqa0+ZkKSwDuotIZxEJAyYCs0445hPgfAARiceznJTpYKZT\nklNQzOpdB+0swXilUSltufvC7sxcuZtXF3jdj4/xMY7dfaSqZSJyBzAHCAbeUNV0EXkMSFPVWZXv\nGyUi64Fy4H5V3edUplM1d30OYEtHxnvd9bPuZOQd4c9fbaRDy0gu6d/O7UjGRzn6BJaqzgZmn/C2\nh6u9rMC9lX+81pz0HLomRNGtdbTbUYypUVCQ8OyEAWQXFHPvB6to1zyCwR1srw/TcPZEcx0OFJaw\nJHO/nSUYrxcRGsz0G4bQJjaCqW+lsWt/kduRjA+yUqjDtxvzKK9QxthTzMYHtIoO542bz6CkrIJf\nvLmMgqO2x7NpGCuFOsxJzyExLoJ+SXFuRzGmXrq1jubVyUPYtreQ2/+xgtJym5Fk6s9KoRZFJWUs\n3JzPqJS2iNi2m8Z3DO8WzxNX9uOHjL08+LHdqmrqz0Z91mL+xnyOlVUwyrbdND7omtT2ZB8s5vlv\nNhMfE8ZvL+7tdiTjA6wUajF7bTbx0WEM7WzbbhrfdNcF3cg/Usy0BZkkRIdzy4gubkcyXs5K4SSO\nlpQzb2MeVw5OIjjIlo6MbxIRHh3Xl/2FJfzxiw3ER4dz+aATR5AZ8292TeEk5m/K42hpOZf0s4eA\njG8LDhKev3YgZ3VpxW8+XM13m7x6GLFxmZXCSXyxNptWUWGc2bml21GMOW3hIcFMv3EIPdrE8Kt3\nVrBy5wG3IxkvZaVQg6Ml5czbkMfovm0JCba/IuMfYiJCefMXZ9A6Npyb3viRdbsL3I5kvJD9i1eD\n72zpyPip1jER/OOWoUSHh3DD35ayOfew25GMl7FSqMEXa7NpGRXGUFs6Mn4ouUUk7946jNDgIK57\nbSmZ+UfcjmS8iJXCCYpLPXcdjU6xpSPjvzrFR/HurUNRVa57bSk799mcJONh/+qd4LtNeRSV2NKR\n8X/dWsfwzi1DKS4r57rXl7D7oG3paawUfuKLtTm0jApjWBdbOjL+r3e7WN7+xVAKjpZy7bTFNlnV\nWClUV1xazrcbchmd0saWjkzA6Jccx7u3DONwcRnXTFvMtr2FbkcyLrJ/+aqZv9GzdDTWlo5MgOmX\nHMeMW4dxrKyCa6ctJiPP7koKVI6WgoiMEZFNIpIhIg/WctzVIqIikupknrp8umoP8dHhnNXFZh2Z\nwNMnMZb3pg6jQmHi9CVsyrFiCESOlYKIBAMvAxcDfYBJItKnhuNigLuApU5lqY+Co6XM25THpf3b\n2dKRCVg92sTw/i+HERwkTJy+mDVZB92OZJqYk//6nQlkqGqmqpYA7wHjazjuceBpoNjBLHWak55D\nSVmFDQszAa9rQjQf/PIsosJDmDR9CT9s2et2JNOEnCyFJGBXtdezKt9WRUQGAe1V9XMHc9TLrFV7\n6NgqkgHJtsOaMR1bRfHxr4bTvmUkP3/zRz5bvcftSKaJOFkKNc2brtr+SUSCgOeB++r8RCJTRSRN\nRNLy8/MbMaJH3qFiFm3dy/gBibbDmjGV2sRG8P4vz2Jg++bc9d5K/r5ou9uRTBNwshSygPbVXk8G\nqv+6EQP0Bb4Tke3AMGBWTRebVXW6qqaqampCQkKjB/18TTYVCuMGJjb65zbGl8U1C+XtKUO5oFcb\nHpmVznNzN9nWnn7OyVJYBnQXkc4iEgZMBGYdf6eqFqhqvKp2UtVOwBJgnKqmOZipRp+u3kNKYizd\nWsc09Zc2xutFhAbz6uTBXJOazEvzMrjvg9UcKyt3O5ZxiGOloKplwB3AHGAD8IGqpovIYyIyzqmv\n21Db9xayetdBxttZgjEnFRIcxJ+v6s+9F/Vg5srd3PD6jxwoLHE7lnGAo9txqupsYPYJb3v4JMee\n52SWk/l01R5E4LIBVgrG1EZEuOuC7nRsFcn9H63hilf+xd9uPoOuCdFuRzONKKBvyK+oUD5ekcVZ\nXVrRLq6Z23GM8QnjByYx49ahHC4u48pXFrF46z63I5lGFNClsGz7fnbuL2JCarLbUYzxKUM6tuST\n288mISacG/62lLcWb7cL0H4ioEvhw+VZRIeHMCbFZh0Z01DtW0Yy89fDGdkjgYc/Tec3H66huNQu\nQPu6gC2FwmNlzF6bzaX929EsLNjtOMb4pNiIUF67MZW7L+zOxyuyuPrVRWQdsPHbvixgS2H22myK\nSsq5eogtHRlzOoKChLsv7MHfbkplx74iLnvpBxuN4cMCthQ+XJ5F5/gohnRs4XYUY/zCBb3bMOuO\nczzXGd5Yyl++3kxZeYXbsUwDBWQp7NhXyI/b9nP1kGQba2FMI+ocH8Unt5/NVYOTefHbLVz32lKy\nC2ybT18SkKXw8fIsROAKm4jw1KSkAAANXElEQVRqTKOLDAvh2QkDeP7aAazbU8DF//M936zPdTuW\nqaeAK4XS8greW7aLc7snkNjcnk0wxilXDErm8zvPIal5M255K40/zErnaIndneTtAq4Uvt2QS97h\nY0we1tHtKMb4vS4J0cz89XBuHt6JNxdt55KXvmfVLtu4x5sFXCn8Y+lOEuMi+Fmv1m5HMSYghIcE\n84dxKbwzZSjFJeVc9ddFPDd3EyVldhHaGwVUKWzfW8j3W/Yy8cwOBAfZBWZjmtI53eP56p5zuXxg\nEi/Ny+Dyl//FhuxDbscyJwioUpjx407P3rNntK/7YGNMo4uNCOW5awYw7YYh5B0u5rKXfuDprzba\nk9BeJGBKobi0nA/SdjGqTxtax0a4HceYgDY6pS1z7xnJ5YOSeOW7rYx+YaE98OYlAqYUvlqXw4Gi\nUq4faheYjfEGLaPCeHbCAN69ZSgCTP7bUu59fxX7jhxzO1pAE1+bbJiamqppaQ3fnK2gqJTZ67K5\nNrU9QXY9wRivUlxazv/Oy+DVBVuJCg/h3ot6cP3QDoQEB8zvrY4TkeWq+pPtjn9yXKCUgjHG+23O\nPcyjn6Xzr4x99GgTzSOXpXB2t3i3Y/mF+paCozUsImNEZJOIZIjIgzW8/14RWS8ia0TkWxGxtR1j\nAliPNjG8M2Uo024YwtHScq5/fSlT30pjx75Ct6MFDMdKQUSCgZeBi4E+wCQR6XPCYSuBVFXtD3wE\nPO1UHmOMbxARRqe05et7RnL/6J78kLGXi/6ykEc/S2evXW9wnJNnCmcCGaqaqaolwHvA+OoHqOp8\nVT0+fH0JYHOsjTEARIQGc/v53Zj/m/O4akgSby3ewcin5/OXrzdzuLjU7Xh+y8lSSAJ2VXs9q/Jt\nJzMF+NLBPMYYH9QmNoInr+zP3HvO5byerXnx2y2c+/R8Xv8+055vcICTpVDTLT41XtUWkclAKvDM\nSd4/VUTSRCQtPz+/ESMaY3xF14RoXr5+MLPuOJu+SXH88YsNnPPn+UxbsJUjx8rcjuc3nCyFLKD6\no8PJwJ4TDxKRC4HfAeNUtcYFQ1WdrqqpqpqakJDgSFhjjG/on9yct6cM5f2pw+jdLoYnv9zIOX+e\nx4vfbqHgqC0rnS7HbkkVkRBgM3ABsBtYBlynqunVjhmE5wLzGFXdUp/Pa7ekGmOqW7nzAC/Pz+Cb\nDXnEhIdw/bCO3DS8I+3ibDR+dV7xnIKIjAVeAIKBN1T1TyLyGJCmqrNE5BugH5Bd+SE7VXVcbZ/T\nSsEYU5P0PQW8Mn8rX67LJkiEsf3aMeWczgxo39ztaF7BK0rBCVYKxpja7NpfxN8Xbef9Zbs4fKyM\nIR1bMOWczlzUpw2hAfyEtJWCMSagHS4u5cO0LN5ctJ2d+4tIiAlnwpBkJp7RgQ6tIt2O1+SsFIwx\nBiivUL7blMeMH3cxb2MuFQojuscz6cwOXNi7DWEhgXH2YKVgjDEnyC44yodpWby/bBe7Dx6lZVQY\nl/Zvx+WDkhjUvjki/jss00rBGGNOorxCWbg5n49WZPHN+lyOlVXQsVUk4wckMn5QEl0Tot2O2Ois\nFIwxph4OF5fy1bocPlm1m0Vb96EK/ZLiGNO3LWP6tvWbgrBSMMaYBso9VMxnq/fw2eo9rM4qAKBb\n62jGpLRldEpb+ibF+uwSk5WCMcachuyCo8xNz+WrdTks3baPCoWk5s0Y2TOB83okMLxbPNHhIW7H\nrDcrBWOMaST7C0v4ZkMuX6/PZVHGXgpLygkJElI7tWBkj9ac1zOBXm1jvPoswkrBGGMcUFJWwfId\nB/hucx4LNuWzMecwAK2iwhjWpRVDu7RkWJdWdG8d7VUlYaVgjDFNIPdQMQs257Nk6z6WZO5jT0Ex\n4CmJoV1aMrRzKwZ3aEGvdjGuPlFtpWCMMU1MVck6cJTFmZ6CWJq5n90HjwIQHhJEv6Q4BnVozsD2\nLRjYoTmJcRFNdjZhpWCMMV4g60ARK3ceZNWug6zceYB1ew5RUlYBQOuYcPolxdEnMZY+7WLp3S6W\nDi0jCQpq/KKobyn4zqVzY4zxQcktIkluEcllAxIBzzWJDdmHqkoifc8hvtucT3mF5xf0qLBgereL\npU+ipyR6tImmW+sY4pqFNkleO1MwxhiXFZeWszn3MOv3HGJD9iHWZx9iQ/bh/9hRrk1sOLeO6MIt\nI7qc0tewMwVjjPEREaHB9E9uTv/kf+/9UFGh7D54lM25h9mce4QteYdJiAl3PIuVgjHGeKGgIKF9\ny0jat4zkgt5tmu7rNtlXMsYY4/UcLQURGSMim0QkQ0QerOH94SLyfuX7l4pIJyfzGGOMqZ1jpSAi\nwcDLwMVAH2CSiPQ54bApwAFV7QY8D/zZqTzGGGPq5uSZwplAhqpmqmoJ8B4w/oRjxgN/r3z5I+AC\n8abnwo0xJsA4WQpJwK5qr2dVvq3GY1S1DCgAWjmYyRhjTC2cLIWafuM/8aGI+hyDiEwVkTQRScvP\nz2+UcMYYY37KyVLIAtpXez0Z2HOyY0QkBIgD9p/4iVR1uqqmqmpqQkKCQ3GNMcY4WQrLgO4i0llE\nwoCJwKwTjpkF3FT58tXAPPW1R6yNMcaPODrmQkTGAi8AwcAbqvonEXkMSFPVWSISAbwNDMJzhjBR\nVTPr+Jz5wA7HQnvEA3sd/hpOsezu8OXs4Nv5LXv9dFTVOpdafG72UVMQkbT6zAjxRpbdHb6cHXw7\nv2VvXPZEszHGmCpWCsYYY6pYKdRsutsBToNld4cvZwffzm/ZG5FdUzDGGFPFzhSMMcZUCehSqMcU\n13tFZL2IrBGRb0Wkoxs5a1JX9mrHXS0iKiJec4dDfbKLyDWVf/fpIvJuU2c8mXp8z3QQkfkisrLy\n+2asGzlrIiJviEieiKw7yftFRF6s/H9bIyKDmzrjydQj+/WVmdeIyCIRGdDUGU+mruzVjjtDRMpF\n5OqmylYjVQ3IP3iendgKdAHCgNVAnxOOOR+IrHz5V8D7bueub/bK42KAhcASINXt3A34e+8OrARa\nVL7e2u3cDcg+HfhV5ct9gO1u566W7VxgMLDuJO8fC3yJZ/zMMGCp25kbkH14te+Xi30pe7XvrXnA\nbOBqN/MG8plCnVNcVXW+qhZVvroEz6gOb1CfCbQAjwNPA8VNGa4O9cl+K/Cyqh4AUNW8Js54MvXJ\nrkBs5ctx/HS0i2tUdSE1jJGpZjzwlnosAZqLSLumSVe7urKr6qLj3y94189qff7eAe4EPgZc/14P\n5FKozxTX6qbg+S3KG9SZXUQGAe1V9fOmDFYP9fl77wH0EJF/icgSERnTZOlqV5/sfwAmi0gWnt/6\n7myaaI2ioT8T3sqbflbrJCJJwBXAq25ngcDeo7leE1oBRGQykAqMdDRR/dWaXUSC8GxadHNTBWqA\n+vy9h+BZQjoPz29834tIX1U96HC2utQn+yTgTVV9TkTOAt6uzF7hfLzTVu+fCW8lIufjKYVz3M7S\nAC8AD6hquTdsJxPIpVCfKa6IyIXA74CRqnqsibLVpa7sMUBf4LvKb7K2wCwRGaeqaU2Wsmb1nZ67\nRFVLgW0isglPSSxrmognVZ/sU4AxAKq6uHK+VzxesCxQD/X6mfBWItIfeB24WFX3uZ2nAVKB9yp/\nVuOBsSJSpqqfuBEmkJeP6pziWrkEMw0Y50Xr2lBHdlUtUNV4Ve2kqp3wrLF6QyFA/abnfoLnIj8i\nEo9nOanWQYlNpD7ZdwIXAIhIbyAC8JVNQGYBN1behTQMKFDVbLdD1YeIdABmAjeo6ma38zSEqnau\n9rP6EfBrtwoBAvhMQVXLROQOYA7/nuKaXn2KK/AMEA18WNniO1V1nGuhK9Uzu1eqZ/Y5wCgRWQ+U\nA/d7w29+9cx+H/CaiNyDZ+nlZq28vcRtIjIDz5JcfOU1j0eAUABVfRXPNZCxQAZQBPzcnaQ/VY/s\nD+PZtfGVyp/VMvWSQXP1yO5V7IlmY4wxVQJ5+cgYY8wJrBSMMcZUsVIwxhhTxUrBGGNMFSsFY4wx\nVawUjDmBiHwnIqNPeNvdIvJKLR9zxPlkxjjPSsGYn5qB58G06iZWvt0Yv2alYMxPfQRcKiLhACLS\nCUgEVlXuq7FCRNaKyE8m04rIeSLyebXX/1dEbq58eYiILBCR5SIyx1smkBpTnZWCMSeofHr6Rypn\nGOE5S3gfOApcoaqD8YzheE7qOcFMREKBl/DMyh8CvAH8qbGzG3O6AnbMhTF1OL6E9Gnlf3+BZ4ro\nEyJyLlCBZ6x0GyCnHp+vJ54hhV9X9kgw4BNzhUxgsVIwpmafAH+p3JKymaquqFwGSgCGqGqpiGzH\nM/CuujL+8wz8+PsFSFfVs5yNbczpseUjY2qgqkeA7/As8xy/wBwH5FUWwvlATXt27wD6iEi4iMRR\nOTEV2AQkVO6xgIiEikiKk/8PxpwKO1Mw5uRm4BnHfPxOpH8An4lIGrAK2HjiB6jqLhH5AFgDbMGz\n1zSqWlK5IfuLlWURgmdzlXTH/y+MaQCbkmqMMaaKLR8ZY4ypYqVgjDGmipWCMcaYKlYKxhhjqlgp\nGGOMqWKlYIwxpoqVgjHGmCpWCsYYY6r8f9GBFz0GM5ooAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fd68816e828>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEKCAYAAADjDHn2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAFDpJREFUeJzt3X+QXXd93vH3g2RhMmBwozUDkhaZ\nZpXiuAngGyegwNjNGBSmtZMZ11WaBEgzVtOMmZBpPbEzndCKaadpQ9pJUCZ1iSchE+M6KbUFTRFM\nYhyPg4hWYGy0rkDIIdqKjoyQQ2jSSDKf/nHP2tfLar8rec/ulfx+zezonnO/995HP46ee77n3HNT\nVUiStJgXrHYASdL4sywkSU2WhSSpybKQJDVZFpKkJstCktRkWUiSmiwLSVKTZSFJalq72gGWy/r1\n62vz5s2rHUOSziv79+//alVNtMZdMGWxefNmpqenVzuGJJ1Xknx5KeOchpIkNVkWkqQmy0KS1GRZ\nSJKaLAtJUpNlIUlqsiwkSU2WhSSpybKQJDVZFpKkJstCktRkWUiSmiwLSVKTZSFJarIsJElNloUk\nqcmykCQ1WRaSpCbLQpLUZFlIkposC0lSk2UhSWqyLCRJTZaFJKmp17JIsi3JwSSHktx2hjE3JZlJ\nciDJXSPrn0rycPezu8+ckqTFre3riZOsAXYB1wGzwL4ku6tqZmTMFHA7sLWqTiS5bOQp/rqqXttX\nPknS0vW5Z3E1cKiqDlfVSeBu4IZ5Y24GdlXVCYCqOtZjHknSOeqzLDYAR0aWZ7t1o7YAW5I8lGRv\nkm0j912cZLpb/8M95pQkNfQ2DQVkgXW1wOtPAdcAG4EHk1xZVU8Ck1V1NMmrgT9K8mhVfelZL5Ds\nAHYATE5OLnd+SVKnzz2LWWDTyPJG4OgCY+6rqlNV9ThwkGF5UFVHu18PA58EXjf/BarqjqoaVNVg\nYmJi+X8HkiSg37LYB0wluTzJOmA7MP+spnuBawGSrGc4LXU4yaVJXjiyfiswgyRpVfQ2DVVVp5Pc\nAuwB1gB3VtWBJDuB6ara3d33liQzwFPArVV1PMkbgf+c5JsMC+3fjZ5FJUlaWamafxjh/DQYDGp6\nenq1Y0jSeSXJ/qoatMb5CW5JUpNlIUlqsiwkSU2WhSSpybKQJDVZFpKkJstCktRkWUiSmiwLSVKT\nZSFJarIsJElNloUkqcmykCQ1WRaSpCbLQpLUZFlIkposC0lSk2UhSWqyLCRJTZaFJKnJspAkNVkW\nkqQmy0KS1NRrWSTZluRgkkNJbjvDmJuSzCQ5kOSuefddkuR/J3l/nzklSYtb29cTJ1kD7AKuA2aB\nfUl2V9XMyJgp4HZga1WdSHLZvKd5L/BAXxklSUvT557F1cChqjpcVSeBu4Eb5o25GdhVVScAqurY\n3B1JrgJeDny8x4ySpCXosyw2AEdGlme7daO2AFuSPJRkb5JtAEleALwPuLXHfJKkJeptGgrIAutq\ngdefAq4BNgIPJrkS+HHgD6rqSLLQ03QvkOwAdgBMTk4uQ2RJ0kL6LItZYNPI8kbg6AJj9lbVKeDx\nJAcZlscbgDcl+RngxcC6JN+oqmcdJK+qO4A7AAaDwfwikiQtkz6nofYBU0kuT7IO2A7snjfmXuBa\ngCTrGU5LHa6qH6uqyaraDPwL4IPzi0KStHJ6K4uqOg3cAuwBHgPuqaoDSXYmub4btgc4nmQGuB+4\ntaqO95VJknRuUnVhzN4MBoOanp5e7RiSdF5Jsr+qBq1xfoJbktRkWUiSmiwLSVKTZSFJarIsJElN\nloUkqcmykCQ1WRaSpCbLQpLUZFlIkposC0lSk2UhSWqyLCRJTZaFJKnJspAkNVkWkqQmy0KS1GRZ\nSJKaLAtJUpNlIUlqsiwkSU2WhSSpybKQJDX1WhZJtiU5mORQktvOMOamJDNJDiS5q1v3qiT7kzzc\nrf/pPnNKkha3tq8nTrIG2AVcB8wC+5LsrqqZkTFTwO3A1qo6keSy7q6vAG+sqr9J8mLg891jj/aV\nV5J0Zn3uWVwNHKqqw1V1ErgbuGHemJuBXVV1AqCqjnW/nqyqv+nGvLDnnJKkhj7/E94AHBlZnu3W\njdoCbEnyUJK9SbbN3ZFkU5JHuuf4pYX2KpLsSDKdZPqJJ57o4bcgSYJ+yyILrKt5y2uBKeAa4EeB\nDyR5GUBVHamq7wa+A3hHkpd/y5NV3VFVg6oaTExMLGt4SdIz+iyLWWDTyPJGYP7ewSxwX1WdqqrH\ngYMMy+Np3R7FAeBNPWaVJC2iz7LYB0wluTzJOmA7sHvemHuBawGSrGc4LXU4ycYkL+rWXwpsZVgk\nkqRV0FtZVNVp4BZgD/AYcE9VHUiyM8n13bA9wPEkM8D9wK1VdRx4DfDpJJ8DHgB+uaoe7SurJGlx\nqZp/GOH8NBgManp6erVjSNJ5Jcn+qhq0xnlKqiSpybKQJDVZFpKkpkXLIsnHR27f3n8cSdI4au1Z\njH7S7R/2GUSSNL5aZXFhnColSXpOWledfXWS3Qwv3TF3+2lVdf3CD5MkXUhaZTF6ldhf7jOIJGl8\nLVoWVfXA3O0kE906L+8qSc8zi5ZFkgC/CLyL4VTUC5KcBn6tqnauQL4V8a8/coCZo19f7RiSdE6u\neOUlvOcffFevr9E6wP1u4AeA762qb6+qS4HvA7Ym+blek0mSxsai14ZK8lnguqr66rz1E8DHq+p1\nPedbMq8NJUlnb7muDXXR/KKAp49bXHSu4SRJ55dWWZw8x/skSReQ1qmz35Pk6zzzFalzc1YBLu4t\nlSRprLROnV2zUkEkSeOrdersxcBPA98BPALc2X0DniTpeaR1zOK3gQHwKPA24H29J5IkjZ3WMYsr\nqurvAiT5TeBP+48kSRo3rT2LU3M3nH6SpOevpZ4NBcMzoF40cnZUVdUlvaaTJI0Fz4aSJDX1+h3c\nSbYlOZjkUJLbzjDmpiQzSQ4kuatb99okn+rWPZLkH/WZU5K0uNY01DlLsgbYBVwHzAL7kuyuqpmR\nMVPA7cDWqjqR5LLurr8C3l5VX0zySmB/kj1V9WRfeSVJZ9bnnsXVwKGqOlxVJ4G7efaXKQHcDOyq\nqhMAVXWs+/ULVfXF7vZR4BjP/j5wSdIK6rMsNgBHRpZnu3WjtgBbkjyUZG+SbfOfJMnVwDrgS70l\nlSQtqrdpKJ65ntSo+ddDXwtMAdcAG4EHk1w5N92U5BXA7wDvqKpvfssLJDuAHQCTk5PLl1yS9Cx9\n7lnMAptGljcCRxcYc19Vnaqqx4GDDMuDJJcA/wP4l1W1d6EXqKo7qmpQVYOJCWepJKkvfZbFPmAq\nyeVJ1gHbgd3zxtwLXAuQZD3DaanD3fj/Dnywqn6vx4ySpCXorSy6T3zfAuwBHgPuqaoDSXYmub4b\ntgc4nmQGuB+4taqOAzcBbwbemeTh7ue1fWWVJC1u0a9VPZ/4taqSdPaW62tVJUmyLCRJbZaFJKnJ\nspAkNVkWkqQmy0KS1GRZSJKaLAtJUpNlIUlqsiwkSU2WhSSpybKQJDVZFpKkJstCktRkWUiSmiwL\nSVKTZSFJarIsJElNloUkqcmykCQ1WRaSpCbLQpLUZFlIkpp6LYsk25IcTHIoyW1nGHNTkpkkB5Lc\nNbL+Y0meTPLRPjNKktrW9vXESdYAu4DrgFlgX5LdVTUzMmYKuB3YWlUnklw28hT/Afg24J/2lVGS\ntDR97llcDRyqqsNVdRK4G7hh3pibgV1VdQKgqo7N3VFVfwj8ZY/5JElL1GdZbACOjCzPdutGbQG2\nJHkoyd4k23rMI0k6R71NQwFZYF0t8PpTwDXARuDBJFdW1ZNLeoFkB7ADYHJy8tyTSpIW1eeexSyw\naWR5I3B0gTH3VdWpqnocOMiwPJakqu6oqkFVDSYmJp5zYEnSwvosi33AVJLLk6wDtgO75425F7gW\nIMl6htNSh3vMJEk6B72VRVWdBm4B9gCPAfdU1YEkO5Nc3w3bAxxPMgPcD9xaVccBkjwI/B7wg0lm\nk7y1r6ySpMWlav5hhPPTYDCo6enp1Y4hSeeVJPuratAa5ye4JUlNloUkqcmykCQ1WRaSpCbLQpLU\nZFlIkposC0lSk2UhSWqyLCRJTZaFJKnJspAkNVkWkqQmy0KS1GRZSJKaLAtJUpNlIUlqsiwkSU2W\nhSSpybKQJDVZFpKkJstCktRkWUiSmnotiyTbkhxMcijJbWcYc1OSmSQHktw1sv4dSb7Y/byjz5yS\npMWt7euJk6wBdgHXAbPAviS7q2pmZMwUcDuwtapOJLmsW/+3gPcAA6CA/d1jT/SVV5J0Zn3uWVwN\nHKqqw1V1ErgbuGHemJuBXXMlUFXHuvVvBT5RVV/r7vsEsK3HrJKkRfRZFhuAIyPLs926UVuALUke\nSrI3ybazeKwkaYX0Ng0FZIF1tcDrTwHXABuBB5NcucTHkmQHsANgcnLyuWSVJC2izz2LWWDTyPJG\n4OgCY+6rqlNV9ThwkGF5LOWxVNUdVTWoqsHExMSyhpckPaPPstgHTCW5PMk6YDuwe96Ye4FrAZKs\nZzgtdRjYA7wlyaVJLgXe0q2TJK2C3qahqup0klsY/ie/Brizqg4k2QlMV9VunimFGeAp4NaqOg6Q\n5L0MCwdgZ1V9ra+skqTFpepbDgWclwaDQU1PT692DEk6ryTZX1WD1jg/wS1JarIsJElNloUkqcmy\nkCQ1WRaSpCbLQpLUZFlIkposC0lSk2UhSWqyLCRJTZaFJKnJspAkNVkWkqQmy0KS1GRZSJKaLAtJ\nUpNlIUlqsiwkSU2WhSSpybKQJDVZFpKkJstCktRkWUiSmiwLSVKTZSFJakpVrXaGZZHkCeDLZ/mw\n9cBXe4jzXI1rLjDbuTLb2RvXXHBhZXtVVU20Bl0wZXEukkxX1WC1c8w3rrnAbOfKbGdvXHPB8zOb\n01CSpCbLQpLU9HwviztWO8AZjGsuMNu5MtvZG9dc8DzM9rw+ZiFJWprn+56FJGkJLviySLItycEk\nh5Lctsi4G5NUkhU7w2Ep2ZLclGQmyYEkd41LtiSTSe5P8tkkjyR52wrlujPJsSSfP8P9SfKrXe5H\nkrx+JXItMduPdZkeSfInSb5nXLKNjPveJE8luXGcsiW5JsnD3XbwwLhkS/LSJB9J8rku20+uUK5N\n3fb3WPe6P7vAmOXdFqrqgv0B1gBfAl4NrAM+B1yxwLiXAH8M7AUG45INmAI+C1zaLV82RtnuAP5Z\nd/sK4M9WKNubgdcDnz/D/W8D/icQ4PuBT6/gv7dWtjeO/F3+0DhlG/l7/yPgD4AbxyUb8DJgBpjs\nlldkO1hitl8Afqm7PQF8DVi3ArleAby+u/0S4AsLbKPLui1c6HsWVwOHqupwVZ0E7gZuWGDce4F/\nD/y/Mct2M7Crqk4AVNWxMcpWwCXd7ZcCR1ciWFX9McMN8kxuAD5YQ3uBlyV5xThkq6o/mfu7ZPjG\nZONK5Opeu/XnBvAu4L8BK/XvDFhStn8MfLiq/rwbv2L5lpCtgJckCfDibuzpFcj1lar6THf7L4HH\ngA3zhi3rtnChl8UG4MjI8izz/kCTvA7YVFUfXclgLCEbsAXYkuShJHuTbBujbP8K+PEkswzfib5r\nZaI1LSX7OPgphu/6xkKSDcCPAL+x2lkWsAW4NMknk+xP8vbVDjTi/cBrGL5ZehT42ar65koGSLIZ\neB3w6Xl3Leu2sPZcH3ieyALrnj79K8kLgP8IvHOlAo1YNFtnLcOpqGsYvgt9MMmVVfXkGGT7UeC3\nqup9Sd4A/E6XbUU3lAUsJfuqSnItw7L4gdXOMuI/AT9fVU8N3ySPlbXAVcAPAi8CPpVkb1V9YXVj\nAfBW4GHg7wF/G/hEkger6usr8eJJXsxwb/DdC7zmsm4LF/qexSywaWR5I8+eLnkJcCXwySR/xnBe\nb/cKHeRuZZsbc19Vnaqqx4GDDMtjHLL9FHAPQFV9CriY4TVpVttSsq+aJN8NfAC4oaqOr3aeEQPg\n7m47uBH49SQ/vLqRnjYLfKyq/m9VfZXh8cUVOzmg4ScZTpFVVR0CHgf+zkq8cJKLGBbF71bVhxcY\nsqzbwoVeFvuAqSSXJ1kHbAd2z91ZVX9RVeuranNVbWY4j3x9VU2vdrbOvcC1AEnWM9wdPzwm2f6c\n4Ts9kryGYVk8sQLZWnYDb+/OBPl+4C+q6iurHQqGZ5ABHwZ+YkzeFT+tqi4f2Q5+H/iZqrp3lWPN\nuQ94U5K1Sb4N+D6Gc/TjYHQ7eDnwnazANtodI/lN4LGq+pUzDFvWbeGCnoaqqtNJbgH2MDzT486q\nOpBkJzBdVfP/Axy3bHuAtySZAZ4Cbl2Jd6NLzPbPgf+S5OcY7tq+s7pTMPqU5EMMp+XWd8dL3gNc\n1OX+DYbHT94GHAL+iuE7vxWxhGy/CHw7w3ftAKdrhS5Gt4Rsq6aVraoeS/Ix4BHgm8AHqmrRU4BX\nKhvDk2N+K8mjDKd9fr7b++nbVuAngEeTPNyt+wVgciTbsm4LfoJbktR0oU9DSZKWgWUhSWqyLCRJ\nTZaFJKnJspAkNVkW0lnoLjnx1nnr3p3k1xd5zDf6Tyb1y7KQzs6HGH5IcdT2br10wbIspLPz+8Df\nT/JCePoibq8EHk7yh0k+k+TRJN9ydePuOxk+OrL8/iTv7G5fleSB7kJ5e1bqSrnSUlkW0lnoPkH/\np8DcFYC3A/8V+GvgR6rq9Qwv0fK+LPGKfN01fn6N4XdIXAXcCfyb5c4uPRcX9OU+pJ7MTUXd1/36\nTxhe6uHfJnkzw0tSbABeDvyfJTzfdzK8oOUnun5ZA4zF9aykOZaFdPbuBX6l+5rKF1XVZ7rppAng\nqqo61V299eJ5jzvNs/fm5+4PcKCq3tBvbOncOQ0lnaWq+gbwSYbTRXMHtl8KHOuK4lrgVQs89MvA\nFUlemOSldFcrZXjp+Ynue0FIclGS7+rz9yCdLfcspHPzIYaXG587M+p3gY8kmWb4ZTj/a/4DqupI\nknsYXj31iwy/X52qOpnkRuBXuxJZy/DLiA70/ruQlsirzkqSmpyGkiQ1WRaSpCbLQpLUZFlIkpos\nC0lSk2UhSWqyLCRJTZaFJKnp/wMjjgQDZ5nPJgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fd6590424a8>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "gamma = Gamma(alpha=5, beta=8)\n",
    "gamma.plot()\n",
    "\n",
    "uniform = Uniform(lower=0.4, upper=2.0)\n",
    "uniform.plot()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "They can also represent a set of coefficients. For example, the Gaussian channel weighting module takes the set of input channels, weights them using a Gaussian function, and then sums the weighted channels to produce a single output channel. The Gaussian is specified by a mean, `mu`, and standard deviation, `sd`.\n",
    "\n",
    "The channel weighting module can produce multiple outputs (each weighted using a Gaussian function). However, we would expect that the output channels are somewhat evenly distributed across the rang eof input channels. Assuming that the input channels are normalized to the range [0, 1] where 0 is the lowest-frequency input channel and 1 is the highest-frequency input channel, we can define a set of priors that indicate that we want `mu` for the first input channel to fall at the lower end of the range and `mu` for the last input channel to fall at the upper end of the range."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{'mu': Beta(α=array(3), β=array(3)), 'sd': HalfNormal(σ=array(3))}\n"
     ]
    }
   ],
   "source": [
    "# Initialize the module with three outputs\n",
    "wc = WeightChannelsGaussian(3)\n",
    "\n",
    "# For this particular module, the generation of priors does not\n",
    "# require any actual data, so we can get away with providing fake\n",
    "# data.\n",
    "wc_priors = wc.get_priors(signals)\n",
    "print(wc_priors)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "As shown above, the weight channels module returned a set of two priors (one for `mu` and one for `sd`). Each set of priors is an array of shape `(3,)` (i.e., a 1D array of 3 elements). Each entry in the array corresponds to one of the output channels."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x7fd658fd3320>"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEKCAYAAAD9xUlFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzs3XdcleX/x/HXxd64QEEJHLi3qLgq\nV5o5SnOmLcv8plm/yrJvaZbfdmmlZlpWpmlmmVlquUpN08Qt4sCtoCAIyIZzrt8fN5kZLuRwnwOf\n5+NxHpxxnXO/FeXDfV9Laa0RQgghAJzMDiCEEMJ+SFEQQghxkRQFIYQQF0lREEIIcZEUBSGEEBdJ\nURBCCHGRFAUhhBAXSVEQQghxkRQFIYQQF7mYHeBGVapUSYeFhZkdQwghHMq2bdvOaa0DrtXO4YpC\nWFgYUVFRZscQQgiHopQ6fj3t5PKREEKIi6QoCCGEuEiKghBCiItsXhSUUs5KqR1KqZ8Kec1dKbVQ\nKRWrlNqilAqzdR4hhBBXVhJnCk8CMVd4bThwXmtdC5gCvFUCeYQQQlyBTYuCUqoacBfw6RWa9AHm\nFNz/FuislFK2zCSEEOLKbH2m8D7wHGC9wutVgZMAWut8IBWoaONMQgghrsBmRUEp1RNI0Fpvu1qz\nQp771/6gSqkRSqkopVRUYmJikfKcTDvJW3++RZ41r0jvF0IIM32w+hDbjifb/Di2PFNoB/RWSh0D\nvgY6KaXmXdbmFBACoJRyAfyBf/2ptdaztNYRWuuIgIBrTsgr1JHUI8yLmceS2CVFer8QQphl96kU\npqw+yO+Hkmx+LJsVBa31C1rralrrMGAQsFZrPfSyZkuBBwru31vQ5l9nCsXh1mq30iSgCR/v/Jjs\n/GxbHEIIIWzi3ZUHKe/lysPtw2x+rBKfp6CUelUp1bvg4WygolIqFngaGGfD4/Jk8ydJyEpg4YGF\ntjqMEEIUqz+PJrP+YCIjb6uJr4erzY9XImsfaa1/A34ruD/hkuezgf4lkQGgZZWWtA1uy6d7PqVf\neD983HxK6tBCCHHDtNa8+8sBAnzdub9NWIkcs8zNaB7TfAwpOSl8ue9Ls6MIIcRVrT90jj+PJfNE\np1p4ujmXyDHLXFFoULEBXUO7Mid6Duezz5sdRwghCmW1at7+eT/VynsysGVIiR23zBUFgNFNR5Nt\nyWb2ntlmRxFCiEIt2xNPdFwaT3etjbtLyZwlQBktCjXK1aBXjV4s2L+AMxlnzI4jhBD/kGex8t7K\nA9St4kufplVL9NhlsigA/Kfpf7BiZebumWZHEUKIf1i49STHkjIZ260Ozk4lu/JPmS0KVX2q0r92\nf74/9D0n0k6YHUcIIQDIzM3ngzWHaBlWnk51A0v8+GW2KACMaDwCN2c3pu2cZnYUIYQA4PONx0i8\nkMPz3etixvqgZbooVPKsxH317mPF0RXsT95vdhwhRBmXkpnLx+sO06VeIBFhFUzJUKaLAsBDDR/C\nz82P97e9b3YUIUQZN+O3w6Tn5DO2W13TMpT5ouDn5seIxiPYGLeRTXGbzI4jhCij4lOz+GLTMe5p\nVpU6VXxNy1HmiwLA4LqDqepTlSnbpmDVV9r6QQghbOeD1YfQGv6vS21Tc0hRANyc3Xii2RPsT97P\nsiPLzI4jhChjYhPS+SbqJPdF3kJIBS9Ts0hRKHBn9TupV6EeU3dMJceSY3YcIUQZ8vbP+/F0dWZ0\nx1pmR5Gi8Bcn5cQzEc8QnxHPgpgFZscRQpQRm48ksXLfWR7vWIuKPu5mx5GicKnWQa1pX7U9s/bM\nIjUn1ew4QohSzmrVvL48hiB/Dx5uV93sOIAUhX/5vxb/R3puOp/s/sTsKEKIUu7H3XHsPpXK2G51\nSmxp7GuxWVFQSnkopf5USu1SSkUrpV4ppM2DSqlEpdTOgtsjtspzvWqXr02fWn2Yv38+p9NPmx1H\nCFFKZedZePvnAzQI9uPuEl707mpseaaQA3TSWjcBmgLdlVKRhbRbqLVuWnD71IZ5rtuopqNwUk5M\n3THV7ChCiFLq843HOJ2SxYt31cOphBe9uxqbFQVtSC946Fpw07Y6XnGq4l2FYfWHsezIMqKTos2O\nI4QoZZLSc/jo11i61Aukbc1KZsf5B5v2KSilnJVSO4EEYJXWekshzfoppXYrpb5VSpXc9kLX8HDD\nh6ngUYG3/3wbrR2ilgkhHMQHaw6RmWdh3J3mLWdxJTYtClpri9a6KVANaKWUanhZkx+BMK11Y2A1\nMKewz1FKjVBKRSmlohITE20Z+SJfN19GNR3F9oTtrDq+qkSOKYQo/Q4npvPVlhMMbhVCrUDzlrO4\nkhIZfaS1TgF+A7pf9nyS1vqvmWKfAC2u8P5ZWusIrXVEQECATbNeqm94X8LLhzN522SZ0CaEKBZv\nLDcmqj1l8nIWV2LL0UcBSqlyBfc9gS7A/svaBF3ysDcQY6s8ReHi5MJzLZ/jdPpp5u2bZ3YcIYSD\n++NwEqtjzvKf22tSyQ4mqhXGlmcKQcCvSqndwFaMPoWflFKvKqV6F7QZUzBcdRcwBnjQhnmKJDIo\nkttDbueTPZ9wLuuc2XGEEA7KYtVM+mkfwf4eDG9vHxPVCmPL0Ue7tdbNtNaNtdYNtdavFjw/QWu9\ntOD+C1rrBlrrJlrrjlpru9zp5tmIZ8mx5DBth+zQJoQomm+iTrIvPo0XetTDw9U+JqoVRmY0X4dQ\nv1CG1B3C4kOLZYc2IcQNS83K451fDtAyrDw9Gwdd+w0mkqJwnR5r8hj+7v689edbMkRVCHFDPlxz\niPOZubzcq4Ep+y7fCCkK18nPzY/RTUcTdTaKtSfWmh1HCOEgYhPSmbPpGINahtCwqr/Zca5JisIN\n6Fe7H7XK1eLdqHfJteSaHUcIYee0NjqXPV2deeaOOmbHuS5SFG6Ai5MLYyPGcir9FHP3zTU7jhDC\nzv16IIF1BxN5sku43Q5BvZwUhRvUtmpbbg+5nZm7Z3I246zZcYQQdio338qkn2KoEeDN/W3CzI5z\n3aQoFMFzLZ/DYrXw3rb3zI4ihLBTX2w6ytFzGYzvWR83F8f5Ues4Se1IiG8IwxsNZ8XRFWw9s9Xs\nOEIIO5N4IYcP18TSsU4AHesEmh3nhkhRKKKHGz5MVZ+qvL7ldfKseWbHEULYkbd/3k92noXxPeub\nHeWGSVEoIg8XD8a2HEtsSiwL9y80O44Qwk5sO57Mom2nGN6+OjUCfMyOc8OkKNyETiGdaBfcjuk7\np8u6SEII8i1WXloSTRU/D8Z0Djc7TpFIUbgJSinGtRpHtiWbKdummB1HCGGyuZuPExOfxvie9fF2\ndzE7TpFIUbhJYf5hPFD/AZYeXsrOhJ1mxxFCmCQhLZvJKw/SIbwSPRpVMTtOkUlRKAYjGo8g0CuQ\n17e8jsVqMTuOEMIEry+PISffyiu97X99o6uRolAMvFy9GBsxlpjkGBYdXGR2HCFECfvjcBJLdsYx\n4tYaDtm5fCkpCsWkW1g3Wge15sPtH0qnsxBlSJ7FyoQf9lK1nCejOtYyO85Ns+V2nB5KqT+VUrsK\ndld7pZA27kqphUqpWKXUFqVUmK3y2JpSipdav0S2JZu3/3zb7DhCiBLy+cajHEpIZ2LvBni62e/m\nOdfLlmcKOUAnrXUToCnQXSkVeVmb4cB5rXUtYArwlg3z2FyYfxiPNHqEFcdWsOn0JrPjCCFsLD41\ni/dXH6Jz3UC61q9sdpxiYcvtOLXWOr3goWvB7fLdafoAcwrufwt0Vo7cQwMMbzScUL9QJm2eRHZ+\nttlxhBA2NOmnfVismom9G5gdpdjYtE9BKeWslNoJJACrtNZbLmtSFTgJoLXOB1KBirbMZGvuzu6M\njxzPqfRTzNo9y+w4QggbWRNzluV7zjC6Yy1CKniZHafY2LQoaK0tWuumQDWglVKq4WVNCjsr+Nde\nl0qpEUqpKKVUVGJioi2iFqvWQa3pWaMnn0d/zpGUI2bHEUIUs4ycfMYv2Uvtyj48dltNs+MUqxIZ\nfaS1TgF+A7pf9tIpIARAKeUC+APJhbx/ltY6QmsdERAQYOO0xePZiGfxcvHi1c2vyp7OQpQy7608\nSFxqNm/0beRQy2JfD1uOPgpQSpUruO8JdAH2X9ZsKfBAwf17gbW6lPwErehZkf9r8X9sO7uNJbFL\nzI4jhCgmu06m8MWmowyNvIUWoRXMjlPsbFnigoBflVK7ga0YfQo/KaVeVUr1LmgzG6iolIoFngbG\n2TBPiesb3pdmgc2YvG0y57PPmx1HCHGT8ixWxi3eQyUfd57rXtfsODZhy9FHu7XWzbTWjbXWDbXW\nrxY8P0FrvbTgfrbWur/WupbWupXWulRdgHdSToyPHE96bjrvRr1rdhwhxE367PejxMSn8WqfBvh5\nuJodxyZK18UwOxRePpyHGj7E0sNLZe6CEA7sRFImU1YfpGv9ynRr4LgL3l2LFIUS8FiTx6juX51X\n/niFzLxMs+MIIW6Q1poXl+zBWSmHX/DuWqQolAB3Z3deafsK8RnxfLD9A7PjCCFu0NJdcWw4dI6x\n3eoQXM7T7Dg2JUWhhDQLbMbguoNZsH8BOxJ2mB1HCHGdktJzeOXHfTQJKcewNmFmx7E5KQol6Mnm\nTxLkHcTLm14mx5JjdhwhxHWYsDSaC9l5vN2vMc5Opfey0V+kKJQgL1cvXm7zMkdTjzJz10yz4wgh\nruHnvfEs2x3PmE7h1Knia3acEiFFoYS1rdqWPjX78Pnez9mffPlcPiGEvTifkctLS/bSINiPkbeX\nrqUsrkaKggnGthyLv7s/EzZOIN+ab3YcIUQhXvkxmpTMPN65twmuzmXnR2XZ+ZPaEX93f16KfImY\n5BjmRM+59huEECVq9b6zLNkZx6iOtagf7Gd2nBIlRcEkXUK70DW0K9N3Tif2fKzZcYQQBVIz8/jv\n93uoW8W3VGyveaOkKJjopciX8HXz5b+//5c8a57ZcYQQwKs/7SMpI5d3+zcpdSugXo+y9ye2IxU8\nKjChzQRikmP4dPenZscRosz79UAC320/xX9uq0nDqv5mxzGFFAWTdb6lMz1r9GTW7llEJ0WbHUeI\nMis1K48XvttD7co+PNG57F02+osUBTswrtU4KnhU4KXfX5JJbUKY5OUf9nIuPYd3+zfB3cXZ7Dim\nkaJgB/zd/Xml3SvEpsQyfed0s+MIUeYs2x3Pkp1xPNEpnMbVypkdx1RSFOxE+6rtubf2vXyx9wt2\nJuw0O44QZUZCWjYvLtlDk5ByjOpYdiapXYktt+MMUUr9qpSKUUpFK6WeLKTN7UqpVKXUzoLbBFvl\ncQTPRjxLsE8wL/7+oiyxLUQJ0Frz3He7yc6zMHlAE1zK0CS1K7Hl30A+8IzWuh4QCYxSStUvpN0G\nrXXTgturNsxj97xdvZnUbhInLpxgyrYpZscRotT7assJfjuQyAt31qNmgI/ZceyCLbfjjNdaby+4\nfwGIAara6nilRcsqLRlWfxhfH/ia9afWmx1HiFLr6LkMXlsWQ4fwSgyLDDU7jt0okXMlpVQY0AzY\nUsjLbZRSu5RSK5RSDa7w/hFKqSilVFRiYqINk9qHJ5s/SXj5cMZvHE9SVpLZcYQodfItVp7+Zieu\nzop37m2CUxlYEvt62bwoKKV8gO+Ap7TWaZe9vB0I1Vo3AaYCSwr7DK31LK11hNY6IiAgwLaB7YC7\nsztvdXiL9Nx0JmyagNba7EhClCofrzvMjhMpTLq7IVX8PcyOY1dsWhSUUq4YBeErrfXiy1/XWqdp\nrdML7i8HXJVSlWyZyVGElw/n6YinWX9qPQsPLDQ7jhClxo4T55my+hC9mgTTp6lc0b6cLUcfKWA2\nEKO1nnyFNlUK2qGUalWQR66XFBhSdwgdqnbg3ah3ZdE8IYrBhew8xny9gyp+Hvzv7oZmx7FLtjxT\naAcMAzpdMuS0h1JqpFJqZEGbe4G9SqldwIfAIC3XSi5SSvFqu1fxdvXm+Q3Pk2vJNTuSEA5t/JK9\nxKVk8+Hgpvh7upodxy652OqDtda/A1ftvdFaTwOm2SpDaVDJsxKT2k1i1JpRfLD9A8a2HGt2JCEc\n0uLtp1iyM46nu9amRWgFs+PYLZmp4QBurXYrg+oM4st9X7IpbpPZcYRwOMfOZTB+yV5aVa9QJvdI\nuBFSFBzEMxHPUMO/Bi/+/qIMUxXiBuTmW3ny6x24ODvx/sCmOMvw06uSouAgPFw8eOe2d7iQe4H/\n/v5frNpqdiQhHMLkVQfZdSqVt/o1Iricp9lx7J4UBQdSu3xtxrUax6a4TczeM9vsOELYvd8PnWPm\n+sMMbnUL3RsGmR3HIUhRcDD9wvtxZ/U7mbZzGlFnosyOI4TdSkjL5qmFO6gV4MOEnoUtuyYKI0XB\nwSileLnNy4T4hvD8+udJzk42O5IQdiffYuWJBTvIyLHw0X3N8XQru5vm3CgpCg7I29Wbd297l5Sc\nFOlfEKIQ768+xJajybx2T0PCK/uaHcehSFFwUHUr1OW5ls+x8fRGPt/7udlxhLAbvx5IYNqvsQyM\nCKFv82pmx3E4UhQc2IA6A7gj9A6m7pjKjoQdZscRwnRxKVk8vXAndav48kqfQhddFtdgsxnNwvaU\nUkxsO5F9SfsYu24s3/T6hgoeMlPTprLT4OxeSD4K54/ChTOQkwbZqaCt4OQKzq7gVRF8AsE3GALq\nQGB98Cn9K/yaKc9iZfT87eTmW/novuZ4uEo/QlFIUXBwvm6+vHf7ewxbPozn1j/Hx10+xsVJvq3F\nJus8HF4LsWvh1FY4dxAoWJ5LOYF3IHj4g4cfKGewZoIlB+J3Q0YCWPP//izfYAhta9xqdYHysrFL\ncXr75/1sP5HC1MHNqCG7qBWZ/PQoBepXrM9LkS8xYdMEpu2YxlMtnjI7kmPLSoHoxbDnWzixGbQF\nPMtDSCQ0uheCmkLFmlDuFuOs4EqsVkg/C4n7ISEGTkfBsd9h77fG60FNoV4vaNRfCsRN+iX6DJ9s\nOMqwyFB6NQk2O45Dk6JQStwTfg97zu1h9t7ZNKrUiM6hnc2O5HhO/glbPoaYn4zf9gPqQvv/g9rd\noGoLcLrByxFOTuAXZNxqdjSe0xqSYuHActi3FNZOgrX/M84cIh6C8G7gLP8tb0RsQjrPfLOLJtX8\nefGuembHcXjqaitVK6VWaq3vKLj/gtb6jRJLdgURERE6KkombRUm15LLgz8/yJHUI8y/az41/GuY\nHcn+Wa1wYBlsmgontxiXghoPhCaDIbgZKBuvk5NyEnbMg+1z4EI8lA8zClGTweDibttjlwIXsvPo\nM30jqZl5/PhEe1nG4iqUUtu01hHXbHeNorBDa92s4P52rXXzYsxYJFIUru5MxhkG/DiA8h7lWXDX\nArxcvcyOZJ+0htg1sGYinNkD5UKhzShoeh+4m3A92pJvnD38PgXitoNfVejwDDR/QM4crsBq1Yyc\nt401+xOYN7w1bWpWNDuSXbveonCtIamy4Y2DqeJdhXdue4djaccYv3G87O9cmPhd8EVP+KqfMZro\nnlnwxHZo/Zg5BQGMH/z1e8Oja2HoYvAPgWVPw8ft4OBKo4iJf5ix7jAr953lhTvrSkEoRtcqCjWU\nUkuVUj9ecv/i7WpvVEqFKKV+VUrFKKWilVJPFtJGKaU+VErFKqV2K6VMPxMpDVoHtebJ5k+y8vhK\nvtz3pdlx7Ed2Kqx4HmbdbnT+9ngXRkdBk4H289u4UlCrMzz8Mwz8Ciy5ML8/fNUfzh83O53d+O1A\nAu+uPECfpsEMb1/d7DilyrX+J/S55P67N/jZ+cAzWuvtSilfYJtSapXWet8lbe4EwgturYEZBV/F\nTXqowUPsSdzDlG1TqF2+Nm2C25gdyVzR3xsFIT0BWg6HTi8ZI4rslVJQryeE3wFbP4G1r8FHkUbu\n1iNvvNO7FDmRlMmTX++kTmVf3uzbGGXrfp8y5qp9Cv9oqFQAgNY6sUgHUuoHYJrWetUlz80EftNa\nLyh4fAC4XWsdf6XPkT6F65eRl8HQ5UNJyExg/l3zCfUrg8Mes87DsmeNYaBBTaHnFKjqgCekKSeN\ny0mHVkJwc+j7CVQqezuIZebm0/ejTcSnZvPj6PbcUlH6zK5XsfQpFFzeeVkpdQ7YDxxUSiUqpSbc\nYJgwoBmw5bKXqgInL3l8quC5y98/QikVpZSKSkwsUk0qk7xdvfmw04c4KSfGrB3DhdwLZkcqWYfX\nwkdtYd8S6PgiPLLGMQsCQLkQGPIN9JttzKSe2QG2zy1TfQ1Wq+bphbs4ePYCHwxqKgXBRq7Vp/AU\n0B5oqbWuqLUuj3F5p51S6v+u5wBKKR/gO+AprXXa5S8X8pZ//SvXWs/SWkdorSMCAmSpgBsR4hvC\n5NsncyLtBM+vfx6L1WJ2JNuzWmDNJJh7D7j7wvBVcNtz9tNvUFRKGZPn/rPJmDexdDQsesA4GyoD\npqw+yM/RZ/hvj3rcXifQ7Dil1rWKwv3AYK310b+e0FofAYYWvHZVSilXjILwldZ6cSFNTgEhlzyu\nBsRd63PFjWlZpSXjWo1jw+kNfLD9A7Pj2FZ6Asy9Gza8C82GwWPrHPfs4Er8guH+H6DLRNi/DGZ1\nhLP7rvUuh/bDztNMXWusfCody7Z1raLgqrU+d/mTBf0KV5nfb1x6AmYDMVrryVdothS4v+AyVSSQ\nerX+BFF0A+sOZGCdgXwe/Tk/Hv7R7Di2cWILzLzVmJncZzr0mQaupXQyk5OzMcntweWQlwmfdoa9\nhf3e5fh2nkxh7Le7aRVWgUl3N5SOZRu7VlHILeJrAO2AYUAnpdTOglsPpdRIpdTIgjbLgSNALPAJ\n8Pj1hBZF83yr52lZpSUTN01kd+Jus+MUr10LYU5PcPEwLhc1G2p2opJxS2t4bD1UaQTfPgQrxxuX\nz0qJM6nZjPgyikBfd2YMbY6bi6z2b2vXmtFsATL4+9r/X40V4KG1vurZgi3I6KObcz77PIOXDSY7\nP5v5d80n2MfBFw/TGn59Hda/DWEdYMCX4FUGlw/Pz4Wfx0HUbKjb0xid5ObYHbFZuRYGzPyDI4np\nLH68HXWqyA5qN6NYRh9prZ211n5aa9+Cm98lj0u8IIibV96jPNM7TyfXksuoNaMce0RSXhZ8+7BR\nEJoNNWYCl8WCAODiBj0nQ/e3jH6GOb0g3XFH6lmtmmcX7WJvXCofDm4mBaEEXWtIqodS6iml1LSC\nYaEOPnxDANQsV5PJHSdzLPUYT//2NHnWPLMj3bisFGN0UfRio8O19zTjB2NZFzkSBs6Ds9Ewuwuc\nO2R2oiJ565f9LNsTzwt31qVzvcpmxylTrnWBbg4QAewBegDv2TyRKBGRQZG83PZlNsdvZtIfkxxr\njaQLZ421i05Fwb2fGx2u0vn4t3o94cGfICcdZneF09vMTnRD5m4+zsx1RxgWGcqjHWSl35J2raJQ\nX2s9VGs9E7gX6FACmUQJubvW3YxoPILvY7/n0z2fmh3n+pw/Bp91g+TDMGQhNOxrdiL7VC0CHlkF\n7n4wpw8c22h2ouuyJuYsL/+wl851A3m5V30ZaWSCaxWFi9cVtNb5V2soHNPopqPpUb0HH+74kOVH\nlpsd5+oSYuCz7sZkrft/MBaOE1dWoYaxsJ5fEMzrC4dWm53oqnafSmH0/B00CPZn6pBmuDjLSCMz\nXOtvvYlSKq3gdgFo/Nd9pdTls5OFA1JKMandJJoHNueljS+x/ex2syMV7swe+LyHMdrooRUQ0srs\nRI7BL9j4+6pUGxYMgn0/mJ2oUCeTM3n4iygqeLsx+8EIvNyk+9Is1zv66K8RRy6X3PcrqZDCttyc\n3fig4wdU9anKE2ufIPZ8rNmR/unMHpjT25iI9vAKqFzf7ESOxbsSPPCjMbN70UPGirF2JDUzj4e+\n2EpuvoU5D7ck0NfD7EhlmpyfCQDKeZRjRpcZuDm7MXL1SM5knDE7kuHM3r8LwoM/GZdExI3zLGcM\n2Q1pBd8ON/aHtgPZeRYenRvFiaRMZt0fQa1AGXpqNikK4qJqvtX4uMvHZORl8Niqx0jJTjE30Nlo\n+LK3MUv5gR+lINwsdx+4b5GxmN63DxnzGUyUb7Eyev4Oth5L5r0BTYisIbun2QMpCuIf6lSow9RO\nUzl14RSj1o4iMy/TnCBn9xkTsJzdjTOEijXNyVHauPvC0G8hqAl88wAc/MWUGFprXli8h9UxZ3m1\ndwN6NXHwmfWliBQF8S8RVSJ4+7a32XtuL8+ue7bkJ7clHzVWOnV2k4JgCx7+xqWkyg1g4VCILflR\nSW+u2M+ibad4snM4w9qElfjxxZVJURCF6nxLZ16KfIkNpzcwcdPEkpvclhYPX/Yx9iYetkQKgq14\nloNh30NAHVg4zFhZtoTMXHeYmeuPcH+bUJ7qEl5ixxXXR4qCuKL+tfszqukolh5eyjtR79i+MGQm\nG+PpM5Pgvu8gsK5tj1fWeVUwzhh8q8BX/UtkT4ZFUSd5Y8V+ejYOYmKvBjI5zQ5JURBX9Vjjxxha\nbyhz981l+s7ptjtQTjrMHwBJsTBoPlRrYbtjib/5BBpnZK6exlpS54/Z7FAro88wbvEeOoRXYvKA\npjg5SUGwR1IUxFUppXiu5XP0C+/HzN0zmb1ndvEfJD/HuLZ9epuxllGN24r/GOLKyocaZwz52UZh\nSE8o9kP8diCB0fN30KiqPx8PbSH7Itgx+c6Ia1JKMT5yPD2q9+D97e/zVcxXxffhVit8/xgc+dVY\n6bRez+L7bHH9Ktc3hqteOGNcwstOLbaP3hR7jsfmbiO8sg9zHm6Ft7vMVrZnNisKSqnPlFIJSqm9\nV3j9dqVU6iW7sk2wVRZx85ydnPlf+//RKaQTb/75Jt8fKqZZsavGGzNsu06CZvcVz2eKoglpBQPn\nGmtMLRxqbNxzk7YeS2b4nCjCKnozd3hr/D1lGxZ7Z8szhS+A7tdos0Fr3bTg9qoNs4hi4Orkyju3\nvUO74Ha8vOllVhxdcXMfuGUm/DENWj0GbZ8onpDi5tTqYpyxHV0PP44x1poqop0nU3jo860ElfNg\n3iOtqeAt+104ApsVBa31eiDZVp8vzOHm7MaUjlNoXrk5L2x4gVXHVxXtg2J+ghXPQ527oPsbsh+C\nPWk6GDq+CLsWwG9vFOkjouPm2f0sAAAgAElEQVRSuX/2Fip4uzH/kUgCfN2LOaSwFbP7FNoopXYp\npVYopRpcqVHBrm9RSqmoxETH3WKwtPB08WR65+k0qtSIsevGsvLYyhv7gJNb4bvhxnIL/T4FJ2fb\nBBVFd+tYY4vTdW/Bjnk39Nb9Z9IYNvtPfNxdmP9oa6r4ywJ3jsTMorAdCNVaNwGmAkuu1FBrPUtr\nHaG1jggICCixgOLKvF29+bjrxzQOaMxz65/jl2PXuVxC8hFYMBB8g4xNchx8c/lSSyno+T7U6Ag/\nPgmH117X2/bFpTF41mbcnJ2Y/2gk1crL99fRmFYUtNZpWuv0gvvLAVelVCWz8ogb5+3qzYwuM2gS\n0ITn1z/Pz0d/vvobMpJg3r3Gdeqh3xlLOgv75ewKA76EgLqw8H5jxdqr2Hs6lSGfbsbT1ZmFj0US\nVsm7hIKK4mRaUVBKVVEF0xmVUq0KsiSZlUcUzaWFYdyGcVfufM7Pga+HQNpp4wxBlq9wDB5+MOQb\nYyG9r/pDWlyhzXadTGHIJ5vxdnNh4WNtCK0oBcFR2XJI6gLgD6COUuqUUmq4UmqkUmpkQZN7gb1K\nqV3Ah8Ag7VC7x4u/eLl6MaPLDJoGNmXchnEsO3LZksxaG5cgTm6Gez6WXdMcjX9VYw5DThosGAy5\n/1w5d/uJ8wz9dAv+Xq4sfCySkApyyciRKUf7ORwREaGjoqLMjiEKkZmXyag1o9iesJ2JbSZyT/g9\nxgu/T4HVE40RLbc9Z2pGcRMO/Gxs6Vm/N9z7BTg5EXUsmQc/30olHzfmPxpJcDlPs1OKK1BKbdNa\nR1yrndmjj0Qp4uXqxUddPiIyKJIJmyYwd99cYyOX1a9Aw37GiBbhuOp0hzsmGfs8r3uT3w+d4/7P\n/iTQ152vR7SRglBKyHxzUaw8XTyZ2mkqz69/nre3vk16aiYjg5ui+kyXuQilQZvRkLgf1r3Fd/mZ\nhFS8g7nDWxHoJ8NOSws5UxDFzs3ZjXdaPEfvbM1H/l680+B2tIv80CgVlOLboGfYYq3LWy4f820v\nVykIpYwUBVH88rJx+eZ+JiWeY0hIV+bGfsfEPyZisVrMTiZu0qcbjvDs4hi+qDYJl3LB+C55AFJP\nmR1LFCMpCqJ4aW2smXPqT5z6zmRcx/cY0XgEiw8tZuz6seRYcsxOKIpAa817Kw/wv2Ux9GhUhfcf\n7ozTkG8gL8vofM7NMDuiKCZSFETx+n0K7F4IHV+C+n1QSvFEsyd4NuJZVh1fxchVI0nLTTM7pbgB\nFqtmwg/RTF0by8CIEKYObo67i7OxM969n8HZaFg8wlgGXTg8KQqi+MT8CGtegYb3wq3P/uOlBxo8\nwJsd3mRn4k4eWPEAZzLOmBRS3IisXAv/mbeNuZuP89itNXizXyOcL90xLbwrdHsd9v8Ev/7PvKCi\n2EhREMUjfrfx22LVCOgzrdCRRnfVuIsZXWYQnxHP0OVDiT0fa0JQcb2S0nMY8ulmVsWcZWKv+rzQ\no17heyq3HgktHoQN78HuRSWeUxQvKQri5qUnGDNdPcsb+yu7Xnm8emRQJHO6z8Gqrdz/8/1EnZGJ\niPbo2LkM+s3YxL64NGbc14IH21W/cmOl4M53ILQ9/DAKTsn31JFJURA3Jz8XFg6DzCSjIPhWvuZb\n6lSow7we86jkWYkRq0bw87FrLKQnStSOE+fpO2MTqVl5zH+0Nd0bVrn2m1zcjMXz/IKMNa5ST9s+\nqLAJKQqi6LSGZU8baxrd/REEN73utwb7BDP3zrkX92SYuWsmjrbkSmm0MvoMgz/ZjI+7C9/9py0t\nQitc/5u9K8Lgr421kb7+9xpJwjFIURBFt2Um7JhrLF/RsO8Nv93f3Z9P7viEXjV6MW3nNF74/QUZ\nsmoSrTUf/RbLY/O2UaeyL4sfb0uNAJ8b/6DAenDvbKOPacl/ZESSA5KiIIrm8Fr45QWo2xNu/2+R\nP8bN2Y3X2r/GmGZjWHZkGcN/GU5SlqygXpKy8yw8880u3v75AD0bB7PwsTZU8rmJ7TNrd4Our8K+\nJbD+7eILKkqEFAVx45IOw6IHIaAe3DMTnG7un5FSikcbP8rk2ydzIPkAQ5YN4dD5Q8WTVVxVwoVs\nBn+ymcU7TvN019p8OKgpHq7FsD1q2yeg6X3GHs/R39/854kSI0VB3JjsVGMGq3KGwfPBvQiXGK6g\na2hXvrjzC/Kt+QxdPpTfTv5WbJ8t/m3v6VTunraR/fEXmHFfc8Z0Di98yGlRKAU9p0BIa/j+PxC3\ns3g+V9icLTfZ+UwplaCUKnQPP2X4UCkVq5TarZRqbqssophYLfDdI8Y+ywO+hPJhxX6IBhUbMP+u\n+YT5h/HE2ieYvnM6Vi3XpYvbst3x9P/4DzSwaGQb7mwUVPwHcXGHgV8Z264uGAwXZMKiI7DlmcIX\nQPervH4nEF5wGwHMsGEWURxWT4RDK+HOt6F6B5sdprJ3ZeZ0n0Ofmn34eNfHPLH2CVkao5jkW6y8\ntmwfo+Zvp26QLz+MakfDqv62O6BPAAxeYJxhfj3EWCtJ2DWbFQWt9Xog+SpN+gBfasNmoJxSyga/\nrohisetr2PQhtHwEWg63+eE8XDyY1G4SL7Z+kU2nNzH4p8EyA/omJV7IYejsLXyy4Sj3twll4Yg2\nJbPsdZVG0HcWnN4GS8cYQ5mF3TKzT6EqcPKSx6cKnhP25lSU8Z85rAN0f7PEDquUYlDdQczuNpvM\n/EyGLB/CL8d+KbHjlybbT5yn19Tf2XEihff6N+HVPg1xcynB//71ekKn8bDnG/h9cskdV9wwM4tC\nYT1ahf4KoZQaoZSKUkpFJSYm2jiW+Ie0OPj6PvCtAv3ngLNriUdoXrk5C3supHb52jy77lne+vMt\n8ix5JZ7DEWmtmbv5OANn/oGri2Lx423p16KaOWE6PAON+sOaV41tWoVdMrMonAJCLnlcDYgrrKHW\nepbWOkJrHREQEFAi4QTG9d+vh0BuujFT1buiaVECvQL5rNtn3FfvPubFzGPYimGcvHDy2m8swy5k\n5zHm652MX7KX9rUq8dPoDjQItmH/wbUoBb2nQtUW8N2jcKbQMSjCZGYWhaXA/QWjkCKBVK11vIl5\nxKW0hqVPGEMJ+34CleubnQg3ZzfGtRrH+7e/z4kLJxjw4wBWHV9ldiy7tPtUCj2n/s7yPfE8e0dt\nZj/QEn+vkj/L+xdXT2ONLA9/Y2hzupz52xtbDkldAPwB1FFKnVJKDVdKjVRKjSxoshw4AsQCnwCP\n2yqLKILfJ8OeRdDpJajbw+w0/9A5tDOLei2iun91nv7taV7b/Josj1FAa82nG47Qb8Ym8vKtfD0i\nktGdwnFyKqb5B8XBt4oxxyXjHCwcCvnyvbMnytEWIYuIiNBRUbI0r01FL4FFDxib5fT7tNC9EexB\nniWP97e/z5f7vqRehXq8detbVPe/yhLPpVxyRi5jF+1izf4EutSrzLv9G1POy83sWFcW/b0xM77p\n0CvuwSGKj1Jqm9Y64lrtZEaz+KfT2+D7kVCtFfSZbtf/UV2dXRnbcixTO00lLiOOAT8OYOH+hWVy\ntdUNhxLp8cEGNhw6x8u96vPJ/S3suyAANLgHbhsHO+fBH9PNTiMKSFEQf0s9Zcw89Qko2CynBMaw\nF4PbQ25nce/FtKjcgv9t+R+Pr3mcc1nnzI5VIrJyLbz8w16Gzf4Tb3dnFj/elofaVS++5Sps7bbn\noX4fWDUeDq40O41AioL4S046zB9kjDga8o1RGBxIoFcgM7rM4IVWL7D1zFbu+eEe1hxfY3Ysm9p1\nMoW7pm5gzh/HebBtGMvGdLDt7GRbcHKCu2dA5Ybw7cOQsN/sRGWeFAVRsKbRcEiIhv6fG2viOyCl\nFEPqDeGbnt8Q5B3EU789xfiN40nPTTc7WrHKs1h5f/VB+s7YRFauhXnDWzOxd4PiWd3UDG7expBn\nNy9YMBAyr7YQgrA1KQoCVo6Hgz8baxrV6mJ2mptWo1wNvurxFY82epSlh5dy9w93s/7UerNjFYuD\nZy9w74xNvL/6EL0aB/HzU7fSPryS2bFunn9V45JlWjx8c7+xzaswhRSFsi7qM9g8HVo9Bq0eNTtN\nsXF1dmVM8zHMvXMuPq4+jFoziv9u+C+pOalmRyuS3Hzj7OCuDzdw8nwW04Y04/1BzfD3tIO5B8Wl\nWoQxCunYBlgxVtZIMomL2QGEiQ6thmXPQq2u0O11s9PYROOAxnzT6xtm7p7JZ3s+Y1PcJl6MfJGu\noV3Njnbddpw4z/Pf7ebg2XT6NA1mQs/6VLyZndHsWeMBkBBjzJMpXx3aP2V2ojJH5imUVXE74PO7\noEINeGg5ePiZncjm9ifvZ8LGCcQkx9A1tCv/bf1fKnna76WXzNx83v3lIJ9vOkoVPw9eu6chnepW\nNjuW7VmtsPgR2PudMZu+8QCzE5UK1ztPQYpCWZR8FGZ3BRdPGL4S/MrOiuV51jzmRM9hxs4ZuDm7\nMbrZaAbWGYiLk32dNK/df5YJP0Rz6nwWwyJDea57HXw9StGlomvJz4F5/eDEZhj6LdS43exEDk+K\ngihcRpJREDKTjIIQUMfsRKY4nnac17e8zqa4TdSrUI8XI1+kSUATs2NxMjmTV37cx+qYs9QM8OaN\nvo1pVb2C2bHMkZUCn98JKSfh4RXGvgyiyKQoiH/LzYQve0P8bnhgKdwSaXYiU2mtWXl8JW9vfZuE\nzAT6hffjqeZPUc6jXIlnyc6zMGv9Eab/Gouzk2JM53Aeble9ZPc8sEepp+DTroCG4augXMg13yIK\nJ0VB/JMlH74ZBgdWGPsr1+9tdiK7kZGXwYydM5gXMw9fN1/GNB9D31p9cXYqmXH/vx1IYOLSaI4l\nZXJXoyBevKseweU8S+TYDuFsNHzWHXyD4OGfwauMnjndJCkK4m9aw49PwvY5cOc70HqE2Yns0sHz\nB3lt82tsT9hOePlwxkaMpU1wG5sd79DZC7y+PIZfDyRSo5I3r/RpQIdwx5pJXmKOrjf6GIKbw7Dv\njYlu4oZIURAGrY11ZTZNhfZPQ5eXzU5k17TWrDq+isnbJnM6/TS3VbuNZyKeKdbVV8+l5/D+6oMs\n+PMkXq7OjOpUi4faheHu4qAzkkvK3sXGUhi1OsOgBeBi5wv+2RkpCsKw7h349X/Q8hHo8a5dr3pq\nT3IsOXwV8xWzds8iJz+HgXUHMrLxyJvqb8jOs/DZxqN89OthsvIs3Nf6Fp7sHF565xzYwrY58OMY\nqNcb7v0cnO1r1Jg9k6IgYPPH8PPz0HiQseiYUxnvtCyCc1nnmL5zOosPLcbLxYsHGzzIsPrD8HK9\n/ssXFqtmyY7TTF51kNMpWXSpF8i4O+tRK9DHhslLsU3TYOWL0PQ+6D1N/l1fJ7soCkqp7sAHgDPw\nqdb6zctefxB4Bzhd8NQ0rfWnV/tMKQrXacdX8MPjULcn9J8jv1HdpEPnDzF1x1R+PfkrFTwq8Gij\nRxlQZwBuzle+hGG1an6OPsPkVQeJTUinQbAfL/aoR9ta9jthzmH8+gasexNaj4Tub8oZ8HUwvSgo\npZyBg0BX4BSwFRistd53SZsHgQit9ejr/VwpCtch+nvj2mv122DIQnCRyxPFZVfiLj7c/iF/nvmT\nIO8g/tPkP/Sq2esfk9+01vx2IJF3Vx4gOi6NWoE+PNO1Nt0aVLGvbTEdmdbwy4vGul0dnoFO46Uw\nXMP1FgVb/vrYCojVWh8pCPQ10AfYd9V3iZuzdzF894ixc9qgr6QgFLMmAU2Y3W02m+M38+H2D5mw\naQKz987mkUaP0KN6D7YcSeX91YfYdvw8t1TwYvKAJvRpWhVnKQbFSyno9hrkpsOG9wBl7CcuheGm\n2bIoVAVOXvL4FNC6kHb9lFK3YpxV/J/W+mQhbcT12PMtLB4BIa3gvkXGOvXCJiKDImndozVrT67l\n410fM37jeF7ZMIX0s7dSwdqe1+5pyICIEFyd5Xq3zSgFPd837m94F7QVOk+QwnCTbPkvtrDvzOXX\nqn4EwrTWjYHVwJxCP0ipEUqpKKVUVGJiYjHHLCV2L4LFj0JIa7jvW3D3NTtRqWexatKT65IaO5rM\nEw+i8/3wCFqCd623sPquJ1/nmB2x9HNyMgpDi4eMlVVXT5Qlt2+SLfsU2gATtdbdCh6/AKC1fuMK\n7Z2BZK31VfcTlD6FQuxaCEtGwi1tjT4EdxnVYkuZufl8t/00n244wvGkTMIDfRjVsRZ3NarCtsSt\nzNo9i61ntlLOvRwD6wxkUN1Bdr0aa6lgtcLyZyFqNrQdA11flTOGy9hDn8JWIFwpVR1jdNEgYMil\nDZRSQVrr+IKHvYEYG+Ypnf78BJaPhbD2RkGQS0Y2cyY1mzl/HGP+lhOkZuXRuJo/Hw9twR31K1/s\nQI4MiiQyKJLtZ7fzefTnzNo9i8/2fsZdNe7i/vr3E14+3Nw/RGnl5AR3vQdOzrDpQ8jNgB7vGI/F\nDbFZUdBa5yulRgO/YAxJ/UxrHa2UehWI0lovBcYopXoD+UAy8KCt8pQ6WsO6t+G316FOD7j3M3CV\n9XJsYfepFGb/fpRlu+Oxas0d9aswvEN1IkLLo67w22jzys1pXrk5x1KPMS9mHj/E/sCS2CW0DW7L\n/fXvp01wG5yU9DcUK6WMLWVdvWDj+8ZKwH1nyWCLGyST1xyR1Qo/j4M/Z0KTIdB7qsxDKGbZeRZW\n7I3nq80niDp+Hh93FwZEhPBQuzBCKtz4ujsp2SksOriI+fvncy7rHKF+ofSv3Z8+NfuYsiprqbdp\nKqx8yRiWPegr6WPDDuYp2EqZLwp52caktL3fQZvR0HWSzOgsRocT01mw5QTfbj9FSmYeYRW9GBoZ\nysCWIcWyyU2uJZdfjv3CooOL2JGwAzcnN7qFdWNAnQE0CWhyxTMPUQQ7F8APo4x9GO77FnzK9mKD\nUhRKo/REWHgfnNwCXSZCu6ekM60YZOdZWLnvLPO3HGfzkWRcnBTdGlRhSOtbaFOjos0mnB1IPsCi\ng4v46chPZORlUKd8HfqG96VH9R5y9lBcDv4C3zwAPoFGn1tgPbMTmUaKQmmTsB/m94f0BLhnJjS4\n2+xEDk1rzbbj5/lu+2mW7Y4jLTufkAqeDG51C/1bhBDgW3LXoTPyMlh+dDmLDiwiJjkGFycXbq92\nO71r9qZ9tfa4OpWhbTht4dQ2+HqwscnUvZ9B7TvMTmQKKQqlSexqWPSQ0ZE8eAFUbWF2Iod1PCmD\nxdtP8/2O05xIzsTT1ZnuDavQt3lV2tWsZPoyFAeSD/DD4R9YdmQZydnJVPCoQI/qPehdszd1K9SV\ny0tFlXoaFgyCs3vhjv9B5ONl7ixbikJpYLXC+nfgtzegcgMY/LVsR1gEcSlZrNh7hmW749h+IgWl\noG3NivRtVo3uDavg7W5/nfR51jw2nt7I0sNL+e3kb+RZ8wjzC+OOsDvoFtaN8HLhUiBuVG4GfP8Y\nxPwITQYbQ1jL0BBuKQqOLjPZWLIidhU0Hgg9p5Spf8A369T5TFbsOcPyvfHsOJECQN0qvvRpWpW7\nmwUT5O84w3dTc1JZeXwlvxz7ha1ntmLVVqr7V6dbWDe6hXajVvlaZkd0HFYrrHvLuAXUNbamDaht\ndqoSIUXBkZ3801jlNP2ssSxwxMNl7lT3RmmtOXg2nbX7E/h5bzy7TqUC0CDYjx6NgrizYRVqBDj+\nTO+krCTWnFjDL8d+IepsFFZtJcwvjI4hHbkt5DaaBDT5x4qt4goOr4XvHoW8LOj1ATTub3Yim5Oi\n4Ijyc4014n+fAv7VoP8X0n9wFdl5FjYfSWLt/gTWxCRwOiULgEZV/S8WgrBKpffs6lzWOVYdX8Wv\nJ35l69mt5Fvz8Xf3p0PVDtwWchvtgtvh6ybj868oLc745evEH9BoAPR4GzzLm53KZqQoOJqE/caC\ndmd2Q9Oh0P0N8PAzO5Vd0Vpz5FwGm2LPse7gOTbGniMrz4KHqxPtawXQuV4gHesEUsXfw+yoJS49\nN52NcRtZd3Id60+vJzUnFRcnF5oGNCUyKJI2wW1oULEBzrLswz9Z8mDDZFj/NngHGDu5hXcxO5VN\nSFFwFHnZxpT8DZONhex6fQj1epqdym6cSc1mY+w5Nh4+x6bYJM6kZQNQtZynUQTqBtKmRkU8XOWH\n3V/yrfnsStzFulPr+CPuD/Yn7wfA182X1lVaXywSIb4h0ln9l7gd8P1/IDEGmg0zFtTzqmB2qmIl\nRcERHF4Ly56B5CPQ8F7j7MAn0OxUptFacyI5k6hj54k6fp4/jyZxODEDgPJerrStWYm2tSrSrmYl\nQit6yQ+065ScncyW+C1sjt/MH3F/EJ9hrEEZ6BVI80Bjjabmgc2pVa5W2T6TyMs2Rvptmgoe/sYE\n0WbDSs2KAVIU7FniAVj9ChxYBhVqGEPjanYyO1WJy823Eh2Xyrbj5y8WgnPpxh4Evh4utAgtT7uC\nQlCvip/pcwhKA601Jy6cYHPcZrad3ca2hG0kZCYA4OvqS5PAJjQPbE7jgMbUr1i/bPZJnN1nLMN9\nfKPRp3fH/yC0rdmpbpoUBXuUFge/vQk75oKrN7R7Eto+Aa6l/xp4Tr6Fg2fS2XM6lT2nU9l7OpUD\nZy6Qa7ECcEsFL1qElqdFaHkiwspTO9BXikAJ0FoTlxHH9rPb2Z6wnR1nd3A49fDF18P8wmhQqQEN\nKzakQaUG1K1QF08XxxnOW2Raw+5vYNUESD8Dtboau7oFNTY7WZFJUbAn52Jh0wew62vjH1vL4XDr\nWPAufRuvaK1JTM/h0Nl0Dp69wIEzF9hzOpWDZy+QZzH+rfl5uNComj8Nq/rTtFo5WoSWJ9Cv9BdG\nR5GSnUJ0UjTRSdHsPbeX6KToi2cTTsqJGv41CC8XTnh541arXC2CfYJL51LguZnw5yxjRGB2CtTt\naWzic0thOwvbNykKZrNa4eg6YyeomJ+MNd2b3gftxkD5MLPT3TSrVZNwIYej5zI4lHCBg2cvcPBs\nOofOXuB8Zt7FduW8XGkYbBSARgW3kAqe0h/gYBIyE4g+F83epL0cTD7IoZRDnE4/ffF1LxcvapWv\nRXi5cKr7VyfML4xb/G6hmk81XJ1LwdpNWSnwx3TY+glknYdqraDN41DnLnBxMzvddZGiYJbzx2HP\nItj+JaQcN8Y9RzwMrUc6XCdynsXKmdRsjiVlcDwpk+NJGRxLyuREUibHkzPIzrNebOvr4ULtyr7U\nruxDeKDvxfsBvu5SAEqp9Nx0YlNiOZRyiEPnjVtsSiwpOSkX2zgrZ6r6VCXUL5RQv1Bu8buFqj5V\nCfYOJsgnCG9XB5tHkpsBO+fDH9Pg/DHwrGCsONDsPqjc0K4nmdpFUVBKdQc+wNh57VOt9ZuXve4O\nfAm0AJKAgVrrY1f7TLsrClarMYztwHJjTZX4XcbzYR2gxYPG6aYd9hnk5FtISs8lPjWb+NQszqRm\nE5di3P/rucQLOVgv+efh7uJEaEUvQit6E1bRi1sKvoYH+lLZT374C0NKdgrHLxzneNpxjqUe43ja\ncU5cOMHxtONk5Wf9o62fmx/BPsEEeQdd/FrZqzIVPStSybMSlTwr4ePqY3//tqwWOPwr7JwH+5eB\nJRfKV4e6dxk7IYa0truNr0wvCkopZ+Ag0BU4hbFn82Ct9b5L2jwONNZaj1RKDQLu0VoPvNrnml4U\nctIhYR/E7YTjv8Ox341t/8A4pazXC+r3LtFLRPkWK2nZ+aRl5ZGWnUdaVj6pWXkkZeRwLj2XpPQc\nktJzOZeeQ1KG8fVCdv6/PsfLzZkgfw+C/D2Nr+U8qVrOo6AIeBPo6y6dv6LItNYkZCYQnxFPfEY8\ncelx/7qfkZfxr/d5OHv8o0hU8qxEOfdy+Lv74+/uTzn3cvi5+V187OfmV7JLfWQmw74lRnE4ut4o\nEG4+RmEIbWt8rdzA9HkP9lAU2gATtdbdCh6/AKC1fuOSNr8UtPlDKeUCnAEC9FVC2bQoWPIhO9Xo\nUMpKMUYdpJz4+5a4H5IOAwXx/KpB9Q7GWUHNjuAXfPGjrFZNntVKnkWTbzG+5lms5Fs0uRYr+VYr\neflGm9x8K1m5FrLyLGQWfM3KzScr10pmXj7ZuZc+b7n4g9/4mkdGruWKfySloLyXG5V83Kjo7U5F\nHzcq+bgbj33cqeLnQVA5D4L8PPHzdLG/38hEmaG15kLeBRIzEzmXde6qt7TcNKzaesXP8nH1wdfN\nF29Xb7xcvPB09cTbxRsvV6+Lz3m5el386u7sjpuzG25Obrg7u+Pq7Pr3c5c87+bshquTK07KCWfl\njLOT8z872LPTjPlHxzbA8U3GL5B/8Q2GyvWhQk3wr2osZeNXzRhw4u5nrGBgw/2kr7co2LKcVgVO\nXvL4FHB5l/3FNlrrfKVUKlAROFfcYXb/+i3+G17GWVtwwoozFuOmrThhwZV8vMgu9L1ZuBNHACdU\nMPtVCw4Qyn7CiE+rhN4Feheg92DVu8m3avKtGou1eIqtm4sTXm7OeLo641nw1c/DleqVvPHzdMHP\nwxU/T1f8PFwKvhY89nShorc75b1ccXEuhaNCRKmjlMLPzQ8/Nz9qlqt51bZWbeVC7gXSctJIzU0l\nNce4peSkkJqbSlpOGmm5aWTlZ5GRl0FmXibJ2clk5mUat/xMciw5xZbdWRnFwcXJBSflZNyv4IpT\nhabGzxxrPk6WfFTeIYiPgTiNKvjl8q9fwxSgUKCcCp5TF/solPEKfStH8sBds4otd2FsWRQK+5Xz\n8p+U19MGpdQIYATALbfcUqQwbj7lOOdVC6sySoJVuaCVM9aLNxeynH3IdvYl29mXLBdfMpzLc949\niCwXf+MbpcBJgS+KVqrgm1jwTXNSCqXAxUnh6uyEi7Px1dX5r8dOuDkrXJyccHVxwvWSdsYPfhc8\nXZ3xcnPG45IC4CyXa4T4FyfldPFyUQhF22Mkz5pHVn4WmXmZ5FpyybXkkmPNIc+SR44lhxzL3/dz\nrQWvW3LIt+Zj0Ras2nkQl2AAAAa0SURBVIpFW7BY/75v1Vbyrfl/v/ZXO6sFjeaviyDakgt5mejc\nDLDkoi15xnPWPOOKBRqtraA1GqtxH6joXbm4/gqvyJZF4RT847tVDYi7QptTBZeP/IHkyz9Iaz0L\nmAXG5aOihKnbsgu0LJ0LXQkhbpyrkyuubq74ucnCk5ey5XWFrUC4Uqq6UsoNGAQsvazNUuCBgvv3\nAmuv1p8ghBDCtmx2plDQRzAa+AVjSOpnWutopdSrQJTWeikwG5irlIrFOEMYZKs8Qgghrs2m47a0\n1suB5Zc9N+GS+9lA6d/ySAghHIQMSxFCCHGRFAUhhBAXSVEQQghxkRQFIYQQF0lREEIIcZHDLZ2t\nlEoEjhfx7ZWwwRIadk7+zGXD/7d3vyFSVWEcx7+/0rRILNKiP6YFGplvUgl9Y4pRIqEEEhtIWdIL\no6A/BEEviqJeFFbYH6xIKijTJHQTQqI0I1prUzNXksxMhcKIEsRIracX5+60zf67uvfONLO/D4j3\nzj2z8zx7Z/eZc+7dc5zz4HGqeY+NiNH9NWq4ojAQktrzTAjVTJzz4OCcB4+y8/bwkZmZVbgomJlZ\nxWArCuXOOfv/5JwHB+c8eJSa96C6pmBmZn0bbD0FMzPrQ9MVBUlzJO2WtEfSQz0cHyZpVXZ8i6Rx\ntY+yeDnyvl/SLkk7JH0kaWw94ixSfzl3abdAUkhq+DtV8uQs6ebsXHdIervWMRYtx3v7UkkbJW3L\n3t9z6xFnkSStkHRI0s5ejkvSsux7skPS5MJePCKa5h9piu7vgcuBM4CvgYlVbe4ClmfbLcCqesdd\no7xnAWdl20saPe88OWftRgCbgTZgar3jrsF5Hg9sA87N9s+vd9w1yPkVYEm2PRHYV++4C8h7BjAZ\n2NnL8bnAB6QFIKcBW4p67WbrKVwD7ImIvRFxDHgHmF/VZj7wRra9Bpitxl+tvt+8I2JjRBzNdttI\nK+E1sjznGuBx4CnoZQHuxpIn5zuBFyPiN4CIOFTjGIuWJ+cAOpdPG0n3FR4bTkRspodVKLuYD7wZ\nSRtwjqQLi3jtZisKFwMHuuwfzB7rsU1EnAAOA+fVJLry5Mm7q8WkTxmNrN+cJV0NjImI9bUMrER5\nzvMEYIKkzyS1SZpTs+jKkSfnR4GFkg6S1m+5pzah1dXJ/sznVuoiO3XQ0yf+6tur8rRpNLlzkrQQ\nmApcW2pE5eszZ0mnAc8Ci2oVUA3kOc9DSENIM0m9wU8lTYqI30uOrSx5cr4FeD0ilkqaTlrNcVJ0\nrnbfnEr7PdZsPYWDwJgu+5fQvStZaSNpCKm72Vc3rRHkyRtJ1wEPA/Mi4s8axVaW/nIeAUwCNkna\nRxp3bW3wi81539/rIuJ4RPwA7CYViUaVJ+fFwGqAiPgcGE6aH6iZ5fqZPxXNVhS+BMZLukzSGaQL\nya1VbVqB27LtBcDHkV25aWD95p0NpbxMKgiNPs4M/eQcEYcjYlREjIuIcaTrKPMior0+4RYiz/t7\nLemmAiSNIg0n7a1plMXKk/N+YDaApCtJReGXmkZZe63ArdldSNOAwxHxUxFfuKmGjyLihKS7gQ2k\nuxZWRESHpMeA9ohoBV4jdS/3kHoILfWLuBg5834aOBt4N7uuvj8i5tUt6AHKmXNTyZnzBuB6SbuA\nv4AHI+LX+kU9MDlzfgB4VdJ9pCGURY3+QU/SStIQ4KjsWskjwFCAiFhOunYyF9gDHAVuL+y1G/x7\nZ2ZmBWq24SMzMxsAFwUzM6twUTAzswoXBTMzq3BRMDOzChcFsyqSNkm6oeqxeyW91MdzjpQfmVn5\nXBTMultJ979fackeN2tqLgpm3a0BbpQ0DCBbc+MiYHu2FsVWSd9I6jYrq6SZktZ32X9B0qJse4qk\nTyR9JWlDUbNamhXJRcGsSvYXwF8AnTOMtgCrgD+AmyJiMmkqiaV5p12XNBR4HlgQEVOAFcATRcdu\nNlBNNc2FWYE6h5DWZf/fQZqZ8klJM4C/SVMVXwD8nOPrXUGaoO/DrI6cDhQyV41ZkVwUzHq2Fngm\nW+bwzIjYmg0DjQamRMTxbPbV4VXPO8F/e+CdxwV0RMT0csM2GxgPH5n1ICKOAJtIwzydF5hHAoey\ngjAL6Gmd6x+BiUprgY8km72TNIX16Gy+fyQNlXRVmTmYnQr3FMx6txJ4j3/vRHoLeF9SO7Ad+Lb6\nCRFxQNJqYAfwHWm9ZCLimKQFwLKsWAwBngM6Ss/C7CR4llQzM6vw8JGZmVW4KJiZWYWLgpmZVbgo\nmJlZhYuCmZlVuCiYmVmFi4KZmVW4KJiZWcU/ldsF3BH9f8oAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fd658f75438>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "wc_priors['mu'].plot()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "As shown in the plot above, the probabilities for each channel are weighted such that they are distributed fairly well along the normalized range of input channels. However, we don't require that channel 1 have a mean value that is less than channel 2 or 3."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can also create priors with multidimensional coefficients. This is useful in the case of the n-tap FIR filter which takes multiple input channels. The number of coefficients for the FIR filter is the number of inputs times the number of taps. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "n_inputs = 3\n",
    "n_taps = 7\n",
    "mu = np.zeros((n_inputs, n_taps))\n",
    "sd = np.full((n_inputs, n_taps), 5)\n",
    "coefficient_prior = Normal(mu, sd)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can get the mean (i.e., expected) value of each coefficient or draw a random sample from each distribution. This can be useful for initializing the fitter."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[ 0.,  0.,  0.,  0.,  0.,  0.,  0.],\n",
       "       [ 0.,  0.,  0.,  0.,  0.,  0.,  0.],\n",
       "       [ 0.,  0.,  0.,  0.,  0.,  0.,  0.]])"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "coefficient_prior.mean()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[-2.63, -1.69, -0.38, -4.98, -1.09, -5.09,  2.67],\n",
       "       [ 3.03,  4.17,  2.97,  6.34, -2.88,  2.82,  0.05],\n",
       "       [ 1.12, -0.68, -1.42,  5.18, -2.9 ,  8.62, -2.85]])"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "coefficient_prior.sample()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Note that these methods return data in the same shape as the set of coefficients specified by the prior. This makes it easy to generate initial data in the correct shape required by the model and module evaluation system.\n",
    "\n",
    "We can also draw more than one sample. Let's go back to the priors for the weight channels module."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(10000, 3)"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "samples = wc_priors['mu'].sample(10000)\n",
    "samples.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x7fd658f75e10>"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEKCAYAAAARnO4WAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzsnXd8W+X1/9+Ptrwlj8QrcXZC4uxF\nCCMhjBAIo5RSWkppC5R+6Y+2lNkCAQoNlAKllLa0paWUVcoIEHYgezrTieNsx4nt2I4tD0nWvr8/\nFMt7S5YtPe/Xy6/4Xt17nyNH+txzz3Oec4SiKEgkEokk8lGF2wCJRCKR9A9S8CUSiSRKkIIvkUgk\nUYIUfIlEIokSpOBLJBJJlCAFXyKRSKIEKfgSiUQSJUjBl0gkkihBCr5EIpFECZpwG9CclJQUJScn\nJ9xmSCQSyaBh+/btpxVFSe3OsQNK8HNycsjLywu3GRKJRDJoEEIc7+6xMqQjkUgkUYIUfIlEIokS\npOBLJBJJlCAFXyKRSKIEKfgSiUQSJUjBl0gkkihBCr5EIpFECVLwJRKJJEoYUAuvekvuK7mB3/Nv\nyg+jJRKJRDJwkR6+RCKRhJGc+1aSc9/KfhlLCr5EIpFECVLwJRKJJEqQgi+RSCRRghR8iUQiiRKk\n4EskEkmUIAVfIpFIogQp+BKJRBIlSMGXSCSSKEEKvkQikUQJUvAlEokkTPTXCttGpOBLJBJJlCAF\nXyKRSKIEKfgSiUQSJUjBl0gkkighogU/95XcFrXyJRKJJJqJaMGXSCQSSRMRJ/jSo5dIJJL2iTjB\nb6S58MvQjkQikUSw4EskEomkJVLwJRKJJErQhNsAiUQiiQSal0koWr4kjJZ0jPTwJRKJJEqQgi+R\nSCRRghR8iUQiiRJkDF8ikUj6kf4uidwc6eFLJBJJlCAFXyKRSKIEKfgSiUQSZHLuWxnW0E1HhFzw\nhRBqIcROIcRHoR5LIpFIJB3THx7+ncD+fhhHIpFIJJ0Q0iwdIUQWsAR4HPhFKMfSu1T4hIJbq/To\nvOZF1fJvyg+2WRKJJIrpavVtvMuGV6ixaw39Yk+oPfzngHsAX0cHCCFuFULkCSHyKisrezyAz+dl\nbHEcV6/JYNqhpD6YKpFIJP3LLfkf8rcvn+TC4jwUX4cyGTRCJvhCiMuBCkVRtnd2nKIoLymKMlNR\nlJmpqak9Huf4nl3M25uMwa1mfFE8pjptb02WSCSSfmN8dREXncjD7KznlzvepGHX7pCPGUoP/xxg\nqRCiCHgTWCiE+E+wB8mZMp3S5AYAVAjm7jOjKD0L60gkEkl/olJ8/GT3e4HtDemTiJk+LfTjhurC\niqLcryhKlqIoOcD1wFeKonw32OMIIdgysRqf8Iv8EIuB/etXB3sYiUQiCRqXHN/KmNoSAJwqDS9N\nWtov40ZEHn5tnIeCnLrA9tr/vIzWLcJokUQikbSkcQI3zmXnpoKPA/vfHrOAilhzv9jQL4KvKMpq\nRVEuD+UYu8bUYtN7ALDVWJgqJ3AlEskA5Mb9n5LosgNQbjTx9tiF/TZ2xBRP82gU8iZYOH+Xf+J3\nwvF4DmVbqYl3h9kyiUQykAhno5KGvftYcmxTYPul3KW41P2XaBIRIZ1GjqXbOWV2AKBSBHP2mUHO\n30okkgGAUHycevRR1GdEaXvaWDamT+pXGyLGwwdAwOaJ1Vy9IQvF5yO92sCIshiOZdjbHNp8wZVE\nIpG0JthPApce34pjzx4A3Co1L06+GkT/zjVGlIcPUBPvZu+wmsD2zP0mNB45gSuRSMJHotPKzfua\nbiBvj1lAaVzP1x31lYgTfIBdY2po0HkBiHVqmHI4EZBevUQiCQ837/uYeLd/vVBZjJm3xl4YFjsi\nUvDdWoVtEyyB7YnHEki0Rlb0SiKRDA7OqjrGJcVbA9t/nnxVv07UNiciBR/gaIaNclPTBO7Ze5Pl\nBK5EIulXVD4v/7f73cD2xvRJbBt6VvjsCdvIoebMBG7jCtyh1QZGn4wNs1ESiSSaWHp0AyPrygBw\nqLX8Nbd/VtR2ROQKPmBJcDPnimsD2zMLTdjrasNokUQiiRbMDbXcWPhZYPv1cRdREdO0orZo+ZJ+\nXwcQ0YIPcPY3vo3V6F+Ba3CrWfufl8NskUQiiQZu2fshMR4nAMVxabw3+rwwWxQFgq81GNg0sSqw\nvW/NKor37gmjRRKJZKAQqt6zM8v3c0HJrsD2i1OuxqMKf+JIxAs+QEmag2NDbYHtL//+J1TeMBok\nkUgiFr3HyR27miZqv8qazu7UMWG0qImoEHyArWdZcGn8HWUsZSXkHk0Ms0USiWSw094Two2FnzOk\nwZ8WXqeN4aUwT9Q2J/zPGP1Eg8HLjrEW5hYkAzD5SCLH0m3UxXnCbJlEIhlsdBQGGl1zkqsOrw1s\n/y33Cmr1cW2O6+/J2kaixsMHODDcytBR/kcrtU9w9r5k9nxPxvMlEknfUfm83Lnz7UBxtF0po/ky\ne2aYrWpJVAm+IuCiW3+KUPnfdnqVgYK1X4V0zNxXcgM/EokkcrnqyDpGN+ti9fzUa/u9OFpXRJXg\nA6TljGT6ZVcGtlf/++8YnFH3Z5BIJEHEdfJky5z78RdRFpcSRovaJ+qULveVXOZ984ZAbr7DWs+c\ngv5pLyaRSAYujROwPU7TVBROPfIoBq+/2dKxhHTeGX1B8A0MAlEn+AAz3prNxklNufkjymIZdsoY\nRoskEslg5fySXdjWrQPAh+APU6/Fq1KH2ar2iWjBz78pv8PXSlMdHMq0Brbn7jPjsFk7PF4ikUha\nk+i08uM97we2Pxx5DgfMw8NoUedEtOB3xbYJFmIS/c3OY5wa7nzocjnBKpFI2tBRqOfHe94nyeVf\n1FlpTOTfEy7tb9N6RFQLvkvnY9EPfxLYHnsyjvTThk7PkTcEiUQCcHbp3hblE56fei12bef6EW6i\nZuFVR4yZM4+ioTZyTvlLJ8/LT2bFuaVhtkoikQxEGr38OJedO3a/E9j/RfZM8oZMCJdZ3SaqPfxG\ntkysxqn1F9eJb9Aw/WBSj68hPX+JJHq4LX8FZmc9AFWGhAFVPqEzIkrwO5uk7YwGvY+tzVoiTiiK\nZ+FzM6WISyRRSmepmbNOFbDoxPbA9h+nfAOrLiYs9e17SkQJfkfk35Tf5c3gSKaNk6n+JsMCwfw9\nyai9Ha+SkzcCiST6iHU18P92/S+wnXD55WxJnxhGi3pG1MfwAwjYNKmKq9ZmoPWqSLRpmX4giW1n\nWaS4SyQSwN/UJMVRB4BFH8eYXz0AT27s8Pii5UsCTwsDwfuPCg+/u9iM3hahnbOK4hlSrQ+jRRKJ\nZKAwvfwAlxRvDWz/afI1aEymLs8bSKEe6eG34lC2leHlMWRVGv2hnd3JrDi3DI9GCbdpEokkTMS6\nGrhz19uB7bUZk9mQObndYweKuLdHVHn43YnlI2BjbhVOTWPWjpaZhV3fxSUSSeTy4/z3SWuoAaBW\nF8OLU64Os0W9I6oEv7vYDV62TGwK7YwvjiejcmAvqJBIJKFhXml+i6ycF6Z8g1p9fBgt6j1S8Dvg\naIaN40Psge1z8pPRuQdWbWuJRBJaTI46/l+zUM5XWdNZnzkljBb1DSn4HXEma8eh84d2Yh0aZssy\nyhJJ9KAo3LnzbRJdfsfvtCGRFycPzlBOI1LwO8Gh97Uoozy6JI7s8s7LKMsOVxJJZHDx8a3MKd8f\n2P799G9h0w3uMuohE3whhEEIsVUIsVsIsU8I8UioxgolxUMbOJLRVDZ5Xn6y7JAlkUQ4Q2xV3Lb3\ng8D2ipHz2ZU2NowWBYdQpmU6gYWKoliFEFpgvRDiE0VRNodwzJCwZaKFoVUGYp0ajC415+xJZtXM\nyh5fp7nX39syEBKJJLSoFB937XiTGI8TgBNxqfzzrMvCbFVwCJmrqvhpdI21Z34GZTK7S+tj/eSm\n0E52ZQzjj8eF0SKJRBIqrj68ltyqYwB4hYqnZ3wbp0YXZquCQ0gXXgkh1MB2YDTwJ0VRtoRyvO7S\nG++6LNXBvpw6JhYlADCz0ETVyWKSs4YF2zyJRBImRtaUcNP+TwLbb4y9kIOmyPmOhzQYrSiKV1GU\nqUAWMFsIMan1MUKIW4UQeUKIvMrKnodJ+pMd4yxUx7sA0PhUrHz+d3jc7jBbJZFENr1qLN4L9B4n\n9+a9htbnz8w7mJTFm+MWhXzc/qRfZh8VRakBVgNt+n8pivKSoigzFUWZmZqa2h/m9BqvGtZOPY1X\n5Y9MVR4/xi0PLgyzVRKJJBjcuvdDhlkrAGhQ63hq5nfwqtTt1sIZSPVxekIos3RShRBJZ343AouA\nwlCN11/UxLvZNr5pFe6kY4kcz9/VyRkSiaQ/6c0TQd3nn3NZUVM+yZ8nX0VJ3MB2QHtDKD38dOBr\nIcQeYBvwhaIoH4VwvH6jcHh9oHY+wKcvPoveJVM1JZLBSIq9hsK77gtsr8mcwhfDZoXRotARyiyd\nPYqiTFMUZbKiKJMURXk0VGP1OwI2TD4dWIVrra7i7HxzmxwkufhKIhnYqBQfd29/nXi334ErNybx\nwpRvgIjMMirSLe0lDXofG3KbUjVzymMZVyxTNSWSwcQ3D37N5KqjAHgR/G7GDVh1MWG2KnRIwe8D\nJ4Y0UDisPrA9e78Zc502jBZJJJLuMq76ODcWfhbYfnPcIvaljAyjRaFHNkDpI9smVJNm0WOu16H2\nCc7fmcqH58iGKRLJQCbG3cC9ea+jVnwA7DPn8HovUzAHU7ZOxHn4jU1OQlW6oPV1vWpYPa0St9r/\nwUm0aTl7b9t4fmtkgTWJJEwoCj/f+V/S7f6QrE1j4KkZN+BTqcNsWOiJOMEPB3VxHjZNqg5sjyqN\nY3RJbBgtknTKskT/T6SOJ+mUpUc3ML+0yXH7w7RrqYiNjtLnUvAJzlPB0UwbhzKbqmrO2WcmsV7G\n8yWSgURDfj4/2vthYPvDEfNYlzm1W+cOptBNR8gYfhDZMrGa1Bo9STYtWq+KC3am8NE5p/CqO47v\nNIZ1ZPXMbtDoJS+rDa8dkkFHzn0riXPZeX/3X9Eq/nTqQ4mZ/G3S0jBb1r9IwQ8iHo3C6mmVXLsl\nB4/bhcmqY3aBiU251V2f3A16Wl5ZlmMehDQP/cgbW/A4E7d3l5UA/rj9b2ffiFsdXRIoQzpBpibB\nzYKbbw1sjzsRz6iTXcfz5QSuRBJ8GsssXHVkHfPK9gb2PzvtOspiU8JoWXiQgt+KYMTzcxdewtF0\nW2D77L1mTDI/XyLpE43i3dM6OeOqj/PDfU1VXd4fOZ8NmZPbPTYS4vSd0angCyE+b/b7/aE3JzIQ\nQrAxtwpLXFMp5QU7UtG5Q79cW6Z7RhAys6fPxLts3L/tP2jO5NsfSMrmH5MuD7NV4aMrD795ubhv\nhtKQSGPnD/ewenpTfn6CXcv83Skt8vNlXD0MSBGNSNrz/hWvl3vzXmNIg7+6bb3WyBOzbsSjiq64\nfXO6Eny5XLQP1MZ5WrRGHFYRQ+6RhDBaJJFED5XP/5EZFQcD27+ffn3U5Nt3RFe3upFCiA8A0ez3\nAIqiRFdOUw/Jvymf3Fdy2VtTx6RjfqGfdjCJ00kuylIcYbZOIolc5pXmU/X+K4Ht18ctYkv6xDBa\nNDDoSvCvbPb706E0JJLZPs5CSo2OoRYDKgTn70rhg3PKwm1WZDCY0hjlOoIWNA+/BHOyNKu+grt2\nvBnY3pY2jtfGXxy06w9mOhV8RVHWNP4uhEg9s29gN54dgCgqWDPtNFdsGEqMU4PBpWbBztTo7oc7\nGMVvMNocZRjdDh7c8i9iPE4AymLMPDXzO/iETEiELgRfCCGAh4Cf4g/rqIQQHuCPEdXQpBv0dYK1\nweBl9bTTXLplCCpFkFqj58u//wmM+P+yzWidZdPTRVYDnsEycToQnh6WJcobTHdRFH6x461AX1qH\nWstjc74f0fXte0pXt72fAfOBWYqiJCuKYgLmAOcIIX4ecusijAqzk7xm/XD3rf6S8cfjw2iRRBI5\nfPPQ18wva3KOnp96LccSM9o9NtLz7Tuiqxj+94CLFEU53bhDUZSjQojvAp8Dz4bSuIFKo8fdG6+6\nIKcec52O0SX+7liz95uoiXNzSk7i9j/BCNEMhKeAKKX5HMCM8kJuKvgksL1i5Hy+zp4RDrMGNF15\n+NrmYt/ImTi+XDraGwRsmlRFZaI/xqhSBBfsTCHOHtzc4EEV4gkmzUsRy7LEUUFWfQX3b/sP6jNZ\n5PvMOfw9ihdXdUZXKuPq5WsRQagWRnnV8PWMSm7ekYutxoLBrWbh9lSm/WOy7JQVbgbDDUI+VQS8\n+ziXnWWbXybW439CrjAm8fjs70X14qrO6MrDnyKEqBNC1J/5qWvcBqLUhQwOdoOXpXc9gFflF3hz\nvY75u5PbXeoWtd56TxkMYi0JGiqfl/u3vUqmzR+EcKi1PDrnZiyG0C1uLFq+ZFDH/7tKy4z8nl9h\nJGPsBDZNrGJ+vr9qX055LFMOu9k9Jjq9tg7pDyHvagx5MwkLnRVKu3Xvh0yvPBTY/v306zmSlNnm\nuMEs0MGmq+JpBiHEz4QQLwghbhVCyOekIHM420ZBTl1ge9qhJHLK2qaRSS+/D4Qzli9vFCFh8bFN\nXHl0fWD71fEXsz5zShgtGhx0JeCvAG5gHXAZMBG4M9RGRTrN5wbyb8pn8j9zSarXklFlBODc3SnY\nDKeoNEX8NIlE0mP2XZ3K0e+/F9hemzGZ18ddFJKxIu3poKsY/lmKonxXUZS/AtcC5/aDTVGHooLV\n005TG+tfeav2CRZuTyPOLiNqvSacHr3MDgoZ6dbTlNx5Z6Dc8aHETJ6Zfj2I0JcejwS6EvzA2n9F\nUTwhtiWqcel8fDmzAofW32/T6FJzYV4a2n6ooT/okGIalcS7bDy66e94a2oAqNbH8+icm3FqdGG2\nbPDQ3Sydxsycyc2yduq6OFfSQ+pjPXw1oxKv8KfqmKw6LtiZivCF2TBJ7+jPG1OEP1lovW4e2vIv\nss5k5DhVGh6b831OxySF2bLBRaeCryiKWlGUhDM/8YqiaJr9Lgu7d4OetkusMDvZ0KyGfuZpI3P3\nmaOnM8EAEi1FAY+ixa3ocPkMeBWZsxAOhOLjFzveYlLVsTM7BL+bcQOF5uHhNWwQIj/BA5CjmTYS\nbBqmHvZ7L+NOxFMb66ZgZH2L42TmTt/xKmpOF9VRZb8QiyeLGm8m9d4UGnyJOHwJ+Fp9RdQ40f/k\nfxhVtSTmziax/nskqUtI0x7G/LAJlYzA9Yqc+1Z2OEH6vf2fcUHJrsB22t13s+FgWshsibSJ2uZI\nwe8nerpqd9eYWhJsWkaWxQIwq9CEzejleLo9FOZFDYoClZ5RFDlmUuqeSLl7DJ7lecAd3Trfix67\nT4/dZ6ZqZyVwdeA1jWhgiPYw2bpdDNfvIFlTJOcS+8glRVu4/uCqwPaHI+Zx983fh/s/Dp9Rgxgp\n+H2guYj3xdtu7IzVAgEbJp8m1qFmiMWAQHDe7hQ+15dTbnb2eqzu0NyWAd13t5uhH0WBcvcYDrxx\ngGOVf8PmS+n2ECrcCHwI4cOnaNt4/M3xKEZKXLmUuHLZbL2RWNVpRhs2MM64lhTN0fCIf5jKMHTU\n3KSzhVStz59ecYBHd78T2LdlyAT+knsl93TzD9mVp94buwY7UvAHII03AK8avppRyWWbhpJo06L2\nCS7MS+Pjs0+F28RBQYMvnsKGhey3X4jFmw1rSoC2Yh+fbCDNtgqzphiT+iSJmnKMqhqMqjo0oqlJ\njT+mr8epxGHzmqn1DqXGm06VewTl7tFtbiQ2Xwq77Vey234lJvUJJsV8ynjj17TJKelL1c4IrZc/\nsqaEX239N+oz6ZeHEzNZPuu7+FShSVWO5DBOc6TgD3CcOh9fzKrgsk1DiHFq0HlUXLQtjbrTUd54\nrBPvvtYzlF32pRTaF+JB3+Z1vagnR5/HcP0O0nX7iXvsKCxb3OWQQoBWONHiJE5dxRAOtXjd6jVz\n0jWZYud0ip1TcSpNvQ4s3mzW1d/CZut3mPDWQSYvzCYx1diDNxw9uIqLeWzT3wNdqzTp6Tw8+Qc4\nNE3/l40CHS2eebCQgj8IsMZ4+HJmBZduGYrOoyLWoeHd3z6MboIKl7ZnOZuN4ZoBHarpJTWeDLZa\nr+ewYx4KLT1BjWhgjGED4wyrSdftR9U817V12KOXWUJx6mrGG1cz3rgan6LihGsKBxvO46hzDh7F\nL+5uJYY9X59k75oSxp09lJmLc5Dpbk14Tp+m+Ee3YHb6ExSsGgO5f/kL1f881MWZku4QMsEXQmQD\n/waGAj7gJUVR/hCq8SKd6kQ3X0+vZNG2NNSKoOpkMQttqXwxuxxvJCzI7UMqptWbzDbrdexvuLCN\n0KdojpIb8zGjDRvQqfqvyYxK+Biu38lw/U5cPgMHHeezx3Y5Fm8WAD6fwv4NZRzYdIoJ+h8zO+4N\nYgZIOmq48FqtFN96K+7iYsCfa79s7g9YOW4sIAU/GITSw/cAdymKskMIEQ9sF0J8oShKQQjHjGjK\nUhysn3ya83enAjDUYuDcXSmsmX4aZTBng/RS6DyKlh22q9lhvQZvq9DNMN0Opsa+T5YuP/STpY0x\n9A7eh07lYFLMZ0w0fsZJ1xS2pz1HyQH/alGfT2FfwyUcdJzLrLj/MjlmJWrRx0Xtg7DZutbr5uT/\n3YGzYD8AXqHit7NuZF/KyDBbFlmETPAVRSkDys78Xi+E2A9kAlLw+8CxTDsxzmpmFZoBf0llV76P\njbnVbZqhRzJFzumsq7uFOu/QFvuzdLuZG/caQ3RdeIS9Dd30IeQjBGTrd5P98+mUHLCw9YX/UOqe\nCPhDPRvrv88++8Wcn/AS2frdvRpjMKJSfNyz/XXspS370W5JnxjY7umkarRMwvaUrkorBAUhRA4w\nDdjSH+NFOvtG1rN3RFNli7En45m139TpatxIWaRl8ybxieVeVloebCH2qZojLDU9zJXmZV2LfbhZ\nlkjmOBNXmX/N5abHMKlPBl6q9WbwgWUZq2rvwOGLC6ORvaf1RGqnE6uKwk92v8v8ZmL/8lmX8fnw\n2aEyL6oJ+aStECIOeAf4maIobervCCFuBW4FGDZsWKjNGTR01Sg9b7wFvUvFmDPN0CcWJeDS+iK6\necphxzzW1N6Go1lVD72oZ278a5xl/KLlRGxf6CJEEyyEgOH6HWSl7GavfTHbrN/Cqfj/PwsbLuS4\nczrnJfyN0SG1IowoCt8v+IQlRZsDu94ddR5vj1kQRqMim5AKvhBCi1/sX1MU5d32jlEU5SXgJYCZ\nM2cO2ooxfc166fH5AjbmVqH1CnJO+VfjTjuUhFvjo2BEfRcntyXoTwB9XfDT7HyHL541dbdw2NGy\nOvd44yrmxf8bo6qXdfz6EJ7p83jNxlULL1NiP2KMYR3r6n/EYcd8ABp8Jj6ruYfjrxRw7rfGojMM\nnKS6jhZWdfccgOsPruJbh74KbK/Kmu5vPt6DSZdoXDzVF0IW0hFCCOAfwH5FUZ4J1TjRjKKCtVNO\nU5LSENg3e7+Z0Sdiw2hVcCl1ncWbp59tIfZxqkqWmh7mwsQXei/23aGfJz1j1LVckvR7Fif9llhV\nUwG9wk2n+O/j2ygvipwCtVcdXstN+z8NbG8eehbPTv8WiuiXKHPUEkqX4RzgRiBfCNFY+egBRVFk\nEYwg4lPD19MruWhbGkMsBgDm5Sfj1iht6u409+K7Wz4hqHn7PfCmFUWw3XYNW63fbpFqOd64ivnx\nL6NXDfCaQn14ehhp2Eqmbi9r627loON8AGorG3j3qe3MuXIk0xSBEM0ehgdZOufiY5u4be8Hge0d\nqWN5YtaNeEO0ilbSRMhup4qirFcURSiKMllRlKlnfqTYhwCPRuHLmRVUJfhXJqoQnL8rhaxy48BY\nYNXDkscNvgQ+svyaLdbvBsTeEKflsqTHuTDxhfCJfT96/HqVnYuSnmNR4nNoDf6/gc+nsOm9I3xS\ncw8uXy9W6Q6AG8PC4u3csbspurs3eQSPzrkJt1obRquih4ETFJS0S1eTt424tQqfz6pg8eahJNm0\nqBTBgh2pHNm+tT/M7D2tRKjCPZJPLPdjbVaXJl1bwMW/+hFxf8jr4bVDPPna0Q0giDeGccY1DP3l\nm3zxcgHlx/whnWPOubxdlc1i03LMmpNdXKF/6E78fH7Jbn6x401UZ9LJDiRl8/DcH+LUtC1/IQkN\nMmAWQTj1Pu55+jXqYs70xlUEHz7zBJkVhjBb1j0ONczn3aonWoj99Nh3uMr8IHGmfngP3RXqZbX9\n6u0npsZw9S+nM2VhdmBfjTeTt6ue4rDj7H6zoy/Uf/UV9+a9hvqM2B9NSOfX827Bru34/7Vo+RKZ\nTx9kpIcfYcSbU/h0TjmXbhlCgl2L1+Nh4Y40vppRQUlq/5UW6AmKIthivYHttmsD+3TCxkVJz5Cj\n3xFGy8JIq4we9bJa5hfMJC1xPl/X3oEHPR7FyGc191AT9xozYv8X9BXFHXntPRXhs0v3cvLO1wKN\nx0/EpfLAObdi1cX02UZJz5AefgRiN3r5bE459cYznr5PsHB7GhmVQfKS24vJ9yRO38xDdvv0fFJz\nbwuxT1KXcG3yPdEt9h0w1riebyTfS4K6LLBvi/U7fFX70wHZgnFeaT4PbPs3uP2fxdLYZB445zZq\n9fFdnCkJBQPvEyLpNu02TjmDzejl07nl/HDvdOoqK86IfiqrZlZSljIwPH27N5GVll9R4RkT2DdM\nt4OLk37fcmK2q/BJ69fb3IzCtBitr3n+HZyfoj3ON5Pv5tOaeyhxTQag0LGQOm8ai5ffjuF33W/w\nEkrml+zmvrzXAjXtT8amcN/826kytnxPRcuXBCWHXoZ/ukZ6+BGMzejluoeewGrwF+PS+FRcmJfa\nxtPv06KrXjYdr6mw8071b1uI/ZSYFSwxPT7wUy7DQau/sUFl4wrTY0wwfhnYV+qexDtPbafOE7p+\nr92ljdjHpbYr9pL+RQp+hJOYNpRP55a3FP3taWSXdy+tL/9YMfnHioNq06ljtX5h8qYDIPByfsJf\nmJ/wr+CUR+jnSdVwoRYeFiT8iblxrwb21ZT7b6RV7tCVKenKGz/v5K4WYn8iLpV7pNgPCKTgRxAd\n5dxbYzx8NqdJ9NU+f8rmiNJNZ/LGAAAgAElEQVSWk2b9kbNf9I0jrHhmJw6rP6arwcmlSU8xKeaz\nlgdGumB396bUxTFCwIy4d7kk6XeocQFg95l5r/pxSh+Y2871/KKbc9/KwE9faX6NBSe2c08zsdeN\nHMm982/HYpBtXgYCUvCjhPpYD5+cfSqQsqlSBOftSmlThiGUon+oYT6f/Dkfj9svBgZRx5Xmhxlp\nGOBrBQYBow0bucL8KLozi7ScShwfVD9MkWNG24NDtC5hydEN3LP9jUDq5fH4IQz/9yudir2Mu/cv\nctJ2kNMTgbYZvXwyt5yLt6ZhsuoQCObnp6DxCgpzrCG0EgoW7efr1wr9ncCBhBQDV/B/JGlKQzpu\nnxlETxqZun1c9ePpfLj8Cxp8SXjR83HN/SxMfIHxxtWhG1hRuO7QV9xc8Elg17GEdB6Ydyu7UgbG\nBLLEj/Two4wGg5dP55RTleAK7JtbkMzEo6F75N5tu5yv/1MYqNdvSo/lml/O6FrsoyQW32M6+Zuk\nDovnGvP9JKhPAaCgZlXtney1XxIaWxSFHxSsbCH2haZh3DP/dmoMMvVyoCEFPwpx6n18NqecyiRn\nYN+sQhPr3ngFRemgQnUPwgCNWT+KAnnWa1lf/8PAa6nD4rn6rmnEJkXJcvrGm1Y/3riSNKe4xvwA\nyZqiwL41dT9mj+0ychyvk+N4PSjjqBQfd+x+h28eWh3YtytlNPefc1uPF1XJVbX9gxT8KMWl9fHZ\n7HJOmZty8re+/zaf//V5upso01k6p6LAZuuNbLF+J7CvLP4oTw75CcY4Xa/tlnSPWLWFq8y/Jk3b\n1P1rXf0tzHQEJ4qr9nn55fY3WjQviVu4kIfO/iEOWRtnwCJj+FGMR6PwxawKLtiZQnaF3yPb+/UX\nLEhLZc2003jVrbz9bnj5ua/kggKbrDey03ZNYH/WeBN/T/gzHrWrk7Pp/xBOuENGvS3w1o1FXQaV\njaWmZXxkeZBT7vEALHBoUQFbDb1vlG7wOHlg67+ZVXEgsO+rrOn85A/P4X7w815fF+QkbqiRHn6U\n41UrfDW9kkOZTZO2wypiuHhrGjp3Lz4eCswuvryF2Ofot7Lk/yZ3LfaDnQE456BX2bnC9Ajp2oLA\nvvMdWs7upaef5KjnyfV/biH2K3PO5ukZ1yO0ssTxQEcK/iAk/6b8oKZPKirYMLmKWVc21bMZYjGw\neNMQ6t0dh19aL8hSFJh9YgnTSy8K7MvRb+XSpN+h0crmFr2mhzeRxjh9Y6xep3JwhelRMnVNn5n5\nDi3zGnom+hnWSp5Z+wJja5pKMr8+dhEvTLlGdqoaJMj/JYkfAefd8H0u+N4tgV0mq443iqZQ5ex6\nAk5RYKv120wvuTiwr1Hs1aL34YOII0xPAVqVkyWm31Ck8Qb2nePUMqcDT7/1gqyxlmKeWfsC6XZ/\n60UvguenfINXz7q0qQftAGiwIukcGcOPQPri/c9YciUxiYl89MLvUCmCeo+BN4qmcEXWfobH1nR4\n3jbrt8izXRfYHq7P84v9I1UdntNjmgulFJeWdOMmohUu3o11cbVNxwiP/4nrPIeWXV8WM3VRx6UY\nZp8q4P5tr2Lw+hftOdRals/8LlvSJwbHdkm/IT18SRsmzL+AL2dW4Fb703WcPg3vFk8kv2ZIu8dv\ntV7HNtv1ge3h+jwWJz0pPfsBQPPQTo7jdbwC3o91cbyZp7/hf4fJv/v6dgvhLT62iYe2/Csg9rW6\nGO4/5zYp9oMUKfhRSldPAaWpDj45+xRxGn+uvg8Vn5eNZV1FDs1T9XdYr2Kb9duB7WG67Vya9JQU\n+wGMR8B7sS5OqJtEf239rVxR+1bg5qBSfNy2ZwX/b/c7gbo4p2JM3HXuHRSac9pcs7vZNTILJ7zI\nkM4gIRiTtK2v0bxfbv6x4mYdlvzhgeoENzfk7OK9ExOpdMYBsLUqmxqXgUszDnLAsYhN1psC1xum\n28Fi05NohLvPtnbJAMuG6TP9/H7cAt6Jc3GdVU+G1+/3XdKgxSfAa7WybPPLzCovDBx/KDGTh8/+\noSyCNsiRgh8F9OVmEa91cf3wPawsHcdRazIAB+tTqSxKwW74Lo3JGZm6/P4Te0lQcAv4X6yTx+wl\nVHpGIxBcatey6QcPtxD79em5PD3jetlsPAKQgi/pEp3ay5VZBXxdPopdlgwALE4BrjfQxV3JUGMt\nlyU9gUZEeJ79QKPVU0FvSh07VfBQTCbfsvpI86lQIdhjupTc5GJSqvby5tgL+feES1CEiiLDDf5x\neliaoXUYR4Z1wocUfEm3UAm4cOgRhCqJnVVn0jQVK676NxgXdxSdamC0TZT0HIcKymu3kKPKwR6b\njiLU5E/8EYUNhXyQOTrc5kmCiJy0jXaWJXa7o9Up1xgOK3ehjbsaROPjvY8vynJYUz4CXwd111rT\np5aKkqCi9nm5ffd73Ln9P0zb/UcMDZUAKCotI2InMdQjwmyhJJhIwZd0iyr3MD6yPIhbMaLWjsBj\nvpwkbVO1zbzqLN47MRGHt28PjfJm0H+YHHX8dsNfWHpsAwB6Vy2ph98gJt6fo69DcK1NT7JXin6k\nIAVf0iU1nqF8YHkYp+Kvb24QdXyU+xrfHbGdkXFNC6uKbGZeOzaV046uV+ZKYe8hfSyz3DruPr76\nOM+vfo7cqmOBfesyJnPP2d9j6c9n0CD8j2tGRfBNq57EZqJfZLgh8CMZXEjBl3SK1ZvMB9WPYPeZ\nAdAJG1eYH6Emphy92stVWQXMSW4KCdW4jbxWNJWC2jTAL+z9Ke79Pd5gZHHRZp5a9yIpjjrAXybh\n5bMu44lZN+LQ6Jnx/Br+F+vEdaZjTbwiuM6mw+Y1hdNsSRCQk7bRTBflCbYeqeG96sep9/nFW42T\nJabHSdMeBfxL8YWA+WnHSTPY+KR0LB5FjUdR80npOJ5TZ6D2VuPtoG5a7gj/Nbo7hxCtNM++aZ7h\n0tOsHL3Hxe173ueS4qYewnXaGJ6c9R12pI1rcewpjcK7sS6utenQIEjyqfjAsoyrzb/GoKpve3FZ\n6mJQIAVf0i46j4EPqh/C4s0CQIWbxaanyNDtb/f4sQmnMevtfHhyAtUuf0hn3Il4Ump1fD39dI/H\nb/TSQ9lUPZoYVneK+7e9Sk59eWDfkYQMHptzE+Wxye2ec0LrY0Wsi6tsOtQIqj3D+NDya640Pdxl\nVpZMvRyYyJCOpA1ul5fFhbdy2jMKAB8+Lkp6juH6HeSOGBbwzFuTorfznRE7OZpuC+xLrtOzdH06\nh7dtbvccSei5sDiPP6z5QwuxX5U1nbvOuyMg9kXLl7Qblz+q9fFJjBvlTHinwj2Wj2vux6PI2veD\nEenhS1rgVTR8+stnSXdND+xbO/Itfmrf2K3zdSofa6eeptzsYHaBGbUi0HlUrHj6N8zKMbF9nAVf\nL0vjdxSbj8SngN4somqN3uPkJ3ve4+LivMA+odfzzISlfD58dlNZ4y7Yr/OiV+CiBn9vhBLXZD6v\nuYtLkn7XZxsl/Yv08CUBfIqKL2p+TnEzsT8n/mUKh/TQOxdwYLiVj88+hdXYVERtYlECl29MJ9Ha\nfT9DTsD2jmF1p/jDmudbiL0uwU3Of//L5zlzui32jezSe5kb92pg+5hzDl/V3gHdXHshGRhIDz/S\n6WgyrVV6n6IIVtfdzhHnvMC+mbFvMTX2QxonaKFjb7o9Ya5KcvHBOWWcuyc50DPXXK/jivXpbD3L\nwsFsKwQhxTsS4v0dTcz2FKH4uOLoBn64byU6X9PNdlXWdH4yZyWqcWOBQx1fAFqEdZqnc86IexeX\nEsuOM+0rDzouYJHPw5fGduonLUuMvAJ3EYD08CUoCmyo/z77GxYF9k2O+ZDZcW/2+dounY9VMypZ\n+IMf41X53UGNT8W8vcks2JGK3iU/gsHC3FDLY5v+zu35KwJi71BreXbadTw949uotH13x+fGvcpE\n42eB7WkuDef2sj+upP8J2bdNCPGyEKJCCLE3VGNIgsM267fYbV8a2B5vXMX8+H8Gnvr77DkL+N6p\n+/lwXhmWuKYCa8PLY1i6Pp3jtqROT5e59V0zrzSfF7/+PTMqDgb2HUnI4M7z7+xRvL4F7XjoQsB5\nCS8xxrAmsG+uU8t26zVtjpUMPEJ5a/4X8ALw7xCOIemMbuRG77Jd0aJb1Sj9RhYkvIgQwQ/O1iS4\n+eicU8wsTGLCcX9d9ViHhv8V5zLFVMp5ace6uEJk0pcJWqPbwW35H7TIrfcheGf0+bw64VLc6uB/\nxVXCx4WJf8StGClyzgZgs/VGdCo7RTFy9e1AJmSCryjKWiFETqiuL+k7BfZFbKj/QWB7mG4HFyU9\ni0r4QjamV62wZaKFkhQH8/ckY3D7U3Z2WzJY50wlzRTEHrgRTM59K5laeYif7fgvQxosgf2VxkSe\nnv5t9qS2rHLZPC4fjBx5tfBySdLTPFD1BsPP9MddW3cbOtHAOOMZ779VQx1J+Al7AFUIcasQIk8I\nkVdZWRluc6KGQ3nlfF13e2A7XVvApabg9KHtTgjo5JAGVpxXSnGaPbAvwa5l8eYh/ODe81DLgl0d\nctZd7/DTXf/jtxv+2kLsV2dO5ScL7gqIfagXP2mEm/diXZSqmxyEVbU/5ahjdkjHlfSesM+2KIry\nEvASwMyZM2WSV7DoJJxT5JjBly8X0Hi/Tx0Wz5Kf/xjtk7/q9XC9ifM36H089/wqrn7ibOYUmNF5\nVAgEE4sSyKo0sn5yFZUmZ9cXGqD0tiRCZ69fc/Mz/HnX26Q11AT21WuNvDj5alZnT+/wvFDhFvBO\nrJNl9nKqPDkoqPms5pdcbvoN2fo9/oNae/rNP5vS++9Xwu7hS/qXEtdEPq25G9+Z4vUm9Umu+OkU\n9MYu7v09rNTYnRtA/k35CCE4kmXj/XNLKUlpCLyWaNNy2aYhzN1rRuuW3n6sq4Gf7fgvj2/6Wwux\n35g+idsuvLv/xL6dz4BDBVeYHiFRXQqADy0f19zPKde4NsdKwkvYPXxJ/1HuHs1Ky6/w4m9eEq8u\nZ6n5YYzx32v/hH70vuxGL1/MqmDMiThm7zeh9fq9/fHF8WSXG9kysZrioQ1dXyjSUBTOK9nNbfkr\nMDubipbV6mL48+SrWZM5NZCBU7R8SVBW6DanZVio489DrLqGpeZlvFf1BFZfCh7FwEeWX3OV+UFS\ntEVBtUnSe0Im+EKIN4ALgBQhxEngYUVR/hGq8aKCPjwKV7mH8WH1Q7gVIwAxiTqu1C4jTl0dTAv7\nhoBDw6yUpjYwd6+Z7Er/Yq1Yp4aFO9I4PsTOlonV2A3eMBvalmAtnGpOuu00P9n9HjMrDrTYvzZj\nMi9OuZpafXxQxgkWCepKlpof5r3qJ2jwJeJU4vjA8jDXmH9FkqZUVtQcAIQyS+fbobq2pGe018Bk\n6Z2LSHzpVJgtax+b0cuqmZXknIphzj4zRpc/C2R4eQzpVQZ2jq2hcFg9ShADkj1drdso8KGYGNV6\nPXzj8GquP/Al+marZasMCfwl90rWZ04J7BtoVSlNmlKuMC3jffsfcTV4aPAlsaJ6GdckP0C8up2q\nqXJFbr8iY/iRxrLEFp5U7f+dYkX1Y4EGJlph5wrzoyS/lBkuC7uHgKJ0O++dX8qB7KZQhs6jYk6B\nmaXr0xlape/kAgOTrkIukysP88LXz3DT/k8DYu9DsGLkfG698O4WYj9QSdUWcfkdU9AIfwllqy+V\nFdWPyAYqAwAZw49g6rypvP/sDqy+FAA0OFlieoI07ZGuTx4gXpdL62NTbjVHM23My08m0eYvy2uy\n6rh0y1COpdvYNt6C3egN6WrcUNfrGWKr4kf7PmJ+acvrH0rM5I9Tr+WQKbtX1y1aviSkoZSO2hym\nj0rksqTlfGT5FT601HozWFH9KFeaHyRWXdPuOZLQIwU/Qqn3JrOi+lGsXn9aoxonl5meIFO3L8yW\n9Y5ys5MV80s5qyiBKYcT0Xr9D6cjymLJqjCyZ1Qt+0bUBUov90Sgw1m2weBxct3Br7ihaB2Kq6ns\nhF2j55UJi/lo5Dx8op8fxIN0s8/W7+aSpKf5tOYeFNRYvFmssDzCVaaHiFE3G0Mu0Oo3ZEgnArF6\nzayofpQ671DA363qMtPyprzoQYpPDXtH1fHeeaUtmqxovSpmHDRx1boMhpfFBK1kb859K1uEYOr3\nLw9aFoxQfCw4sZ2/ffkk3z64qoXYr8qazq0X3sMHo+b3v9gHmZGGrVyc9AwC/0S7xTOMFZZHafAl\ntH9Cq5CkJLhIDz9SOPMlsXmTWFH9KLXeDABUasHi+CcZpt8VTusCBCMksuXH/vey6JmZzN1nxmT1\nN+ZIsGtZsDOViiQneeMtVJidLbz31mO39uybC3p3J0PbuwHk3LeyzfnNj5t0+gg/3LeS8ZaWvXwP\nJGXzl8lXUmjOaZNi2V172j0uzJ7zaMNGlETBF7U/R0FNtWcYK6of4UrzQxjb648rCRlS8CMIuzeR\nFZZHqPH6J2RVeLjklmnkvL89zJaFhvJkJx/ML2NccTzTDiaiP1PTJa1Gz2Wbh3J8iJ3t4yzUxTVl\nuoSzdn5ObRnfL/iYOeUt+wKrU1NI+8VdXLZZjRJOjz6EN4Yxxg0oCL6s/RkKaqo8OayofoSrzA+3\n3xRdEhIG9/OiJIDVa+b96t9g8fiblQi8XJz0e0ZOTQ2bTfk35fdYWHt6vKKCwpx63rmglL0j6gI1\n98GfxnnVugzm7DNhdPbsox7MuH6avZq7tr/Bn75+poXYu1Vqkm+5hVGffErS1VeFV+z7SnsrsVtt\njzWu58KbcwNNb6o8I1hRvQyHL66fjJRIDz8CqK928H71b6jzpgN+sb8o8VlGGaKncbhL5yNvgoXC\n4fVMP5DEyLJYAFSKYMLxBMaciGNNzMvonSqc+tBVA21OotPKdQe/4opjG9D6mhaL+RB8lT2dV8df\nwta7buwXWwYK4+YMRVEUVv1rH6DitGck71c/ylLTI8SE27goQAr+QKOj1bQdTGTVVtpZ8ewu6s+I\nvQoPFyU9w2jDplBa2a/k35TfbY/bGuNh7bTT7BtRx6xCE0OrDYC/y1beh+9yrTqT/Tn1jLv7HZzq\n5QDET7gvqPYmOuspf+p3/OvzVzF4W7b/2zJkAv86azFFiRltzmsdfx9oi6ra0JMQULNjx89NR3n3\ndr6quwNQUeUZwXvVv+FKiwPp64cWKfiDlWWJWDwZrHD/DVutP8NDhZtLk55ihCGv4y/jstqIyILw\nT2gu71Csq5JcfDqnnKxKI9MOJpJc51+kpfWqmHwkkfHiP+xKnMzOhK4XMuXct5L4CV3blOSo59rD\nq1lybCPVXjeGZq/tNw3n5YmXsTdlVDfe3SCmmzeBCTFfoxI+VtX+FAU1Nd4s3n16B1eKISRqyltd\nU6ZtBgsp+IOUKvcwVliW0eDzi70/z355+9k43f2iDJAvVG8nVNs8CQg4mdbAydQGhpUbmXYwKZDR\no1PczK7ZztTaPRzQJFEwoo7eTh2aHHV849BqlhRtauPRH0nI4LXxF7MpfWLv2gxGMOOMa9AIJ5/X\n/AIfWuqrHLynepyl5mWYw21chCIFfxBS6prAx5YHcCr+B2CNaGBJ0hNk6WX74PaoL1zOPmBfikLu\nmMeYejCJpDMrdnWKm9yjiZxVlEBB7Bp2JE6jTttBjngrMusr+MbhNSw6kdciRg9wJDGD18ZdzOb0\nswb3ZGyIGWXYzGLTcj613IMXPTZfMu9V/YalxfWkDmtVHK513R1ZV7/HSMEfZBxzzOKzmrsCJY61\nws4VpsdI1xWG2bLw0q2nAiEoSrdzfKid9LwfMLNmB2a3v2OU2ifIrS9gYv1+DsWOpjBT26HHb9+x\nkwe3/Iu5ZftQtVrldTgxk9fGX8TmoZ179AM+Pt8ZQRbXHP0OLjf9hpX2J/A4vTiURN57ZgeLf5xL\n9njp6wcTKfiDiAL7hayuux0Ff765UWXhCtNjpGojv/l3TyZuu0IRcCBuLAdixzDSfoyZNTsZ4qoA\nQIXCONshxq3P4IThA3Yn5FIUMxxbwRPMPbWXaw+t5vj7Rcxrdc1CUzZvjl3ElqFntRD6UNSoj0Sy\n9Hu5UtzNh64HcSmxuB1ePvrjbi68aQJjw21cBCEFfxCgKAo7rNew2dqUwpegLmOp6VESNQOzxHEo\nCPpiKSE4GjuSozEjGJv2CBP2zGCY42Tg5WxHCdmOEjyKlpzTNUw4VY7W19Kj3zJkAv8bcwF7k0fK\nGH0fGao7wDXmB/jQ8hA2XzI+r8IXLxdgjb+KaTHvyz9vEJCCH066kS3jfTiZNXW3sb+hSexTNEe4\nwvRYywJUgWv2LlUulIRjVWtH1O9vJ7NHCMpSHBxMv4I0ZwXTa3Yy2n4EcWaFkEa4OZkaS5k5h+zq\nOjKqrGwdkss7oy+gOGFoGN5F5JKsLeYbyffxoeXBwCLCTfU3YfMmc078P1FFQIZZOJGCP4Bx2Nx8\nanmIEldTKCNTt4fLkpajU/Wh3Z+c4GpD/f7l6LxuJlacZNGxQ4ysP8nxlAROmBNwa/whNK9aRVFq\nEkWpSZw0GDGo6lD7UvCq5NcomMSrT3ON+QE+ttxPmXsiAHvsl1PvTWVR4nPoVI4wWzh4kZ/UAUqN\nJ52VT22nppnYjzN8zYLEF1ELTydnDj66egLoLAbe+tz6/csDv3eUo9/8GIBRNSe55PhWLji5k3h3\n0410fFk1o8st7B4WT3FiNnqahCbLUUqWo5QG1XoK48axL34CFl3bCUb/vMPyDvZ3/d47K/7WG4J9\nvV7TWbbNskQMKhtLzY/wZc3POOL0z5gcc87hnerlLEl6ggRNRT8bHBlIwR+AnHBO5rOaX+JU7IF9\nc+NeZXrsuzKOGSTiXHYWntjBxcVbGVVb2uZ1jwrWZUzjg5HnUDLvz9Tv/z7ZjpPk1u1jhL0okJ1j\n9DmZVreHaXV7KNUPpTBuHIdjR+FU97wbVzgLuw1ENMLNxUm/Z2N9JbvtVwJQ7RnO21W/41LTU4O2\nt0M4kYLfX3QjZ1hRBNtt17DV+u1AJo4aJ4uS/tC9UgkyVNMpRx9ZxHdveZYFJ3cw+9R+dL62T0pl\nMcl8PnwWnw2fjcXgz8ePFwKE4IQxmxPGbGI9NiZYCzmrfj+JnqbkzQznKTKcpzi/ah3HYoZzIG4c\nVb0o29NX4R8wXnwQUAkf8xP+RbLmOKvrbseHFoeSwAfVy5if8A8mGT/1O0FdxfbldwOQgj9gcPpi\n+LL2ToqcswP7YlTVXJa0nCG6Q2G0bGDTXu35Rur3L0fl8zL99AHOP7mTQ+cs40Grtc1xTpWG9RmT\n+Wz4bPamjOxyoZRNE0te0gzyEqeT3XCSifUFjGo4gkrxP36p8THafozR9mM4VmVxyLCGA3FjKdMP\nbZPJk/tKboei3F4aavObwaB4IgiS0E6I+ZokTSmf1NxLg8+EDw1r626jzHUWFyS8KOP63UQK/gCg\nwj2Sz2vuCjQtAcjQ7uPipKdl/89eoPJ5mVR1jHNK8zm3dDcmp1/kWzvbB5Oy+Gz4bPIWrmDrba9x\nB103GW+BEJyIyeZETDapIx9gRGkso0piSalrCucY3Gpy3QXk1hdgVcdwJHYkh2NGYVX86wEk7dBB\nfD9dd4BvJt/NJ5b7qfT4axIdcpxLpXsklyY9RbK2uPWVJK2Qgh9GFEWw03YlW6w34EMb2D8lZgVn\nx7+KWng7Obt/6E3XpXBc3+dwYNu4kZ/veJM5pwpIdNnbPa4kNoXVWdN48Nmfsfglf236eMMHvR63\nEYfex/4R9ewfUY9m972Msx5knPUgCd6mJ4o4r50pdXuZUreXhlVZFA+xU5Rux+f1olKrez12JIVw\nuiJeXcU1yQ+wbuw6Ctb5515qvJm8XfUU8xNeZqLxcznP1QlS8MOE1eLgS8sySlyTA/u0ooEFCS8w\nxrixexeJwLhkT1oMJjhtzKg4wMk7v6Diy68xel1c3M5xVYYE1mZO4efLf8ri/xwHIfjNyJGAX/CD\nLZIWnYnN5jlsNs1m9NBlZO9bwCjbMYy+prCD0aVm3Il4xp2I58+3fIfhU6Yz0hZLSWoDTl3/1Osf\nrGiEiwXfGU/G6CRW/2s7HsWIFz1r6m6nyDmLBQkvEqu2hNvMAYkU/H5GUeDAvVexvu5mnEqT2A/R\nHuCih24iMXVJxxNQEVLauLeoFB9jLCeYWV7IzIpCxlpOokKhHjC2Ova0IYFN6ZPYkJFLfsoofELF\nA7m5IPr+2N+Y7pl/U36bEFDrEhDlwOGUC1idfB4ZjjIm6F5jeLmRGGfTV89hs3Jg41rOIwUfCpUm\nJydTGziZ1oAl3h3oENWb0hKdzREMCjrpCTFuzlBSP7yHz2rupvrMIq3jzpm8cfoPnJ/wl+47TlGE\nFPze0Mv63HWeNFbX3c4J19Rme33MiH2HWXFvoU69vZMxB5c339Nm4B3hLivjwuI8ZlQcYHrFgQ5D\nNQAn41LZmD6JjemTOGjKbnfytbk9vbWtuYAWLV/SLSFWhIoSYyZ1E6p56cldXPSHOeScimHYqRji\nHE1fQxWCIRYDQywGZhw0Ydd7KEt2UJbioDTZgd3YeZgvmK0ZBwNmzUm+mXw3m+u/w277UgCcSjyf\n197NYccmzk34B3Hqqp5dNILr70vB7wtdpVqeed2raNhtv4pt1m/hUZraYsSry1mU+DwZuoIOrj84\nP3C9LRaWc99KkhtqmXz6MM9e8l+mVB4h3V7FLzs43oug0DycvCHj2Zg+ieL4Ie3Ws+mtsDeel/vK\nfS28+o7orietUqmpMDupMDvZOsHCVws+4M6/f5esSiOpFj0qmt5DjFPDqNI4RpX6S2HXxbgpS3Zw\nYNN69C5Vt8I/wfby27uphPMpQiNczE/4Jzn6bawST2GtdgJw1Hk2J05PZXbcG0x+2IRKtPpbDdLv\nV1+Qgh9MWnkGigJFzplsqL+5RQaOwMvkmJXMiXsdrcrZ72YGy/vuK4rPR1Z9BROqi5hQfZzJp4+Q\naTvd6TnV+njyhowjb0sh1loAABBJSURBVMh4dqaOxaoLfyfUPomdgIWrl8JoyB9dh96lIrPSSFaF\nkczTBvTulpO5CXYtCXYtHz23nG+TTU2smwqTgwqTkwqTk7pYD7QzaRkNnn+Wfi/X3zuHDQ88xP6G\nRQC4FSMb6n9AYcMCzon/J9n6QRzeCgJS8EPBskQq3CPZXP9wq/ANmDXHWZjwp97n1gfRK+ksh73x\n9e5ep5HW12v+2pFfnUfDnj007NpFw67dnNqcx9/cndcEcqi1FJhHsDt1FHlp4zmamOGva99snGCW\nHw7FTTB+wn3Nbgr+63dU7tmp83E008bRTBtCAXOtjvQqA+lVBoZU69H4Woapkmxakmxaxp70Nwtp\n0HmpNDmpSHJyOslJVYILt1ZpM05zuuOxd3XDaHyK6FXGUBA7sumNGhY++wTjDlpY8/oBLKf8IcAq\nzwg+sDzKMN0Ozo7/Nyna490bMzB2ZDRbkYIfLM58ICrdI9hqvb7FAioAvbAyK+4tJsV80pRu2fjB\naT0RO4g/UI3Eu2yMqilhVG0Jo2pLGVVTwsEVd/sfexqPaec8p0pDQfIILr5hMd/b4eWgKRtPD4qT\nRVL9eUX4e/NWJbnYO6oOlRfSavT+G8BpA8m1etStkvmNLjXDymMYVt705FMb46Yq0UVVopPTiS6q\nz9wE+sPr79fFYc2+N5ljTXzr17PZdd8PybNeh+dMw6Bi13SKq6YyxrCOGSVWkjOjq226FPye0EGG\njKJAiSuX3fbL2wi9wMtE4+fMjn8Do6q3XVN7Rqhy59sLBSkuF8PqTpFdX8GIujJevuCfvFJbQlpD\n9xaM1epi2W8ezhXXX8z385wcMA3HrdZww+1LKOijcIc7ZNVTuvKqfWr44hd55L6Sy86xtai9gpRa\nHY/n3MMbq/5CWo2+TQgIINGuJdGuZWRZLAAKCnUxHmri3VjiXf5/49zUxbpRVC3H7otQhzuMpNao\nmBH3LuOMa9hqvZ7ChgVnSpaoOOQ4n0OPbSVHv5UZse8wVHcwrLb2F1Lw+4Dbp+eg4zz22JdQ7Rne\n5vXRhvXMinsLs+Zk25MHiBffVVinEb3HSZa1kmH1FVQ8dwjXkSM4jxzFVVzMXz3dq97pFSqOJaSz\n3zyc/ebhFJqGUxabTNGTlwOwN0ie+WAT+t7iVSuUm53Mufo6flT3GCiQaNWSVqMn1aIjuU6PqV4b\nKPnQiEAEbgLDmz0JeIVCbZybmjM/9bFuyo8eJmloetBt70/PP05dxcLEPzEl5kM2WW/kuHNm4LUi\n52yKnLNJ0x5iovEzxhjWh2Verb+Qgt9DfIqKEtckDjRcwBHnXDxK6wxwGKXfyKy4t9ou9e6GyAfb\nO28U9K7CHIqiMONnb5Bur2aorYp0WxW/tFWRbq9iqK0as7Pp6aRqe9fjulQajiWkczQxgyNJmTx9\n/zeZ+p9jODW6Xr2Pzt5DqAW+t6LUk/N6kv3TXvy/vfNzX8lF7QVTvY7kWh3JtXpSanUkWdveBADU\nisBcr8Nc3/R/9J9dPwPgW7os6mLd1Md4Av9ajR5sRi92vRdE5x59R5lCncX8uzUf0APHKVlbzOWm\nxyl3jWG77RqOOecGXqtwj6HCPYYN9TczxrCe0YYNZLSX2TPIkYLfHq2ybVwNHk7sr6ao9g6OO2fQ\n4Etqc4pGr2a8+iMmx6zEpCnp0XA5jtcBKOqluV2KudeLuaGWFEctKQ21JJ/5t/H3w5c8j+dUOa85\ne+fZlBtNFMenURw/hKOJmRxJyuBEXBo+VVN44U9Tp+J8s2d/l8FKsLzWvl6n8fzcV3I5neQC/GUe\n1F5BolVDklWHqV6LqV5HUr22xXqA1hhdaowuNUPaWcDqFQp2gwer0YvN2HQjaNA3+9GFv0xII0N0\nh7hM9yRV7mx22q7mkGN+oLSJS4llX8Ml7Gu4BKOqlpH6zeTo88jQ7Y2IAm0hFXwhxKXAHwA18HdF\nUdp2ghhgKIqC3Wui3D2Wsv8douxILZXF9fi8CnBhm+NN6hOcFfMlEx55Gf1TS9u/aBDDN4qioDQ0\n4K2vx1dfj9di4Zz/397Zx8hRl3H88+z73u7dXm2vjZZWfCmkvCi2DbTBII28yEsgmKol8qYIikqi\nGCO+JBqN/iFBSZEXESpgACnEwIE1BATEiEVqaY+WQPoiXEto2qNly93t7s3uPv4xcy+93t7N3e3c\ndW6fTzKZmZ3f7Dzfndlnfq/P7+0OWvu6yZW6yZV6yPX1kCu97617eP2JH/BApfYfzhm9JyQAZYmw\nNzObzuxcOlvmsTs7j86WeezJtlGMjRz7vVbOe6INqzOtqmYyDn0ypYeT7zuZAzmHAznnsDRL7v7E\nwEugpSdGS0+clt4Yzb1xotXaAWqiKjQX4jQX4jXTKMrtL36ZTOssMq2z+PS7sykmKhQTVUrxKqVE\nhd3bOkhlm70lC8qIXUzrxez4bs5qXcPp1T/yRmEl23rP4b3K/IHjhWpuwPlHKDM3vp35ia3M27Kf\ntoUtZFoTSMgC9wTm8EUkCtwGnA3sAV4WkXZVrTHKKECGdakqOxWK3WW6DxbpPlii+2CR/PpbOeAs\n5EB5AUVd66Z9ZveIX5duSbBo2Vx3aPfaS9yxPk2DD7tW4bjCfcSqFaJaZXNXF9ViCS0WuPCmZ0hW\n+khWHG/p47yKQ7LcR9ede7jitVdJl/toKhe4e+W9ZJwCGadI1inSVHa3Y3p4MfMndfqZDsWb2Jv5\nAO9kZvP5C04jvuAYEgsWsPLB7XSlW6mOETbYCDdO3A3rsH/W4SU9UcgUojR7L4CWnjjZQoxMIUq2\nGCPVN3bgN0EoHMpTOJSnq/NNPs6RvWPWbfrRYfuXRxZSilcoxavc+OLZONEqTkxxYlWcqFKOefvR\nKmXv83JEqUSVSsRdHik2EYtUiUqVmKi7jlSJoANj9NKR9zkl084nm9p5x1nMjuLp7CyuoLc6OINZ\nlRh7ncXsdRbDHa965x1kduwtcssvINeWJteWJjsrSSobJ92cIJ6YeEC8oAgyh38qsENVdwGIyJ+B\ni4G6O/xvPXslQoJ1Lz6AIqiKu/aWat8tVCSJE22i/PWnqcpIN+K8Ua/R3NvJ3EMdzM130JKMw6Yy\n+dsc3iufjJbL0L4ULX7Q7XWoQjs/HDh3+/rB71kzyjX2b4VLx6XcP/lEE12pVrrSObrSOd5N5bzt\nVt5NtdCVztEbH2yP+MYNg7npfY+Nc2h6DWZaDn2mMbQkMLT+XAW6myp0N1V4Z4TzNn3pZVbes4JM\nIUamGCNbiJIpxEiXogNLqi8yMCm8X6JVoakUo2kSbaj3U7vBOSYVIqIDzt9dKxHZgrCZOHHKmnYX\nUrjFjQh4GZ8SQh7gr0PH1AxqlIgQiQgSiSARQQRv7U6oI176a3pK9EVS9I/RCJIgHf58YGgWeQ9w\nWhAX0sSJVKJJ9vd/MPy5Ss0Z1/dFy0WyPW+Ty+8il99J7tAuEk7PwHGn5plTU7wrRWL0xNP0xlMc\nSjRx6ikf45GdPeSTGfKJrLtOZge3E1mcaCzwPuqTdejDz7cXRDCMd4BUrRfB0M/z2TL5bO3eWlsu\n28zye5aQLkYPexEknQgJx10nnQjJvggJJ0LKiY5ajVQPyhqF0cek4f7bHeDQiEfHOt1Py0UKSErG\nR8rJI6pjKp7YF4t8AThXVb/m7V8OnKqq1w9Ldy1wrbd7PPDGBC85B/BRGz2jMM0zn0bTC6Z5vHxY\nVdv8JAwyh78HWDBk/xjgiNmiVfUu4K7JXkxENqrqsrFTzhxM88yn0fSCaQ6SIFvhXgYWichHRCQB\nrAYmP7WQYRiGMSECy+GrallEvg08hdstc62qbgvqeoZhGMboBNoPX1XXA+vHTFgfJl0tFEJM88yn\n0fSCaQ6MwBptDcMwjKMLG0ljGIbRIITO4YvI50TkDRHZISI3jnA8KSIPe8dfEpFjp97K+uFD7w0i\n8pqIdIjI30XkyLCdIWMszUPSrRIRFZHQ9+jwo1lEvujd620i8uBU21hvfDzbC0XkORF5xXu+z58O\nO+uFiKwVkX0isrXGcRGRNd7v0SEiS+puhKqGZsFt/N0JfBRIAFuAE4al+SZwp7e9Gnh4uu0OWO9K\noMnbvi7Mev1q9tI1Ay8AG4Bl0233FNznRcArwCxvf+502z0Fmu8CrvO2TwDenG67J6n5DGAJsLXG\n8fOBv+GO4FwOvFRvG8KWwx8I16CqfUB/uIahXAzc520/CnxWwhbhaJAx9arqc6ra6+1uwB3vEGb8\n3GOAXwC/BsIfwtCf5muA21T1IICq7ptiG+uNH80KtHjbOUYYxxMmVPUF4MAoSS4G7leXDUCriNR1\nMoKwOfyRwjXMr5VGVctAHpg9JdbVHz96h3I1bg4hzIypWUQ+BSxQ1Sen0rAA8XOfjwOOE5F/icgG\nLxJtmPGj+WfAZSKyB7e33/XMbMb7fx83YYuHP1JOfXg3Iz9pwoJvLSJyGbAM+EygFgXPqJpFJAL8\nFrhqqgyaAvzc5xhutc6ZuKW4f4rISarqby7Jow8/mi8F7lXVm0VkBfAnT/PMmpVkkMB9V9hy+H7C\nNQykEZEYblFwtGLU0Yyv8BQichbwY+AiVQ37/GxjaW4GTgKeF5E3ces620PecOv3uX5cVR1V/R9u\nzKlFU2RfEPjRfDWwDkBV/40bZ2x8kRDDha//+2QIm8P3E66hHbjS214FPKtei0gIGVOvV73xe1xn\nH/Z6XRhDs6rmVXWOqh6rqsfitltcpKobp8fcuuDnuX4Mt4EeEZmDW8Wza0qtrC9+NHfizTokIotx\nHf5+Zi7twBVeb53lQF5VR4pIPWFCVaWjNcI1iMjPgY2q2g7cg1v024Gbs189fRZPDp96bwKywCNe\n23SnqtaYeuvox6fmGYVPzU8B54jIa7hRd7+vqvWZqGAa8Kn5e8AfROS7uFUbV4U484aIPIRbJTfH\na5f4KbhzK6rqnbjtFOcDO4Be4Ct1tyHEv59hGIYxDsJWpWMYhmFMEHP4hmEYDYI5fMMwjAbBHL5h\nGEaDYA7fMAyjQTCHbzQUIvK8iJw77LPviMjto5zTHbxlhhE85vCNRuMhjhybsdr73DBmNObwjUbj\nUeBCEUkCePMlfAjY7M0nsElEXhWRIyJ0isiZIvLkkP3fichV3vZSEfmHiPxXRJ6qd5RDw6gH5vCN\nhsIbnfofoD/a5GrgYaAAXKKqS3BDGNzsN6y2iMSBW4FVqroUWAv8st62G8ZkCVVoBcOoE/3VOo97\n66/iRir8lYicAVRxw9LOA/b6+L7jcQO6Pe29I6JAXWOgGEY9MIdvNCKPAb/xppBLq+omr2qmDViq\nqo4XiTM17Lwyh5eK+48LsE1VVwRrtmFMDqvSMRoOVe0GnseteulvrM0B+zxnvxIYaW7gt4ATxJ03\nOYcXyRE3VHGbF7MdEYmLyIlBajCMiWA5fKNReQj4C4M9dh4AnhCRjcBm4PXhJ6jqbhFZB3QA23Hn\nmEVV+0RkFbDGexHEgFuAbYGrMIxxYNEyDcMwGgSr0jEMw2gQzOEbhmE0CObwDcMwGgRz+IZhGA2C\nOXzDMIwGwRy+YRhGg2AO3zAMo0Ewh28YhtEg/B9oAQ3n+JEV5QAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fd659042400>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "figure, ax = pl.subplots(1, 1)\n",
    "ax.hist(samples, bins=500, normed=1, histtype='step', lw=2)\n",
    "wc_priors['mu'].plot(ax=ax, lw=3)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# How distributions work with the NEMS model/module system\n",
    "\n",
    "We have discussed how to specify distributions for a coefficient. However, how might we use these in the actual NEMS fit process? For Bayesian fitting algorithms, the priors would become part of the actual fit itself. For standard descent algorithms (e.g., those provided by `scipy.optimize`), we can't use the distributions directly. \n",
    "\n",
    "Let's build a rudimentary linear STRF model. Note that the FIR module knows that it will be getting two input channels because the `model.get_priors` method asks each module to transform the input data appropriately. The channel weighting collapses the 18 channels in the input data to two channels, which are then fed to the FIR filter so it can generate an appropriate set of priors.\n",
    "\n",
    "This approach (i.e., allowing the modules to inspect the data before generating a set of priors) allows the modules to specify some reasonable starting points for the fitting process."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[{'mu': Beta(α=array(2), β=array(2)), 'sd': HalfNormal(σ=array(2))},\n",
       " {'coefficients': Normal(μ=array(2, 10), σ=array(2, 10))}]"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model = Model()\n",
    "model.append(WeightChannelsGaussian(2))\n",
    "model.append(FIR(n_taps=10))\n",
    "\n",
    "priors = model.get_priors(signals)\n",
    "priors"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Here, `model.get_priors` returns a list of dictionaries. Each entry in the list matches the corresponding module. The dictionary is the set of coefficients for that module. If a module has no coefficients, then the dictionary will be empty.\n",
    "\n",
    "We can use this to generate some reasonable constraints for our model fitting system. For example, let's compute the starting value of `phi` for fitting (the actual coefficient values)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Initial value for FIR coefficients: \n",
      " [[ 0.1  0.1  0.1  0.1  0.1  0.1  0.1  0.1  0.1  0.1]\n",
      " [ 0.1  0.1  0.1  0.1  0.1  0.1  0.1  0.1  0.1  0.1]]\n"
     ]
    }
   ],
   "source": [
    "phi = []\n",
    "for module_priors in priors:\n",
    "    module_phi = {}\n",
    "    for name, prior in module_priors.items():\n",
    "        module_phi[name] = prior.mean()\n",
    "    phi.append(module_phi)\n",
    "    \n",
    "# Starting values for the FIR coefficients\n",
    "print('Initial value for FIR coefficients: \\n', phi[1]['coefficients'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now, let's get the lower bound for each coefficient using the percentile."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Lower bound for weight channels mu: \n",
      " [ 0.  0.]\n",
      "Lower bound for FIR coefficients: \n",
      " [[-inf -inf -inf -inf -inf -inf -inf -inf -inf -inf]\n",
      " [-inf -inf -inf -inf -inf -inf -inf -inf -inf -inf]]\n"
     ]
    }
   ],
   "source": [
    "phi_lower = []\n",
    "for module_priors in priors:\n",
    "    module_phi = {}\n",
    "    for name, prior in module_priors.items():\n",
    "        module_phi[name] = prior.percentile(0)\n",
    "    phi_lower.append(module_phi)\n",
    "    \n",
    "print('Lower bound for weight channels mu: \\n', phi_lower[0]['mu'])\n",
    "print('Lower bound for FIR coefficients: \\n', phi_lower[1]['coefficients'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This is a bit broad. It may be better to use a percentile of 0.001 instead."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Lower bound for weight channels mu: \n",
      " [ 0.1  0. ]\n",
      "Lower bound for FIR coefficients: \n",
      " [[-2.99 -2.99 -2.99 -2.99 -2.99 -2.99 -2.99 -2.99 -2.99 -2.99]\n",
      " [-2.99 -2.99 -2.99 -2.99 -2.99 -2.99 -2.99 -2.99 -2.99 -2.99]]\n"
     ]
    }
   ],
   "source": [
    "phi_lower = []\n",
    "for module_priors in priors:\n",
    "    module_phi = {}\n",
    "    for name, prior in module_priors.items():\n",
    "        module_phi[name] = prior.percentile(0.001)\n",
    "    phi_lower.append(module_phi)\n",
    "    \n",
    "print('Lower bound for weight channels mu: \\n', phi_lower[0]['mu'])\n",
    "print('Lower bound for FIR coefficients: \\n', phi_lower[1]['coefficients'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This information can be used by fit algorithms that support bounds (such as `scipy.optimize`). The code above has been wrapped up into some utility functions available via `nems.fitters.util`. Some basic fitting algorithms have already been implemented in `nems.fitters.scipy`."
   ]
  }
 ],
 "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.6.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}