3. Landmark Detection and Tracking.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Project 3:  Implement SLAM \n",
    "\n",
    "---\n",
    "\n",
    "## Project Overview\n",
    "\n",
    "In this project, you'll implement SLAM for robot that moves and senses in a 2 dimensional, grid world!\n",
    "\n",
    "SLAM gives us a way to both localize a robot and build up a map of its environment as a robot moves and senses in real-time. This is an active area of research in the fields of robotics and autonomous systems. Since this localization and map-building relies on the visual sensing of landmarks, this is a computer vision problem. \n",
    "\n",
    "Using what you've learned about robot motion, representations of uncertainty in motion and sensing, and localization techniques, you will be tasked with defining a function, `slam`, which takes in six parameters as input and returns the vector `mu`. \n",
    "> `mu` contains the (x,y) coordinate locations of the robot as it moves, and the positions of landmarks that it senses in the world\n",
    "\n",
    "You can implement helper functions as you see fit, but your function must return `mu`. The vector, `mu`, should have (x, y) coordinates interlaced, for example, if there were 2 poses and 2 landmarks, `mu` will look like the following, where `P` is the robot position and `L` the landmark position:\n",
    "```\n",
    "mu =  matrix([[Px0],\n",
    "              [Py0],\n",
    "              [Px1],\n",
    "              [Py1],\n",
    "              [Lx0],\n",
    "              [Ly0],\n",
    "              [Lx1],\n",
    "              [Ly1]])\n",
    "```\n",
    "\n",
    "You can see that `mu` holds the poses first `(x0, y0), (x1, y1), ...,` then the landmark locations at the end of the matrix; we consider a `nx1` matrix to be a vector.\n",
    "\n",
    "## Generating an environment\n",
    "\n",
    "In a real SLAM problem, you may be given a map that contains information about landmark locations, and in this example, we will make our own data using the `make_data` function, which generates a world grid with landmarks in it and then generates data by placing a robot in that world and moving and sensing over some numer of time steps. The `make_data` function relies on a correct implementation of robot move/sense functions, which, at this point, should be complete and in the `robot_class.py` file. The data is collected as an instantiated robot moves and senses in a world. Your SLAM function will take in this data as input. So, let's first create this data and explore how it represents the movement and sensor measurements that our robot takes.\n",
    "\n",
    "---"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Create the world\n",
    "\n",
    "Use the code below to generate a world of a specified size with randomly generated landmark locations. You can change these parameters and see how your implementation of SLAM responds! \n",
    "\n",
    "`data` holds the sensors measurements and motion of your robot over time. It stores the measurements as `data[i][0]` and the motion as `data[i][1]`.\n",
    "\n",
    "#### Helper functions\n",
    "\n",
    "You will be working with the `robot` class that may look familiar from the first notebook, \n",
    "\n",
    "In fact, in the `helpers.py` file, you can read the details of how data is made with the `make_data` function. It should look very similar to the robot move/sense cycle you've seen in the first notebook."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      " \n",
      "Landmarks:  [[34, 84], [19, 52], [81, 26], [55, 48], [96, 20]]\n",
      "Robot: [x=24.39911 y=81.05220]\n"
     ]
    }
   ],
   "source": [
    "import numpy as np\n",
    "from helpers import make_data\n",
    "\n",
    "# your implementation of slam should work with the following inputs\n",
    "# feel free to change these input values and see how it responds!\n",
    "\n",
    "# world parameters\n",
    "num_landmarks      = 5        # number of landmarks\n",
    "N                  = 20       # time steps\n",
    "world_size         = 100.0    # size of world (square)\n",
    "\n",
    "# robot parameters\n",
    "measurement_range  = 50.0     # range at which we can sense landmarks\n",
    "motion_noise       = 2.0      # noise in robot motion\n",
    "measurement_noise  = 2.0      # noise in the measurements\n",
    "distance           = 20.0     # distance by which robot (intends to) move each iteratation \n",
    "\n",
    "\n",
    "# make_data instantiates a robot, AND generates random landmarks for a given world size and number of landmarks\n",
    "data = make_data(N, num_landmarks, world_size, measurement_range, motion_noise, measurement_noise, distance)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### A note on `make_data`\n",
    "\n",
    "The function above, `make_data`, takes in so many world and robot motion/sensor parameters because it is responsible for:\n",
    "1. Instantiating a robot (using the robot class)\n",
    "2. Creating a grid world with landmarks in it\n",
    "\n",
    "**This function also prints out the true location of landmarks and the *final* robot location, which you should refer back to when you test your implementation of SLAM.**\n",
    "\n",
    "The `data` this returns is an array that holds information about **robot sensor measurements** and **robot motion** `(dx, dy)` that is collected over a number of time steps, `N`. You will have to use *only* these readings about motion and measurements to track a robot over time and find the determine the location of the landmarks using SLAM. We only print out the true landmark locations for comparison, later.\n",
    "\n",
    "\n",
    "In `data` the measurement and motion data can be accessed from the first and second index in the columns of the data array. See the following code for an example, where `i` is the time step:\n",
    "```\n",
    "measurement = data[i][0]\n",
    "motion = data[i][1]\n",
    "```\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Example measurements: \n",
      " [[0, -14.79420387999314, 32.01170947997674], [1, -32.821820115075965, 2.2298155261087915], [2, 29.423424704604315, -22.330698649091364], [3, 4.886080294200587, -1.6954373289649842], [4, 45.34924418067613, -28.0654896536633]]\n",
      "\n",
      "\n",
      "Example motion: \n",
      " [-8.463450506808663, -18.120982465606005]\n"
     ]
    }
   ],
   "source": [
    "# print out some stats about the data\n",
    "time_step = 0\n",
    "\n",
    "print('Example measurements: \\n', data[time_step][0])\n",
    "print('\\n')\n",
    "print('Example motion: \\n', data[time_step][1])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Try changing the value of `time_step`, you should see that the list of measurements varies based on what in the world the robot sees after it moves. As you know from the first notebook, the robot can only sense so far and with a certain amount of accuracy in the measure of distance between its location and the location of landmarks. The motion of the robot always is a vector with two values: one for x and one for y displacement. This structure will be useful to keep in mind as you traverse this data in your implementation of slam."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Initialize Constraints\n",
    "\n",
    "One of the most challenging tasks here will be to create and modify the constraint matrix and vector: omega and xi. In the second notebook, you saw an example of how omega and xi could hold all the values the define the relationships between robot poses `xi` and landmark positions `Li` in a 1D world, as seen below, where omega is the blue matrix and xi is the pink vector.\n",
    "\n",
    "<img src='images/motion_constraint.png' width=50% height=50% />\n",
    "\n",
    "\n",
    "In *this* project, you are tasked with implementing constraints for a 2D world. We are referring to robot poses as `Px, Py` and landmark positions as `Lx, Ly`, and one way to approach this challenge is to add *both* x and y locations in the constraint matrices.\n",
    "\n",
    "<img src='images/constraints2D.png' width=50% height=50% />\n",
    "\n",
    "You may also choose to create two of each omega and xi (one for x and one for y positions)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### TODO: Write a function that initializes omega and xi\n",
    "\n",
    "Complete the function `initialize_constraints` so that it returns `omega` and `xi` constraints for the starting position of the robot. Any values that we do not yet know should be initialized with the value `0`. You may assume that our robot starts out in exactly the middle of the world with 100% confidence (no motion or measurement noise at this point). The inputs `N` time steps, `num_landmarks`, and `world_size` should give you all the information you need to construct intial constraints of the correct size and starting values.\n",
    "\n",
    "*Depending on your approach you may choose to return one omega and one xi that hold all (x,y) positions *or* two of each (one for x values and one for y); choose whichever makes most sense to you!*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "def initialize_constraints(N, num_landmarks, world_size):\n",
    "    ''' This function takes in a number of time steps N, number of landmarks, and a world_size,\n",
    "        and returns initialized constraint matrices, omega and xi.'''\n",
    "    \n",
    "    ## Recommended: Define and store the size (rows/cols) of the constraint matrix in a variable\n",
    "    \n",
    "    ## TODO: Define the constraint matrix, Omega, with two initial \"strength\" values\n",
    "    ## for the initial x, y location of our robot\n",
    "    total_size = N+num_landmarks\n",
    "    omega = np.zeros((2, total_size, total_size))\n",
    "    \n",
    "    ## TODO: Define the constraint *vector*, xi\n",
    "    ## you can assume that the robot starts out in the middle of the world with 100% confidence\n",
    "    xi = np.zeros((2, total_size, 1))\n",
    "    \n",
    "    for i in range(2):\n",
    "        omega[i, 0, 0] = 1\n",
    "        xi[i, 0, 0] = world_size/2.0\n",
    "    \n",
    "    return omega, xi\n",
    "    "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Test as you go\n",
    "\n",
    "It's good practice to test out your code, as you go. Since `slam` relies on creating and updating constraint matrices, `omega` and `xi` to account for robot sensor measurements and motion, let's check that they initialize as expected for any given parameters.\n",
    "\n",
    "Below, you'll find some test code that allows you to visualize the results of your function `initialize_constraints`. We are using the [seaborn](https://seaborn.pydata.org/) library for visualization.\n",
    "\n",
    "**Please change the test values of N, landmarks, and world_size and see the results**. Be careful not to use these values as input into your final smal function.\n",
    "\n",
    "This code assumes that you have created one of each constraint: `omega` and `xi`, but you can change and add to this code, accordingly. The constraints should vary in size with the number of time steps and landmarks as these values affect the number of poses a robot will take `(Px0,Py0,...Pxn,Pyn)` and landmark locations `(Lx0,Ly0,...Lxn,Lyn)` whose relationships should be tracked in the constraint matrices. Recall that `omega` holds the weights of each variable and `xi` holds the value of the sum of these variables, as seen in Notebook 2. You'll need the `world_size` to determine the starting pose of the robot in the world and fill in the initial values for `xi`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "# import data viz resources\n",
    "import matplotlib.pyplot as plt\n",
    "from pandas import DataFrame\n",
    "import seaborn as sns\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "# define a small N and world_size (small for ease of visualization)\n",
    "N_test = 5\n",
    "num_landmarks_test = 2\n",
    "small_world = 10\n",
    "\n",
    "# initialize the constraints\n",
    "initial_omega, initial_xi = initialize_constraints(N_test, num_landmarks_test, small_world)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABFsAAAHiCAYAAAA3XcMSAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzt3X/UpnddH/j3ZzKJoYLYlc4kS0agm7iWH2dFs7hbuvJDgYFwIBRtwVqrFaeHFavHHlsoXVzjAWt7Wn8cQBgC/lqVtbrYSLLBlsKCVJQoyC+lO7A0mUKeYLVYimnI5Lt/5Jmc8cpkfiTXfd/P9blfL859+tzPc0+e7ztOruvd7/W57rvGGAEAAABgHvs2vQAAAACATmy2AAAAAMzIZgsAAADAjGy2AAAAAMzIZgsAAADAjGy2AAAAAMzIZgsAAADAjGy2wIZU1YOr6pNV9c2nfO8hVXVzVX3jJtcGANBVVf18Vb1p8r0nVdV/rKpLN7UuoJcaY2x6DbC1qurpSX4+yaPHGJ+pqp9McnCM8Vc3vDQAgJaq6suSfCTJ3xxj/KuqujjJB5O8aozx0xtdHNCGzRbYsKr66SRflOT1SX4lyWPHGJ/e6KIAABqrqm9K8k+SPDbJP0ryVWOMZ252VUAnNltgw6rqzyf5aJILk3z/GOOnNrwkAID2quqXk1yU5IlJHj/GuHnDSwIasdkCe0BV/eskfznJpWOMz256PQAA3VXVwSQfT/LyMcaPb3o9QC/eIBc2rKq+Jckjk/zrJD+y2dUAAGyHMcZOkj/M3e/fAjCr/ZteAGyzqjqQ5EeT/LUkf5DkI1X1C2OMd212ZQAAANxfJltgs16d5FfHGO/YfVPcv5/kDVX1RRteFwAAAPeTzRbYkKq6OslfSfL9J783xrg2yfEkr9jUugAAAHhgvEEuAAAAwIxMtgAAAADMyGYLAOwRVfWmqrqtqj58Hz+vqvqJqjpWVR+sqq9e9xoBALpZRQez2QIAe8dPJzl8hp8/M8kVu48jSX5yDWsCAOjupzNzB7PZAgB7xO7Hvv/RGV7y3CQ/O+723iRfWlWXrmd1AAA9raKD2WwBgOV4eJJbTnl+fPd7AACsznl3sP0rXc7dfNwRAHOqVf+CBz3+JSs5d93+gdf8ndw9enrS0THG0fP4R5wuu/Ms98XfDQDmpIP9WWdc6zo2W/Kgx79kHb9mY/70/a9Oktx+54YXsmIX7/5tkbMHOXvZtpxLtXtSP58T+9TxJIdOeX5Zkk89oEXRmg7Ww7Yd4+XsQc5edLDz72BuIwKAqdq3mscDd12Sb919R/z/KclnxxifnuMfDACwcY062ML3pwCgj6r6xSRPTvKwqjqe5AeSXJgkY4zXJbkhybOSHEvy+STfvpmVAgD0sYoOZrMFAKZq5bckn9YY44Vn+flI8l1rWg4AwHo16mBuIwIAAACYkckWAJia595eAADOR6MOZrMFAKY2NMIKALDVGnWwPttGAAAAAHuAyRYAmGo0wgoAsBiNOlifJAAAAAB7gMkWAJhqdL8wAMBiNOpgNlsAYKrRCCsAwGI06mB9kgAAAADsASZbAGCq0QgrAMBiNOpgJlsAAAAAZmSyBQCmGt0vDACwGI06mM0WAJhqNMIKALAYjTpYn20jAAAAgD3AZAsATDUaYQUAWIxGHaxPEgAAAIA9wGQLAEw1ul8YAGAxGnUwky0AAAAAMzLZAgBTje4XBgBYjEYdzGYLAEw1OtEDACxGow7WJwkAAADAHtBqs+V1P/A38u/f/sO56V/8w00vZeXe8+535TlXPSPPPvy0vPENRze9nJWRsxc5e2mdc1+t5gFNbUsHa33cO4WcvcjZS/ucjTpYq82Wn/u19+a53/WaTS9j5U6cOJFXvfKavPZ11+Yt112fG294az5+7NimlzU7OXuRs5dtyQmcm23oYNty3JOzFzl72ZacXbTabHnP7348f/TZz296GSv34Q99MIcOPSKXHTqUCy+6KIefdVXe+Y63b3pZs5OzFzl7aZ+z9q3mAU1tQwdrf9zbJWcvcvayFTkbdbCz/taq+sqq+gdV9RNV9eO7X/+ldSyO07ttZyeXXHrJPc8PHDyYnZ2dDa5oNeTsRc5e2uesWs0DzoMOtre0P+7tkrMXOXvZipyNOtgZN1uq6h8keXOSSvLbSd63+/UvVtVLV788Tmdk3Ot71bDEy9mLnL1sS07YFB1s79mW456cvcjZy7bk7OJsH/38HUkeM8b4wqnfrKp/nuQjSf7x6f5QVR1JciRJXv/618+wTE518OAlufXTt97z/LadnRw4cGCDK1oNOXuRs5f2Od3yw+bpYHtM++PeLjl7kbOXrcjZqIOdLcldSf7b03z/0t2fndYY4+gY48oxxpVHjhx5IOvjNB7z2Mfl5ps/mePHb8kX7rgjN95wfZ70lKduelmzk7MXOXvZlpywQTrYHrMtxz05e5Gzl23J2cXZJlu+N8nbq+r/TXLL7ve+PMnlSV6yyoXdHz/zw9+W/+VrrsjDvvTBOXbjD+WHXndDfuZXf3PTy5rd/v3787KXvyIvPvKi3HXXiVz9vOfn8suv2PSyZidnL3L20j6nkVw2TwfbY9of93bJ2YucvWxFzkYdrMa4931ff+YFVfuSPCHJw3P3vcLHk7xvjHHiHH/HeNDj91wnmNWfvv/VSZLb79zwQlbs4t2tOTl7kLOXLcu58rPwg57+T898cryf/vTXv79Pg2DldLCz08F6kbMXOXvRwc7f2SZbMsa4K8l717AWAAB26WAAsFxn3WwBgK3TaIQVAGAxGnWwPm/1CwAAALAHmGwBgKlGHzsIALAYjTqYzRYAmGo0wgoAsBiNOlifbSMAAACAPcBkCwBMNRphBQBYjEYdrE8SAAAAgD3AZAsATDW6XxgAYDEadTCTLQAAAAAzMtkCAFON7hcGAFiMRh3MZgsATDU60QMALEajDtYnCQAAAMAeYLIFAKYavTkbAMBiNOpgJlsAAAAAZmSyBQCmGt0vDACwGI06mM0WAJhqNMIKALAYjTpYn20jAAAAgD3AZAsATDUaYQUAWIxGHaxPEgAAAIA9wGQLAEw1ul8YAGAxGnUwmy0AMFGNTvQAAEvRqYO5jQgAAABgRiZbAGCi01UVAICl6NTBTLYAAAAAzMhkCwBM9bmoAgCwHI062Fo2W/70/a9ex6/ZuIu3ZOtKzl7k7GVbcgLnRgfrRc5e5OxlW3Jy7vyVAICJTvcLAwAsRacOtpbNltvvXMdv2ZyTu5hy9iBnL3L2sq6rRp1O9Gy3bTkmyNmDnL3I2YsOdv68QS4AAADAjNxGBAATna6qAAAsRacOZrIFAAAAYEYmWwBgotNVFQCApejUwWy2AMBUn/M8AMByNOpgbiMCAAAAmJHJFgCY6DTCCgCwFJ06mMkWAAAAgBmZbAGAiU5XVQAAlqJTB7PZAgATnU70AABL0amDuY0IAAAAYEYmWwBgotNVFQCApejUwUy2AAAAAMzIZAsATPW5qAIAsByNOpjJFgAAAIAZmWwBgIlO9wsDACxFpw5mswUAJjqd6AEAlqJTB3MbEQAAAMCMTLYAwESnqyoAAEvRqYOZbAGAPaSqDlfVx6rqWFW99DQ///KqekdVvb+qPlhVz9rEOgEAulhF/7LZAgBTtaLH2X5t1QVJXpPkmUkeneSFVfXoycv+UZJfGmM8PskLkrz2fucEANhLNtDBVtW/3EYEABMbHGF9QpJjY4xP7K7jzUmem+Sjp7xmJPmS3a8fmuRTa10hAMCKbKiDraR/2WwBgL3j4UluOeX58SRfO3nN/57k16vqu5N8cZJvWM/SAABaWkn/chsRAExU1aoeR6rqplMeR6a/+jTLGZPnL0zy02OMy5I8K8nPVZXzOQCweBvqYCvpX+3K2Xve/a4856pn5NmHn5Y3vuHoppezMnL2ImcvcnJfxhhHxxhXnvKY/os7nuTQKc8vy73HVL8jyS/t/vN+M8nFSR62qjXDudiW44GcvcjZi5ycyVk62Er6V6vNlhMnTuRVr7wmr33dtXnLddfnxhvemo8fO7bpZc1Ozl7k7EXOHlZ1VeUcvC/JFVX1qKq6KHe/Adt1k9fcnOTrd9f5l3L3yf4zM8aH89L9eHCSnL3I2YucfWyog62kf7XabPnwhz6YQ4cekcsOHcqFF12Uw8+6Ku98x9s3vazZydmLnL3I2cOmNlvGGHcmeUmStyX5/dz9rvcfqaprquo5uy/7e0m+s6p+L8kvJvm2McZ01BXWpvvx4CQ5e5GzFzn72EQHW1X/arXZctvOTi659JJ7nh84eDA7OzsbXNFqyNmLnL3IyQM1xrhhjPEVY4z/bozxyt3vvWKMcd3u1x8dYzxxjPE/jDG+aozx65tdMdtuW44HcvYiZy9y8kCton/d782Wqvr2+/tnV2Xc6z1sNvrxnSsjZy9y9iJnE7WiB8xgr3Ww9seDXXL2ImcvcjbSqIM9kMmWH7yvH5z6Tr9Hj67vTXsOHrwkt3761nue37azkwMHDqzt96+LnL3I2YucwBrsqQ62LccDOXuRsxc52YvOuNlSVR+8j8eHkhy8rz936jv9Hjky/VTL1XnMYx+Xm2/+ZI4fvyVfuOOO3HjD9XnSU566tt+/LnL2ImcvcvawqfdsgZOW1MG6Hw9OkrMXOXuRs49OHWz/WX5+MMkzkvzx5PuV5N+uZEUPwP79+/Oyl78iLz7yotx114lc/bzn5/LLr9j0smYnZy9y9iInMJPFdLBtOR7I2YucvcjJXlRnegPdqnpjkp8aY/zGaX72C2OMbz6H3zFuv/MBrHABLt7dspKzBzl7kbOX3Zwrvzxx2f/6qyv5dJ/jr73aeAvnRAc7N1t27JOzCTl72bKcOth5OONkyxjjO87ws3M5yQPA4rjlh03TwQDYRp06WKuPfgYAAADYtLO9ZwsAbJ8+F1UAAJajUQcz2QIAAAAwI5MtADDR6X5hAICl6NTBbLYAwESnEz0AwFJ06mBuIwIAAACYkckWAJjodFUFAGApOnUwky0AAAAAMzLZAgATna6qAAAsRacOZrMFAKb6nOcBAJajUQdzGxEAAADAjEy2AMBEpxFWAICl6NTBTLYAAAAAzMhkCwBMdLqqAgCwFJ06mMkWAAAAgBmZbAGAiUYXVQAAFqNTB7PZAgATnUZYAQCWolMHcxsRAAAAwIxMtgDARKOLKgAAi9Gpg5lsAQAAAJiRyRYAmOh0vzAAwFJ06mA2WwBgotF5HgBgMTp1MLcRAQAAAMzIZAsATOzb1+iyCgDAQnTqYCZbAAAAAGa0lsmWi7dkfkbOXuTsRU7OR6f7hdlu23JMkLMXOXuRk/PRqYP5KwEAE53eCR8AYCk6dbC1bLbcfuc6fsvmnNzFlLMHOXuRsxdXjeD8bMsxQc4e5OxFzl50sPPnXxkATDS6qAIAsBidOpg3yAUAAACYkckWAJjodL8wAMBSdOpgJlsAAAAAZmSyBQAmOl1VAQBYik4dzGYLAEw0Os8DACxGpw7mNiIAAACAGZlsAYCJTiOsAABL0amDmWwBAAAAmJHJFgCYaHRRBQBgMTp1MJstADDRaYQVAGApOnUwtxEBAAAAzMhkCwBMNLqoAgCwGJ06mMkWAAAAgBmZbAGAiU73CwMALEWnDmazBQAmGp3nAQAWo1MHcxsRAAAAwIxMtgDARKcRVgCApejUwUy2AAAAAMzIZAsATDS6qAIAsBidOpjJFgAAAIAZmWwBgIlO9wsDACxFpw5mswUAJhqd5wEAFqNTB3MbEQAAAMCMTLYAwESnEVYAgKXo1MHaTba8593vynOuekaeffhpeeMbjm56OSsjZy9y9iInsG225XggZy9y9iIne02rzZYTJ07kVa+8Jq993bV5y3XX58Yb3pqPHzu26WXNTs5e5OxFzh6qVvOAjrofD06Ssxc5e5Gzj04d7KybLVX1lVX19VX14Mn3D69uWffPhz/0wRw69IhcduhQLrzoohx+1lV55zvevullzU7OXuTsRc4eqmolDzgfS+lg3Y8HJ8nZi5y9yNlHpw52xs2Wqvq7Sf5lku9O8uGqeu4pP37VKhd2f9y2s5NLLr3knucHDh7Mzs7OBle0GnL2ImcvcgJzWFIH25bjgZy9yNmLnOxFZ5ts+c4kXzPGuDrJk5P8b1X1Pbs/u8/toao6UlU3VdVNR4+u7z6ykXG6tazt96+LnL3I2YucPXS6qsJiLaaDdT8enCRnL3L2ImcfnTrY2T6N6IIxxueSZIzxyap6cpJfrqpH5Awn+jHG0SQnz/Dj9jvnWOrZHTx4SW799K33PL9tZycHDhxYzy9fIzl7kbMXOYGZLKaDbcvxQM5e5OxFTvais0223FpVX3Xyye5J/9lJHpbkcatc2P3xmMc+Ljff/MkcP35LvnDHHbnxhuvzpKc8ddPLmp2cvcjZi5w9dHpzNhZrMR2s+/HgJDl7kbMXOfvo1MHONtnyrUn+zDWRMcadSb61ql6/slXdT/v378/LXv6KvPjIi3LXXSdy9fOen8svv2LTy5qdnL3I2YucPXQbyWWRFtPBuh8PTpKzFzl7kbOPTh2sxrj3fV8zW9ttRJty8e6WlZw9yNmLnL3s5lz5WfjJP/ZvV3JyfOf3/uU+DYIl0MGakLMXOXvZspw62Hk422QLAGydRhdVAAAWo1MHO9t7tgAAAABwHky2AMBEp/uFAQCWolMHs9kCABONzvMAAIvRqYO5jQgAAABgRiZbAGBiX6fLKgAAC9Gpg5lsAQAAAJiRyRYAmGh0UQUAYDE6dTCTLQAAAAAzMtkCABOdPnYQAGApOnUwky0AMLGvVvM4F1V1uKo+VlXHquql9/Gav1ZVH62qj1TVL8yZHQBgUzbVwVbRv0y2AMAeUVUXJHlNkqclOZ7kfVV13Rjjo6e85ookL0vyxDHGH1fVgc2sFgBg+VbVv2y2AMDEBkdYn5Dk2BjjE7vreHOS5yb56Cmv+c4krxlj/HGSjDFuW/sqAQBWYEMdbCX9y21EALB3PDzJLac8P777vVN9RZKvqKr3VNV7q+rw2lYHANDPSvqXyRYAmFjVRZWqOpLkyCnfOjrGOHrqS07zx8bk+f4kVyR5cpLLkry7qh47xvhPc64VAGDdNtTBVtK/bLYAwESd9pz7wO2e1I+e4SXHkxw65fllST51mte8d4zxhST/X1V9LHef/N8351oBANZtQx1sJf3LbUQAsHe8L8kVVfWoqrooyQuSXDd5za8meUqSVNXDcvdY6yfWukoAgD5W0r9MtgDAxLl+TPPcxhh3VtVLkrwtyQVJ3jTG+EhVXZPkpjHGdbs/e3pVfTTJiSTfP8b4j5tZMQDAfDbRwVbVv2y2AMAeMsa4IckNk++94pSvR5Lv230AAPAAraJ/2WwBgIkNfvQzAMDW6tTBbLYAwESj8zwAwGJ06mDeIBcAAABgRiZbAGBiX6fLKgAAC9Gpg5lsAQAAAJiRyRYAmGh0UQUAYDE6dTCTLQAAAAAzWstky8VbMj8jZy9y9iIn56PTxw6y3bblmCBnL3L2Iifno1MHW8tfidvvXMdv2ZyT/2HJ2YOcvcjZy7qKTKPzPFtuW44JcvYgZy9y9qKDnT+3EQEAAADMyLATAEx0+thBAICl6NTBTLYAAAAAzMhkCwBM9LmmAgCwHJ06mM0WAJjo9E74AABL0amDuY0IAAAAYEYmWwBgYl+fiyoAAIvRqYOZbAEAAACYkckWAJjodL8wAMBSdOpgNlsAYKLReR4AYDE6dTC3EQEAAADMyGQLAEx0GmEFAFiKTh3MZAsAAADAjEy2AMBEp48dBABYik4dzGQLAAAAwIxMtgDARKf7hQEAlqJTB7PZAgATfU7zAADL0amDuY0IAAAAYEYmWwBgYl+jEVYAgKXo1MFMtgAAAADMyGQLAEw0uqgCALAYnTqYzRYAmOj0TvgAAEvRqYO5jQgAAABgRiZbAGCi0UUVAIDF6NTBTLYAAAAAzKjdZst73v2uPOeqZ+TZh5+WN77h6KaXszJy9iJnL3Iu376qlTygq87Hg1PJ2YucvcjZQ6cO1mqz5cSJE3nVK6/Ja193bd5y3fW58Ya35uPHjm16WbOTsxc5e5Gzh6rVPKCj7seDk+TsRc5e5OyjUwc762ZLVT2hqv7H3a8fXVXfV1XPWv3Szt+HP/TBHDr0iFx26FAuvOiiHH7WVXnnO96+6WXNTs5e5OxFTmAuS+lg23I8kLMXOXuRk73ojJstVfUDSX4iyU9W1Q8neXWSByd5aVW9fA3rOy+37ezkkksvuef5gYMHs7Ozs8EVrYacvcjZi5w9VNVKHnCultTBuh8PTpKzFzl7kbOPTh3sbJMt35jkiUm+Lsl3Jbl6jHFNkmck+ev39Yeq6khV3VRVNx09ur77yEbG6daytt+/LnL2ImcvcgIzWUwH25bjgZy9yNmLnOxFZ/vo5zvHGCeSfL6qPj7G+JMkGWP8aVXddV9/aIxxNMnJM/y4/c55Fns2Bw9ekls/fes9z2/b2cmBAwfW88vXSM5e5OxFzh5avaEZS7WYDtb9eHCSnL3I2YucfXTqYGfLckdV/bndr7/m5Der6qFJ7vNEvymPeezjcvPNn8zx47fkC3fckRtvuD5PespTN72s2cnZi5y9yAnMZDEdbFuOB3L2ImcvcrIXnW2y5evGGP81ScYYp57YL0zyt1a2qvtp//79ednLX5EXH3lR7rrrRK5+3vNz+eVXbHpZs5OzFzl7kbMHI7nsAYvpYN2PByfJ2YucvcjZR6cOVmPc+76vma3tNqJNuXh3y0rOHuTsRc5ednOu/Cz8vf/yD1Zycvyx535lnwbBEuhgTcjZi5y9bFlOHew8dLolCgAAAGDjznYbEQBsnX3mTwAA1q5TBzPZAgAAADAjky0AMNHpzdkAAJaiUwez2QIAE51GWAEAlqJTB3MbEQAAAMCMTLYAwESjCVYAgMXo1MFMtgAAAADMyGQLAEzs63RZBQBgITp1MJstADBh7BMAYP06dbBOWQAAAAA2zmQLAEw0mmAFAFiMTh3MZAsAAADAjEy2AMBEpzdnAwBYik4dzGQLAAAAwIxMtgDARKOLKgAAi9Gpg9lsAYCJfY1O9AAAS9Gpg7mNCAAAAGBGJlsAYKLTm7MBACxFpw5msgUAAABgRiZbAGCi0UUVAIDF6NTBbLYAwESnN2cDAFiKTh3MbUQAAAAAMzLZAgATlUaXVQAAFqJTBzPZAgAAADAjky0AMNHpfmEAgKXo1MHWstly8ZZs6cjZi5y9yMn56HSiZ7ttyzFBzl7k7EVOzkenDraWvxK337mO37I5J//DkrMHOXuRsxdFBs7PthwT5OxBzl7k7EUHO3/+lQHARFWjyyoAAAvRqYN5g1wAAACAGZlsAYCJTvcLAwAsRacOZrIFAAAAYEYmWwBgotHtwgAAi9Gpg9lsAYCJfZ3O9AAAC9Gpg7mNCAAAAGBGNlsAYGJfreZxLqrqcFV9rKqOVdVLz/C6b6yqUVVXzpUbAGCTNtXBVtG/bLYAwB5RVRckeU2SZyZ5dJIXVtWjT/O6hyT5u0l+a70rBADoZVX9y2YLAExUreZxDp6Q5NgY4xNjjDuSvDnJc0/zuh9K8k+S3D5baACADdtQB1tJ/7LZAgAT+1IreZyDhye55ZTnx3e/d4+qenySQ2OMt86XGABg8zbUwVbSv2y2AMCaVNWRqrrplMeR6UtO88fGKX9+X5IfTfL3VrlOAIBOztLBVtK/fPQzAEys6lMHxxhHkxw9w0uOJzl0yvPLknzqlOcPSfLYJO+suxd5SZLrquo5Y4ybZl4uAMBabaiDraR/mWwBgL3jfUmuqKpHVdVFSV6Q5LqTPxxjfHaM8bAxxiPHGI9M8t4kNloAAO6/lfQvky0AMHGuH9M8tzHGnVX1kiRvS3JBkjeNMT5SVdckuWmMcd2Z/wkAAMu1iQ62qv5lswUAJvataob1HIwxbkhyw+R7r7iP1z55HWsCAFiHTXWwVfQvtxEBAAAAzMhkCwBMbHCwBQBga3XqYCZbAAAAAGZksgUAJjb5ni0AANuqUwcz2QIAAAAwI5MtADDR6KIKAMBidOpgNlsAYMLYJwDA+nXqYJ2yAAAAAGycyRYAmKhOM6wAAAvRqYO1m2x5z7vfledc9Yw8+/DT8sY3HN30clZGzl7k7EVOYNtsy/FAzl7k7EVO9ppWmy0nTpzIq155TV77umvzluuuz403vDUfP3Zs08uanZy9yNmLnD3Uih7QUffjwUly9iJnL3L20amDnfdmS1X97CoWMocPf+iDOXToEbns0KFceNFFOfysq/LOd7x908uanZy9yNmLnD3sq1rJA+4v/Wvz5OxFzl7k7KNTBzvjZktVXTd5/FqSv3ry+ZrWeM5u29nJJZdecs/zAwcPZmdnZ4MrWg05e5GzFzmBB0r/2pvk7EXOXuRkLzrbG+ReluSjSa5NMnL3BM6VSf7Zmf5QVR1JciRJXv/61+db//aRB77SczAyTreWtfzudZKzFzl7kbOHPklYqPvVv5LNdLDux4OT5OxFzl7k7KNTmrNttlyZ5HuSvDzJ948xPlBVfzrG+H/O9IfGGEeTnHy3nnH7nQ98oefi4MFLcuunb73n+W07Ozlw4MB6fvkaydmLnL3ICczgfvWvZDMdbFuOB3L2ImcvcrIXnfE2ojHGXWOMH03y7UleXlWvzh7+uOjHPPZxufnmT+b48VvyhTvuyI03XJ8nPeWpm17W7OTsRc5e5OyhajUPOBf6194kZy9y9iJnH5062DmduMcYx5N8U1VdleRPVruk+2///v152ctfkRcfeVHuuutErn7e83P55Vdselmzk7MXOXuRs4duI7ksk/61t8jZi5y9yNlHpw5WY9z7vq+Zre02ok25eHfLSs4e5OxFzl52c678LPyL7/8PKzk5vvDxD+/TIFgCHawJOXuRs5cty6mDnYc9O5IKAJtyxntsAQBYiU4drFMWAAAAgI0z2QIAE53uFwYAWIpOHcxkCwAAAMCMTLYAwESfayoAAMvRqYPZbAGAiU4jrAAAS9Gpg7mNCAAAAGBGJlsAYMKVCACA9evUwTplAQAAANg4ky0AMNFOeO+xAAASHElEQVTpfmEAgKXo1MFstgDARJ/TPADAcnTqYG4jAgAAAJiRyRYAmGg0wQoAsBidOpjJFgAAAIAZmWwBgIl9re4YBgBYhk4dzGYLAEx0GmEFAFiKTh3MbUQAAAAAMzLZAgAT1WiEFQBgKTp1MJMtAAAAADMy2QIAE53uFwYAWIpOHcxmCwBMdHonfACApejUwdxGBAAAADAjky0AMNFphBUAYCk6dTCTLQAAAAAzMtkCABOdrqoAACxFpw5msgUAAABgRmuZbLl4S+Zn5OxFzl7k5HxUo3fCZ7ttyzFBzl7k7EVOzkenDraWvxK337mO37I5J//DkrMHOXuRs5d1FZl9fc7zbLltOSbI2YOcvcjZiw52/txGBAAAADAjw04AMNFphBUAYCk6dTCTLQAAAAAzMtkCABOdPnYQAGApOnUwmy0AMNFphBUAYCk6dTC3EQEAAADMyGQLAEx0+thBAICl6NTBTLYAAAAAzMhkCwBMdLpfGABgKTp1MJstADDR6Z3wAQCWolMHcxsRAAAAwIxMtgDARKOLKgAAi9Gpg5lsAQAAAJiRyRYAmNjX6YZhAICF6NTBTLYAAAAAzMhkCwBM9LmmAgCwHJ06mM0WAJjqdKYHAFiKRh3MbUQAAAAAMzLZAgAT1emyCgDAQnTqYCZbAAAAAGZksgUAJhp96iAAwGJ06mA2WwBgotF5HgBgMTp1MLcRAQAAAMzIZAsATHW6rAIAsBSNOpjJFgAAAIAZtdtsec+735XnXPWMPPvw0/LGNxzd9HJWRs5e5OxFzuWrFf0Puup8PDiVnL3I2YucPXTqYK02W06cOJFXvfKavPZ11+Yt112fG294az5+7NimlzU7OXuRsxc5e6hazQM66n48OEnOXuTsRc4+OnWw89psqaq/UlXfV1VPX9WCHogPf+iDOXToEbns0KFceNFFOfysq/LOd7x908uanZy9yNmLnMAq7OUOti3HAzl7kbMXOdmLzrjZUlW/fcrX35nk1UkekuQHquqlK17bebttZyeXXHrJPc8PHDyYnZ2dDa5oNeTsRc5e5OyhVvSAc7WkDtb9eHCSnL3I2YucfXTqYGebbLnwlK+PJHnaGOMHkzw9yd9Y2arup5Fxr+9Vw7ltOXuRsxc5gZkspoNty/FAzl7k7EVO9qKzbbbsq6o/X1VflqTGGJ9JkjHGf0ly5339oao6UlU3VdVNR4+u7017Dh68JLd++tZ7nt+2s5MDBw6s7fevi5y9yNmLnE10uqzCUi2mg7U/HuySsxc5e5GzkUYd7GybLQ9N8jtJbkry31TVJUlSVQ/OGZY8xjg6xrhyjHHlkSNHZlvs2TzmsY/LzTd/MseP35Iv3HFHbrzh+jzpKU9d2+9fFzl7kbMXOYGZLKaDbcvxQM5e5OxFTvai/Wf64Rjjkffxo7uSPG/21TxA+/fvz8te/oq8+MiLctddJ3L1856fyy+/YtPLmp2cvcjZi5w9bPJjmqvqcJIfT3JBkmvHGP948vPvS/Ki3D3d8Jkkf3uM8e/XvlBWakkdrPvx4CQ5e5GzFzn72FQHW0X/qjHufd/XzMbt9zns2sPFu1tWcvYgZy9y9rKbc+Vn4Q/c/J9XcnL8qi9/yBnXXlUXJPl3SZ6W5HiS9yV54Rjjo6e85ilJfmuM8fmqenGSJ48x/voq1svi6WBNyNmLnL1sWc6WHWxV/eu8PvoZAFipJyQ5Nsb4xBjjjiRvTvLcU18wxnjHGOPzu0/fm+SyNa8RAKCTlfQvmy0AMLHB92Z7eJJbTnl+fPd79+U7kvzf5/aPBgDY2zbUwVbSv874ni0AwHyq6kju/hjfk46OMU79yJjT9YHTjtNW1bckuTLJk+ZbIQBAP2fpYCvpXzZbAGBqRXck757Uz/R5vMeTHDrl+WVJPjV9UVV9Q5KXJ3nSGOO/zrpIAIBN2UwHW0n/stkCABMb/DSi9yW5oqoeleQ/JHlBkm8+9QVV9fgkr09yeIxx2/qXCACwGhvqYCvpX96zBQD2iDHGnUlekuRtSX4/yS+NMT5SVddU1XN2X/ZPkzw4yb+oqg9U1XUbWi4AwOKtqn/56OcZbNnHfcnZhJy9bFnOlV/y+NDxz63k5Pi4yx68sZEZtpIO1oScvcjZy5bl1MHOg8kWAAAAgBl5zxYAmDB+AgCwfp06mM0WAJjqdKYHAFiKRh3MbUQAAAAAMzLZAgATG/zoZwCArdWpg5lsAQAAAJiRyRYAmKg+F1UAABajUwcz2QIAAAAwI5MtADDR6KIKAMBidOpgNlsAYKrTmR4AYCkadTC3EQEAAADMyGQLAEx0+thBAICl6NTBTLYAAAAAzMhkCwBMdPrYQQCApejUwWy2AMBEo/M8AMBidOpgbiMCAAAAmJHJFgCY6nRZBQBgKRp1MJMtAAAAADMy2QIAE50+dhAAYCk6dbC1bLZcvCVbOnL2ImcvcnI+Or0TPtttW44JcvYiZy9ycj46dbC1/JW4/c51/JbNOfkflpw9yNmLnL0oMnB+tuWYIGcPcvYiZy862PnzrwwAJhpdVAEAWIxOHcwb5AIAAADMyGQLAEx1uqwCALAUjTqYyRYAAACAGZlsAYCJTh87CACwFJ06mM0WAJjo9LGDAABL0amDuY0IAAAAYEYmWwBgotFFFQCAxejUwUy2AAAAAMzIZAsATHW6rAIAsBSNOpjNFgCY6PRO+AAAS9Gpg7mNCAAAAGBGJlsAYKLTxw4CACxFpw5msgUAAABgRiZbAGCi0UUVAIDF6NTBbLYAwESnEVYAgKXo1MHcRgQAAAAwI5MtAHAvjS6rAAAsRp8OZrIFAAAAYEYmWwBgotP9wgAAS9Gpg5lsAQAAAJiRyRYAmGh0UQUAYDE6dTCbLQAw0WmEFQBgKTp1MLcRAQAAAMzIZAsATFSrIVYAgGXo1MHaTba8593vynOuekaeffhpeeMbjm56OSsjZy9y9iInsG225XggZy9y9iIne02rzZYTJ07kVa+8Jq993bV5y3XX58Yb3pqPHzu26WXNTs5e5OxFziZqRQ9oqP3xYJecvcjZi5yNNOpgZ9xsqaqvraov2f36QVX1g1X1a1X1I1X10PUs8dx9+EMfzKFDj8hlhw7lwosuyuFnXZV3vuPtm17W7OTsRc5e5Oyh0XmehVpSB+t+PDhJzl7k7EXOPjp1sLNNtrwpyed3v/7xJA9N8iO73/upFa7rfrltZyeXXHrJPc8PHDyYnZ2dDa5oNeTsRc5e5ARmspgOti3HAzl7kbMXOdmLzvYGufvGGHfufn3lGOOrd7/+jar6wArXdb+MjHt9rzp9dtQuOXuRsxc5e2gUheVaTAfrfjw4Sc5e5OxFzj46xTnbZMuHq+rbd7/+vaq6Mkmq6iuSfOG+/lBVHamqm6rqpqNH1/emPQcPXpJbP33rPc9v29nJgQMH1vb710XOXuTsRU5gJovpYNtyPJCzFzl7kZO96GybLS9K8qSq+niSRyf5zar6RJI37P7stMYYR8cYV44xrjxy5Mh8qz2Lxzz2cbn55k/m+PFb8oU77siNN1yfJz3lqWv7/esiZy9y9iJnD7Wi/8F5WEwH6348OEnOXuTsRc4+OnWwM95GNMb4bJJvq6qHJPmLu68/PsbYkzeG7d+/Py97+Svy4iMvyl13ncjVz3t+Lr/8ik0va3Zy9iJnL3I2YV+EDVtSB2t/PNglZy9y9iJnI406WI1x7/u+ZjZuv/PsL1qyi3e3rOTsQc5e5OxlN+fKT8Of+dydKzk5/oUH729UIVgAHawJOXuRs5cty6mDnYezvUEuAGwdOyIAAOvXqYOd7T1bAAAAADgPJlsAYKLTxw4CACxFpw5msgUAAABgRiZbAGDCxzQDAKxfpw5mswUAJjqNsAIALEWnDuY2IgAAAIAZ2WwBAAAAmJHNFgAAAIAZec8WAJjodL8wAMBSdOpgNlsAYKLTO+EDACxFpw7mNiIAAACAGZlsAYCJTiOsAABL0amDmWwBAAAAmJHJFgCYaHRRBQBgMTp1MJstADDV6UwPALAUjTqY24gAAAAAZmSyBQAmOn3sIADAUnTqYCZbAAAAAGZksgUAJjp97CAAwFJ06mAmWwAAAABmZLIFACYaXVQBAFiMTh3MZgsATHU60wMALEWjDuY2IgDYQ6rqcFV9rKqOVdVLT/PzL6qq/3P3579VVY9c/yoBAPpYRf+y2QIAE7Wi/53191ZdkOQ1SZ6Z5NFJXlhVj5687DuS/PEY4/IkP5rkR2aODwCwEZvoYKvqXzZbAGDveEKSY2OMT4wx7kjy5iTPnbzmuUl+ZvfrX07y9VWd3rsfAGCtVtK/1vKeLRdvyTvDyNmLnL3IyfnY4NbFw5Pccsrz40m+9r5eM8a4s6o+m+TLkvzhWlbIomzLMUHOXuTsRU7Ox4Y62Er61zomW2rdj6r6O5v4vXLKKaeccq7lsXIX70+t4lFVR6rqplMeRya/+nT5xv14DSTbc0yQU0455ZRzPY+V21AHO122B9y/ut5GNC2vXcnZi5y9yMm9jDGOjjGuPOVxdPKS40kOnfL8siSfuq/XVNX+JA9N8kerWjOcp205JsjZi5y9yMm9nKWDraR/dd1sAYAlel+SK6rqUVV1UZIXJLlu8prrkvyt3a+/Mcm/GWOYbAEAuH9W0r/cWQYAe8TuPcAvSfK2JBckedMY4yNVdU2Sm8YY1yV5Y5Kfq6pjufuKygs2t2IAgGVbVf/qutkyHcvuSs5e5OxFTu6XMcYNSW6YfO8Vp3x9e5JvWve64BxtyzFBzl7k7EVOztsq+leZPAYAAACYj/dsAQAAAJhRu82WqjpcVR+rqmNV9dJNr2cVqupNVXVbVX1402tZpao6VFXvqKrfr6qPVNX3bHpNq1BVF1fVb1fV7+3m/MFNr2lVquqCqnp/Vb1102tZpar6ZFV9qKo+UFU3bXo9q1JVX1pVv1xVf7D73+n/vOk1AZuxDf0r2Y4Opn/1tA0dTP9ir2l1G1FVXZDk3yV5Wu7+aKb3JXnhGOOjG13YzKrq65J8LsnPjjEeu+n1rEpVXZrk0jHG71bVQ5L8TpKrG/7fs5J88Rjjc1V1YZLfSPI9Y4z3bnhps6uq70tyZZIvGWM8e9PrWZWq+mSSK8cYf7jptaxSVf1MknePMa7dfef2PzfG+E+bXhewXtvSv5Lt6GD6V7/+lWxHB9O/2Gu6TbY8IcmxMcYnxhh3JHlzkudueE2zG2O8K2f5TO8OxhifHmP87u7X/znJ7yd5+GZXNb9xt8/tPr1w99FnF3RXVV2W5Kok1256LTxwVfUlSb4ud78ze8YYdzjRw9baiv6VbEcH07/60cH60L+Wpdtmy8OT3HLK8+NpeHLYRlX1yCSPT/Jbm13JauyOdn4gyW1J/tUYo2POH0vy95PctemFrMFI8utV9TtVdWTTi1mRv5jkM0l+ancs+dqq+uJNLwrYCP2rKf2rjW3pYPoXe0q3zZY6zfda7lBvk6p6cJJfSfK9Y4w/2fR6VmGMcWKM8VVJLkvyhKpqNZpcVc9OctsY43c2vZY1eeIY46uTPDPJd+2OnXezP8lXJ/nJMcbjk/yXJG3fpwE4I/2rIf2rhy3rYPoXe0q3zZbjSQ6d8vyyJJ/a0FqYwe49tL+S5OfHGP/XptezartjgO9McnjDS5nbE5M8Z/de2jcneWpV/R+bXdLqjDE+tfv/3pbkLbl7xL6b40mOn3IV8Jdz98kf2D76VzP6Vytb08H0L/aabpst70tyRVU9avfNgl6Q5LoNr4n7afeNy96Y5PfHGP980+tZlar6C1X1pbtfPyjJNyT5g82ual5jjJeNMS4bYzwyd/93+W/GGN+y4WWtRFV98e4bCmZ3rPPpSdp9asUY49Ykt1TVf7/7ra9P0urNE4Fzpn81on/1si0dTP9iL9q/6QXMaYxxZ1W9JMnbklyQ5E1jjI9seFmzq6pfTPLkJA+rquNJfmCM8cbNrmolnpjkbyb50O79tEnyD8cYN2xwTatwaZKf2f00h31JfmmM0fZj+bbAwSRvuburZn+SXxhj3LjZJa3Mdyf5+d3/z9Unknz7htcDbMC29K9kazqY/sUS6V/sOa0++hkAAABg07rdRgQAAACwUTZbAAAAAGZkswUAAABgRjZbAAAAAGZkswUAAABgRjZbAAAAAGZkswUAAABgRjZbAAAAAGb0/wN501a0zVtMdQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f11fbac8630>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "titles = ['X', 'Y']\n",
    "\n",
    "# display omega\n",
    "f, axes = plt.subplots(1, 2, figsize=(20, 8))\n",
    "for i in range(2):\n",
    "    sns.heatmap(DataFrame(initial_omega[0]), cmap='Blues', annot=True, linewidths=.5, \n",
    "                ax=axes[i])\n",
    "    axes[i].set_title(titles[i])\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAdoAAAE/CAYAAADhbQKeAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAHYtJREFUeJzt3X+Q3PV93/HX+7jDsQLFjqU7uToVDCHVGNPag4I7YCoCkWUHDVKKQ7Ebp01Nrh0mrQ2dBGfSgmHagNtR0swo03ijnB2oEk9jmtSVxiJIjSygwbFoXVVYpJGMCkLoTnId21gZg1bv/nEnRXxOd3v7/ex3P+/dfT5mdrgfut3PDXrq/f2x+11zdwEAgHoMlV4AAAD9jEELAECNGLQAANSIQQsAQI0YtAAA1IhBCwBAjRi0AADUiEEblJldZGaHzewj53ztYjN70cw+VHJtABbPzLaa2WTytTVm9k0ze3updaF7jAtWxGVm75e0VdI73f24mf0HSWPu/vcKLw3AIpnZ2yQ9J+mj7v6Emf2ApH2SfsXdP1d0cegKBm1wZvY5SW+S9BlJj0l6l7u/UnRRANpiZj8l6d9Kepekfynp3e7+wbKrQrcwaIMzs7dK+rqkEUm/4O6fLbwkABWY2RckXSjpeknvcfcXCy8JXcKg7QFmtlPSdZLe7u7fLr0eAO0zszFJhyT9srv/eun1oHt4MlRwZvbTki6TtFPSp8uuBkBV7j4l6YRmztdigAyXXgDmZ2ajkn5N0u2Snpf0nJn9rrvvKbsyAMBisUcb22ZJf+jufzz7BKhflPRbZvamwusCACwSgzYoM9so6X2SfuHM19x9i6Qjku4rtS4AQHt4MhQAADVijxYAgBrxZCigS8zssKTvSmpKOuXuq8uuCEBV7fTMoAW668fc/UTpRQDoiEX1zKFjAABqxKAFuscl/ZGZPWtmE6UXAyDLonvuxqFjntaMKGwxf+hTq0ba/jv7wJ+d+ieSzo2t4e6N5I9d7+5HZy9E8oSZPd+DFx+hZ0RQW8tS53vuyjnaT60a6cbDAPP61POv13r/sxGmIaZ/5ujsf6fN7A8kXSuppwYtLaO0uluWOt8zh46BhFW4tbxPsx80s4vPfCzp/ZL2d3zxAM6q0nIdPfOsYyBhizoo1bYxSX9gM3c+LOl33X1HLY8EQFJtLUtt9sygBRJ1HOZx929I+ts13DWAedR1yLbdnhm0QKLGrWAAXRSlZQYtkAjSJoBMUVpm0AKJKFvBAPJEaZlBCyR4Kj7QH6K0zKAFElG2ggHkidIygxZIBGkTQKYoLTNogUSUrWAAeaK0zKAFEkHaBJApSssMWiAxFKVOAFmitMygBRJB2gSQKUrLDNrgPrHrz/X9770qbzZ1unlKjQ/9ndJL6ntRzuugv9By90VpmUHbA37nZ35cJ//im6WXMTCCtIk+RMvdFaXlloPWzFZJ2iBphWbe9PmopC+6+4Ga1wYUMWT9+d7mtIxBE6XlBS+cYWb3Svq8ZjYM/lTSV2c//j0z+2T9y4O766O//SVNPPYVXXP7naWXgx5Fy+XR8uBqtUf7MUlXufsb3tLezH5V0nOSHj7fD5nZhKQJSfrMZz7TgWUOrsmPrNF3p1/RD/7QMn10codOfON5/d+9T5VeVl+Lcripwyq1PPtnzva8fmxI17wlyoXtegstd1+UllsVc1rSXz/P198++73zcveGu69299UTExM56xt4351+RZL0vf93XM/v/EOt+Fs/WnhF/c8q3HpApZalN/bMkK2OlruvSst19Nxqj/YTknaZ2Z9Lemn2a39D0g9L+vka1oNzjLx5iWxoSK9971WNvHmJrrh+rb78G/+69LL6XpRnKnYYLRdEy2VEaXnBQevuO8zsRyRdq5knUJikI5K+6u7NLqxvoF30tjH9/c1fkCQNXXCB/ve2z+vgU39UeFX9L0ibHUXLZdFyGVFabvmsY3c/LemZLqwFiW8deUG/ufGa0ssYOFGuJtNptFwOLZcRpWVeRwskgrQJIFOUlhm0QCLKeR0AeaK0zKAFEkHaBJApSssMWiAR5bwOgDxRWmbQAokgbQLIFKVlBi2QiHJeB0CeKC0zaIFEkDYBZIrSMoMWSETZCgaQJ0rLDFogwdV8gf4QpWUGLZCIshUMIE+Ulhm0QCJImwAyRWk5yp41AAB9iT1aIBHlRe4A8kRpmUELJIK0CSBTlJa7Mmg/9fzr3XgYoCOibAVHRMvoJVFaZo8WSPDEBaA/RGm5O4P25ImuPAwwryVLF/1Ho7wkICRaRmk92DJ7tEAiylYwgDxRWmbQAokoW8EA8kRpmUELJIbMSy8BQAdEaZlBCySiHG4CkCdKywxaIBHlcBOAPFFaZtACiShbwQDyRGmZQQskomwFA8gTpWUGLZCIshUMIE+Ulhm0QCLKZdsA5InSMoMWSEQ53AQgT5SWo+xZA2EMVbgtlpldYGb/08y2dXTRAOao0vJie26nZQYt0F0fl3Sg9CIAZFt0ywza4PY8/YzWbbxDa2+9XY3JR0svZyCYtX9b3P3auKRbJG2pc/2IiZa7r0rLi+m53ZYZtIE1m009+PAmbdm8Sdsf26ptO3bq4KEXSi+r71U51GRmE2a295zbxHnu+t9L+kVJp+v/LRAJLZdR9dDxInpuq2WeDBXYvv0HdOnKca0cXyFJumXdzdq1+0n98BXvKLyy/lblmYru3pDUmO/7ZrZe0rS7P2tmN1ZeHHoSLZdR9VnHC/VcpWX2aAObmj6u5WOjZz8fGxvV1PHjBVc0GKzCbRGul3SrmR2W9HlJN5nZf+zsyhEVLZdRpeVF9Nx2y5UHrZn9bNWfxeK45r7zhC32n3VUNmTt31px919y93F3v0zSHZL+m7v/dM2/yqLRc71ouYwqLbfquUrLOXu0D8z3jXOPbzca8x5NQwvLR0d1bGr67OdTU9MaXba04IoGQ017tNEtrufJR7q5pr5By2XUtEfbtgXP0ZrZvvm+JWlsvp9Ljm+7Tp6otroBd/VVq3T4xSN66eWjGhtdpu2P79Kmh+4vvay+V/fVZNx9t6Td9T7KXB3p+eSJGG/w2WNouYwoLbd6MtSYpHWSvpV83ST99yoLw+INDw/rvnvv1p133aPm6aZu27BeV15xeell9b0obxZdA3ouhJbLiNJyq0G7TdJF7v619BtmtruWFeEN1txwndbccF3pZQyUPjkUfD70XBAtd1+UlhcctO7+sQW+95HOLwcoL8qFyDuNnjFoorTM62iBRJA2AWSK0jKDFkhE2QoGkCdKywxaIMFVXID+EKVlBi2QiPIelgDyRGmZQQskohxuApAnSssMWiARpE0AmaK0HOUQNgAAfYk9WiBhUU7sAMgSpWUGLZAI0iaATFFaZtACqSh1AsgTpGUGLZAI0iaATFFaZtACiSjndQDkidIygxZIRIkTQJ4oLTNogRQvegP6Q5CWuzNolyztysMAnRBlKzgkWkYPidJydwbtyRNdeRhgXm0MiCBtxkTLKK0HW+bQMZCIshUMIE+Ulhm0QCpGmwByBWmZQQskomwFA8gTpWUGLZAI0iaATFFaZtACiShbwQDyRGmZQQukgsQJIFOQlhm0QCJImwAyRWmZQQskohxuApAnSstBLlAFAEB/Yo8WSATZCAaQKUrLDFogFaVOAHmCtMygBRJB2gSQKUrLDFogEeUJFADyRGmZJ0MFt+fpZ7Ru4x1ae+vtakw+Wno5A8HM2r4BrdBy91VpuY6eGbSBNZtNPfjwJm3ZvEnbH9uqbTt26uChF0ovq++ZtX8DFkLLZVRpuY6eGbSB7dt/QJeuHNfK8RW6cGREt6y7Wbt2P1l6Wf0vQpnoK7RcSJBJ23LQmtkqM7vZzC5Kvv6Bjq8GbzA1fVzLx0bPfj42Nqqp48cLrmgwBOiyFrRcDi2XEWTOLjxozeyfS/ovkv6ZpP1mtuGcb/9K55eDc7l8ztcsyhss9rEI53Q6jZbLouUyeuUc7c9JusbdN0q6UdK/MrOPn/kd5vshM5sws71mtrfRaHRmpQNo+eiojk1Nn/18ampao8uWFlzRYIiwBVyDSi1LSc+Tj9S8zP5Ey2VE2aNt9fKeC9z9VUly98NmdqOkL5jZpVogTndvSDozYV0nT3RirQPn6qtW6fCLR/TSy0c1NrpM2x/fpU0P3V96Wf2vRyZnmyq1PPvn/6rnkyfm7pqhJVouJEjLrQbtMTN7t7t/TZLc/VUzWy9pUtLVta9uwA0PD+u+e+/WnXfdo+bppm7bsF5XXnF56WX1vV44FFwBLRdEy2VEadnc599ANbNxSafc/dh5vne9uz+9iMdgjxblLVkqtdhzO+Pb//Rtbe+1XfKb34xR9Dw61DJ7tCiv5palzve84B6tux9Z4HuLCxPoMVG2gjuJljGIorTMJRiBVIw2AeQK0jKDFkjYENdxAfpBlJZjrAIAgD7FHi2QquG8jpn9gKQ9kt6kme6+4O68vgOoU5CWGbRAqp4nUHxf0k2zL6sZkfSUmX3J3Z+p48EAKEzLDFogYdb5Myo+8zq6V2c/HZm98VIZoEZRWuYcLZCq6ZptZnaBmX1N0rSkJ9z9K7X+HsCgq+kajO22zKAFUhXCPPd6wLO3ifRu3b3p7u+WNC7pWjN7V/d/OWCAVBy0rXput2UOHQOJKi9yT67v3erP/oWZ7Zb0AUn7234wAItS9YIVi+15sS2zRwukbKj9W6u7NFtmZm+Z/fjNkn5c0vM1/ybAYKvScoueq7TMHi2QsKFanqn4dkm/Y2YXaGYD9z+5+7Y6HgjAjCgtM2iBVA0vCXD3fZLe0/E7BjC/IC0zaIFUDS8JAFBAkJYZtEAiyjt+AMgTpeXuDNqZ9w8EekOQOEOiZfSSIC13Z9Dyxu8orZ0BESTOkGgZpfVgyxw6BhJ1XLYNQPdFaZlBC6SCbAUDyBSkZQYtkKjptXcAuixKyzH2qwEA6FPs0QKpIOd1AGQK0jKDFkgFOa8DIFOQlhm0QCLKi9wB5InSMoMWSAWJE0CmIC0zaIFUkPM6ADIFaZlBCySiHG4CkCdKywxaIBXktXcAMgVpmUELJKJctg1AnigtM2iBVJDDTQAyBWk5xrjHvPY8/YzWbbxDa2+9XY3JR0svZzCYtX8DWqDlAqq0XEPPDNrAms2mHnx4k7Zs3qTtj23Vth07dfDQC6WX1ffMrO0bsBBaLqNKy3X03HLQmtm1Zvajsx+/08zuMbOf6PhKMMe+/Qd06cpxrRxfoQtHRnTLupu1a/eTpZfV/2yo/VsPoOVyaLmQKi3X0POC52jN7H5JH5Q0bGZPSHqvpN2SPmlm73H3f9PxFeGsqenjWj42evbzsbFR7dv/XMEVDYg+3EOl5bJouZAgLbd6MtSHJL1b0pskHZM07u7fMbN/J+krkoizRi6f8zVTjL84/axPDwXTckG0XEaUllvtI59y96a7n5R0yN2/I0nu/peSTs/3Q2Y2YWZ7zWxvo9Ho4HIHy/LRUR2bmj77+dTUtEaXLS24ogExNNT+Lb5KLUtJz5OPdGOtfYeWC6nScg09t9qjfc3MlszGec2ZL5rZJVogTndvSDozYV0nT2QvdBBdfdUqHX7xiF56+ajGRpdp++O7tOmh+0svq/8F2QrusEotS0nPJ0/M3TVDS7RcSJCWWw3av+vu35ckdz83xhFJ/7C2VUGSNDw8rPvuvVt33nWPmqebum3Del15xeWll4XeRMsF0fJgW3DQngnzPF8/IYnd1C5Yc8N1WnPDdaWXMVh65FnE7aDl8mi5gCAtc2UoIBXkcBOATEFaZtACqSBbwQAyBWmZQQukgmwFA8gUpGUGLZAKshUMIFOQlhm0QCrIVjCATEFaZtACqSBbwQAyBWmZQQukgmwFA8gUpGUGLZAKEieATEFaZtACqSCHmwBkCtIygxZIBdkKBpApSMsMWiAVZCsYQKYgLTNogVSQrWAAmYK0zKAFUkG2ggFkCtIygxZIBdkKBpApSMsMWiAVZCsYQKYgLXdn0C5Z2pWHAVAzWgba1p1Be5L3lUZh7QyIIIebQqJllNaDLXPoGEgFOdwEIFOQlhm0QCrIVjCATEFajjHugUhsqP1bq7s0W2lmf2xmB8zsOTP7eBd+E2CwVWm5hp7ZowVS9WwFn5L0L9z9f5jZxZKeNbMn3P3rdTwYANW5R9tWzwxaIFXDeR13f0XSK7Mff9fMDkhaIYlBC9SlpnO07fbMoAVSNZ/XMbPLJL1H0ldqfSBg0HXhHO1ieuYcLZCqcE7HzCbMbO85t4nz3rXZRZIek/QJd/9Od38xYMBUPEfb6Z7ZowVSFbaC3b0hqbHw3dqIZqLc6u7/udriACxaxT3aTvfMoAVSNZzXMTOT9NuSDrj7r3b8AQDMVdM52nZ75tAxkBqy9m+tXS/po5JuMrOvzd5+ot5fBBhwVVquoWf2aIFUDU+gcPenJMV49TwwKGp6MlS7PTNogVSQy7YByBSkZQYtkApy2TYAmYK0HGPcY157nn5G6zbeobW33q7G5KOllzMYarhkG0DLBdR0CcZ28S9EYM1mUw8+vElbNm/S9se2atuOnTp46IXSy+p/AcJEf6HlQnp10JrZIx1fBc5r3/4DunTluFaOr9CFIyO6Zd3N2rX7ydLL6n8BwuwGWu4eWi4kyKBd8BytmX0x/ZKkHzOzt0iSu9/a8RXhrKnp41o+Nnr287GxUe3b/1zBFaFX0XJZtDzYWj0ZalwzF0neIsk1E+dqSZtqXhckuXzO14xXiNQvyBMoOoyWC6LlQoK03GofebWkZyX9sqRvu/tuSX/p7l929y/P90PnXiey0VjwKlZYwPLRUR2bmj77+dTUtEaXLS24ogER4FBTDSq1LCU9T3K0uQpaLqQXDh27+2lJv2Zmvz/736lWPzP7c+deJ9J18kT2QgfR1Vet0uEXj+ill49qbHSZtj++S5seur/0svpfbwzOtlRtefZn/6rnkyfm7pqhJVouJEjLiw3tiKSfMrNbJPGOI10yPDys++69W3fedY+ap5u6bcN6XXnF5aWX1f+CHG6qAy2XQcuFBGnZ3GvfQGWPFuUtWSppcSfFTn/pl9qOYuiDD8Uoum7s0aK0mluWOt8zV4YCUkEONwHIFKRlBi2QChIngExBWmbQAqkg53UAZArSMoMWSAXZCgaQKUjLDFogFSROAJmCtMygBVJDMeIEkClIywxaIBXkvA6ATEFaZtACqSCHmwBkCtIygxZIBYkTQKYgLTNogVSQw00AMgVpmUELpIJsBQPIFKTlGKsAAKBPsUcLpIJsBQPIFKRlBi2QChIngExBWu7OoJ15WyOgNwR5AkVItIxeEqRl9miBVJCtYACZgrTcnUHLG7+jtHb2xILEGRIto7QebJk9WiAV5HATgExBWmbQAqkgW8EAMgVpmUELpILECSBTkJYZtEAqSJwAMgVpmUELpIZinNcBkClIywxaIBVkKxhApiAtM2iBVJA4AWQK0jKDFkgFiRNApiAtM2iBVJDX3gHIFKRlBi0wR4w4AeSK0TKDFkgFOdwEIFOQlmOsAgCAPsUeLZAKcl4HQKYgLbNHG9yep5/Ruo13aO2tt6sx+Wjp5QyIoQo3YGG0XEKVljvfM/9CBNZsNvXgw5u0ZfMmbX9sq7bt2KmDh14ovaz+Z9b+DVgALRdSpeUaem5r0JrZ+8zsHjN7f8dXgjn27T+gS1eOa+X4Cl04MqJb1t2sXbufLL2s/hcgzG6g5+6h5UJ6YdCa2Z+e8/HPSdos6WJJ95vZJzu+GrzB1PRxLR8bPfv52Niopo4fL7iiQVH+UFMd6LkcWi6lNw4dj5zz8YSkte7+gKT3S/oHHV8N3sDlc75mQV4X1tcCbAHXhJ4LoeVCemGPVtKQmb3VzN4mydz9uCS5+/cknZrvh8xswsz2mtneRqPRweUOluWjozo2NX3286mpaY0uW1pwRQOihjDNbNLMps1sfxd+g/nk9zz5SLfW2ldouZCaBm27PbcatJdIelbSXkk/ZGbLZx/kIi1wyQ13b7j7andfPTExsZh14DyuvmqVDr94RC+9fFSvvf66tj++Szfd+L7SyxoAtRxq+pykD3R+rW3J7/kf/0x3VtpnaLmU2g4df05t9Lzg62jd/bJ5vnVa0k8u9kFQzfDwsO67927dedc9ap5u6rYN63XlFZeXXlb/q+HQkbvvMbPLOn7H7a1hvsen55rRciE1ndZpt+dKF6xw95OSeG56F6y54TqtueG60ssYLBUu22ZmE5o573lGw9174rwJPXcHLRdQ8RKMne6ZK0MBc7S/FTwbYU8MVmBwVNuj7XTPDFog1TvPIgawkCAtM2iBVJB3/ACQKUjLMVYBBGJmbd8WcZ+/J+lPJP1NMztiZh+r/RcBBlyVluvomT1aYI7Ob3+6+4c7fqcAWqhnX7Ldnhm0QCrIeR0AmYK0zKFjAABqxB4tkAqyFQwgU5CWGbTAHBzoAfpDjJYZtEAqyFYwgExBWmbQAqkgcQLIFKRlBi0wR4zDTQByxWiZQQukgmwFA8gUpGUGLZAKctk2AJmCtNydQbtkaVceBuiMGFvBIdEyekqMltmjBVJBDjcByBSk5e4M2pMnuvIwwLza2RMLcrgpJFpGaT3YMnu0QCrIVjCATEFaZtACc8SIE0CuGC0zaIFUkMNNADIFaZlBC8wRYysYQK4YLTNogVSQ8zoAMgVpmUELzBEjTgC5YrQc4wA2AAB9ij1aIBXkcBOATEFaZo8WAIAasUcLpIJsBQPIFKRlBi0wR4w4AeSK0TKDFkgF2QoGkClIywxaYI4YcQLIFaNlngwV3J6nn9G6jXdo7a23qzH5aOnlDAaz9m9AC7RcQJWWa+iZQRtYs9nUgw9v0pbNm7T9sa3atmOnDh56ofSyBoBVuAHzo+VSqrTc5UFrZu81s782+/GbzewBM/uvZvZpM7uk46vBG+zbf0CXrhzXyvEVunBkRLesu1m7dj9Zeln9L8AWcKfRclm0XEiP7NFOSjo5+/GvS7pE0qdnv/bZjq8GbzA1fVzLx0bPfj42Nqqp48cLrmhQlN8CrgEtF0TLpfTAHq2kIXc/Nfvxanf/hLs/5e4PSLp8vh8yswkz22tmexuNRscWO2hcPudr1hv/qPe2AFvANajUspT0PPlI/SvtQ7RcSJA92lbPOt5vZj/r7p+V9L/MbLW77zWzH5H0+nw/5O4NSWcmrOvkiQ4td7AsHx3Vsanps59PTU1rdNnSgisaFH35D2CllqWk55Mn5k4MtETLpcRoudUe7Z2S1pjZIUnvlPQnZvYNSb81+z3U6OqrVunwi0f00stH9drrr2v747t0043vK70s9CZaLoiWB9uCe7Tu/m1J/8jMLtbM4aVhSUfcfaobixt0w8PDuu/eu3XnXfeoebqp2zas15VXLHiUDx1gvXEouC20XBYtlxGlZXOv/UgQh45R3pKl0mKPI33rhfajeOs7YhRdNw4do7S6W5Y63jNXhgJSQbaCAWQK0jKDFpgjRpwAcsVomStDAQBQI/ZogVSQw00AMgVpmUELzBEjTgC5YrTMoAVSQbaCAWQK0jKDFpgjRpwAcsVomUELpIJsBQPIFKRlBi0wR4w4AeSK0TKDFkjFaBNAriAtM2iBOYLUCSBTjJa5YAWQqun9K83sA2b2Z2Z20Mw+WfNvAaCm96Ntt2UGLTCHVbi1uEezCyT9hqQPauZt6j5sZu+sYfEAzqrS8sI9V2mZQQuk6tmjvVbSQXf/hru/JunzkjbU+nsAg66ePdq2W+7OOdqZtzUCekQt53VWSHrpnM+PSHpvHQ9UK1pGT4nRcjcGbYyz0T3MzCbcvVF6HQNjydK2/86a2YSkiXO+1Ej+n53vPnvxvV3pORM9d1GFlqWWPbfdMs867g0TkggzsNkIF/p/dETSynM+H5d0tNZFISp6Dq5Fz223zDlaoDu+KulKM3uHmV0o6Q5JXyy8JgDta7tl9miBLnD3U2b285Iel3SBpEl3f67wsgC0qUrL5t6Lp4kGC+d0gP5Bz4OHQQsAQI04RwsAQI0YtMFx2T6gP9Dy4OLQcWCzl/r6P5LWauYp5V+V9GF3/3rRhQFoCy0PNvZoY+OyfUB/oOUBxqCN7XyX+lpRaC0AqqPlAcagja1fLtsHDDpaHmAM2ti4bB/QH2h5gDFoY+OyfUB/oOUBxiUYA+OyfUB/oOXBxst7AACoEYeOAQCoEYMWAIAaMWgBAKgRgxYAgBoxaAEAqBGDFgCAGjFoAQCoEYMWAIAa/X/YlD7ZNKaijAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f11fbac8390>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# display xi\n",
    "f, axes = plt.subplots(1, 2, figsize=(8, 5))\n",
    "for i in range(2):\n",
    "    sns.heatmap(DataFrame(initial_xi[i]), cmap='Oranges', annot=True, linewidths=.5,\n",
    "                ax=axes[i])\n",
    "    axes[i].set_title(titles[i])\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---\n",
    "## SLAM inputs \n",
    "\n",
    "In addition to `data`, your slam function will also take in:\n",
    "* N -   The number of time steps that a robot will be moving and sensing\n",
    "* num_landmarks - The number of landmarks in the world\n",
    "* world_size - The size (w/h) of your world\n",
    "* motion_noise - The noise associated with motion; the update confidence for motion should be `1.0/motion_noise`\n",
    "* measurement_noise - The noise associated with measurement/sensing; the update weight for measurement should be `1.0/measurement_noise`\n",
    "\n",
    "#### A note on noise\n",
    "\n",
    "Recall that `omega` holds the relative \"strengths\" or weights for each position variable, and you can update these weights by accessing the correct index in omega `omega[row][col]` and *adding/subtracting* `1.0/noise` where `noise` is measurement or motion noise. `Xi` holds actual position values, and so to update `xi` you'll do a similar addition process only using the actual value of a motion or measurement. So for a vector index `xi[row][0]` you will end up adding/subtracting one measurement or motion divided by their respective `noise`.\n",
    "\n",
    "### TODO: Implement Graph SLAM\n",
    "\n",
    "Follow the TODO's below to help you complete this slam implementation (these TODO's are in the recommended order), then test out your implementation! \n",
    "\n",
    "#### Updating with motion and measurements\n",
    "\n",
    "With a 2D omega and xi structure as shown above (in earlier cells), you'll have to be mindful about how you update the values in these constraint matrices to account for motion and measurement constraints in the x and y directions. Recall that the solution to these matrices (which holds all values for robot poses `P` and landmark locations `L`) is the vector, `mu`, which can be computed at the end of the construction of omega and xi as the inverse of omega times xi: $\\mu = \\Omega^{-1}\\xi$\n",
    "\n",
    "**You may also choose to return the values of `omega` and `xi` if you want to visualize their final state!**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "def get_motion_matrix(motion):\n",
    "    return np.array([-motion, motion])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "def print_omega(omega):\n",
    "    print('----------------------------------------')\n",
    "    for row in omega:\n",
    "        print('[', end=' ')\n",
    "        for cell in row:\n",
    "            print(cell, end=' ')\n",
    "        print(']')\n",
    "    print('----------------------------------------')\n",
    "                \n",
    "\n",
    "def print_xi(xi):\n",
    "    print('----------------------------------------')\n",
    "    print('[', end=' ')\n",
    "    for row in xi:\n",
    "        print(row[0], end=' ')\n",
    "    print(']')\n",
    "    print('----------------------------------------')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "import random\n",
    "import pdb\n",
    "## TODO: Complete the code to implement SLAM\n",
    "\n",
    "## slam takes in 6 arguments and returns mu, \n",
    "## mu is the entire path traversed by a robot (all x,y poses) *and* all landmarks locations\n",
    "def slam(data, N, num_landmarks, world_size, motion_noise, measurement_noise):\n",
    "    \n",
    "    ## TODO: Use your initilization to create constraint matrices, omega and xi\n",
    "    omega, xi = initialize_constraints(N, num_landmarks, int(world_size))\n",
    "    \n",
    "    x, y = xi[0, 0, 0], xi[1, 0, 0]\n",
    "    ## TODO: Iterate through each time step in the data\n",
    "    ## get all the motion and measurement data as you iterate\n",
    "    dt_matrix = [[1., -1.], [-1., 1.]]\n",
    "    \n",
    "    for time_step in range(N-1):\n",
    "        measurements, motion = data[time_step]\n",
    "        # Loop to process both omega and xi matrix for x and y axis\n",
    "        for ax in range(2):\n",
    "            omega[ax, time_step:time_step+2, time_step:time_step+2] += np.array(dt_matrix) / motion_noise\n",
    "            xi[ax, time_step:time_step+2, 0] += get_motion_matrix(motion[ax]) / motion_noise\n",
    "            \n",
    "            for measurement in measurements:\n",
    "                landmark_index = measurement[0]\n",
    "                # For omega\n",
    "                omega[ax, time_step, time_step] += 1. / measurement_noise\n",
    "                omega[ax, time_step, N+landmark_index] += -1. / measurement_noise\n",
    "                omega[ax, N+landmark_index, time_step] += -1. / measurement_noise\n",
    "                omega[ax, N+landmark_index, N+landmark_index] += 1. / measurement_noise\n",
    "\n",
    "                # For xi\n",
    "                xi[ax, time_step, 0] -= measurement[1+ax] / measurement_noise\n",
    "                xi[ax, N+landmark_index, 0] += measurement[1+ax] / measurement_noise\n",
    "    \n",
    "            \n",
    "    ## TODO: update the constraint matrix/vector to account for all *measurements*\n",
    "    ## this should be a series of additions that take into account the measurement noise\n",
    "            \n",
    "    ## TODO: update the constraint matrix/vector to account for all *motion* and motion noise\n",
    "    \n",
    "    ## TODO: After iterating through all the data\n",
    "    ## Compute the best estimate of poses and landmark positions\n",
    "    ## using the formula, omega_inverse * Xi\n",
    "    \n",
    "    omega_x_inv = np.linalg.inv(omega[0])\n",
    "    omega_y_inv = np.linalg.inv(omega[1])\n",
    "    mu = [np.matmul(omega_x_inv, xi[0]), np.matmul(omega_y_inv, xi[1])]\n",
    "    mu = np.array(mu)\n",
    "    \n",
    "    return mu # return `mu`\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Helper functions\n",
    "\n",
    "To check that your implementation of SLAM works for various inputs, we have provided two helper functions that will help display the estimated pose and landmark locations that your function has produced. First, given a result `mu` and number of time steps, `N`, we define a function that extracts the poses and landmarks locations and returns those as their own, separate lists. \n",
    "\n",
    "Then, we define a function that nicely print out these lists; both of these we will call, in the next step.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "# a helper function that creates a list of poses and of landmarks for ease of printing\n",
    "# this only works for the suggested constraint architecture of interlaced x,y poses\n",
    "def get_poses_landmarks(mu, N):\n",
    "    # create a list of poses\n",
    "    poses = []\n",
    "    for i in range(N):\n",
    "        poses.append((mu[0, i].item(), mu[1, i].item()))\n",
    "\n",
    "    # create a list of landmarks\n",
    "    landmarks = []\n",
    "    for i in range(num_landmarks):\n",
    "        landmarks.append((mu[0, N+i].item(), mu[1, N+i].item()))\n",
    "\n",
    "    # return completed lists\n",
    "    return poses, landmarks\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "def print_all(poses, landmarks):\n",
    "    print('\\n')\n",
    "    print('Estimated Poses:')\n",
    "    for i in range(len(poses)):\n",
    "        print('['+', '.join('%.3f'%p for p in poses[i])+']')\n",
    "    print('\\n')\n",
    "    print('Estimated Landmarks:')\n",
    "    for i in range(len(landmarks)):\n",
    "        print('['+', '.join('%.3f'%l for l in landmarks[i])+']')\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Run SLAM\n",
    "\n",
    "Once you've completed your implementation of `slam`, see what `mu` it returns for different world sizes and different landmarks!\n",
    "\n",
    "### What to Expect\n",
    "\n",
    "The `data` that is generated is random, but you did specify the number, `N`, or time steps that the robot was expected to move and the `num_landmarks` in the world (which your implementation of `slam` should see and estimate a position for. Your robot should also start with an estimated pose in the very center of your square world, whose size is defined by `world_size`.\n",
    "\n",
    "With these values in mind, you should expect to see a result that displays two lists:\n",
    "1. **Estimated poses**, a list of (x, y) pairs that is exactly `N` in length since this is how many motions your robot has taken. The very first pose should be the center of your world, i.e. `[50.000, 50.000]` for a world that is 100.0 in square size.\n",
    "2. **Estimated landmarks**, a list of landmark positions (x, y) that is exactly `num_landmarks` in length. \n",
    "\n",
    "#### Landmark Locations\n",
    "\n",
    "If you refer back to the printout of *exact* landmark locations when this data was created, you should see values that are very similar to those coordinates, but not quite (since `slam` must account for noise in motion and measurement)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "\n",
      "Estimated Poses:\n",
      "[50.000, 50.000]\n",
      "[40.551, 34.009]\n",
      "[33.306, 15.187]\n",
      "[21.221, 27.682]\n",
      "[9.129, 40.578]\n",
      "[27.691, 47.809]\n",
      "[45.022, 56.613]\n",
      "[64.162, 66.494]\n",
      "[83.182, 74.112]\n",
      "[94.942, 89.339]\n",
      "[80.576, 75.181]\n",
      "[67.245, 61.822]\n",
      "[53.749, 46.700]\n",
      "[40.319, 32.123]\n",
      "[26.275, 17.393]\n",
      "[12.266, 2.467]\n",
      "[15.589, 20.866]\n",
      "[17.271, 40.476]\n",
      "[21.293, 60.518]\n",
      "[24.285, 80.293]\n",
      "\n",
      "\n",
      "Estimated Landmarks:\n",
      "[33.851, 83.343]\n",
      "[18.713, 52.043]\n",
      "[80.792, 26.007]\n",
      "[54.592, 47.998]\n",
      "[95.081, 20.628]\n"
     ]
    }
   ],
   "source": [
    "# call your implementation of slam, passing in the necessary parameters\n",
    "mu = slam(data, N, num_landmarks, world_size, motion_noise, measurement_noise)\n",
    "\n",
    "# print out the resulting landmarks and poses\n",
    "if(mu is not None):\n",
    "    # get the lists of poses and landmarks\n",
    "    # and print them out\n",
    "    poses, landmarks = get_poses_landmarks(mu, N)\n",
    "    print_all(poses, landmarks)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Visualize the constructed world\n",
    "\n",
    "Finally, using the `display_world` code from the `helpers.py` file (which was also used in the first notebook), we can actually visualize what you have coded with `slam`: the final position of the robot and the positon of landmarks, created from only motion and measurement data!\n",
    "\n",
    "**Note that these should be very similar to the printed *true* landmark locations and final pose from our call to `make_data` early in this notebook.**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Last pose:  (24.285446291620133, 80.29295521414194)\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABIoAAARfCAYAAABz4ghpAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzs3XuwbnddJvjnTXZuhDeYhHU4gcCoRGmCNkKbYEBuoRLA4EzCgOgUpaYhOoWdNKSpHimkgBHEUScmTA0FaZwi2NU2GiAgNHLIGSHcLxokzWXUyVgNaHJeMDnZueecrPkjiyCFIXvnPYv1/a18Pv8YTp783getb5X11H53Fn3f9wEAAADgfu+gqQsAAAAAUIOhCAAAAIAkhiIAAAAABoYiAAAAAJIYigAAAAAYGIoAAAAASLKFoegVr3hFTjnllDznOc+5+8+uv/76nH322Tn99NNz9tlnZ+/evUmSvu/zute9Lqeddlp+9md/Nl/84hfHaw4AAADAAXWvQ9Fzn/vcvPWtb/2OP7v44otzyimnZNeuXTnllFNy8cUXJ0muuOKK/N3f/V127dqV3/zN38xrXvOaUUoDAAAAcODd61B00kkn5UEPetB3/Nnu3btz5plnJknOPPPMXH755d/x54vFIj/xEz+RG264IXv27BmhNgAAAAAH2n36HUXf/OY3s2PHjiTJjh078o//+I9JkmuvvTY7d+68O7dz585ce+21B6AmAAAAAGM7oL/Muu/77/qzxWJxID8CAAAAgJFs3Jd/6Nhjj82ePXuyY8eO7NmzJ8ccc0ySu36C6Jprrrk7d80119z9k0f/nNVqc0uf13XL2War9GgtW6VHhex28l23TDLf26vSo0K2So/WsmO93eLtVelRIVulR4VslR6t3d6Yb7eWrdKjQrZKjwq3N1aPuWer9KiQrdKjSnZd9+knik499dRcdtllSZLLLrssz3jGM77jz/u+z+c///ksl8vvORQBAAAAUMe9/kTR+eefn8985jO57rrr8pSnPCXnnntufuVXfiUvfelLc+mll+a4447LRRddlCR56lOfmo985CM57bTTcsQRR+S3fuu3Rv8vAAAAAMCBca9D0QUXXPDP/vkll1zyXX+2WCzy6le/ev1WAAAAAHzfHdBfZg0AAABAuwxFAAAAACQxFAEAAAAwMBQBAAAAkMRQBAAAAMDAUAQAAABAkmTR930/dQkAAAAAprcx5YevVptbynXdcrbZKj1ay1bpUSG7nXzXLZPM9/aq9KiQrdKjtexYb7d4e1V6VMhW6VEhW6VHa7c35tutZav0qJCt0qPC7Y3VY+7ZKj0qZKv0qJJdl6+eAQAAAJDEUAQAAADAwFAEAAAAQBJDEQAAAAADQxEAAAAASQxFAAAAAAwMRQAAAAAkMRQBAAAAMDAUAQAAAJAkWfR9309dAgAAAIDpbUz54avV5pZyXbecbbZKj9ayVXpUyG4n33XLJPO9vSo9KmSr9GgtO9bbLd5elR4VslV6VMhW6dHa7Y35dmvZKj0qZKv0qHB7Y/WYe7ZKjwrZKj2qZNflq2cAAAAAJDEUAQAAADAwFAEAAACQxFAEAAAAwMBQBAAAAEASQxEAAAAAA0MRAAAAAEkMRQAAAAAMDEUAAAAAJEkWfd/3U5cAAAAAYHobU374arW5pVzXLWebrdKjtWyVHhWy28l33TLJfG+vSo8K2So9WsuO9XaLt1elR4VslR4VslV6tHZ7Y77dWrZKjwrZKj0q3N5YPeaerdKjQrZKjyrZdfnqGQAAAABJDEUAAAAADAxFAAAAACQxFAEAAAAwMBQBAAAAkMRQBAAAAMDAUAQAAABAEkMRAAAAAANDEQAAAABJkkXf9/3UJQAAAACY3saUH75abW4p13XL2War9GgtW6VHhex28l23TDLf26vSo0K2So/WsmO93eLtVelRIVulR4VslR6t3d6Yb7eWrdKjQrZKjwq3N1aPuWer9KiQrdKjSnZdvnoGAAAAQBJDEQAAAAADQxEAAAAASQxFAAAAAAwMRQAAAAAkMRQBAAAAMDAUAQAAAJDEUAQAAADAYNH3fT91CQAAAACmtzHlh69Wm1vKdd1yttkqPVrLVulRIbudfNctk8z39qr0qJCt0qO17Fhvt3h7VXpUyFbpUSFbpUdrtzfm261lq/SokK3So8LtjdVj7tkqPSpkq/Sokl2Xr54BAAAAkMRQBAAAAMDAUAQAAABAEkMRAAAAAANDEQAAAABJDEUAAAAADAxFAAAAACQxFAEAAAAwMBQBAAAAkCRZ9H3fT10CAAAAgOltTPnhq9XmlnJdt5xttkqP1rJVelTIbiffdcsk8729Kj0qZKv0aC071tst3l6VHhWyVXpUyFbp0drtjfl2a9kqPSpkq/SocHtj9Zh7tkqPCtkqPapk1+WrZwAAAAAkMRQBAAAAMDAUAQAAAJDEUAQAAADAwFAEAAAAQBJDEQAAAAADQxEAAAAASQxFAAAAAAwMRQAAAAAkSRZ93/dTlwAAAABgehtTfvhqtbmlXNctZ5ut0qO1bJUeFbLbyXfdMsl8b69KjwrZKj1ay471dou3V6VHhWyVHhWyVXq0dntjvt1atkqPCtkqPSrc3lg95p6t0qNCtkqPKtl1+eoZAAAAAEkMRQAAAAAMDEUAAAAAJDEUAQAAADAwFAEAAACQxFAEAAAAwMBQBAAAAEASQxEAAAAAA0MRAAAAAEmSRd/3/dQlAAAAAJjexpQfvlptbinXdcvZZqv0aC1bpUeF7HbyXbdMMt/bq9KjQrZKj9ayY73d4u1V6VEhW6VHhWyVHq3d3phvt5at0qNCtkqPCrc3Vo+5Z6v0qJCt0qNKdl2+egYAAABAEkMRAAAAAANDEQAAAABJDEUAAAAADAxFAAAAACQxFAEAAAAwMBQBAAAAkMRQBAAAAMDAUAQAAABAkmTR930/dQkAAAAAprcx5YevVptbynXdcrbZKj1ay1bpUSG7nXzXLZPM9/aq9KiQrdKjtexYb7d4e1V6VMhW6VEhW6VHa7c35tutZav0qJCt0qPC7Y3VY+7ZKj0qZKv0qJJdl6+eAQAAAJDEUAQAAADAwFAEAAAAQBJDEQAAAAADQxEAAAAASQxFAAAAAAwMRQAAAAAkMRQBAAAAMDAUAQAAAJAkWfR9309dAgAAAIDpbUz54avV5pZyXbecbbZKj9ayVXpUyG4n33XLJPO9vSo9KmSr9GgtO9bbLd5elR4VslV6VMhW6dHa7Y35dmvZKj0qZKv0qHB7Y/WYe7ZKjwrZKj2qZNflq2cAAAAAJDEUAQAAADAwFAEAAACQxFAEAAAAwMBQBAAAAEASQxEAAAAAA0MRAAAAAEkMRQAAAAAMDEUAAAAAJEkWfd/3U5cAAAAAYHobU374arW5pVzXLWebrdKjtWyVHhWy28l33TLJfG+vSo8K2So9WsuO9XaLt1elR4VslR4VslV6tHZ7Y77dWrZKjwrZKj0q3N5YPeaerdKjQrZKjyrZdfnqGQAAAABJDEUAAAAADAxFAAAAACQxFAEAAAAwMBQBAAAAkMRQBAAAAMDAUAQAAABAEkMRAAAAAANDEQAAAABJkkXf9/3UJQAAAACY3saUH75abW4p13XL2War9GgtW6VHhex28l23TDLf26vSo0K2So/WsmO93eLtVelRIVulR4VslR6t3d6Yb7eWrdKjQrZKjwq3N1aPuWer9KiQrdKjSnZdvnoGAAAAQBJDEQAAAAADQxEAAAAASQxFAAAAAAwMRQAAAAAkMRQBAAAAMDAUAQAAAJDEUAQAAADAYNH3fT91CQAAAACmtzHlh69Wm1vKdd1yttkqPVrLVulRIbudfNctk8z39qr0qJCt0qO17Fhvt3h7VXpUyFbpUSFbpUdrtzfm261lq/SokK3So8LtjdVj7tkqPSpkq/Sokl2Xr54BAAAAkMRQBAAAAMDAUAQAAABAEkMRAAAAAANDEQAAAABJDEUAAAAADAxFAAAAACQxFAEAAAAwMBQBAAAAkCRZ9H3fT10CAAAAgOltTPnhq9XmlnJdt5xttkqP1rJVelTIbiffdcsk8729Kj0qZKv0aC071tst3l6VHhWyVXpUyFbp0drtjfl2a9kqPSpkq/SocHtj9Zh7tkqPCtkqPapk1+WrZwAAAAAkMRQBAAAAMDAUAQAAAJDEUAQAAADAwFAEAAAAQBJDEQAAAAADQxEAAAAASQxFAAAAAAwMRQAAAAAkSRZ93/dTlwAAAABgehtTfvhqtbmlXNctZ5ut0qO1bJUeFbLbyXfdMsl8b69KjwrZKj1ay471dou3V6VHhWyVHhWyVXq0dntjvt1atkqPCtkqPSrc3lg95p6t0qNCtkqPKtl1+eoZAAAAAEkMRQAAAAAMDEUAAAAAJDEUAQAAADAwFAEAAACQxFAEAAAAwMBQBAAAAEASQxEAAAAAA0MRAAAAAEmSRd/3/dQlAAAAAJjexpQfvlptbinXdcvZZqv0aC1bpUeF7HbyXbdMMt/bq9KjQrZKj9ayY73d4u1V6VEhW6VHhWyVHq3d3phvt5at0qNCtkqPCrc3Vo+5Z6v0qJCt0qNKdl2+egbA2j7wi+/Jm3ZckC+89crv+nuf/u2P5007Lsifv2zXBM0AAIDtMBQBsLanX3R6Hnj8Mp987RVZXbXn7j+/evfV+csLP5OjH3Vsfvr1T5+wIQAAsBWGIgDWdvjRR+S0N5+Rfn+fXee8L3fceHtuXt2cd7/w3Tn4sIPzzP9wRg55wCFT1wQAAO7FpL+jCID5OO7kh+bkX39iPvW6j+XDL788t3zzltx4zY152gWn5Zh/8eCp6wEAAFtgKALggHncuSfl6x//av7mXV9JkvzYL/xYTnzhj0/cCgAA2CpfPQPggFksFvnhn/mRu//zT730pyZsAwAAbJehCIAD5vqrr8snXvORHPYDh2Vx0CLvffF7s+/WfVPXAgAAtshQBMABsf+2fdl1zvtzx8135LS3nJHH/9uTs+eqPfnYb3x46moAAMAW+R1FABwQn3jNFfnGVXvyuHNPyiOe/oM5/imPyJ7P/H2+9PYv5PgnPzwn/A+PmroiAABwLxZ93/dTlwCgbV+57Ct5x1nvyMOe8LD864/96xy0cdcPrO796t685Sfekjv33ZlfvfJXc/QPHz1xUwAA4HuZdCharTa3lOu65WyzVXq0lq3So0J2O/muWyaZ7+1V6VEh+/3ssfm1G/LHp/5h+juTn9v9whz13z3oO7KfefuV+cAvvTc7HveQnPWnP5+DDz140r5TvN3i7VXpUSFbpUeFbJUerd3emG+3lq3So0K2So8KtzdWj7lnq/SokK3So0p2Xb56BsBalscflRf99a/d49//oWefkJfsOf/72AgAALiv/DJrAAAAAJIYigAAAAAYGIoAAAAASGIoAgAAAGBgKAIAAAAgiaEIAAAAgIGhCAAAAIAkhiIAAAAABou+7/upSwAAAAAwvY0pP3y12txSruuWs81W6dFatkqPCtnt5LtumWS+t1elR4VslR6tZcd6u8Xbq9KjQrZKjwrZKj1au70x324tW6VHhWyVHhVub6wec89W6VEhW6VHley6fPUMAAAAgCSGIgAAAAAGhiIAAAAAkhiKAAAAABgYigAAAABIYigCAAAAYGAoAgAAACCJoQgAAACAwcbUBQC+py9/OUf+/htz6Ec/koO+/vUs7rg9dz64yx2P/8ncdtbzcvsZPzt1QwAAgNlY9H3fT10C4J/16lcnr399sn//PWee+tTkne9Mjj32+9cLAABgpiYdilarzS3lum4522yVHq1lq/SokN1OvuuWSdq4vSNf95o84I0XJEn6gw/ObWf+j7njyU9Nf/jh2fjyl3L4f/rDHLTakyTZd+KP5boP7E6OOKLE/00qZKv0aC071tst3V61HhWyVXpUyFbp0drtjfl2a9kqPSpkq/SocHtj9Zh7tkqPCtkqPapk1+WrZ0A5G5/9dI74P37/rv9w5JHZ+58uzR2nPOnuv39bkptfcm4e9IKzcsjnr8zGl/5rjvzfXp+bXvO6aQoDAADMhF9mDZTzgP/zjVl864cdf/d3v2Mk+pb+6GNywx/8YfoHPCBJcsTb/iCLvdd/P2sCAADMjqEIqOW223Lo7l1JkjuPOSZ50YvuMXrnwx+RW896XpJkcfNNOfTD//f3oyEAAMBsGYqAUja+eFUWt92WJLnjiU9ODj30e+bveNqp3/5n/+Jzo3YDAACYO0MRUMpB115791/vf+QJ95rf98Pfzhy055pROgEAANxfrPXLrN/2trflT/7kT7JYLPKjP/qjecMb3pA9e/bk/PPPz969e3PiiSfmd37nd3LovfxEAMC3LG789m/z/9bvH/pe+iOP/Cf/7I2jdAIAALi/uM8/UXTttdfm7W9/e975znfmfe97X/bv35/3v//9+b3f+7388i//cnbt2pWjjjoql1566YHsC8xc/8Bv/+scFzfffK/5xU03/ZN/9oGjdAIAALi/WOurZ/v378+tt96affv25dZbb03XdfnUpz6VZz7zmUmSs846K7t37z4gRYH7hzsf8pC7//rgq//fe80f/P99O3PnQ44bpRMAAMD9xaLvv/XvoN6+Sy65JBdeeGEOO+ywPOlJT8orX/nKvOAFL8iHPvShJMk//MM/5Jxzzsn73ve+A1YYmLnbbkse9KC7/ueDH5z8/d8nhxxyz/lzzkne+ta7/vod70h+7ue+Pz0BAABm6D7/RNHevXuze/fu7N69Ox/96Edzyy235Iorrviu3GKxWKsgcD9z2GHJGWfc9dff+Ebytrfdc/arX03+6I/u+usjj0xOP330egAAAHN2n3+Z9Sc+8Ykcf/zxOeaYY5Ikp59+eq688srccMMN2bdvXzY2NnLNNddkx44d9/jGarV5j3/vn+q65WyzVXq0lq3So0J2O/muu+v3/0zd+d6yGy9+SX7gssuyuPPO5N/9u1x33A9m3xN+6jsyi+uvy4N+/rk5ZPgdRTf/0oty0x0Hp0v9/37fj2yVHq1lx3q7ldur2KNCtkqPCtkqPVq7vTHfbi1bpUeFbJUeFW5vrB5zz1bpUSFbpUeV7Lru808UPfShD81f/dVf5ZZbbknf9/nkJz+ZE044IU94whPywQ9+MEny7ne/O6eeeuraJYH7l30/eXJuOfdld/2Hzc38wJnPzvIl5+SwP/qPOezdl+YBb/hfc8yTTsohf/kXd+VP/LHc9L+8csLGAAAA83Cff6LosY99bJ75zGfmrLPOysbGRh796EfnBS94QZ72tKflZS97WS688MI8+tGPzvOf//wD2Re4n7jpla9Of/DBOfKi/z2L/ftz+KXvyOGXvuO7crc/8adzw//1h8kRR0zQEgAAYF7u81CUJOedd17OO++87/izhz/84bn00kvXKgWQJDf/+m/kyBf9Um7+/Tfm0I9+JAd9/etZ3HF77jz2wdn3+J/Mrc99fm5/zn8/dU0AAIDZWGsoAhjdiSfmpt/63dw0dQ8AAID7gfv8O4oAAAAAmBdDEQAAAABJDEUAAAAADAxFAAAAACRJFn3f91OXAAAAAGB6k/5bz1arzS3lum4522yVHq1lq/SokN1OvuuWSeZ7e1V6VMhW6dFadqy3W7y9Kj0qZKv0qJCt0qO12xvz7dayVXpUyFbpUeH2xuox92yVHhWyVXpUya7LV88AAAAASGIoAgAAAGBgKAIAAAAgiaEIAAAAgIGhCAAAAIAkhiIAAAAABoYiAAAAAJIYigAAAAAYGIoAAAAASJIs+r7vpy4BAAAAwPQ2pvzw1WpzS7muW842W6VHa9kqPSpkt5PvumWS+d5elR4VslV6tJYd6+0Wb69KjwrZKj0qZKv0aO32xny7tWyVHhWyVXpUuL2xesw9W6VHhWyVHlWy6/LVMwAAAACSGIoAAAAAGBiKAAAAAEhiKAIAAABgYCgCAAAAIImhCAAAAICBoQgAAACAJIYiAAAAAAaGIgAAAACSJIu+7/upSwAAAAAwvY0pP3y12txSruuWs81W6dFatkqPCtnt5LtumWS+t1elR4VslR6tZcd6u8Xbq9KjQrZKjwrZKj1au70x324tW6VHhWyVHhVub6wec89W6VEhW6VHley6fPUMAAAAgCSGIgAAAAAGhiIAAAAAkhiKAAAAABgYigAAAABIYigCAAAAYGAoAgAAACCJoQgAAACAgaEIAAAAgCTJou/7fuoSAAAAAExvY8oPX602t5TruuVss1V6tJat0qNCdjv5rlsmme/tVelRIVulR2vZsd5u8faq9KiQrdKjQrZKj9Zub8y3W8tW6VEhW6VHhdsbq8fcs1V6VMhW6VEluy5fPQMAAAAgiaEIAAAAgIGhCAAAAIAkhiIAAAAABoYiAAAAAJIYigAAAAAYGIoAAAAASGIoAgAAAGBgKAIAAAAgSbLo+76fugQAAAAA09uY8sNXq80t5bpuOdtslR6tZav0qJDdTr7rlknme3tVelTIVunRWnast1u8vSo9KmSr9KiQrdKjtdsb8+3WslV6VMhW6VHh9sbqMfdslR4VslV6VMmuy1fPAAAAAEhiKAIAAABgYCgCAAAAIImhCAAAAICBoQgAAACAJIYiAAAAAAaGIgAAAACSGIoAAAAAGBiKAAAAAEiSLPq+76cuAQAAAMD0Nqb88NVqc0u5rlvONlulR2vZKj0qZLeT77plkvneXpUeFbJVerSWHevtFm+vSo8K2So9KmSr9Gjt9sZ8u7VslR4VslV6VLi9sXrMPVulR4VslR5Vsuvy1TMAAAAAkhiKAAAAABgYigAAAABIYigCAAAAYGAoAgAAACCJoQgAAACAgaEIAAAAgCSGIgAAAAAGhiIAAAAAkiSLvu/7qUsAAAAAML2NKT98tdrcUq7rlrPNVunRWrZKjwrZ7eS7bplkvrdXpUeFbJUerWXHervF26vSo0K2So8K2So9Wru9Md9uLVulR4VslR4Vbm+sHnPPVulRIVulR5Xsunz1DAAAAIAkhiIAAAAABoYiAAAAAJIYigAAAAAYGIoAAAAASGIoAgAAAGBgKAIAAAAgiaEIAAAAgIGhCAAAAIAkyaLv+37qEgAAAABMb2PKD1+tNreU67rlbLNVerSWrdKjQnY7+a5bJpnv7VXpUSFbpUdr2bHebvH2qvSokK3So0K2So/Wbm/Mt1vLVulRIVulR4XbG6vH3LNVelTIVulRJbsuXz0DAAAAIImhCAAAAICBoQgAAACAJIYiAAAAAAaGIgAAAACSGIoAAAAAGBiKAAAAAEhiKAIAAABgsOj7vp+6BAAAAADT25jyw1erzS3lum4522yVHq1lq/SokN1OvuuWSeZ7e1V6VMhW6dFadqy3W7y9Kj0qZKv0qJCt0qO12xvz7dayVXpUyFbpUeH2xuox92yVHhWyVXpUya7LV88AAAAASGIoAgAAAGBgKAIAAAAgiaEIAAAAgIGhCAAAAIAkhiIAAAAABoYiAAAAAJIYigAAAAAYGIoAAAAASJIs+r7vpy4BAAAAwPQ2pvzw1WpzS7muW842W6VHa9kqPSpkt5PvumWS+d5elR4VslV6tJYd6+0Wb69KjwrZKj0qZKv0aO32xny7tWyVHhWyVXpUuL2xesw9W6VHhWyVHlWy6/LVMwAAAACSGIoAAAAAGBiKAAAAAEhiKAIAAABgYCgCAAAAIImhCAAAAICBoQgAAACAJIYiAAAAAAaGIgAAAACSJIu+7/upSwAAAAAwvY0pP3y12txSruuWs81W6dFatkqPCtnt5LtumWS+t1elR4VslR6tZcd6u8Xbq9KjQrZKjwrZKj1au70x324tW6VHhWyVHhVub6wec89W6VEhW6VHley6fPUMAAAAgCSGIgAAAAAGhiIAAAAAkhiKAAAAABgYigAAAABIYigCAAAAYGAoAgAAACCJoQgAAACAgaEIAAAAgCTJou/7fuoSAAAAAExvY8oPX602t5TruuVss1V6tJat0qNCdjv5rlsmme/tVelRIVulR2vZsd5u8faq9KiQrdKjQrZKj9Zub8y3W8tW6VEhW6VHhdsbq8fcs1V6VMhW6VEluy5fPQMAAAAgiaEIAAAAgIGhCAAAAIAkhiIAAAAABoYiAAAAAJIYigAAAAAYGIoAAAAASGIoAgAAAGBgKAIAAAAgSbLo+76fugQAAAAA09uY8sNXq80t5bpuOdtslR6tZav0qJDdTr7rlknme3tVelTIVunRWnast1u8vSo9KmSr9KiQrdKjtdsb8+3WslV6VMhW6VHh9sbqMfdslR4VslV6VMmuy1fPAAAAAEhiKAIAAABgYCgCAAAAIImhCAAAAICBoQgAAACAJIYiAAAAAAaGIgAAAACSGIoAAAAAGBiKAAAAAEiSLPq+76cuAQAAAMD0Nqb88NVqc0u5rlvONlulR2vZKj0qZLeT77plkvneXpUeFbJVerSWHevtFm+vSo8K2So9KmSr9Gjt9sZ8u7VslR4VslV6VLi9sXrMPVulR4VslR5Vsuvy1TMAAAAAkhiKAAAAABgYigAAAABIYigCAAAAYGAoAgAAACCJoQgAAACAgaEIAAAAgCSGIgAAAAAGhiIAAAAAkiSLvu/7qUsAAAAAML2NKT98tdrcUq7rlrPNVunRWrZKjwrZ7eS7bplkvrdXpUeFbJUerWXHervF26vSo0K2So8K2So9Wru9Md9uLVulR4VslR4Vbm+sHnPPVulRIVulR5Xsunz1DAAAAIAkhiIAAAAABoYiAAAAAJIYigAAAAAYGIoAAAAASGIoAgAAAGBgKAIAAAAgiaEIAAAAgMGi7/t+6hIAAAAATG9jyg9frTa3lOu65WyzVXq0lq3So0J2O/muWyaZ7+1V6VEhW6VHa9mx3m7x9qr0qJCt0qNCtkqP1m5vzLdby1bpUSFbpUeF2xurx9yzVXpUyFbpUSW7Ll89AwAAACCJoQgAAACAgaEIAAAAgCSGIgAAAAAGhiIAAAAAkhiKAAAAABgYigAAAABIYigCAAAAYGAoAgAAACBJsuj7vp+6BAAAAADT25jyw1erzS3lum4522yVHq1lq/SokN1OvuuWSeatkb51AAAgAElEQVR7e1V6VMhW6dFadqy3W7y9Kj0qZKv0qJCt0qO12xvz7dayVXpUyFbpUeH2xuox92yVHhWyVXpUya7LV88AAAAASGIoAgAAAGBgKAIAAAAgiaEIAAAAgIGhCAAAAIAkhiIAAAAABoYiAAAAAJIYigAAAAAYGIoAAAAASJIs+r7vpy4BAAAAwPQ2pvzw1WpzS7muW842W6VHa9kqPSpkt5PvumWS+d5elR4VslV6tJYd6+0Wb69KjwrZKj0qZKv0aO32xny7tWyVHhWyVXpUuL2xesw9W6VHhWyVHlWy6/LVMwAAAACSGIoAAAAAGBiKAAAAAEhiKAIAAABgYCgCAAAAIImhCAAAAICBoQgAAACAJIYiAAAAAAaGIgAAAACSJIu+7/upSwAAAAAwvY0pP3y12txSruuWs81W6dFatkqPCtnt5LtumWS+t1elR4VslR6tZcd6u8Xbq9KjQrZKjwrZKj1au70x324tW6VHhWyVHhVub6wec89W6VEhW6VHley6fPUMAAAAgCSGIgAAAAAGhiIAAAAAkhiKAAAAABgYigAAAABIsuZQdMMNN+S8887Ls571rDz72c/OlVdemeuvvz5nn312Tj/99Jx99tnZu3fvgeoKAAAAwIjWGope//rX58lPfnL+7M/+LO95z3vyyEc+MhdffHFOOeWU7Nq1K6ecckouvvjiA9UVAAAAgBHd56HoxhtvzGc/+9k873nPS5IceuihOeqoo7J79+6ceeaZSZIzzzwzl19++YFpCgAAAMCoFn3f9/flH/zyl7+cV73qVTnhhBPyla98JY95zGPyyle+Mk95ylPyuc997u7cSSedlM9+9rMHrDAAAAAA47jPP1G0b9++fOlLX8ov/MIv5LLLLssRRxzha2YAAAAADdu4r//gzp07s3Pnzjz2sY9NkjzrWc/KxRdfnGOPPTZ79uzJjh07smfPnhxzzDH3+MZqtbmlz+q65WyzVXq0lq3So0J2O/muWyaZ7+1V6VEhW6VHa9mx3m7x9qr0qJCt0qNCtkqP1m5vzLdby1bpUSFbpUeF2xurx9yzVXpUyFbpUSW7rvv8E0Vd12Xnzp25+uqrkySf/OQn88hHPjKnnnpqLrvssiTJZZddlmc84xlrlwQAAABgfPf5J4qS5FWvelVe/vKX54477sjDH/7wvOENb8idd96Zl770pbn00ktz3HHH5aKLLjpQXQEAAAAY0VpD0aMf/ei8613v+q4/v+SSS9Z5FgAAAIAJ3OevngEAAAAwL4YiAAAAAJIYigAAAAAYGIoAAAAASGIoAgAAAGBgKAIAAAAgSbLo+76fugQAAAAA09uY8sNXq80t5bpuOdtslR6tZav0qJDdTr7rlknme3tVelTIVunRWnast1u8vSo9KmSr9KiQrdKjtdsb8+3WslV6VMhW6VHh9sbqMfdslR4VslV6VMmuy1fPAAAAAEhiKAIAAABgYCgCAAAAIImhCAAAAICBoQgAAACAJIYiAAAAAAaGIgAAAACSGIoAAAAAGBiKAAAAAEiSLPq+76cuAQAAAMD0Nqb88NVqc0u5rlvONlulR2vZKj0qZLeT77plkvneXpUeFbJVerSWHevtFm+vSo8K2So9KmSr9Gjt9sZ8u7VslR4VslV6VLi9sXrMPVulR4VslR5Vsuvy1TMAAAAAkhiKAAAAABgYigAAAABIYigCAAAAYGAoAgAAACCJoQgAAACAgaEIAAAAgCSGIgAAAAAGhiIAAAAAkiSLvu/7qUsAAAAAML2NKT98tdrcUq7rlrPNVunRWrZKjwrZ7eS7bplkvrdXpUeFbJUerWXHervF26vSo0K2So8K2So9Wru9Md9uLVulR4VslR4Vbm+sHnPPVulRIVulR5Xsunz1DAAAAIAkhiIAAAAABoYiAAAAAJIYigAAAAAYGIoAAAAASGIoAgAAAGBgKAIAAAAgiaEIAAAAgMGi7/t+6hIAAAAATG9jyg9frTa3lOu65WyzVXq0lq3So0J2O/muWyaZ7+1V6VEhW6VHa9mx3m7x9qr0qJCt0qNCtkqP1m5vzLdby1bpUSFbpUeF2xurx9yzVXpUyFbpUSW7Ll89AwAAACCJoQgAAACAgaEIAAAAgCSGIgAAAAAGhiIAAAAAkhiKAAAAABgYigAAAABIYigCAAAAYGAoAgAAACBJsuj7vp+6BAAAAADT25jyw1erzS3lum4522yVHq1lq/SokN1OvuuWSeZ7e1V6VMhW6dFadqy3W7y9Kj0qZKv0qJCt0qO12xvz7dayVXpUyFbpUeH2xuox92yVHhWyVXpUya7LV88AAAAASGIoAgAAAGBgKAIAAAAgiaEIAAAAgIGhCAAAAIAkhiIAAAAABoYiAAAAAJIYigAAAAAYGIoAAAAASJIs+r7vpy4BAAAAwPQ2pvzw1WpzS7muW842W6VHa9kqPSpkt5PvumWS+d5elR4VslV6tJYd6+0Wb69KjwrZKj0qZKv0aO32xny7tWyVHhWyVXpUuL2xesw9W6VHhWyVHlWy6/LVMwAAAACSGIoAAAAAGBiKAAAAAEhiKAIAAABgYCgCAAAAIImhCAAAAICBoQgAAACAJIYiAAAAAAaGIgAAAACSJIu+7/upSwAAAAAwvY0pP3y12txSruuWs81W6dFatkqPCtnt5LtumWS+t1elR4VslR6tZcd6u8Xbq9KjQrZKjwrZKj1au70x324tW6VHhWyVHhVub6wec89W6VEhW6VHley6fPUMAAAAgCSGIgAAAAAGhiIAAAAAkhiKoAkf+MX35E07LsgX3nrld/29T//2x/PaxWvz5y/bNUEzAAAA5sRQBA14+kWn54HHL/PJ116R1VV77v7zr13x3/KXF34m3Yldfvr1T5+wIQAAAHNgKIIGHH70ETntzWek399n1znvyx033p6bVzfn8pd8IAcfdnCe98fPyyEPOGTqmgAAADTOUASNOO7kh+bkX39i9l59fT788stz+Us+kJv33JSffv3Ts+MxO6auBwAAwAxsTF0A2LrHnXtSvv7xr+Zv3vWVJMmPPPdROfGFPz5xKwAAAObCTxRBQxaLRX74Z37k7v/8L3/l8RO2AQAAYG4Wfd/3U5cAtuabf/PNXPz4i3PQIQfltr23pXtMl3M+c042DvfDgQAAAKxv0qFotdrcUq7rlrPNVunRWrZKj+9ndv9t+/LOn/nP+cZ/3ZPn/Ofn5h8+9fX8xe9/Oif+4r/M8y85a0tvd90yyXxvr0qPCtkqPVrLjvV2i7dXpUeFbJUeFbJVerR2e2O+3Vq2So8K2So9KtzeWD3mnq3So0K2So8q2XX56hk04hOvuSLfuGpPHvdvTsojnv6DOenfn5KdJz80X3r7F/LFP/7i1PUAAACYAUMRNODq//K3ueoPPp+H/KudecIrnpQkOejgg3LaW87IYUcfnj8950+z9++un7glAAAArTMUQXGbX7shf/7SD+bQow7LaW8+IwdtfPtslw9b5tQLT89tN9yWD/3q+7P/9v0TNgUAAKB1fgMuFLc8/qi86K9/7R7//g89+4S8un/1tr6/CwAAAP8cP1EEAAAAQBJDEQAAAAADQxEAAAAASQxFAAAAAAwMRQAAAAAkSRZ93/dTlwAAAABgehtTfvhW/3XeXbecbbZKj9ayVXpUyG4n33XLJPO9vSo9KmSr9GgtO9bbLd5elR4VslV6VMhW6dHa7Y35dmvZKj0qZKv0qHB7Y/WYe7ZKjwrZKj2qZNflq2cAAAAAJDEUAQAAADAwFAEAAACQxFAEAAAAwMBQBAAAAEASQxEAAAAAA0MRAAAAAEkMRQAAAAAMDEUAAAAAJEkWfd/3U5cAAAAAYHobU374arW5pVzXLWebrdKjtWyVHhWy28l33TLJfG+vSo8K2So9WsuO9XaLt1elR4VslR4VslV6tHZ7Y77dWrZKjwrZKj0q3N5YPeaerdKjQrZKjyrZdfnqGQAAAABJDEUAAAAADAxFAAAAACQxFAEAAAAwMBQBAAAAkMRQBAAAAMDAUAQAAABAEkMRAAAAAANDEQAAAABJkkXf9/3UJQAAAACY3saUH75abW4p13XL2War9GgtW6VHhex28l23TDLf26vSo0K2So/WsmO93eLtVelRIVulR4VslR6t3d6Yb7eWrdKjQrZKjwq3N1aPuWer9KiQrdKjSnZdvnoGAAAAQBJDEQAAAAADQxEAAAAASQxFAACz9YFffE/etOOCfOGtV37X3/v0b388b9pxQd774vdO0AwAqMpQBAAwU0+/6PQ88PhlPvnaK7K6as/df/61K/5b/vLCz+ToRx2bZ7/x2RM2BACqMRQBAMzU4UcfkdPefEb6/X12nfO+3HHj7bl5dXMuf8kHcvBhB+eZ/+GMHPKAQ6auCQAUsjF1AQAAxnPcyQ/Nyb/+xHzqdR/Lh19+eW755i25ec9NedoFp+WYf/HgqesBAMUYigAAZu5x556Ur3/8q/mbd30lSfIjz31UTnzhj0/cCgCoaNH3fT91CQAAxvW5t3wu7/+f358kefGnX5yHnfywiRsBABVNOhStVptbynXdcrbZKj1ay1bpUSG7nXzXLZPM9/aq9KiQrdKjtexYb7d4e1V6VMhW6bFO9vqrr8ufPOM/5qBDDsrtN9yeox91bJ73wf8pG4dvzPp/F1Vub8y3W8tW6VEhW6VHhdsbq8fcs1V6VMhW6VEluy6/zBoAYMb237Yvu855f+64+Y6c9pYz8vh/e3L+8cvfyMd+48NTVwMACjIUAQDM2Cdec0W+cdWePO7fnJRHPP0Hc9K/PyU7T35ovvT2L+Rv3/P/TF0PACjGUAQAMFNX/5e/zVV/8Pk85F/tzBNe8aQkyUEHH5TT3nJGDjv68Hz4/A/luquvm7glAFCJoQgAYIY2v3ZD/vylH8yhRx2W0958Rg7a+Pb/27d82DKnXnh6bt+8PZf+/KXZf/v+CZsCAJVsTF0AAIADb3n8UXnRX//aPf79H3r2CXnJnvO3/ctCAYB58xNFAAAAACQxFAEAAAAwMBQBAAAAkMRQBAAAAMBg0fd9P3UJAAAAAKY36b/1bKv/ho3t/Ns4WstW6dFatkqPCtnt5LtumWS+t1elR4VslR6tZcd6u8Xbq9KjQrZKjwrZKj1au70x324tW6VHhWyVHhVub6wec89W6VEhW6VHley6fPUMAAAAgCSGIgAAAAAGhiIAAAAAkhiKAAAAABgYigAAAABIYigCAAAAYGAoAgAAACCJoQgAAACAgaEIAAAAgCTJou/7fuoSAAAAAExvY8oPX602t5TruuVss1V6tJat0qNCdjv5rlsmme/tVelRIVulR2vZsd5u8faq9KiQrdKjQrZKj9Zub8y3W8tW6VEhW6VHhdsbq8fcs1V6VMhW6VEluy5fPQMAAAAgiaEIAAD+//buP1izgj4T/PM2VyJTvG3SeJpGJU5AUYNodILATLIbMWBqWMeOEauySY2aRE2yJQu4GJWs5ZJpTKWQ1VQqJT2YGHbKGhOyNlW0SGsDNv7+HayKmwxlUdOwdPeL8uOCSEN79o+c9MwmktzLvafP9xw+n/+6++Gc559vveXjfbsBgI6hCAAAAIAkhiIAAAAAOoYiAAAAAJIYigAAAADoGIoAAAAASGIoAgAAAKBjKAIAAAAgSTJr27YdugQAAAAAw1sa8uWLxfKKck0zn2y2So+xZav0qJBdTb5p5kmme3tVelTIVukxtmxfzx7j7VXpUSFbpUeFbJUeY7u9Pp89tmyVHhWyVXpUuL2+ekw9W6VHhWyVHlWya+WrZwAAAAAkMRQBAAAA0DEUAQAAAJDEUAQAAABAx1AEAAAAQBJDEQAAAAAdQxEAAAAASQxFAAAAAHQMRQAAAAAkSWZt27ZDlwAAAABgeEtDvnyxWF5Rrmnmk81W6TG2bJUeFbKryTfNPMl0b69KjwrZKj3Glu3r2WO8vSo9KmSr9KiQrdJjbLfX57PHlq3So0K2So8Kt9dXj6lnq/SokK3So0p2rXz1DAAAAIAkhiIAAAAAOoYiAAAAAJIYigAAAADoGIoAAAAASGIoAgAAAKBjKAIAAAAgiaEIAAAAgI6hCAAAAIAkyaxt23boEgAAAAAMb2nIly8WyyvKNc18stkqPcaWrdKjQnY1+aaZJ5nu7VXpUSFbpcfYsn09e4y3V6VHhWyVHhWyVXqM7fb6fPbYslV6VMhW6VHh9vrqMfVslR4VslV6VMmula+eAQAAAJDEUAQAAABAx1AEAAAAQBJDEQAAAAAdQxEAAAAASQxFAAAAAHQMRQAAAAAkMRQBAAAA0DEUAQAAAJAkmbVt2w5dAgAAAIDhLQ358sVieUW5pplPNlulx9iyVXpUyK4m3zTzJNO9vSo9KmSr9Bhbtq9nj/H2qvSokK3So0K2So+x3V6fzx5btkqPCtkqPSrcXl89pp6t0qNCtkqPKtm18tUzAAAAAJIYigAAAADoGIoAAAAASGIoAgAAAKBjKAIAAAAgiaEIAAAAgI6hCAAAAIAkhiIAAAAAOoYiAAAAAJIks7Zt26FLAAAAADC8pSFfvlgsryjXNPPJZqv0GFu2So8K2dXkm2aeZLq3V6VHhWyVHmPL9vXsMd5elR4VslV6VMhW6TG22+vz2WPLVulRIVulR4Xb66vH1LNVelTIVulRJbtWvnoGAAAAQBJDEQAAAAAdQxEAAAAASQxFAAAAAHQMRQAAAAAkWYeh6NChQ9m6dWve8pa3JEn27t2b888/P+eee24uvPDCHDx4cM0lAQAAAOjfmoeia665JieffPLhX19xxRV5wxvekF27dmXjxo259tpr1/oKAAAAAI6AWdu27RP9j/ft25ff+Z3fyW/+5m/mwx/+cD74wQ/mzDPPzGc/+9ksLS3l61//ev7oj/4oH/rQh9azMwAAAAA9WNNPFF1++eW55JJLsmHD3z3m3nvvzcaNG7O0tJQk2bJlS/bv37/2lgAAAAD0bumJ/oc333xzNm3alBe+8IX54he/+Li52Wz2uH+2WCyv6F1NM59stkqPsWWr9KiQXU2+aeZJpnt7VXpUyFbpMbZsX88e4+1V6VEhW6VHhWyVHmO7vT6fPbZslR4VslV6VLi9vnpMPVulR4VslR5Vsmv1hIeir33ta7npppuyZ8+ePPLII3nwwQezbdu2PPDAA3nssceytLSUffv2ZfPmzWsuCQAAAED/nvBXz972trdlz549uemmm3LllVfmzDPPzPve976cccYZufHGG5MkH/vYx3L22WevW1kAAAAA+rPmf/XsH7rkkkvyp3/6pznnnHNy33335fzzz1/vVwAAAADQgyf81bP/3hlnnJEzzjgjSXLiiSfm2muvXY/HAgAAAHAErftPFAEAAAAwToYiAAAAAJIYigAAAADoGIoAAAAASGIoAgAAAKAza9u2HboEAAAAAMNbGvLli8XyinJNM59stkqPsWWr9KiQXU2+aeZJpnt7VXpUyFbpMbZsX88e4+1V6VEhW6VHhWyVHmO7vT6fPbZslR4VslV6VLi9vnpMPVulR4VslR5Vsmvlq2cAAAAAJDEUAQAAANAxFAEAAACQxFAEAAAAQMdQBAAAAEASQxEAAAAAHUMRAAAAAEkMRQAAAAB0DEUAAAAAJElmbdu2Q5cAAAAAYHhLQ758sVheUa5p5pPNVukxtmyVHhWyq8k3zTzJdG+vSo8K2So9xpbt69ljvL0qPSpkq/SokK3SY2y31+ezx5at0qNCtkqPCrfXV4+pZ6v0qJCt0qNKdq189QwAAACAJIYiAAAAADqGIgAAAACSGIoAAAAA6BiKAAAAAEhiKAIAAACgYygCAAAAIImhCAAAAICOoQgAAACAJMmsbdt26BIAAAAADG9pyJcvFssryjXNfLLZKj3Glq3So0J2NfmmmSeZ7u1V6VEhW6XH2LJ9PXuMt1elR4VslR4VslV6jO32+nz22LJVelTIVulR4fb66jH1bJUeFbJVelTJrpWvngEAAACQxFAEAAAAQMdQBAAAAEASQxEAAAAAHUMRAAAAAEkMRQAAAAB0DEUAAAAAJDEUAQAAANAxFAEAAACQJJm1bdsOXQIAAACA4S0N+fLFYnlFuaaZTzZbpcfYslV6VMiuJt808yTTvb0qPSpkq/QYW7avZ4/x9qr0qJCt0qNCtkqPsd1en88eW7ZKjwrZKj0q3F5fPaaerdKjQrZKjyrZtfLVMwAAAACSGIoAAAAA6BiKAAAAAEhiKAIAAACgYygCAAAAIImhCAAAAICOoQgAAACAJIYiAAAAADqGIgAAAACSJLO2bduhSwAAAAAwvKUhX75YLK8o1zTzyWar9BhbtkqPCtnV5JtmnmS6t1elR4VslR5jy/b17DHeXpUeFbJVelTIVukxttvr89ljy1bpUSFbpUeF2+urx9SzVXpUyFbpUSW7Vr56BgAAAEASQxEAAAAAHUMRAAAAAEkMRQAAAAB0DEUAAAAAJDEUAQAAANAxFAEAAACQxFAEAAAAQMdQBAAAAECSZNa2bTt0CQAAAACGtzTkyxeL5RXlmmY+2WyVHmPLVulRIbuafNPMk0z39qr0qJCt0mNs2b6ePcbbq9KjQrZKjwrZKj3Gdnt9Pnts2So9KmSr9Khwe331mHq2So8K2So9qmTXylfPAAAAAEhiKAIAAACgYygCAAAAIImhCAAAAICOoQgAAACAJIYiAAAAADqGIgAAAACSGIoAAAAA6BiKAAAAAEiSzNq2bYcuAQAAAMDwloZ8+WKxvKJc08wnm63SY2zZKj0qZFeTb5p5kuneXpUeFbJVeowt29ezx3h7VXpUyFbpUSFbpcfYbq/PZ48tW6VHhWyVHhVur68eU89W6VEhW6VHlexa+eoZAAAAAEkMRQAAAAB0DEUAAAAAJDEUAQAAANAxFAEAAACQxFAEAAAAQMdQBAAAAEASQxEAAAAAnVnbtu3QJQAAAAAY3tKQL18slleUa5r5ZLNVeowtW6VHhexq8k0zTzLd26vSo0K2So+xZft69hhvr0qPCtkqPSpkq/QY2+31+eyxZav0qJCt0qPC7fXVY+rZKj0qZKv0qJJdK189AwAAACCJoQgAAACAjqEIAAAAgCSGIgAAAAA6hiIAAAAAkhiKAAAAAOgYigAAAABIYigCAAAAoGMoAgAAACBJMmvbth26BAAAAADDWxry5YvF8opyTTOfbLZKj7Flq/SokF1NvmnmSaZ7e1V6VMhW6TG2bF/PHuPtVelRIVulR4VslR5ju70+nz22bJUeFbJVelS4vb56TD1bpUeFbJUeVbJr5atnAAAAACQxFAEAAADQMRQBAAAAkMRQBAAAAEDHUAQAAABAEkMRAAAAAB1DEQAAAABJDEUAAAAAdAxFAAAAACRJZm3btkOXAAAAAGB4S0O+fLFYXlGuaeaTzVbpMbZslR4VsqvJN808yXRvr0qPCtkqPcaW7evZY7y9Kj0qZKv0qJCt0mNst9fns8eWrdKjQrZKjwq311ePqWer9KiQrdKjSnatfPUMAAAAgCSGIgAAAAA6hiIAAAAAkhiKAAAAAOgYigAAAABIYigCAAAAoGMoAgAAACCJoQgAAACAjqEIAAAAgCTJrG3bdugSAAAAAAxvaciXLxbLK8o1zXyy2So9xpat0qNCdjX5ppknme7tVelRIVulx9iyfT17jLdXpUeFbJUeFbJVeozt9vp89tiyVXpUyFbpUeH2+uox9WyVHhWyVXpUya6Vr54BAAAAkMRQBAAAAEDHUAQAAABAEkMRAAAAAB1DEQAAAABJDEUAAAAAdAxFAAAAACQxFAEAAADQMRQBAAAAkCSZtW3bDl0CAAAAgOEtDfnyxWJ5RbmmmU82W6XH2LJVelTIribfNPMk0729Kj0qZKv0GFu2r2eP8faq9KiQrdKjQrZKj7HdXp/PHlu2So8K2So9KtxeXz2mnq3So0K2So8q2bXy1TMAAAAAkhiKAAAAAOgYigAAAABIYigCAAAAoGMoAgAAACCJoQgAAACAjqEIAAAAgCSGIgAAAAA6hiIAAAAAkiSztm3boUsAAAAAMLylIV++WCyvKNc088lmq/QYW7ZKjwrZ1eSbZp5kurdXpUeFbJUeY8v29ewx3l6VHhWyVXpUyFbpMbbb6/PZY8tW6VEhW6VHhdvrq8fUs1V6VMhW6VElu1a+egYAAABAEkMRAAAAAB1DEQAAAABJDEUAAAAcATf8++vyx5uvzG1Xf/0f/dlN//tN+ePNV+bmi3YN0Az47xmKAAAA6N3LP3Bujn3WPJ//P/Zk8c0Dh3//zj3/NZ+5/DP5secdl5/Z9vIBGwKJoQgAAIAj4Kk/dkzO+eB5aQ+12fWm6/PogwfzvcX38qnfviFH/chReeV/PC9P+RdPGbomPOktDV0AAACAJ4cTXvaMvOwd/zpf+A+fyS3/26fy8HcezvcOPJRX/cdXZdPznz50PSCGIgAAAI6gl7z19Nz12b35L//3/5Mkee5rnpeX/sZLs1gsD9wMSHz1DAAAgCNoNpvlpH/73MO/ftGbXzpgG+AfmrVt2w5dAgAAgCeH7/yX72T7S7dnw1M25JH7H0lzapM3felNWXqqL7xABYNe4kp/tLBp5pPNVukxtmyVHhWyq8k3zTzJdG+vSo8K2So9xpbt69ljvL0qPSpkq/SokK3SY2y31+ezx5at0qNCtkqPI317hx55LH/5S3+egw8dzP/0n1+Tu79wV776f34xn7jwEznj9/7HwTqPMVulR4VslR5VsmtlsgUAAOCI+Nx79uSebx7IS956en785f8yz/offjx3fXZvvnrVV3PcT2/Jc179vKErwpOev6MIAACA3n3747fnmx/6Ro7/V1tyxjv/TZJkw1Ebcs5V5+WYTcfklos/mfvvuG/gloChCAAAgF4t3/lAbr7wxhy98UdyzgfPy4al/5y28GsAABpbSURBVPY/RefPnOff/cm/y8Hlg/nkW3bm0MFDAzYFfPUMAACAXs2ftTG//rf/y+P++fNf/fz89oGLj2Aj4PH4iSIAAAAAkhiKAAAAAOgYigAAAABIYigCAAAAoGMoAgAAACBJMmvbth26BAAAAADDWxry5YvF8opyTTOfbLZKj7Flq/SokF1NvmnmSaZ7e1V6VMhW6TG2bF/PHuPtVelRIVulR4VslR5ju70+nz22bJUeFbJVelS4vb56TD1bpUeFbJUeVbJr5atnAAAAACQxFAEAAADQMRQBAAAAkMRQBAAAAEDHUAQAAABAEkMRAAAAAB1DEQAAAABJDEUAAAAAdGZt27ZDlwAAAABgeEtDvnyxWF5Rrmnmk81W6TG2bJUeFbKryTfNPMl0b69KjwrZKj3Glu3r2WO8vSo9KmSr9KiQrdJjbLfX57PHlq3So0K2So8Kt9dXj6lnq/SokK3So0p2rXz1DAAAAIAkhiIAAAAAOoYiAAAAAJIYigAAAADoGIoAAAAASGIoAgAAAKBjKAIAAAAgiaEIAAAAgI6hCAAAAIAkyaxt23boEgAAAAAMb2nIly8WyyvKNc18stkqPcaWrdKjQnY1+aaZJ5nu7VXpUSFbpcfYsn09e4y3V6VHhWyVHhWyVXqM7fb6fPbYslV6VMhW6VHh9vrqMfVslR4VslV6VMmula+eAQAAAJDEUAQAAABAx1AEAAAAQBJDEQAAAAAdQxEAAAAASQxFAAAAAHQMRQAAAAAkMRQBAAAA0DEUAQAAAJAkmbVt2w5dAgAAAIDhLQ358sVieUW5pplPNlulx9iyVXpUyK4m3zTzJNO9vSo9KmSr9Bhbtq9nj/H2qvSokK3So0K2So+x3V6fzx5btkqPCtkqPSrcXl89pp6t0qNCtkqPKtm18tUzAAAAAJIYigAAAADoGIoAAAAASGIoAgAAAKBjKAIAAAAgiaEIAAAAgI6hCAAAAIAkhiIAAAAAOoYiAAAAYLJu+PfX5Y83X5nbrv76P/qzL/7+Z/PHm6/MzRftGqBZTbO2bduhSwAAAAD04eHvPpyrXnJVHtz/YH7987+eE15yQpLk27u/nf907n/K05//9Lzpy2/KU/7FUwZuWsOgQ9FisbyiXNPMJ5ut0mNs2So9KmRXk2+aeZLp3l6VHhWyVXqMLdvXs8d4e1V6VMhW6VEhW6XH2G6vz2ePLVulR4VslR4Vbq+vHlPPVulRIbuS/N1f+n9z3dY/z/zHN+a3/+q3cvd/vTd//vL/KweXH8lrb/yfs+n5Tz+infvMrtXSmp8AAAAAUNgJL3tGXvaOf50v/IfP5Pq3XJ9773og3zvwUH7uynMedyR6sjIUAQAAAJP3kreenrs+uzff/Mg3kyTPfc3z8pO/etrArerxl1kDAAAAkzebzXLSv33u4V+/6M0vHbBNXYYiAAAAYPLu+/a9+dx7Pp2n/thTM9swy80XfTKPff+xoWuVYygCAAAAJu3QI49l15t25tHvPZrX/ufX5qX/68vy3W/dk8/87i1DVyvnCf8dRXfffXfe/va355577smGDRvyute9Lq9//etz33335aKLLspdd92VZz7zmXn/+9+fpz3taevZGQAAAGDFPveePbnnmwfykreenpPPPTnHvujpueuze/PX19yWZ/3siXnOq583dMUynvBPFB111FF5xzvekRtuuCEf/ehH85GPfCS33357tm/fnrPOOiu7du3KWWedle3bt69nXwAAAIAV+/bHb883P/SNHP+vtuSMd/6bJMmGozbknKvOy4/82FNzy8WfzP133Ddwyzqe8FC0efPmnHrqqUmSY489NieddFL279+f3bt3Z+vWrUmSrVu35lOf+tT6NAUAAABYheU7H8jNF96Yozf+SM754HnZsPTfZpD5M+c5+/3n5uDywXzyLTtz6OChAZvWMWvbtl3rQ+6888786q/+aq6//vr83M/9XL7yla8c/rPTTz89X/7yl9f6CgAAAAB6tua/zPqhhx7KBRdckHe961059thj16MTAAAAAAN4wn+ZdZI8+uijueCCC/KqV70q5557bpLkuOOOy4EDB7J58+YcOHAgmzZtetz/frFYXtF7mmY+2WyVHmPLVulRIbuafNPMk0z39qr0qJCt0mNs2b6ePcbbq9KjQrZKjwrZKj3Gdnt9Pnts2So9KmSr9Khwe331mHq2So8K2So9qmTX6gn/RFHbtrn00ktz0kkn5Y1vfOPh3z/77LOzY8eOJMmOHTvyile8Ys0lAQAAAOjfE/6Joq9+9au57rrrcsopp+TVr351kuTiiy/Om9/85lx44YW59tprc8IJJ+QDH/jAupUFAAAAoD9PeCj66Z/+6fzN3/zND/2zP/uzP3vChQAAAAAYxpr/MmsAAAAApsFQBAAAAEASQxEAAAAAHUMRAAAAAEkMRQAAAAB0DEUAAAAAJElmbdu2Q5cAAAAAYHhLQ758sVheUa5p5pPNVukxtmyVHhWyq8k3zTzJdG+vSo8K2So9xpbt69ljvL0qPSpkq/SokK3SY2y31+ezx5at0qNCtkqPCrfXV4+pZ6v0qJCt0qNKdq189QwAAACAJIYiAAAAADqGIgAAAACSGIoAAAAA6BiKAAAAAEhiKAIAAACgYygCAAAAIImhCAAAAICOoQgAAACAJMmsbdt26BIAAAAADG9pyJcvFssryjXNfLLZKj3Glq3So0J2NfmmmSeZ7u1V6VEhW6XH2LJ9PXuMt1elR4VslR4VslV6jO32+nz22LJVelTIVulR4fb66jH1bJUeFbJVelTJrpWvngEAAACQxFAEAAAAQMdQBAAAAEASQxEAAAAAHUMRAAAAAEkMRQAAAAB0DEUAAAAAJDEUAQAAANCZtW3bDl0CAAAAgOEtDfnyxWJ5RbmmmU82W6XH2LJVelTIribfNPMk0729Kj0qZKv0GFu2r2eP8faq9KiQrdKjQrZKj7HdXp/PHlu2So8K2So9KtxeXz2mnq3So0K2So8q2bXy1TMAAAAAkhiKAAAAAOgYigAAAABIYigCAAAAoGMoAgAAACCJoQgAAACAjqEIAAAAgCSGIgAAAAA6hiIAAAAAkiSztm3boUsAAAAAMLylIV++WCyvKNc088lmq/QYW7ZKjwrZ1eSbZp5kurdXpUeFbJUeY8v29ewx3l6VHhWyVXpUyFbpMbbb6/PZY8tW6VEhW6VHhdvrq8fUs1V6VMhW6VElu1a+egYAAABAEkMRAAAAAB1DEQAAAABJDEUAAAAAdAxFAAAAACQxFAEAAADQMRQBAAAAkMRQBAAAAEDHUAQAAABAkmTWtm07dAkAAAAAhrc05MsXi+UV5ZpmPtlslR5jy1bpUSG7mnzTzJNM9/aq9KiQrdJjbNm+nj3G26vSo0K2So8K2So9xnZ7fT57bNkqPSpkq/SocHt99Zh6tkqPCtkqPapk18pXzwAAAABIYigCAAAAoGMoAgAAACCJoQgAAACAjqEIAAAAgCSGIgAAAAA6hiIAAAAAkhiKAAAAAOgYigAAAABIkszatm2HLgEAAADA8JaGfPlisbyiXNPMJ5ut0mNs2So9KmRXk2+aeZLp3l6VHhWyVXqMLdvXs8d4e1V6VMhW6VEhW6XH2G6vz2ePLVulR4VslR4Vbq+vHlPPVulRIVulR5XsWvnqGQAAAABJDEUAAAAAdAxFAAAAACQxFAEAAADQMRQBAAAAkMRQBAAAAEDHUAQAAABAEkMRAAAAAB1DEQAAAABJklnbtu3QJQAAAAAY3tKQL18slleUa5r5ZLNVeowtW6VHhexq8k0zTzLd26vSo0K2So+xZft69hhvr0qPCtkqPSpkq/QY2+31+eyxZav0qJCt0qPC7fXVY+rZKj0qZKv0qJJdK189AwAAACCJoQgAAACAjqEIAAAAgCSGIgAAAAA6hiIAAAAAkhiKAAAAAOgYigAAAABIYigCAAAAoGMoAgAAACBJMmvbth26BAAAAADDWxry5YvF8opyTTOfbLZKj7Flq/SokF1NvmnmSaZ7e1V6VMhW6TG2bF/PHuPtVelRIVulR4VslR5ju70+nz22bJUeFbJVelS4vb56TD1bpUeFbJUeVbJr5atnAAAAACQxFAEAAADQMRQBAAAAkMRQBAAAAEDHUAQAAABAEkMRAAAAAB1DEQAAAABJDEUAAAAAdAxFAAAAACRJZm3btkOXAAAAAGB4S0O+fLFYXlGuaeaTzVbpMbZslR4VsqvJN808yXRvr0qPCtkqPcaW7evZY7y9Kj0qZKv0qJCt0mNst9fns8eWrdKjQrZKjwq311ePqWer9KiQrdKjSnatfPUMAAAAgCSGIgAAAAA6hiIAAAAAkhiKAAAAAOgYigAAAABIYigCAAAAoGMoAgAAACCJoQgAAACAjqEIAAAAgCTJrG3bdugSAAAAAAxvaciXLxbLK8o1zXyy2So9xpat0qNCdjX5ppknme7tVelRIVulx9iyfT17jLdXpUeFbJUeFbJVeozt9vp89tiyVXpUyFbpUeH2+uox9WyVHhWyVXpUya6Vr54BAAAAkMRQBAAAAEDHUAQAAABAEkMRAAAAAB1DEQAAAABJDEUAAAAAdAxFAAAAACQxFAEAAADQmbVt2w5dAgAAAIDhLQ358sVieUW5pplPNlulx9iyVXpUyK4m3zTzJNO9vSo9KmSr9Bhbtq9nj/H2qvSokK3So0K2So+x3V6fzx5btkqPCtkqPSrcXl89pp6t0qNCtkqPKtm18tUzAAAAAJIYigAAAADoGIoAAAAASGIoAgAAAKBjKAIAAAAgiaEIAAAAgI6hCAAAAIAkhiIAAAAAOoYiAAAAAJIks7Zt26FLAAAAADC8pSFfvlgsryjXNPPJZqv0GFu2So8K2dXkm2aeZLq3V6VHhWyVHmPL9vXsMd5elR4VslV6VMhW6TG22+vz2WPLVulRIVulR4Xb66vH1LNVelTIVulRJbtWvnoGAAAAQBJDEQAAAAAdQxEAAAAASQxFAAAAAHQMRQAAAAAkMRQBAAAA0DEUAQAAAJDEUAQAAABAx1AEAAAAQJJk1rZtO3QJAAAAAIa3NOTLF4vlFeWaZj7ZbJUeY8tW6VEhu5p808yTTPf2qvSokK3SY2zZvp49xtur0qNCtkqPCtkqPcZ2e30+e2zZKj0qZKv0qHB7ffWYerZKjwrZKj2qZNfKV88AAAAASGIoAgAAAKBjKAIAAAAgiaEIAAAAgI6hCAAAAIAkhiIAAAAAOoYiAAAAAJIYigAAAADoGIoAAAAASJLM2rZthy4BAAAAwPCWhnz5YrG8olzTzCebrdJjbNkqPSpkV5NvmnmS6d5elR4VslV6jC3b17PHeHtVelTIVulRIVulx9hur89njy1bpUeFbJUeFW6vrx5Tz1bpUSFbpUeV7Fr56hkAAAAASQxFAAAAAHQMRQAAAAAkMRQBAAAA0DEUAQAAAJDEUAQAAABAx1AEAAAAQBJDEQAAAAAdQxEAAAAASZJZ27bt0CUAAAAAGN7SkC9fLJZXlGua+WSzVXqMLVulR4XsavJNM08y3dur0qNCtkqPsWX7evYYb69KjwrZKj0qZKv0GNvt9fnssWWr9KiQrdKjwu311WPq2So9KmSr9KiSXStfPQMAAAAgiaEIAAAAgI6hCAAAAIAkhiIAAAAAOoYiAAAAAJIYigAAAADoGIoAAAAASGIoAgAAAKBjKAIAAAAgSTJr27YdugQAAAAAw1sa8uWLxfKKck0zn2y2So+xZav0qJBdTb5p5kmme3tVelTIVukxtmxfzx7j7VXpUSFbpUeFbJUeY7u9Pp89tmyVHhWyVXpUuL2+ekw9W6VHhWyVHlWya+WrZwAAAAAkMRQBAAAA0DEUAQAAAJDEUAQAAABAx1AEAAAAQBJDEQAAAAAdQxEAAAAASQxFAAAAAHQMRQAAAAAkSWZt27ZDlwAAAABgeEtDvnyxWF5Rrmnmk81W6TG2bJUeFbKryTfNPMl0b69KjwrZKj3Glu3r2WO8vSo9KmSr9KiQrdJjbLfX57PHlq3So0K2So8Kt9dXj6lnq/SokK3So0p2rXz1DAAAAIAkhiIAAAAAOoYiAAAAAJIYigAAAADoGIoAAAAASGIoAgAAAKBjKAIAAAAgiaEIAAAAgI6hCAAAAIAkyaxt23boEgAAAAAMb2nIly8WyyvKNc18stkqPcaWrdKjQnY1+aaZJ5nu7VXpUSFbpcfYsn09e4y3V6VHhWyVHhWyVXqM7fb6fPbYslV6VMhW6VHh9vrqMfVslR4VslV6VMmula+eAQAAAJDEUAQAAABAx1AEAAAAQBJDEQAAAAAdQxEAAAAASQxFAAAAAHQMRQAAAAAkMRQBAAAA0Jm1bdsOXQIAAACA4S0N+fLFYnlFuaaZTzZbpcfYslV6VMiuJt808yTTvb0qPSpkq/QYW7avZ4/x9qr0qJCt0qNCtkqPsd1en88eW7ZKjwrZKj0q3F5fPaaerdKjQrZKjyrZtfLVMwAAAACSGIoAAAAA6BiKAAAAAEhiKAIAAACgYygCAAAAIElPQ9GePXvyyle+Muecc062b9/exysAAAAAWGfrPhQdOnQol112Wa6++urs3Lkz119/fW6//fb1fg0AAAAA62zdh6Lbbrstz372s3PiiSfm6KOPznnnnZfdu3ev92sAAAAAWGeztm3b9XzgJz7xidx6663Ztm1bkmTHjh257bbb8u53v3s9XwMAAADAOlv3nyj6YbvTbDZb79cAAAAAsM7WfSjasmVL9u3bd/jX+/fvz+bNm9f7NQAAAACss3Ufik477bTccccd2bt3bw4ePJidO3fm7LPPXu/XAAAAALDOltb9gUtLefe7353f+I3fyKFDh/JLv/RLee5zn7verwEAAABgna37X2b9z9mzZ0+2bduWH/zgBzn//PPz5je/+Ui+Hp407r777rz97W/PPffckw0bNuR1r3tdXv/61+e+++7LRRddlLvuuivPfOYz8/73vz9Pe9rThq4Lk/P3/2fJ8ccfn6uuuip79+7NxRdfnPvvvz8/+ZM/mT/4gz/I0UcfPXRNmJwHHnggv/u7v5u//du/zWw2y+WXX56f+Imf8NkHPfvwhz+cv/iLv8hsNsspp5yS9773vTlw4IDPPlhn73znO3PLLbfkuOOOy/XXX58kj/u/8dq2zbZt2/LpT386T33qU/P7v//7OfXUU//Zd6z7V8/+KYcOHcpll12Wq6++Ojt37sz111+f22+//UhWgCeNo446Ku94xztyww035KMf/Wg+8pGP5Pbbb8/27dtz1llnZdeuXTnrrLOyffv2oavCJF1zzTU5+eSTD//6iiuuyBve8Ibs2rUrGzduzLXXXjtgO5iubdu25Wd/9mfziU98Itddd11OPvlkn33Qs/379+eaa67JX/7lX+b666/PoUOHsnPnTp990IPXvOY1ufrqq/9/v/d4n3N79uzJHXfckV27duX3fu/38p73vGdF7ziiQ9Ftt92WZz/72TnxxBNz9NFH57zzzsvu3buPZAV40ti8efPhtfjYY4/NSSedlP3792f37t3ZunVrkmTr1q351Kc+NWRNmKR9+/bllltuyWtf+9okf/cvgn7hC1/IK1/5yiTJL/7iL/r8gx48+OCD+fKXv3z49o4++uhs3LjRZx8cAYcOHcr3v//9PPbYY/n+97+fpml89kEPTj/99H/0U7GP9zn3978/m83yUz/1U3nggQdy4MCBf/YdR3Qo2r9/f7Zs2XL418cff3z2799/JCvAk9Kdd96Zb33rW3nxi1+c73znO4f/JcLNmzfnu9/97sDtYHouv/zyXHLJJdmw4e8+Zu+9995s3LgxS0t/91cDbtmyxecf9GDv3r3ZtGlT3vnOd2br1q259NJL873vfc9nH/Ts+OOPz6/92q/l5S9/eX7mZ34mxx57bE499VSffXCEPN7n3D/cYFZ6h0d0KPphfx3SbDY7khXgSeehhx7KBRdckHe961059thjh64Dk3fzzTdn06ZNeeELX/hP5nz+wfp77LHH8td//df55V/+5ezYsSPHHHOMr5nBEXD//fdn9+7d2b17d2699dY8/PDD2bNnzz/K+eyDI+uJbjDr/q+e/VO2bNmSffv2Hf71/v37D69ewPp79NFHc8EFF+RVr3pVzj333CTJcccdlwMHDmTz5s05cOBANm3aNHBLmJavfe1ruemmm7Jnz5488sgjefDBB7Nt27Y88MADeeyxx7K0tJR9+/b5/IMebNmyJVu2bMmLX/ziJMkv/MIvZPv27T77oGef+9zn8qxnPevwbZ177rn5+te/7rMPjpDH+5z7hxvMSu/wiP5E0WmnnZY77rgje/fuzcGDB7Nz586cffbZR7ICPGm0bZtLL700J510Ut74xjce/v2zzz47O3bsSJLs2LEjr3jFK4aqCJP0tre9LXv27MlNN92UK6+8MmeeeWbe97735YwzzsiNN96YJPnYxz7m8w960DRNtmzZkm9/+9tJks9//vM5+eSTffZBz57xjGfkr/7qr/Lwww+nbdt8/vOfz3Oe8xyffXCEPN7n3N//ftu2+cY3vpH5fL6ioWjW/rCfRerRpz/96Vx++eWH/9ng3/qt3zqSr4cnja985Sv5lV/5lZxyyimH/56Uiy++OC960Yty4YUX5u67784JJ5yQD3zgA/nRH/3RgdvCNH3xi1/Mn/zJn+Sqq67K3r17c9FFF+X+++/PC17wglxxxRX+iWDowbe+9a1ceumlefTRR3PiiSfmve99b37wgx/47IOe/eEf/mE+/vGPZ2lpKS94wQuybdu27N+/32cfrLOLL744X/rSl3LvvffmuOOOy1vf+tb8/M///A/9nGvbNpdddlluvfXWHHPMMbn88stz2mmn/bPvOOJDEQAAAAA1HdGvngEAAABQl6EIAAAAgCSGIgAAAAA6hiIAAAAAkhiKAAAAAOgYigAAAABIYigCAAAAoGMoAgAAACBJ8v8BP8YWSbGNzuwAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f11f90cd198>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# import the helper function\n",
    "from helpers import display_world\n",
    "\n",
    "# Display the final world!\n",
    "\n",
    "# define figure size\n",
    "plt.rcParams[\"figure.figsize\"] = (20,20)\n",
    "\n",
    "# check if poses has been created\n",
    "if 'poses' in locals():\n",
    "    # print out the last pose\n",
    "    print('Last pose: ', poses[-1])\n",
    "    # display the last position of the robot *and* the landmark positions\n",
    "    display_world(int(world_size), poses[-1], landmarks)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Question: How far away is your final pose (as estimated by `slam`) compared to the *true* final pose? Why do you think these poses are different?\n",
    "\n",
    "You can find the true value of the final pose in one of the first cells where `make_data` was called. You may also want to look at the true landmark locations and compare them to those that were estimated by `slam`. Ask yourself: what do you think would happen if we moved and sensed more (increased N)? Or if we had lower/higher noise parameters."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "**Answer**: Difference between points generated by `slam` and generated by `make_data` is very very much less near to negligible. But why there is a difference in generated point and given point, this is due to `measurement_noise` and `motion_noise` captured by robot affecting the points value. \n",
    "<br>*What Happen If*:<br>\n",
    "- **N is increased**: Length of points generated from `slam` will increase to the value of N.\n",
    "- **lower/higher noise parameters**: When we lower the noise parameters the differece b/w `slam` generated points and true points will be decrease as these is less noise captured by robot, and if we higher the value of noise parameters `slam` generated points will step toward the randomness. **We can say accuracy of `slam` generated points are inversely propostional to value of noise paremeters**"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Testing\n",
    "\n",
    "To confirm that your slam code works before submitting your project, it is suggested that you run it on some test data and cases. A few such cases have been provided for you, in the cells below. When you are ready, uncomment the test cases in the next cells (there are two test cases, total); your output should be **close-to or exactly** identical to the given results. If there are minor discrepancies it could be a matter of floating point accuracy or in the calculation of the inverse matrix.\n",
    "\n",
    "### Submit your project\n",
    "\n",
    "If you pass these tests, it is a good indication that your project will pass all the specifications in the project rubric. Follow the submission instructions to officially submit!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "\n",
      "Estimated Poses:\n",
      "[50.000, 50.000]\n",
      "[37.973, 33.652]\n",
      "[26.185, 18.155]\n",
      "[13.745, 2.116]\n",
      "[28.097, 16.783]\n",
      "[42.384, 30.902]\n",
      "[55.831, 44.497]\n",
      "[70.857, 59.699]\n",
      "[85.697, 75.543]\n",
      "[74.011, 92.434]\n",
      "[53.544, 96.454]\n",
      "[34.525, 100.080]\n",
      "[48.623, 83.953]\n",
      "[60.197, 68.107]\n",
      "[73.778, 52.935]\n",
      "[87.132, 38.538]\n",
      "[80.303, 20.508]\n",
      "[72.798, 2.945]\n",
      "[55.245, 13.255]\n",
      "[37.416, 22.317]\n",
      "\n",
      "\n",
      "Estimated Landmarks:\n",
      "[82.956, 13.539]\n",
      "[70.495, 74.141]\n",
      "[36.740, 61.281]\n",
      "[18.698, 66.060]\n",
      "[20.635, 16.875]\n"
     ]
    }
   ],
   "source": [
    "# Here is the data and estimated outputs for test case 1\n",
    "\n",
    "test_data1 = [[[[1, 19.457599255548065, 23.8387362100849], [2, -13.195807561967236, 11.708840328458608], [3, -30.0954905279171, 15.387879242505843]], [-12.2607279422326, -15.801093326936487]], [[[2, -0.4659930049620491, 28.088559771215664], [4, -17.866382374890936, -16.384904503932]], [-12.2607279422326, -15.801093326936487]], [[[4, -6.202512900833806, -1.823403210274639]], [-12.2607279422326, -15.801093326936487]], [[[4, 7.412136480918645, 15.388585962142429]], [14.008259661173426, 14.274756084260822]], [[[4, -7.526138813444998, -0.4563942429717849]], [14.008259661173426, 14.274756084260822]], [[[2, -6.299793150150058, 29.047830407717623], [4, -21.93551130411791, -13.21956810989039]], [14.008259661173426, 14.274756084260822]], [[[1, 15.796300959032276, 30.65769689694247], [2, -18.64370821983482, 17.380022987031367]], [14.008259661173426, 14.274756084260822]], [[[1, 0.40311325410337906, 14.169429532679855], [2, -35.069349468466235, 2.4945558982439957]], [14.008259661173426, 14.274756084260822]], [[[1, -16.71340983241936, -2.777000269543834]], [-11.006096015782283, 16.699276945166858]], [[[1, -3.611096830835776, -17.954019226763958]], [-19.693482634035977, 3.488085684573048]], [[[1, 18.398273354362416, -22.705102332550947]], [-19.693482634035977, 3.488085684573048]], [[[2, 2.789312482883833, -39.73720193121324]], [12.849049222879723, -15.326510824972983]], [[[1, 21.26897046581808, -10.121029799040915], [2, -11.917698965880655, -23.17711662602097], [3, -31.81167947898398, -16.7985673023331]], [12.849049222879723, -15.326510824972983]], [[[1, 10.48157743234859, 5.692957082575485], [2, -22.31488473554935, -5.389184118551409], [3, -40.81803984305378, -2.4703329790238118]], [12.849049222879723, -15.326510824972983]], [[[0, 10.591050242096598, -39.2051798967113], [1, -3.5675572049297553, 22.849456408289125], [2, -38.39251065320351, 7.288990306029511]], [12.849049222879723, -15.326510824972983]], [[[0, -3.6225556479370766, -25.58006865235512]], [-7.8874682868419965, -18.379005523261092]], [[[0, 1.9784503557879374, -6.5025974151499]], [-7.8874682868419965, -18.379005523261092]], [[[0, 10.050665232782423, 11.026385307998742]], [-17.82919359778298, 9.062000642947142]], [[[0, 26.526838150174818, -0.22563393232425621], [4, -33.70303936886652, 2.880339841013677]], [-17.82919359778298, 9.062000642947142]]]\n",
    "\n",
    "##  Test Case 1\n",
    "##\n",
    "# Estimated Pose(s):\n",
    "#     [50.000, 50.000]\n",
    "#     [37.858, 33.921]\n",
    "#     [25.905, 18.268]\n",
    "#     [13.524, 2.224]\n",
    "#     [27.912, 16.886]\n",
    "#     [42.250, 30.994]\n",
    "#     [55.992, 44.886]\n",
    "#     [70.749, 59.867]\n",
    "#     [85.371, 75.230]\n",
    "#     [73.831, 92.354]\n",
    "#     [53.406, 96.465]\n",
    "#     [34.370, 100.134]\n",
    "#     [48.346, 83.952]\n",
    "#     [60.494, 68.338]\n",
    "#     [73.648, 53.082]\n",
    "#     [86.733, 38.197]\n",
    "#     [79.983, 20.324]\n",
    "#     [72.515, 2.837]\n",
    "#     [54.993, 13.221]\n",
    "#     [37.164, 22.283]\n",
    "\n",
    "\n",
    "# Estimated Landmarks:\n",
    "#     [82.679, 13.435]\n",
    "#     [70.417, 74.203]\n",
    "#     [36.688, 61.431]\n",
    "#     [18.705, 66.136]\n",
    "#     [20.437, 16.983]\n",
    "\n",
    "\n",
    "### Uncomment the following three lines for test case 1 and compare the output to the values above ###\n",
    "\n",
    "mu_1 = slam(test_data1, 20, 5, 100.0, 2.0, 2.0)\n",
    "poses, landmarks = get_poses_landmarks(mu_1, 20)\n",
    "print_all(poses, landmarks)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "\n",
      "Estimated Poses:\n",
      "[50.000, 50.000]\n",
      "[69.181, 45.665]\n",
      "[87.743, 39.703]\n",
      "[76.270, 56.311]\n",
      "[64.317, 72.176]\n",
      "[52.257, 88.154]\n",
      "[44.059, 69.401]\n",
      "[37.002, 49.918]\n",
      "[30.924, 30.955]\n",
      "[23.508, 11.419]\n",
      "[34.180, 27.133]\n",
      "[44.155, 43.846]\n",
      "[54.806, 60.920]\n",
      "[65.698, 78.546]\n",
      "[77.468, 95.626]\n",
      "[96.802, 98.821]\n",
      "[75.957, 99.971]\n",
      "[70.200, 81.181]\n",
      "[64.054, 61.723]\n",
      "[58.107, 42.628]\n",
      "\n",
      "\n",
      "Estimated Landmarks:\n",
      "[76.779, 42.887]\n",
      "[85.065, 77.438]\n",
      "[13.548, 95.652]\n",
      "[59.449, 39.595]\n",
      "[69.263, 94.240]\n"
     ]
    }
   ],
   "source": [
    "# Here is the data and estimated outputs for test case 2\n",
    "\n",
    "test_data2 = [[[[0, 26.543274387283322, -6.262538160312672], [3, 9.937396825799755, -9.128540360867689]], [18.92765331253674, -6.460955043986683]], [[[0, 7.706544739722961, -3.758467215445748], [1, 17.03954411948937, 31.705489938553438], [3, -11.61731288777497, -6.64964096716416]], [18.92765331253674, -6.460955043986683]], [[[0, -12.35130507136378, 2.585119104239249], [1, -2.563534536165313, 38.22159657838369], [3, -26.961236804740935, -0.4802312626141525]], [-11.167066095509824, 16.592065417497455]], [[[0, 1.4138633151721272, -13.912454837810632], [1, 8.087721200818589, 20.51845934354381], [3, -17.091723454402302, -16.521500551709707], [4, -7.414211721400232, 38.09191602674439]], [-11.167066095509824, 16.592065417497455]], [[[0, 12.886743222179561, -28.703968411636318], [1, 21.660953298391387, 3.4912891084614914], [3, -6.401401414569506, -32.321583037341625], [4, 5.034079343639034, 23.102207946092893]], [-11.167066095509824, 16.592065417497455]], [[[1, 31.126317672358578, -10.036784369535214], [2, -38.70878528420893, 7.4987265861424595], [4, 17.977218575473767, 6.150889254289742]], [-6.595520680493778, -18.88118393939265]], [[[1, 41.82460922922086, 7.847527392202475], [3, 15.711709540417502, -30.34633659912818]], [-6.595520680493778, -18.88118393939265]], [[[0, 40.18454208294434, -6.710999804403755], [3, 23.019508919299156, -10.12110867290604]], [-6.595520680493778, -18.88118393939265]], [[[3, 27.18579315312821, 8.067219022708391]], [-6.595520680493778, -18.88118393939265]], [[], [11.492663265706092, 16.36822198838621]], [[[3, 24.57154567653098, 13.461499960708197]], [11.492663265706092, 16.36822198838621]], [[[0, 31.61945290413707, 0.4272295085799329], [3, 16.97392299158991, -5.274596836133088]], [11.492663265706092, 16.36822198838621]], [[[0, 22.407381798735177, -18.03500068379259], [1, 29.642444125196995, 17.3794951934614], [3, 4.7969752441371645, -21.07505361639969], [4, 14.726069092569372, 32.75999422300078]], [11.492663265706092, 16.36822198838621]], [[[0, 10.705527984670137, -34.589764174299596], [1, 18.58772336795603, -0.20109708164787765], [3, -4.839806195049413, -39.92208742305105], [4, 4.18824810165454, 14.146847823548889]], [11.492663265706092, 16.36822198838621]], [[[1, 5.878492140223764, -19.955352450942357], [4, -7.059505455306587, -0.9740849280550585]], [19.628527845173146, 3.83678180657467]], [[[1, -11.150789592446378, -22.736641053247872], [4, -28.832815721158255, -3.9462962046291388]], [-19.841703647091965, 2.5113335861604362]], [[[1, 8.64427397916182, -20.286336970889053], [4, -5.036917727942285, -6.311739993868336]], [-5.946642674882207, -19.09548221169787]], [[[0, 7.151866679283043, -39.56103232616369], [1, 16.01535401373368, -3.780995345194027], [4, -3.04801331832137, 13.697362774960865]], [-5.946642674882207, -19.09548221169787]], [[[0, 12.872879480504395, -19.707592098123207], [1, 22.236710716903136, 16.331770792606406], [3, -4.841206109583004, -21.24604435851242], [4, 4.27111163223552, 32.25309748614184]], [-5.946642674882207, -19.09548221169787]]] \n",
    "\n",
    "\n",
    "##  Test Case 2\n",
    "##\n",
    "# Estimated Pose(s):\n",
    "#     [50.000, 50.000]\n",
    "#     [69.035, 45.061]\n",
    "#     [87.655, 38.971]\n",
    "#     [76.084, 55.541]\n",
    "#     [64.283, 71.684]\n",
    "#     [52.396, 87.887]\n",
    "#     [44.674, 68.948]\n",
    "#     [37.532, 49.680]\n",
    "#     [31.392, 30.893]\n",
    "#     [24.796, 12.012]\n",
    "#     [33.641, 26.440]\n",
    "#     [43.858, 43.560]\n",
    "#     [54.735, 60.659]\n",
    "#     [65.884, 77.791]\n",
    "#     [77.413, 94.554]\n",
    "#     [96.740, 98.020]\n",
    "#     [76.149, 99.586]\n",
    "#     [70.211, 80.580]\n",
    "#     [64.130, 61.270]\n",
    "#     [58.183, 42.175]\n",
    "\n",
    "\n",
    "# Estimated Landmarks:\n",
    "#     [76.777, 42.415]\n",
    "#     [85.109, 76.850]\n",
    "#     [13.687, 95.386]\n",
    "#     [59.488, 39.149]\n",
    "#     [69.283, 93.654]\n",
    "\n",
    "\n",
    "### Uncomment the following three lines for test case 2 and compare to the values above ###\n",
    "\n",
    "mu_2 = slam(test_data2, 20, 5, 100.0, 2.0, 2.0)\n",
    "poses, landmarks = get_poses_landmarks(mu_2, 20)\n",
    "print_all(poses, landmarks)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "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
}