getting_started_notebooks/intro_tutorials/01_Introduction_to_RAPIDS.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<a id=\"introduction\"></a>\n",
    "## Introduction to RAPIDS\n",
    "#### By Paul Hendricks\n",
    "-------\n",
    "\n",
    "While the world’s data doubles each year, CPU computing has hit a brick wall with the end of Moore’s law. For the same reasons, scientific computing and deep learning has turned to NVIDIA GPU acceleration, data analytics and machine learning where GPU acceleration is ideal. \n",
    "\n",
    "NVIDIA created RAPIDS – an open-source data analytics and machine learning acceleration platform that leverages GPUs to accelerate computations. RAPIDS is based on Python, has Pandas-like and Scikit-Learn-like interfaces, is built on Apache Arrow in-memory data format, and can scale from 1 to multi-GPU to multi-nodes. RAPIDS integrates easily into the world’s most popular data science Python-based workflows. RAPIDS accelerates data science end-to-end – from data prep, to machine learning, to deep learning. And through Arrow, Spark users can easily move data into the RAPIDS platform for acceleration.\n",
    "\n",
    "In this notebook, we will discuss and show at a high level what each of the packages in the RAPIDS are as well as what they do. Subsequent notebooks will dive deeper into the various areas of data science and machine learning and show how you can use RAPIDS to accelerate your workflow in each of these areas.\n",
    "\n",
    "**Table of Contents**\n",
    "\n",
    "* [Introduction to RAPIDS](#introduction)\n",
    "* [Setup](#setup)\n",
    "* [Pandas](#pandas)\n",
    "* [cuDF](#cudf)\n",
    "* [Scikit-Learn](#scikitlearn)\n",
    "* [cuML](#cuml)\n",
    "* [Dask](#dask)\n",
    "* [Dask cuDF](#daskcudf)\n",
    "* [Dask cuML (coming soon!)](#daskcuml)\n",
    "* [Conclusion](#conclusion)\n",
    "\n",
    "Before going any further, let's make sure we have access to `matplotlib`, a popular Python library for visualizing data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import os\n",
    "\n",
    "try:\n",
    "    import matplotlib\n",
    "except ModuleNotFoundError:\n",
    "    os.system('conda install -y matplotlib')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<a id=\"setup\"></a>\n",
    "## Setup\n",
    "\n",
    "This notebook was tested using the following Docker containers:\n",
    "\n",
    "* `rapidsai/rapidsai-nightly:0.8-cuda10.0-devel-ubuntu18.04-gcc7-py3.7` from [DockerHub - rapidsai/rapidsai-nightly](https://hub.docker.com/r/rapidsai/rapidsai-nightly)\n",
    "\n",
    "This notebook was run on the NVIDIA Tesla V100 GPU. Please be aware that your system may be different and you may need to modify the code or install packages to run the below examples. \n",
    "\n",
    "If you think you have found a bug or an error, please file an issue here: https://github.com/rapidsai/notebooks/issues\n",
    "\n",
    "Before we begin, let's check out our hardware setup by running the `nvidia-smi` command."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Tue Jun 11 01:30:51 2019       \n",
      "+-----------------------------------------------------------------------------+\n",
      "| NVIDIA-SMI 418.67       Driver Version: 418.67       CUDA Version: 10.1     |\n",
      "|-------------------------------+----------------------+----------------------+\n",
      "| GPU  Name        Persistence-M| Bus-Id        Disp.A | Volatile Uncorr. ECC |\n",
      "| Fan  Temp  Perf  Pwr:Usage/Cap|         Memory-Usage | GPU-Util  Compute M. |\n",
      "|===============================+======================+======================|\n",
      "|   0  Tesla V100-SXM2...  On   | 00000000:06:00.0 Off |                    0 |\n",
      "| N/A   39C    P0    44W / 300W |      0MiB / 16130MiB |      0%      Default |\n",
      "+-------------------------------+----------------------+----------------------+\n",
      "|   1  Tesla V100-SXM2...  On   | 00000000:07:00.0 Off |                    0 |\n",
      "| N/A   39C    P0    42W / 300W |      0MiB / 16130MiB |      0%      Default |\n",
      "+-------------------------------+----------------------+----------------------+\n",
      "|   2  Tesla V100-SXM2...  On   | 00000000:0A:00.0 Off |                    0 |\n",
      "| N/A   40C    P0    44W / 300W |      0MiB / 16130MiB |      0%      Default |\n",
      "+-------------------------------+----------------------+----------------------+\n",
      "|   3  Tesla V100-SXM2...  On   | 00000000:0B:00.0 Off |                    0 |\n",
      "| N/A   37C    P0    43W / 300W |      0MiB / 16130MiB |      0%      Default |\n",
      "+-------------------------------+----------------------+----------------------+\n",
      "|   4  Tesla V100-SXM2...  On   | 00000000:85:00.0 Off |                    0 |\n",
      "| N/A   37C    P0    42W / 300W |      0MiB / 16130MiB |      0%      Default |\n",
      "+-------------------------------+----------------------+----------------------+\n",
      "|   5  Tesla V100-SXM2...  On   | 00000000:86:00.0 Off |                    0 |\n",
      "| N/A   38C    P0    43W / 300W |      0MiB / 16130MiB |      0%      Default |\n",
      "+-------------------------------+----------------------+----------------------+\n",
      "|   6  Tesla V100-SXM2...  On   | 00000000:89:00.0 Off |                    0 |\n",
      "| N/A   40C    P0    43W / 300W |      0MiB / 16130MiB |      0%      Default |\n",
      "+-------------------------------+----------------------+----------------------+\n",
      "|   7  Tesla V100-SXM2...  On   | 00000000:8A:00.0 Off |                    0 |\n",
      "| N/A   37C    P0    42W / 300W |      0MiB / 16130MiB |      0%      Default |\n",
      "+-------------------------------+----------------------+----------------------+\n",
      "                                                                               \n",
      "+-----------------------------------------------------------------------------+\n",
      "| Processes:                                                       GPU Memory |\n",
      "|  GPU       PID   Type   Process name                             Usage      |\n",
      "|=============================================================================|\n",
      "|  No running processes found                                                 |\n",
      "+-----------------------------------------------------------------------------+\n"
     ]
    }
   ],
   "source": [
    "!nvidia-smi"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next, let's see what CUDA version we have:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "nvcc: NVIDIA (R) Cuda compiler driver\r\n",
      "Copyright (c) 2005-2018 NVIDIA Corporation\r\n",
      "Built on Sat_Aug_25_21:08:01_CDT_2018\r\n",
      "Cuda compilation tools, release 10.0, V10.0.130\r\n"
     ]
    }
   ],
   "source": [
    "!nvcc --version"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next, let's load some helper functions from `matplotlib` and configure the Jupyter Notebook for visualization."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "from matplotlib.colors import ListedColormap\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<a id=\"pandas\"></a>\n",
    "## Pandas\n",
    "\n",
    "Data scientists typically work with two types of data: unstructured and structured. Unstructured data often comes in the form of text, images, or videos. Structured data - as the name suggests - comes in a structured form, often represented by a table or CSV. We'll focus the majority of these tutorials on working with these types of data.\n",
    "\n",
    "There exist many tools in the Python ecosystem for working with structured, tabular data but few are as widely used as Pandas. Pandas represents data in a table and allows a data scientist to manipulate the data to perform a number of useful operations such as filtering, transforming, aggregating, merging, visualizing and many more. \n",
    "\n",
    "For more information on Pandas, check out the excellent documentation: http://pandas.pydata.org/pandas-docs/stable/\n",
    "\n",
    "Below we show how to create a Pandas DataFrame, an internal object for representing tabular data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Pandas Version: 0.23.4\n",
      "   key  value\n",
      "0    0   10.0\n",
      "1    0   11.0\n",
      "2    2   12.0\n",
      "3    2   13.0\n",
      "4    3   14.0\n"
     ]
    }
   ],
   "source": [
    "import pandas as pd; print('Pandas Version:', pd.__version__)\n",
    "\n",
    "\n",
    "# here we create a Pandas DataFrame with\n",
    "# two columns named \"key\" and \"value\"\n",
    "df = pd.DataFrame()\n",
    "df['key'] = [0, 0, 2, 2, 3]\n",
    "df['value'] = [float(i + 10) for i in range(5)]\n",
    "print(df)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can perform many operations on this data. For example, let's say we wanted to sum all values in the in the `value` column. We could accomplish this using the following syntax:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "60.0\n"
     ]
    }
   ],
   "source": [
    "aggregation = df['value'].sum()\n",
    "print(aggregation)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<a id=\"cudf\"></a>\n",
    "## cuDF\n",
    "\n",
    "Pandas is fantastic for working with small datasets that fit into your system's memory. However, datasets are growing larger and data scientists are working with increasingly complex workloads - the need for accelerated compute arises.\n",
    "\n",
    "cuDF is a package within the RAPIDS ecosystem that allows data scientists to easily migrate their existing Pandas workflows from CPU to GPU, where computations can leverage the immense parallelization that GPUs provide.\n",
    "\n",
    "Below, we show how to create a cuDF DataFrame."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "cuDF Version: 0.7.2+0.g3ebd286.dirty\n",
      "   key  value\n",
      "0    0   10.0\n",
      "1    0   11.0\n",
      "2    2   12.0\n",
      "3    2   13.0\n",
      "4    3   14.0\n"
     ]
    }
   ],
   "source": [
    "import cudf; print('cuDF Version:', cudf.__version__)\n",
    "\n",
    "\n",
    "# here we create a cuDF DataFrame with\n",
    "# two columns named \"key\" and \"value\"\n",
    "df = cudf.DataFrame()\n",
    "df['key'] = [0, 0, 2, 2, 3]\n",
    "df['value'] = [float(i + 10) for i in range(5)]\n",
    "print(df)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "As before, we can take this cuDF DataFrame and perform a `sum` operation over the `value` column. The key difference is that any operations we perform using cuDF use the GPU instead of the CPU."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "60.0\n"
     ]
    }
   ],
   "source": [
    "aggregation = df['value'].sum()\n",
    "print(aggregation)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Note how the syntax for both creating and manipulating a cuDF DataFrame is identical to the syntax necessary to create and manipulate Pandas DataFrames; the cuDF API is based on the Pandas API. This design choice minimizes the cognitive burden of switching from a CPU based workflow to a GPU based workflow and allows data scientists to focus on solving problems while benefitting from the speed of a GPU!"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<a id=\"scikitlearn\"></a>\n",
    "## Scikit-Learn\n",
    "\n",
    "After our data has been preprocessed, we often want to build a model so as to understand the relationships between different variables in our data. Scikit-Learn is an incredibly powerful toolkit that allows data scientists to quickly build models from their data. Below we show a simple example of how to create a Linear Regression model."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "NumPy Version: 1.16.2\n"
     ]
    }
   ],
   "source": [
    "import numpy as np; print('NumPy Version:', np.__version__)\n",
    "\n",
    "\n",
    "# create the relationship: y = 2.0 * x + 1.0\n",
    "n_rows = 100000  # let's use 100 thousand data points\n",
    "w = 2.0\n",
    "x = np.random.normal(loc=0, scale=1, size=(n_rows,))\n",
    "b = 1.0\n",
    "y = w * x + b\n",
    "\n",
    "# add a bit of noise\n",
    "noise = np.random.normal(loc=0, scale=2, size=(n_rows,))\n",
    "y_noisy = y + noise"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can now visualize our data using the `matplotlib` library."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.legend.Legend at 0x7f8e051565f8>"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAD4CAYAAAAJmJb0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO2de2ATVfbHvzdpaNMWSEt5FrDoCioUCpSHVkREHisPK6j4RHwsK+ruwmIFFKGwuFTBhZ+6im9ZqVoBKSCuiAIquiAtTXkIKCivUEqhtNB3mtzfH22GTDozeTTJ5HE+/0DOTGZO0sl37px77jmMcw6CIAgiNNGo7QBBEAThO0jkCYIgQhgSeYIgiBCGRJ4gCCKEIZEnCIIIYSLUdsCehIQEnpSUpLYbBEEQQUV+fv45znlbqW0BJfJJSUnIy8tT2w2CIIiggjF2XG4bhWsIgiBCGBJ5giCIEIZEniAIIoQJqJi8FGazGadOnUJNTY3arhBBQlRUFDp37gydTqe2KwShOgEv8qdOnULLli2RlJQExpja7hABDucc58+fx6lTp9CtWze13SEI1Ql4ka+pqSGBJ1yGMYY2bdqgpKREbVcIFcgtMGHJ5sM4XVaNTgY9Mkb1QHrfRLXdUpWAF3kAJPCEW9D1Ep7kFpgw57N9qDZbAACmsmrM+WwfAIS10NPEK0EQIcGSzYcFgbdRbbZgyebDKnkUGJDIO6GsrAyvv/662m7g2LFj6NWrl9N9PvroI+F1Xl4e/vrXv3rVj6SkJJw7d66JfcOGDcjKyvLquQjCHU6XVbtlDxdI5J2gJPIWi0XS7in19fXNer+jyKempuKVV15prlsuMX78eMyePdsv5yIIKQzR0tlUnQx6P3sSWJDIO2H27Nk4evQoUlJSkJGRge3bt2PYsGG47777kJyc3GSEvXTpUmRmZgIAjh49itGjR6N///4YMmQIDh061OT4mZmZmDp1KkaOHInJkyfDYrEgIyMDAwYMQO/evfHmm282ec+xY8cwZMgQ9OvXD/369cOPP/4o+Pr9998jJSUFy5Ytw/bt2zF27FgAQGlpKdLT09G7d28MHjwYe/fuFc7/yCOP4Oabb8aVV14p3BQqKysxZswY9OnTB7169UJOTo5w/ldffRX9+vVDcnKy8Jk++OADPPXUUwCAKVOm4PHHH8eQIUPQvXt3fP755839MxBhTG6BCWlZW9Ft9iakZW1FboFJcp+KmqaDJJ2WIWNUD3+4GbAExcSrjenTp8NoNHr1mCkpKVi+fLns9qysLOzfv1847/bt2/HTTz9h//796NatG44dOyb73qlTp2LFihW4+uqrsWvXLjzxxBPYunVrk/3y8/OxY8cO6PV6vPXWW2jdujV2796N2tpapKWlYeTIkaLJxHbt2mHLli2IiorCr7/+invvvRd5eXnIysrC0qVLBVHdvn278J758+ejb9++yM3NxdatWzF58mThMx06dAjbtm3DpUuX0KNHD0ybNg1ffvklOnXqhE2bNgEAysvLhWMlJCRgz549eP3117F06VK88847TT7TsWPH8O233+Lo0aMYNmwYjhw5gqioKNnviiCkcHUydcnmwzBbm7YyjWkREdaTrkCQiXygMHDgQKc52BUVFfjxxx9x1113Cbba2lrJfcePHw+9vuGR8quvvsLevXuxZs0aAA3i+uuvv6J79+7C/mazGU899RSMRiO0Wi1++eUXpz7v2LEDa9euBQDccsstOH/+vCDcY8aMQWRkJCIjI9GuXTsUFxcjOTkZTz/9NGbNmoWxY8diyJAhwrEmTJgAAOjfvz8+++wzyfPdfffd0Gg0uPrqq3HllVfi0KFDSElJceonEZp4mtooN5k6PceIJZsPC8eRi7uXV5u94n8wE1QirzTi9icxMTHC/yMiImC1WoXXtpW5VqsVBoPBpScP++NxzvHqq69i1KhRon3snxiWLVuG9u3bo7CwEFar1aURslTDdtvTQWRkpGDTarWor69H9+7dkZ+fjy+++AJz5szByJEjMW/ePNH+tn2lcExjpLTG8KU5qY1Kk6b2x+lk0MMksW+4x+MBisk7pWXLlrh06ZLs9vbt2+Ps2bM4f/48amtrhVBJq1at0K1bN6xevRpAg8gWFhY6Pd+oUaPwxhtvwGxuGIH88ssvqKysFO1TXl6Ojh07QqPR4MMPPxQmgJV8vemmm5CdnQ2gIYyTkJCAVq1ayfpx+vRpREdH44EHHsDTTz+NPXv2OPXdntWrV8NqteLo0aP47bff0KNHeMdFw5nmpDY6E2nbcTJG9YBepxVt0+u0YR+PB4JsJK8Gbdq0QVpaGnr16oU//vGPGDNmjGi7TqfDvHnzMGjQIHTr1g3XXHONsC07OxvTpk3DokWLYDabcc8996BPnz6K53vsscdw7Ngx9OvXD5xztG3bFrm5uaJ9nnjiCUycOBGrV6/GsGHDhCeB3r17IyIiAn369MGUKVPQt29f4T2ZmZl4+OGH0bt3b0RHR2PlypWKfuzbtw8ZGRnQaDTQ6XR44403XPq+bPTo0QNDhw5FcXExVqxYQfH4MEZqhA1cHqUrhXIyRvUQPQXIHce2P612bQqTeoxXi9TUVO7YNOTgwYO49tprVfKI8IQpU6Zg7NixuPPOO1Xzga6bwCC3wIQZOUZIqUxioxA7irhep8XiCcmCQNtuAnI3i0SDHj/MvsUX7gcNjLF8znmq1DYK1xAE4TOWbD4sKfAMDaN0V0I56X0T8cPsW7B8UgqFZDyAwjWE1/nggw/UdoEIEOQmTjkaxHtGjnRigtT7KCTjGSTyBEH4DLmsl8TGCVV3s2LS+yY2EXWqPKkMiTxBhBH2gtharwNjQFmV2WfiKDdxWlVXj9wCk2xM3tUQjK8rT4bCDYREniDCBEdBLLNbKOSrsry2Y2VuOCA634UqM+Z8tg+LJyRj8YRkt4TUXng1jMHikDxii+m7+zkcBX3YNW2xNt8U9KWLKbuGCEnoumlKWtZW2QwVexJ9MGKVO7ezzBj7zBpto6AzQHIy1xEGILqFFlV1FnAAWsZw76AuWJSeLHkexycKufMEYjYPZdcEKDfccINH25TIzMzE0qVLne4XGxuruN1fJZbnzZuHr7/+WnGf7du3C0XYCM9xteSubcQqVQjM2+eWstsKkiXN3oQZOUbh5mAbsbs6LOUAKhsF3vb+VTtPYG7uvib7SmX5yJ0n2EoXk8iriJRw2Vavqi1q/hL5hQsX4tZbb1Xch0RejCtVGaVwZ4m/fRqjp+ezf58crfXi8sC2EbVN2H0RZ/h418kmNneEO9hKJYScyHt6QSqxatUqDBw4ECkpKfjzn/8sCHFsbCxmzZqF/v3749Zbb8VPP/0klOzdsGEDgIZ0wttvvx2jR49Gjx49sGDBAuG4ttG0Y/li+20A8NJLLyE5ORl9+vQRara//fbbGDBgAPr06YOJEyeiqqpK8TP8/vvvuP766zFgwAA8//zzgr2iogLDhw8XSgevX78eQNMSy3L7ORIbG4uZM2eiX79+GD58uNBr1Wg0YvDgwejduzfuuOMOXLhwAUDDwilbMbakpCTMnz9fVMb42LFjWLFiBZYtW4aUlBR8//33WL16NXr16oU+ffrgpptucuVPGDLYiyDH5VH33Nx9Tq97qaX/Spwuq5Y9n7PfleP75MTavqRRboEJMz8tVFzdKoeWMbhaHckxhg/IC7fjMYMxLz+kRN7TC1KJgwcPIicnBz/88INQ9dFWA6ayshI333wz8vPz0bJlS8ydOxdbtmzBunXrhGJeAPDTTz8hOzsbRqMRq1evhuO8g22fF154AT///LPI/t///he5ubnYtWsXCgsL8cwzzwBoqAS5e/duFBYW4tprr8W7776r+Dn+9re/Ydq0adi9ezc6dOgg2KOiorBu3Trs2bMH27Ztw8yZM8E5R1ZWFq666ioYjUYsWbJEdj9HKisr0a9fP+zZswdDhw4VbmqTJ0/Giy++iL179yI5OVl0s7PHVsZ42rRpWLp0KZKSkvD4449jxowZMBqNGDJkCBYuXIjNmzejsLBQuJmGC3KLh7J3nnB63af3TcTiCclINOjBABj0OsTJNNoAAA1jmJ5jlDzfzE8LZX9X7oh1WVXDZOzc3H2YkWOUFGBn6HVavHx3H/yeNQYGvfznsaGVKJYnV/vm/sFdhe8r0aAXrcQNFkIqu0Zp9Zynf5hvvvkG+fn5GDBgQMPxqqvRrl07AECLFi0wevRoAEBycjIiIyOh0+mEZiI2RowYgTZt2gBoEOcdO3YgNVU8RyJXvvjrr7/Gww8/jOjoaABAfHw8AGD//v2YO3cuysrKUFFR0aRqpSM//PCDUGr4wQcfxKxZswA0FE579tln8d1330Gj0cBkMqG4uLjJ++X2s79hAIBGo8GkSZMAAA888AAmTJiA8vJylJWVYejQoQCAhx56SFSC2R5XyhinpaVhypQpuPvuu4X9wwWlxUX2yF33cnnmUmmOSoJr4RzTc4xYsPEA5o/rKRxzbu4+ZO884VbcvO/Cr3ChyvOSwBP7J2LJ5sMNC6tcGM7fO6hLE1soL7QKKZH3RY9HzjkeeughLF68uMk2nU4nlNDVaDRCCV6NRiMqwetK2V37csOO55faf8qUKcjNzUWfPn3wwQcfiBqEyCF1nOzsbJSUlCA/Px86nQ5JSUlCuWRP9nPlnEq4UsZ4xYoV2LVrFzZt2oSUlBQYjUbhJhrqyC0eksLV6z69byLyjpe6Jc42LlSZMT3HiOk5RpezXqSO4Slx0TpRmqMzB2JaaCWzawDpG6A/8HUufkiFa+Tias2ZKBk+fDjWrFmDs2fPAmhoo3f8+HG3jrFlyxaUlpaiuroaubm5SEtLc/m9I0eOxHvvvSfE3EtLSwEAly5dQseOHWE2m4XwkRJpaWn45JNPAEC0f3l5Odq1awedTodt27YJn82xbLHcfo5YrVYhxv7RRx/hxhtvROvWrREXF4fvv/8eAPDhhx8Ko3pXcPTl6NGjGDRoEBYuXIiEhAScPNl0Ii1UkQoryN1G3bnuN+0tavYkp7+TsRmAGrPFrRh+VZ13+zK7itxcoS9CzI54ReQZY+8xxs4yxvbb2eIZY1sYY782/hvnjXMp4Yua0tdddx0WLVqEkSNHonfv3hgxYgSKiorcOsaNN96IBx98ECkpKZg4cWKTUI0So0ePxvjx45GamoqUlBQhPfIf//gHBg0ahBEjRojKG8vxf//3f/j3v/+NAQMGiFr53X///cjLy0Nqaiqys7OFY9mXWM7IyJDdz5GYmBgcOHAA/fv3x9atW4W5iZUrVyIjIwO9e/eG0WgUzVk4Y9y4cVi3bp0w8ZqRkYHk5GT06tULN910k9PyzaGEY1w90aDHDVfFN2uCMLfA1KzRtFpwANVmq9P97NEw5tWkDFdQEvLm1Np3Fa8shmKM3QSgAsB/OOe9Gm0vASjlnGcxxmYDiOOcz1I6jjcWQwXaMuQPPvgAeXl5eO2111TzwZ/ExsaioqJCbTfCZjGU3CKe+wd3lQ1LONLcmHiw4u735ClKC8FONwq/lG+/Z42R2CKN0mIor8TkOeffMcaSHMy3A7i58f8rAWwHoCjy3kCtuBpB+AqlgYvcIp6Pd51E6hXxwj5Kg55gFHhP4//2cACrdp7Aqp0nfLLK14bSXKE/2hb6cuK1Pee8CAA450WMsXZSOzHGpgKYCgBdu3b1oTvqMGXKFEyZMkVtN/xGIIziQwlnBbjkBMTCOTJWFwIMMFu45HttNw9X0QCI0mlQ5WaIxNvEReu8fmNytS6NJ5ECJSFvboE2V1B94pVz/hbnPJVzntq2bVu5ffzsFRHMhNL1IhezteWpK434zFYuCLz9e5dsPtxkZakrWAHVBR4AasxWxfx+T3EWC7fl8rs7Sao0Vyg1x+LtXHxfjuSLGWMdG0fxHQGc9eQgUVFROH/+PNq0aeN2Oh4RfnDOcf78+ZDpKas0Up/z2T5M7J8oTiF0AVNZteTNI1ioNlsQGeHe+FSrYbBYnd/85b7v3AITVu08IemLs3U4Sjn4nHPs3fguatavx+7Nm5GQkODiJ3IdX4r8BgAPAchq/Fd6HbwTOnfujFOnTgnL4wnCGVFRUejcubPL+wfaZL09Snnx1WYLth0qweIJyZj5aaFbq0XdGcEHIvZli13BYnWteqXck1HmhgOy73FlPYLjXKHFYsGjjz6K9957T7BVVFQErsgzxj5GwyRrAmPsFID5aBD3TxljjwI4AUB6iaMTdDqd5EpQgvAGvm460Ry/lJpX2zCVVTdrMVM44ey7kYqF2/4OSjcVdyZJq6qqMH78eHzzzTeCrWfPnvjuu++E1ezeJuDryROEL/G0zrkj3nwakCszIIUtDTDnp5MwuxCOCAcS3VgVbCMyQoO6equo7nzqFfHIWFPYZF7DkeWTUpz+rUtKSnD99dfj6NGjgm3MmDFYs2aNV0KLPk+hJIhgxRulMKSeBmbkGJF3vFQxB1vuxuBOvJwDNIJ34EJlLXRa5lSc7amtvzyhbKs7n/PTCTR3nvnIkSO45pprhMq1APDkk0/ilVdegUbjn7wXEnkirPFGnrJcrnr2zhOyueoAZMNE7tZaIoEX460MIFcPIzXxumvXLgwePFhky8rKEgoD+hMSeSKscTVPWSkco1QZcsHGA6gxW5uIeWSERnY5uztFyAj1sf/7r1+/Hunp6aLtH330Ee69915/uyWgep48QaiFfe0QW41xqTxlZ0WklEb9F6rMkmIuN5FnKqtGUhu9yw0wCPXpZNDj9ddfB2NMJPA9//QykmZ9jteOt/dbnRwpaCRPhCWOcXQL56IFKvY461OQMaoHZuQYvRY2+eFoqZeORDQHDRoWgMnBOUfZtx/g+K61sG9O+X+rv8EbhXWoCJCMLRrJE2GJO9X/nE3OpvdNxP2DuzYZfeu0NB4PZuQEnlvqUbL+RZx4aRwu7mpoxNOhQweYTCZwzpFzBD6vLOkONJInwhJ3smpcmZxNvSIem/YWiWqquJPdQQQ+1toqFOc8j7qiy2Id2akHzh78Ca1atRJsvmhe1BxoJE+EJXJxdA40qTXurE+BLfQTjNUcCefUXzqPk689gJPL7xYEPvqaIej6dC46PPiySOAB3zQvag40kifCEqmsGhuOMVRbHHXBxgOCkNfUWzA9x+h2OQEieKg7dwJF7z4hsrUadCcMQx8S6mglSgi3PypLugOJPBGW2BeNkgrFSBWeqrFLnLbpOgl86FFzYh+KP54jssWPeBwt+40V2eSEO9CagpPIE2GLbZTebfYmycwY+xhqMFdtJFyj8udvcW7jEpGt7YS5iL56MHQahoHd4rDztwuwcA4tY5jYX75BUSA1LyKRJ8IeVyZW1Zo0I3xP+a41KNv+gcjW4YGliEy83Me4RYQGe06UC09uFs6xNt+E1CviA0bM5SCRJ8IeZzHU3AITNIxRaCaE4NyK0i1voqJg02WjRotOj70BXVynJvtX1jV9inOllnwgQCJPhD2O8XktY8IPOO94Kdbmm0jgQwReX4eS3MWoPrpbsEUYOqLDg0uhjW7t9vGC4QmPRJ4IK+Rq0NiE3rFomFQ3ICL4sFRfQvFHs2A+d/nvGXVFH7Sd+Dw0OuVSvzotQ0yLCMlSFGqlRboDiTwRNsiVBJ6eY0SiQY/K2nqaXA0x6suLcfq9p8DrLo+4Y3uPRPyoJ8E0WoV32sE5Msf3DKi0SHcgkSeCGvsOStrGuHmiRMpaboEJMz41wjHqYntJVR9Di9ozR3Bm5XSR7d5pT+NYl9E4XV7j1rHMViDveCkWT0gOmLRId6DOUETQotRBSa/TCtUk5+buo7BLmFD9Wz7Orp4vsrW5bTpik2/FsawxHl8LWsZwdPFt3nLT61BnKCIkUcpdty8IlU0CH/JcKvwKpV++IrK1u3sh9N36AQDionXILTBhbX7Tkr86DaDTahSbjQTzxDuJPBG0OAuxnC6rxoKNB6hzUojCOUf5jo9Q/uPHInvHKa+gRfsrm+yfueGA5KCgXavL/Xy7zdnUJKQHQOg34Axv9vr1FiTyRNCidZK7rtdpqGhYCMKtFpz/7yuo3P+NYGORMej0yKuIaNVO8j1K14F9GuT9g7pKhnMGXxnn1C+piX0168jbIJEnghYlgddpmdd6fRKBgbWuBmfXZKL25H7BpmvXDe3vXQxtVKzHx22t1wn/X5SejN9LKpo0btlzohy5BSZFsXbWXEYtSOSJoCVRoRcq1XIPHSyVF1D0n5mwXDwr2PRXD0bb8bPAInQK73QNx0jMsfOuFaxzJNDqyNsgkSeCFqVywUTwYy414fTbfxbZWvYfj7jhj4Ex77XCcAzleCrWrtRAUgMSeSIocWzCHczZD4SYmlMHUZydIbLFDXsErQZO8Mn5GCAKxXgq1oFWR94GiTzhV9zNPpBa7MQAUcYMCXxoUHX4R5Tk/lNkSxg/CzHXDvHpeTkgpNvarjXHa8wVsQ60OvI2aDEU4TekFi/ZL1qy30/ux0aEHhfzN+LC12+KbO3vy0JUl15+9UOv04quTdu1J7WCOtCgxVCE31AaqbuSfeB4IyCBD0045yjb9h4u7l4nsnd69A3oErr43R9b5VF7QuXaI5EnmoW9qBuidaioqYfZ2vDzcMwTVprQsh+9E6ELt5hxbsMSVP3yo2DTxrZBh4eWISI2XhWfdBomXLNSBEq+u6eQyBMe4zjqllpwYj9Sl5vQim6hxYwcY8iMnIimWGsrUfzJc6g7c0SwRSZeh3Z3ZUITGa2iZwAYYNDrJEsJ2wiEfHdPIZEnPMbVvqe2EbxU9oFOyyS77hChQf2lcyh6/6+wVl8UbDHX3Yw2t00H0waG/JgtHIw1jck74vgkGoglDKQIjG+ZCEpcXeRhSz2Tyj6orK1XHEERwUldyTEUvfeUyNZq8F0w3DQZzMU6MP6krMqMZZNSFEOG9imUgVrCQAoSecJj5MIv9jA0/ADSsrY26cKUW2DC9ByjHzwl/EX18UKc/eQ5kS1+1JNomfJHlTxyjU4GvXBtymWB2adQBmoJAyl8LvKMsWMALgGwAKiXS/Mhgg9XVpzaN+WwH+nYfkhEaFBxYBvOf/6yyNZ24jxE/2GgSh65jqOAu5LvHqglDKTw10h+GOf8nJ/ORfgJ20W/YOOBJpOuUvnt1WYLFmw8gPS+iS7H84nApnznapR9u1Jk6zD5X4js2F0lj9xDLgfe/olTCk9XxaoRx6dwDdEs7B9x7S9euTDOhSoz+i78ikoABzHcakHpljdQYfzyslGrQ6fH3oDO0EE9x9wk0XC5jry7eFLCQK04vj9EngP4ijHGAbzJOX/LfiNjbCqAqQDQtWtXP7hD+AJHsVeCBD44sZprUbLuBdT8vkewRcR3Rof7X4Q2urWKnrmPTsOaVVPGkxIGasXx/SHyaZzz04yxdgC2MMYOcc6/s21sFP23gIayBn7wh/ARSj1XieDFUn0RZ1Y9g/rSU4Itqls/tL3jOWh0kSp65hkMwKSBXZotrM5COo6oFcf3uchzzk83/nuWMbYOwEAA3ym/iwhGKM4eWpjLzqDo3SfB62sFW2zKaMSPmAam0aroWfPgANbmm5B6RcMKW3/FyNUqRexTkWeMxQDQcM4vNf5/JICFvjwnoR6BmFlAuE9t0a84858ZIpvhpsloff3dKnnkGXqdBjX1VsmerbYkgBqz1W8xcrVKEft6JN8ewLrGxQ8RAD7inH+p/BYiWHElb54IXKqO7kbJmgUiW5sxf0dsL88mJ9VEp2Gos3BJgbfhrAyHt1GrFLFPRZ5z/huAPr48B+FflFLAMkb1kK1BY9DrEBMZQTeBAOSS8UuUbn5NZGs3aRH0SSkqedQ8GANioyI8nuD35ROpu3F8b0AplITLOEsBS++bKLuCtazaDOP8kUjL2kpCHwBwzlH2/Spc/F+OyN7x4VfRol03lbzyDvoIDcqcCLxOyxDTIkKypIba7fq8jfcaJRIhj1IKmI1EmR+IrcVaxqge0ARe6ZKwgVstOLdxKU68NE4QeI2+FRKnvY8rZn0e9AIPAFVmKwzRyg2+Y1pEIHN8T+h14gnkQGjX521oJE+4THPqwXMAM3KMuH9wVyiU7iZ8hLWuGmc/nY9a08+CrUWHP6D9PS9AExmjome+gfOG0brZIn2xlVebA7Zdn7chkSdcRm5i1RCtcyk/ngNYtfOEj7wjpLBUXEDRf/4Oy6USwRbd/QYkjM8A0yqPdoOZ8uqGqpIzPjVKTr7aV0YNNVF3hESecBm5FDDOQfnxAYb5/EmcfmeayNZyQDrihj0akKV+vY2tqiQAVdIWAwkSecJl5B5vZ1C54ICh5tQBFGfPEtnihv8JrVJvV8kj/2Mv4uESklGCRD6McacinuO+yyalCPtmbjhAjT9UpvLQDpxbnyWytU1/FtE9blDJI3XQMobFE5JF13E4hGSUIJEPM+wnSO3LAUut9lPad0aOEXnHS7EoPRmVtSTwanFx93pc2Pq2yNb+/pcQ1fk6lTxSDwbg5bv7hLWgS0EiH0Y45rlL1Xu3rfZzti8HkN04iWq2+tZvQgznVlz45h1cyt8gsnd6bAV0bTqr5JX6cARe671AgEQ+jHClgJipMR0yc8MBl7JlPt510oseEkrwejNKNryI6l93CjZtq3boOPllaGPiVPQsMAiD+WSPIJEPI1xdrp2xuhBmF5PZLUrFQQivYK2pQPEnz6Gu+Khgi+zSC+3uzISmRZSKngUWdClKQyIfRrhaQMxVgSd8S/3FEhS9/xdYayoEW0yv4Wjzx78Gdalfwr+QyAcxStkxUttcabxNqE/d2d9Q9P5fRbbWN9yD1jfeHxY57p4S56SUQbhCIh+kKBULAyC5bfGEZCyekCxkzGgZo3BLAFF9zIizOXNFtvjRf0HLPqNU8ih40GkZ5o/rqbYbAQmJfJDirFiY3DZb42Ia0QcOFfu/wflNy0S2dndmQn9VqkoeBReJYbjAyR1I5IMUT/pFmsqq0XfhVyirNtMklcpwzlH+vxyUf79KZO/w0HJEdviDSl4FHwwQBi6ENCTyQYqzfpFyE6yeNlIgvAO3WlC6+d+o2PuVYGO6KHR85DXoDGqKOtsAABsJSURBVB1U9Cw4CbXa776ARD5IkZtEraqrx5jeHbE230ThmADCaq5BydpFqDl+uc6PLqEr2t/3IrT6lip6FrzotCysCo15Col8kCJXN+ZClRlr802Y2D8R2w6VUBcmlbFUlePMqgzUXzgt2PRXDUDb9DlgES1U9Cz4iWkRQXF4FyCRD2LS+yZiyebDTYqDVZst2HaoBD/MvgU9532Jyjoa0fsb84WihlK/1nrBFtv3NsSPeByMUUM2b1BORfFcgkQ+yJGbaLVNspLA+5fa04dx5sOZIpth6BS0HnynSh6FLhSPdw0S+SBHaRUrTbL6j6oju1Cy9h8iW8K4DMRcN1Qlj0KbcGv80RxI5IMcW9MOyohUh0sFX6D0q9dFtvb3/hNRXXur5FHoQ3nx7kEiH+Sk903EdOrM5Fc45yj7diUu7lojsnd85DW0aJukjlMhgl6nFa3KdiTRoKe8eDchkQ8BEl0sPEY0D26px7lN/0LVwe8EmybGgI4PLUdEywQVPQsNHEfo4d6b1VvQNH8IQBe+b7HWVuHMh0/jxNJ0QeBbdOyBLtM/RZenVpHAewFb3bUZOUakZW0FAEzsnwht4wYtY5jYP7zb+HkK4wG0vj01NZXn5eWp7UZQ0nfhVzTR6mXqK0pxZuV0WCpKBVv0NUOQMHYmmJYegn2JTsMABpgtl/XJFsohoW8KYyyfcy5Z7IiuVD/hTtNsT5g/rif+/qkRVAq++dSdO4Gid58Q2VoNnADDzQ9TqV8/IdXTwL49JeE6JPJ+QKkssNwF68lNQcsYrAH0ZBZs1Jzcj+KPZots8SMeR8t+Y1XyiHDE1e5mxGVI5P2AUllgKeF256ZguxnQxKvnVB78Duc2vCSytZ0wF9FXD1bJI0IOWgDlPiTyfsDdssCu3hQcbwaEe5Ss+yeqfvlRZOvwwBJEJl6rkkehS6JCdVSDXofaeqvoOpaLyVOSgfuQyPsBZ2WBHXH1piB1MyCU4VYLTr/9Z9SXnRHZO/3pTejiKdbrC+zFWSotMnN8Q0cnx/CklI3i8e7j8+waxthoAP8HQAvgHc55lty+oZpdIzXiVsoUSMvaKnlTiIvWIbpFhHDRU4jGdax1NTi5rGn9mMTH30VE6/YqeBQe2FpMJpJw+xSl7BqfijxjTAvgFwAjAJwCsBvAvZzzn6X2D1WRB9ybSJW6Kei0DODSWQeEPPUVpTD9e3ITe5fpOdBExqjgUfig0zDR9UopkL5DzRTKgQCOcM5/a3TkEwC3A5AU+VAmvW/ThRxKwh+l0wgib9A3dKF3LClMyFN39ncUvf8XkU0T1RKdn/qQctz9hOOAhFIg1cHXV3sigJN2r08BGGS/A2NsKoCpANC1a1cfuxM4zM3dh+ydJ4TCYqayamSsLsSzn+1Fldkq2reyrl40AUXIU/1bPs6uni+yRXbtjfb3vEA57gEApUD6H1+LvNSvSqRWnPO3ALwFNIRrfOxPQJBbYBIJvA2zlUuGY8wWDsZAzbcVuLRnE0q3vCGytew/DvG3/lklj8Ibg14n+eRJKZD+x9cifwpAF7vXnQGcltk3bFiy+bDbpYFJ4KUp/fpNXMrfKLLFj5iGlv3GqOQR8cDgrki9Ip4KjAUIvhb53QCuZox1A2ACcA+A+3x8zoCHHlmbB+ccZdvfx8WfPhPZ292ZCf1VknNPhJ+Ii9Yhe+cJbDtUIvQZpkwadfGpyHPO6xljTwHYjIYUyvc45wd8ec5ggNIfPYNbLTj/31dRuf9rkb3jw6+iRbtuKnlF2GC43I3MVFaNtfkmyqYJAHyeZsA5/wLAF74+TzCRMaoHrVR1A2tdDc6uXYjaE3sFm65tEtrflwVtVKyKnhH2OEYUKZsmMKBcMhWwXfQLNh6g8sAKWCrLcObDmagvLxZs+j8MRNvbZ4NFtFDRM8IehqYCb4NCk+pDIq8S6X0TsWTzYRJ5CcylJpx++3HYS0fL/uMQN/xPYIz63AQStpWsckXyHLNpfF1ym2gKibwfcbzAKS4vptZ0CGdWPS2yxQ17BK0GTlDJI0IKuUbazrJpPCm5TTQfEnk/IXWBEw1U/fI/lKx7QWRLGP8MYq69SSWPCDnkGmnbRFpplO5uyW3CO5DI+wmqGNmUi/kbceHrN0W29vdlIapLL5U8IpRwluduX7rD9tQ6I8coCL67JbcJ70Ai7ydo5N6AXI57x0dfR4uE8ClrEWwY9Dpkju/p0ohbLixjiNZJzkHRKljfQiLvJ2wlV8MVbjHj3MaXUXV4h2DTxsajw0PLEREbr6JnhCvEREa4HFKRC8tERmig12lpFayfIZH3E+Eq8NbaKhTnPIe6ol8FW2TitWh31wJoIqNV9IxwB3dCKnL7llebsWxSCmXX+BkSeT8RbiP5+kvnUPTB32CtKhds0dcORcKYGVTqNwhxJ6Si1AlNquQ24Vvo1+YlnOX/hovA15UcQ9F7T4lsrQbfBcNNk6nUbxDjTkhFakU3hWXUg0TeC7iS/5sY4nnxNcf3oviTZ0W2+JFPoGXf21TyiPAWBr3OrdG3K+mUhP8gkfcCzvJ/cwtMqKqrV8k731L583ac27hUZGs78XlE/2GQzDuIYMK+0bY7UFgmcCCR9wJK+b9S/VpDgfKdq1H27UqRrcODLyOyEz2SBxNKdWfkVrYSwQWJvBdQmmgKpUVQnFtRumUFKgrsiopqI9Dp0Tegi+uonmOEx3BAMq2RSgSHDlTtyQtkjOoBvU4rstkmmkJhNR+vr0Px6vk48dJ4QeAj4hPR+S/ZuOLpXBL4ICbRoMfiCclINOjB7F6TwIcONJL3AkoTTXLV+YIBS/UlFGc/A/P5y73Yo5L6ou2EudDoIlX0jPAGDMCwa9pS/DzEIZH3EnI/lGBsEFJfXozT7z4Jbq4RbLF9RiF+5BNgGq3CO4lgggNYm29C6hXxJPIhDIm8l1DKk4+M0ASFyNeeOYIzK6eLbIabJqPV4Lsoxz1EoSqQoQ+JvBeQypPPWFOIOZ/tRbXZqrJ3zqk+uhtn1ywQ2dqMmYHYXsNV8ojwJ6Ewb0TIQyLvBaQyaMwWDrMlsFe5XjJ+idLNr4ls7SYtgj4pRSWPiOYSJ1PpUQmqAhnakMg3A1uIJpgmVjnnKN+RjfIfPxHZOz78Klq066aSV4Q30GkZ5o/riek5Rtl9qApk+EEi7yHBtsiJWy04/8VyVB7YJtg0US3R8eFXENGqrYqeEc1ByxisnIvmgeQaxMdF6zB/XE8qNxBmkMg7wX60bqskmWjQo7K2PigE3lpXjbOrM1F76oBga9H+KrS/5wVoomJV9IzwBlbO8XvWGJFt/rieyFhTKAoX2kb5lC4ZfpDIK+A4WrdVknQWnomL1qGiph5mq3oxeUvlBRSt/Dssl0oEm7779Wg77hmwCJ1qfhHeRSqeTgXCCHtI5BXwpCSBrdGxfUqlP6XefP4UTr/zuMjWMvV2xN3yKBijBc6hhFI8nUbshA0SeQXcTS2z/9HZfmS5BSb8/VMjfD2orzn1M4qznxHZ4m75E1oNuN23JyZUQcsYlR8gXIJEXgG5wmNSSP3ocgtMmJFj9OlIvvLwDziXu1hkS7h9NmKuudGHZyXUxsq5pMDnFphEE6/uNOAmQhMSeQfswyyGaB00AFxZzmThHEs2HwYAYQQ/57N9PhP4i3nrceGbt0W29ve/iKjO7tf+JoIPqVh8boGpyYRrWbUZGasLAYCEPkwhkbfDcaLV3UUlprJqTM8xYsHGA+AcXs++4ZzjwtZ3cClvvcje6bEV0LXp7NVzEYGLXCx+yebDkgvwzFZOpQvCGBJ5O7xV+93dm4MzuMWMkg0vofqX/wk2bau26Dj5X9DGxHn1XERgo9TIQ2kOiUoXhC8k8nYE2g/BWlOB4k+eQ13xUcEW2bkn2t2VCU0LWooebsRF6/DD7FtktyvNIVHpgvCFRN4OdyZafUn9xXMoev8vsNZcEmwxPYehzW3TqdRviGPQ61BW3fRJ0LaYSYmMUT2axOQBQKdhVLogjPGZyDPGMgH8CYBtNc6znPMv5N+hPhmjevg8G0aJurO/o+j9v4hsrW+4B61vvJ9K/YYQcn1VDXodjPNHAlAuXS2HbTtl1xD2MM59I2mNIl/BOV/q6ntSU1N5Xl6eT/xxlbm5+5C984ToR2j7USo1PW4O1ceMOJszV2SLH/UUWqaM9sHZCLWJi9ahxmylvqqE12CM5XPOU6W20RJIBxalJ2PZpBRRz0vba28LfMX+rTj+4liRwLe7cz6umPU5CXyQ8sDgrnD2zFVWZaa+qoTf8HVM/inG2GQAeQBmcs4v+Ph8XkFqSfgMhfKt7sA5x8Wdq1H23X9E9g4PLUdkhz945RyEOhj0Omw7VOJ0MNDJoKeyA4TfaFa4hjH2NYAOEpueA7ATwDk0RDj+AaAj5/wRiWNMBTAVALp27dr/+PHjHvvjS9KytjZrUpZbLSj96nVUFG4WbEwXhY6PvAadQeorJIINnZY5bRRjC/kppUIShLsohWt8FpN3cCAJwOec815K+wVCTF4OT+vHW821KPlsEWqOFQg2XZuuaH//i9DqW3rbTUJlbOWopXCc06E4POEtlETel9k1HTnnRY0v7wCw31fn8geO5VsNTtqsWarKcSZ7FupLTwm2qCv7o236s9DoIn3uL6EOFs4luy9FRmiapEZSE23CH/gyJv8SYywFDYOXYwD+7MNz+QXHOOr9b/8PPxwtFe1jLjuD0+9MAyyXf9CxfW9D/IjHqdRvGGALwzimP8rN6QTaAjwi9PCZyHPOH/TVsQMF48ly4f+1Rb/gzH/+LtpuGPoQWg++y99uESphqykjNakq1wuYVqISvoZWvHrI3Nx9qKyzoOrITyhZu1C0rc3YmYjtOUwlzwh/oQHQOlqHsiqz0wVLGaN6NJnToSbahD8gkfeQN1e8hXObXxPZ2t/zT0Rd0Vsljwh/869JKYKo21aozsgxSgo+teQj1IJE3g0453juueeweLG4SUfHR15Di7ZJ6jhFqEbe8VIhDGOfOWMqq8acz/YBQBOhJ1En/E3YibxSTRCpbQDw4hcHsC/7n6g8+K1wHG20AR2mLEdEywRVPgehPvblLxyTJilzhggUwkrkHXPd7UdcAJpsm5n9P5z+ZB5qTAeFfaI6dcc7H+fiYGk9Vu084d8PQAQUzlaYUOYMEQiElchLNQWxjbhs/weA+opSnFk5A5aK88J+0T1uRMK4mWBaHV7/sUio6/3xrpOwcA4tY7h3UBds2lvk9aYhRHBCmTNEIBBWIi83srLZzedO4vS700TbWg2cAMPND4tK/ZrKqpE0exOAy6VcAXGJVyIwMOh1YKx53bpiWmhRVWeRrEwqB2XOEIFCWIm8XFOQ2Au/Yv9bM0S2uFv/jFb9xzk9Zlm1GdO9VLyM8B6OJQM8LUuh12nxwh3JAMSZMcOuaYu1+SbR8aguDRGIhJXIO+YqVx7agXPrs0T7tL3jWUR3v8Fr51SqZUL4Bg1Dk5owjimMGhf+LgzAxP6XM2IcRTv1inhKiSQCnrASedsPcMbcRTj2xQrRth9//BHFUV1FP9rK2nrJVmzuIFXLhPAtWpkuWvYpjK6M7DmAbYdKZLdTSiQRDISNyFutVkyfPh2vvvqqYNNqtTh48CCuvvpqwWb/o/X0Ed8e+1omgdA/NhwwW7nT9EXHkb3cmJ4yZIhgJ+QrZtXW1iI9PR1arVYQ+G7duqG4uBj19fUigXckvW+iqIOPQa9z69y2BsrpfRPxw+xbsHxSCvQ6asTtD5TEObfAhLSsrULRMFvnLykoQ4YIdkJ2JF9WVoabb74ZhYWFgm3YsGHYuHEjYmJiXD6O4yN5yoKvXA7hLLmrj2RcmCZqfY+GMeQWmJqM5uXWSkzsn9hkIpUyZIhQIORG8idOnIDBYEBcXJwg8A8//DDMZjO2bt3qlsBLkTm+p8uj8dV5J5CWtRXdZm9CWtZWQXRiWtBo3tdYOMecz/Yht8Akssutldh2qIT6rhIhSciM5Pft24fevcXFwRYsWIDnn39elOPeXGw/+swNB5yO6O1rzZvKqvH3T43IO16KyjqahPUUDQOsLiYrSZUWUForQROpRCgSMiP5oUOHCv9/9913wTnHvHnzvCrwNtL7JiJzfE9o3Dy0lYNKITQDvU6Lf92dgmNZY7DcLo4ul00DNBV1uRg7xd6JUCVkRH7nzp3YsWMHOOd45JEm/cK9zpLNh10eURLNxxY+ASCaNF0+KQVHF9/m8sRpxqgeTcJtFHsnQpmQCdd0794d3bt3l9ymVHnSUyi1zn3inPTFddweF63D/HE9ZXPb7QvMudqUg+q6E+FGyIi8HErC0JwftlyJBEIeZ/VjCuaNVNyuVGDOVjDOFfGm2DsRToS8yCsJQ3N+6FIjR52WIaZFBMqrzYjSaVBttnp8fKIpzgrMkXgTRFNCXuSdCYOnuPLYPzd3X9BNtGoAqHFrUpo8tSH39ESTpgQhT8iLvCfC4GoM39nIcVF6MtbtMQVVyiR3VkPXR9w7qIvTfagZNkG4T8iLvLvC4EoM3/4m0LqxXnlZlVmynWAwCTwASBVmZEza7g1szVYWpSc73ZcmTQnCfUJe5N0VBmcxfMebgP2CKMcbgq3jVLDDuXthHFcfBhINemHC1FUo7k4Q7hHyIg+4JwzOYvhSNwF7qs0WzPy0EDNyjGpEPXyGFYBOA7gyl8xxuWkHIF+rh9JQCcL3hMxiKG/hbEWkK8Jk4TykBN6GO8lC9k8/VOGRINSDRN4BZysimytMDEDaVfGIi3avbLEctqwU7xdvaD62GyKtMiUI9SCRd8CxhrxjNUIpwXIF27GWTUpB9p+uR8G8kVg+KcWl1EElrJwj0aD36MnBvvaOXqdBtM6zy0HuE9huiM6+U4IgfAfjAdR/NDU1lefl5anthlPksmuU+obKNXfuNntTs0I7iQa9YmcjBuCGq+Kx50S511sQMkC2qbVjI22CIHwHYyyfc54qtS0sJl6lsAm1qaxaaLYtJ8SOyE3kKrULlEvFdKWhtBy2kIdca0EtY3j57obGJXNz9yF75wm3bihKTcgdM2OoqTVBBCZhKfKOYmwTsubWtbFP15QSXalUTDkRjdZpUKUw02nQ65A5/nLxLqm1APYj6W2HStx+YrByjuWTUlwu/EWiThCBR1jG5JXSIG1C7Cm2fq5ycWpnqZhaxrB8Ugp+/scfFXvK1tZfvgG4EvP2pJhaJ4Oe4ukEEeSExEje3VLCztIgvZG/LVdOgaOhHrqc6Fo5v9x9anxP2fBPtdmC6TlGLNl8WPi8Sp9ZKfQihf1onUbpBBG8NGskzxi7izF2gDFmZYylOmybwxg7whg7zBgb1Tw35bGFPUyNk4+2kItjb097nKVBeiN/WykLx1RW7TQjBbg8QlfClc8LQFHgEw16PDC4K43WCSIEae5Ifj+ACQDetDcyxq4DcA+AngA6AfiaMdadc+71Qi6elBKWqmdjw1v5287i8xxNl//LxbrljmHDldLJiTJPFp6UFiAIInho1kiec36Qcy4VwL4dwCec81rO+e8AjgAY2JxzyeFJKWH7ODNweUGRt0ewzuLzvPGczkbPruTmOwsxeXtBUm6BCWlZW9Ft9iakZW11+iRBEIQ6+Comnwhgp93rU422JjDGpgKYCgBdu3Z1+0Se1hj3Z5xZzkdXR9E2PxdsPCDbXam1wiSt/TG8keboq25bBEF4H6cizxj7GkAHiU3Pcc7Xy71NwiYZFOacvwXgLaBhMZQzfxxRs8a4qxO+3vKxRiGlsrKuHrkFJkWR9daNzVfdtgiC8D5ORZ5zfqsHxz0FwL4LRGcApz04jlPUqjHuzmjWGz46q35ptnC/iayvum0RBOF9fBWu2QDgI8bYv9Aw8Xo1gJ98dC5VUvzcHc0210dXBNRfIktt+AgieGhuCuUdjLFTAK4HsIkxthkAOOcHAHwK4GcAXwJ40heZNWri6WjW0wlLVwTUXyJLVSUJInhobnbNOs55Z855JOe8Ped8lN22FzjnV3HOe3DO/9t8VwMLZ3XnpfAkp9+Gswwbf4osrYIliOAhJFa8qoEnk6nNmbB0jOsbonXgHCivbtpb1h/QKliCCA5I5D3Ek8nU5k5YkrASBOEuJPLNwF3RpQlLgiD8TVhWoVQLmrAkCMLf0Ejej6iV008QRPhCIu9nKK5OEIQ/oXANQRBECEMiTxAEEcKQyBMEQYQwJPIEQRAhDIk8QRBECMO4G82dfQ1jrATAcR+fJgHAOR+fI1ih70Ye+m6koe9FHn9+N1dwzttKbQgokfcHjLE8znmq8z3DD/pu5KHvRhr6XuQJlO+GwjUEQRAhDIk8QRBECBOOIv+W2g4EMPTdyEPfjTT0vcgTEN9N2MXkCYIgwolwHMkTBEGEDSTyBEEQIUxYizxj7GnGGGeMJajtS6DAGFvCGDvEGNvLGFvHGDOo7ZOaMMZGM8YOM8aOMMZmq+1PoMAY68IY28YYO8gYO8AY+5vaPgUSjDEtY6yAMfa52r6ErcgzxroAGAHghNq+BBhbAPTinPcG8AuAOSr7oxqMMS2AfwP4I4DrANzLGLtOXa8ChnoAMznn1wIYDOBJ+m5E/A3AQbWdAMJY5AEsA/AMAJp5toNz/hXnvL7x5U4AndX0R2UGAjjCOf+Nc14H4BMAt6vsU0DAOS/inO9p/P8lNAgaNUoAwBjrDGAMgHfU9gUIU5FnjI0HYOKcF6rtS4DzCID/qu2EiiQCOGn3+hRIyJrAGEsC0BfALnU9CRiWo2EAaVXbESCEO0Mxxr4G0EFi03MAngUw0r8eBQ5K3w3nfH3jPs+h4ZE825++BRhMwkZPfnYwxmIBrAUwnXN+UW1/1IYxNhbAWc55PmPsZrX9AUJY5Dnnt0rZGWPJALoBKGSMAQ3hiD2MsYGc8zN+dFE15L4bG4yxhwCMBTCch/dCilMAuti97gzgtEq+BByMMR0aBD6bc/6Z2v4ECGkAxjPGbgMQBaAVY2wV5/wBtRwK+8VQjLFjAFI551RJDw3ZJAD+BWAo57xEbX/UhDEWgYbJ5+EATAB2A7iPc35AVccCANYwQloJoJRzPl1tfwKRxpH805zzsWr6EZYxeUKR1wC0BLCFMWZkjK1Q2yG1aJyAfgrAZjRMLH5KAi+QBuBBALc0XifGxtErEWCE/UieIAgilKGRPEEQRAhDIk8QBBHCkMgTBEGEMCTyBEEQIQyJPEEQRAhDIk8QBBHCkMgTBEGEMP8Pe6YEQgghYUcAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.scatter(x, y_noisy, label='empirical data points')\n",
    "plt.plot(x, y, color='black', label='true relationship')\n",
    "plt.legend()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We'll use the `LinearRegression` class from Scikit-Learn to instantiate a model and fit it to our data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Scikit-Learn Version: 0.21.2\n"
     ]
    }
   ],
   "source": [
    "import sklearn; print('Scikit-Learn Version:', sklearn.__version__)\n",
    "from sklearn.linear_model import LinearRegression\n",
    "\n",
    "\n",
    "# instantiate and fit model\n",
    "linear_regression = LinearRegression()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "CPU times: user 44.4 ms, sys: 19.9 ms, total: 64.2 ms\n",
      "Wall time: 10.2 ms\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "LinearRegression(copy_X=True, fit_intercept=True, n_jobs=None, normalize=False)"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "%%time\n",
    "\n",
    "linear_regression.fit(np.expand_dims(x, 1), y)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "# create new data and perform inference\n",
    "inputs = np.linspace(start=-5, stop=5, num=1000)\n",
    "outputs = linear_regression.predict(np.expand_dims(inputs, 1))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's now visualize our empirical data points, the true relationship of the data, and the relationship estimated by the model. Looks pretty close!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.legend.Legend at 0x7f8dc80a5f60>"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAD4CAYAAAAJmJb0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO2de2DN9f/HH+/NwZnLRqisXCpRjGFJSW6FIi1KuSQJUSrSXLq4lEKUSiGV+EWFLktXEipKwjaXmtxl+rrNhI0dO+/fH8c5nZ2d27Zzdi57Pf7Zzvl8zufz+uzy/Lw/r/fr/XwprTWCIAhCeBIR6AAEQRAE/yEiLwiCEMaIyAuCIIQxIvKCIAhhjIi8IAhCGFMm0AHYU61aNV2nTp1AhyEIghBSbNq06ZjWurqzbUEl8nXq1GHjxo2BDkMQBCGkUErtd7VN0jWCIAhhjIi8IAhCGCMiLwiCEMYEVU7eGSaTiYMHD3L27NlAhyIIXlO+fHkuu+wyDAZDoEMRSjlBL/IHDx6kUqVK1KlTB6VUoMMRBI9orTl+/DgHDx6kbt26gQ5HKOUEvcifPXtWBF4IKZRSXHTRRRw9ejTQoYQ0ySkZTFu+g0NZOdSMMZLUqT6JTWMDHVbIEfQiD4jACyGH/M0Wj+SUDMZ+tpUcUx4AGVk5jP1sK4AIfSGRiVdBEIKOact32ATeSo4pj2nLdwQootBFRN4DWVlZzJo1K9BhsG/fPho1auRxnw8//ND2euPGjTz++OM+jaNOnTocO3aswPvLli1jypQpPj2XUHo5lJVTqPcF14jIe8CdyOfl5Tl9v6icP3++WJ93FPmEhATeeOON4oblFd26dWPMmDElci4hvElOySDCRbqrZoyxhKMJfUTkPTBmzBh2795NfHw8SUlJrFmzhnbt2tG7d2/i4uIKjLCnT5/OhAkTANi9ezedO3emefPmtG7dmvT09ALHnzBhAoMHD6Zjx47069ePvLw8kpKSuO6662jcuDFvv/12gc/s27eP1q1b06xZM5o1a8Yvv/xii/Xnn38mPj6eGTNmsGbNGrp27QpAZmYmiYmJNG7cmJYtW7Jlyxbb+QcMGEDbtm254oorbDeFM2fO0KVLF5o0aUKjRo1YvHix7fwzZ86kWbNmxMXF2a5p/vz5DBs2DID+/fszZMgQWrduzdVXX81XX31V3F+DEEYkp2TQasoq6o75mlZTVpGckpFv29jPtpLnpGOd0RBJUqf6JRlqWBASE69Whg8fTmpqqk+PGR8fz2uvveZy+5QpU9i2bZvtvGvWrGHDhg1s27aNunXrsm/fPpefHTx4MHPmzKFevXr89ttvPPLII6xatarAfps2bWLt2rUYjUbmzp1LdHQ0v//+O+fOnaNVq1Z07Ngx30RejRo1+P777ylfvjw7d+6kV69ebNy4kSlTpjB9+nSbqK5Zs8b2mfHjx9O0aVOSk5NZtWoV/fr1s11Teno6q1ev5tSpU9SvX5+hQ4fy3XffUbNmTb7++msATp48aTtWtWrV2Lx5M7NmzWL69Om8++67Ba5p3759/Pjjj+zevZt27dqxa9cuypcv7/JnJZQOPE2oOsvFA0QqxeTucTLpWgRCSuSDhRYtWnisfz59+jS//PIL99xzj+29c+fOOd23W7duGI2Wx9AVK1awZcsWPvnkE8Airjt37uTqq6+27W8ymRg2bBipqalERkby119/eYx57dq1fPrppwC0b9+e48eP24S7S5culCtXjnLlylGjRg0OHz5MXFwcTz31FKNHj6Zr1660bt3adqzu3bsD0Lx5cz777DOn5+vZsycRERHUq1ePK664gvT0dOLj4z3GKYQWhS1zdDWhOnJJGuA6527WWgS+iISUyLsbcZckFSpUsH1fpkwZzGaz7bV1Za7ZbCYmJsarJw/742mtmTlzJp06dcq3j/0Tw4wZM7j44otJS0vDbDZ7NUJ21rDd+nRQrlw523uRkZGcP3+eq6++mk2bNvHNN98wduxYOnbsyLhx4/Ltb93XGY4lhFJSGH4UpczRlYjnac3Yz7YSbTSQlWMqsF1y8UVHcvIeqFSpEqdOnXK5/eKLL+bIkSMcP36cc+fO2VIllStXpm7duixduhSwiGxaWprH83Xq1InZs2djMln+0P/66y/OnDmTb5+TJ09y6aWXEhERwQcffGCbAHYX680338yiRYsASxqnWrVqVK5c2WUchw4dIioqir59+/LUU0+xefNmj7Hbs3TpUsxmM7t372bPnj3Ury+51HCjKGWO7sQ6x5SHUpbcuz2Siy8eITWSDwQXXXQRrVq1olGjRtx222106dIl33aDwcC4ceO4/vrrqVu3Lg0aNLBtW7RoEUOHDmXSpEmYTCbuu+8+mjRp4vZ8AwcOZN++fTRr1gytNdWrVyc5OTnfPo888gg9evRg6dKltGvXzvYk0LhxY8qUKUOTJk3o378/TZs2tX1mwoQJPPjggzRu3JioqCgWLFjgNo6tW7eSlJREREQEBoOB2bNne/XzslK/fn3atGnD4cOHmTNnjuTjw5AMN2WOrtI4SZ3q5xv9O5KVbWLGvfGy0tWHKGeP8YEiISFBOzYN+fPPP7nmmmsCFJFQFPr370/Xrl25++67Ax1KQAnnv93klAxGLE7FmXrEGA2cO2/OJ+RGQ6Rt4jQ5JYORS9KcVtDExhhZN6a9HyMPT5RSm7TWCc62SbpGEIRCM235DqcCrwClcJvGSWwayys9m0hapoSQdI3gc+bPnx/oEAQ/42oCVWNJuXj6jDX9ImkZ/yMiLwhCoakZY3Sak4+9MLHqbJvjpGti09gCoi7Ok75HRF4QwhB7sYw2GlDKMsL2lXA6m0BVQLsG1UmoXbXANm9SMb52npQbhgUReUEIMxzF0r7u3FeWvYlNY9m4P5NF6w/YcvMa+HRTBgm1qzK5e5xXAmsvxBFKFZiMtebyvYnV/lgxUQZOnz2Pyax9et2hiIi8IIQZrqwBrNivMC2O4K1OP1pg8tUqyuvGtHcp6hOWbXe64MlZtQ1YBLrhuO/Izs1DY7E46HX95UxKjMt3XPsb2wkn8wKFuWGEE1JdEwAqVqwIWBYceSozfO2118jOzi7U8e2NyYqKNUZXOLpzenMthaVt27Y4ltRC0SyUc3JyaNOmjU+dQ3Nzc7n55puL7R7qa7yx47WuMLU3B/PVeRzftxqS1RnzNcMXpzoVeE+cuSDwYIl94foDPJu81bbd043NU8zhjIi8jyiKeNSsWdPmUeOKooi8NxRXmBxF3ptr8RVFsVCeN28e3bt3JzIy0vPOXlK2bFk6dOiQz6HTH7hzbXSGtxYA1pFtYY9v3d/VCptooyHfvmM/2+py4VRx+Oi3v23feyveQWmPkJEBAwfC1Kl+ObyIvAf27dtHgwYNeOCBB2jcuDF33323TXTr1KnD888/z0033cTSpUtdWgvv3buXG264geuuu47nnnsu37GtNsV5eXk89dRTxMXF0bhxY2bOnMkbb7zBoUOHaNeuHe3atQMsBmY33HADzZo145577uH06dMAfPfddzRo0ICbbrrJpWnY/Pnzueeee7jjjjvo2LEjANOmTbPZGo8fP77AZ06fPk2HDh1s1sJffPEFUNCC2f5azp49y4MPPkhcXBxNmzZl9erVtvN3796dzp07U69ePUaNGmW79v79+9OoUSPi4uKYMWOG7fxLly6lRYsWXH311fz8889A/ieVCRMmcP/999O+fXvq1avHO++84/TaFy1axJ133ml7/fLLLxMXF0eTJk1sPvht27Zl+PDh3HjjjTRq1IgNGzbYzjF9+nTbZxs1amTzEkpMTLTZRfgDe5HUWFIXSUvTaPr8CpeinNSpfoEadFdYc9X2x3c3wvdGtK02RdZFT96MsO2J9NLnyD694414B10d/smT8PTTUK8efPAB+GEwB6GWkx8+HHxsNUx8PHgwPtuxYwfvvfcerVq1YsCAAcyaNYunnnoKgPLly7N27VoAOnTo4NRa+IknnmDo0KH069ePt956y+k55s6dy969e0lJSaFMmTJkZmZStWpVXn31VVavXk21atU4duwYkyZNYuXKlVSoUIGpU6fy6quvMmrUKAYNGsSqVau46qqruPfee11ey6+//sqWLVuoWrUqK1asYOfOnWzYsAGtNd26deOnn37i5ptvtu1fvnx5Pv/8cypXrsyxY8do2bIl3bp1K2DBbG+gZr3GrVu3kp6eTseOHW1OmampqaSkpFCuXDnq16/PY489xpEjR8jIyGDbtm2A5SnByvnz59mwYQPffPMNEydOZOXKlQWuacuWLaxfv54zZ87QtGlTunTpQs2aNW3bc3Nz2bNnD3Xq1AHg22+/JTk5md9++42oqCgyMzNt+545c4ZffvmFn376iQEDBthickWjRo34/fff3e5THJylIUxmbcs5O5tQdKxBjzYa+PesCbOLobc7V0j7/LW7lar2ZGWbeDZ5a75JWW+xrox1lbe3x/5m4KzaxxCpqFC2DCdzfFdV5BNyc2H2bHjhBTh+HHr3hkmTwIOzbVEJLZEPEJdffjmtWrUCoG/fvrzxxhs2kbcKqjtr4XXr1tlsfu+//35Gjx5d4BwrV65kyJAhlClj+ZVUrVq1wD7r16/njz/+sMWSm5vLDTfcQHp6OnXr1qVevXq2GOfOnev0Wm699VbbsVesWMGKFStsHjenT59m586d+URea83TTz/NTz/9REREBBkZGRw+fNjtz2vt2rU89thjADRo0IDatWvbRL5Dhw5ER0cDcO2117J//34aNmzInj17eOyxx+jSpYvtKQPy2xq78u6/8847MRqNGI1G2rVrx4YNG0hMTLRtP3bsGDExMbbXK1eu5MEHHyQqKgrI/7Pu1asXYDF0+/fff/PdcJwRGRlJ2bJlOXXqFJUqVXK7b1HwJg3hbELRsQbdcWLSE3laM3xxKhv3ZzIpMa5Qoq2BhesPeHUeR3o0t3jKe5O373X95bbvQ2JxldkMS5bAM8/Anj3QoQO8/DI0a+bX04aWyAfIatidba7VHMyTtbAnq12ttVf73HrrrXz00Uf53k9NTfXaytfR1njs2LE8/PDDLvdftGgRR48eZdOmTRgMBurUqWOzU3YXpyuc2RpXqVKFtLQ0li9fzltvvcWSJUuYN29evv2LY2tsNBrzxezuZ+3sWK7spK2cO3fObwZsrhYdOeLpZuCs5NEbFq4/UGTBLixRhgg+3ZTh1Y2oQtnIfNU14HxxlT8oUv39qlUwahRs2gRNmsB330HHjv/ltvyI5OS94MCBA/z6668AfPTRR9x0000F9nFnLdyqVSs+/vhjAJf5244dOzJnzhybkFlTCPb2wS1btmTdunXs2rULgOzsbP766y8aNGjA3r172b17ty1Gb+jUqRPz5s2z5fUzMjI4cuRIvn1OnjxJjRo1MBgMrF69mv379xeIyxF7W+O//vqLAwcOuLUaPnbsGGazmR49evDCCy8U2tb4iy++4OzZsxw/fpw1a9Zw3XXX5dtepUoV8vLybOLcsWNH5s2bZ5tbsU/XWCdR165dS3R0NNHR0dSpU8cW0+bNm9m7d69t/+PHj1O9enUMhv8mG32Jt/l1b3LSzkoeg4kck9nrJ43sXN/2V3aGswlpZ3MkbquUtm6F22+3jNqPHIH/+z/YvBk6dSoRgQcfibxSap5S6ohSapvde1WVUt8rpXZe+FrFF+cKBNdccw0LFiygcePGZGZmMnToUKf7LVq0iPfee48mTZrQsGFD2yTl66+/zltvvcV1112Xr42ePQMHDqRWrVo0btyYJk2a2BpyDx48mNtuu4127dpRvXp15s+fT69evWy9WtPT0ylfvjxz586lS5cu3HTTTdSuXdur6+rYsSO9e/fmhhtuIC4ujrvvvruAcPfp04eNGzeSkJDAokWLbFbK9hbMSUlJ+T7zyCOPkJeXR1xcHPfeey/z58/PN4J3JCMjg7Zt2xIfH0///v2ZPHmyV/FbadGiBV26dKFly5Y899xz+fLx9tdqnTvp3Lkz3bp1IyEhgfj4+HyTqlWqVOHGG29kyJAhvPfeewD06NGDzMxM4uPjmT17dr4uXatXr+b2228vVLyFIbFpLJO7xxEbY0RhGe064u2Eoj8qXHxJYW5AEUp5XQ1UFFyJ+YRl273z0P/7b3jwQcuo/ddfYdo0+OsvuP9+iCjZsbVPrIaVUjcDp4H/01o3uvDey0Cm1nqKUmoMUEVrXTAZbUcwWg3v27ePrl27epyAEwLDhAkTqFixom2OxBUpKSm8+uqrfPDBBy73adu2LdOnTychwaljq1O6d+/O5MmTnT6p+Ppv11leXQF9WtYqkLpw9tnhi31ctBAkGCIV0+5u4tNUTaspqwp1U1TA3ildICsLpkyB11+35OAffxzGjgUnc2y+xJ3VsE9y8lrrn5RSdRzevhNoe+H7BcAawK3IC4K/aNq0Ke3atSMvL89ntfK5ubkkJib6vOuVq5yvs0obzX/14qvTj7rME7vr1hTqmPIsk8TTlu/w2WRrYRdN1a4YCa++Ci++CCdOQN++luoZL5+q/YnPmoZcEPmv7EbyWVrrGLvtJ7TWBVI2SqnBwGCAWrVqNbfmfK0EeiQvCEWlKH+7zkbr1rJCV006nGH9DFgEPthTNb7CvjmJI4WZMHU1kq8SZeCs3dyB0mbu3vEzE37/mAqH/rZMpk6dainNLkH8PpIvDlrrucBcsKRrAhyOIAQUV31TRy5Jc9nk2hk5pjwmfrk9nyAFIwqIiTI49ZopCq78aRxLQD0ZljmruzcaIhl/R0PA8nuqm/orz/08n/qHdkHTprBgHtxyi0+uw5f4cwbgsFLqUoALX4942F8QSj2u0gR5WnMm9zyGCO8rMk5km4Ja4MGSbvJ1B1Jn3jkLnZSOums67jjhHRtjtD0hRH41j1cm38XCxc9S32CChQth48agFHjw70h+GfAAMOXC1y/8eC5BKDYnsnM5fPIsuXlmykZGcHF0eapElS3RGNzVxZvyNFWiDPybc97jqtNQoiiGZe5wLCedsGy7y33d5d4d6+5//+QT/q/ZPfQFsoBjY8dSbfx4cFM5Fgz4qoTyI+BXoL5S6qBS6iEs4n6rUmoncOuF14IQlJzIziXjRA65eZZFT7l5ZjJO5HAiO7fEYkhOySDzzDm3+5zINtHr+sspmQrr0MO+nNRa5+7uJuLN+oJ9KSlMU4q4e+6hJ/AycCo1lWovvRT0Ag8+EnmtdS+t9aVaa4PW+jKt9Xta6+Na6w5a63oXvmZ6PlLp4sYbbyzSNnc4mmm5orBWwv5i3LhxTv1o7FmzZg2//PKLX+M4fPIsZofRsVlrDp90v7rXyonsXNL/+ZctB7NI/+ffQt8cklMySFqaRo7J7HHfD4vgCVMaUFhSMMMXp3Ltc9/y5JJUjxPO7tYXZP3vf0ysXJnoZs0YCXwEbPvsM8ZoTe0mTXwauz8J+MRracaZcFlL/Pwtap6wivwjjzzi1/M8//zzHvdZs2YNFStWLPKNzxusI3hv37fH+hRgvUnk5pn5OzObMy6E3lmVx7TlO2xdjDzhOaLSif1PL9uLm6Urzufm8kqzZty3fTvjgW8A/dJLPDh2bHFDDAhhZ2tQWG9sb1i4cCEtWrQgPj6ehx9+2OYdX7FiRUaPHk3z5s255ZZb2LBhA23btuWKK65g2bJlgMVe984776Rz587Ur1+fiRMn2o5rHU2vWbOGdu3a0bt3b+Li4vJtA+e2uO+88w7XXXcdTZo0oUePHh49513ZHXtrJexqP0cqVqzIyJEjadasGR06dODo0aOAxV+nZcuWNG7cmLvuuosTJ04A0L9/f5sPfZ06dRg/frztHOnp6ezbt485c+YwY8YM4uPj+fnnn1m6dCmNGjWiSZMm+czUikPZSOf/Cq7et8fZUwDA6XN5xE9c4dWy+NJS4hhs5Jt41Zq5d9/NtnLlGL19O0eBeX36cLvWdAlRgYcwE/lC+0p4wZ9//snixYtZt24dqampREZG2nxZzpw5Q9u2bdm0aROVKlXi2Wef5fvvv+fzzz9n3LhxtmNs2LCBRYsWkZqaytKlS512O9qwYQMvvvgif/zxR7737W1x09LSbB7s3bt35/fffyctLY1rrrnGtgTfFVa7499//51LLrnE9r7VSnjz5s2sXr2akSNHorVmypQpXHnllaSmpjJt2jSX+zly5swZmjVrxubNm2nTpo3tptavXz+mTp3Kli1biIuLy3ezs6datWps3ryZoUOHMn36dOrUqcOQIUMYMWIEqamptG7dmueff57ly5eTlpZmu5kWl4ujyxPh4CUSoRQXR/9nPOYqJeNutJ+VY/JqWby3HuqCb7FOvH4xfjwrIyIY/OmnVAamN29O09xcBixcGNgAfUBYpWtc1RgXp6/jDz/8wKZNm2ymVzk5OdSoUQOwdAbq3LkzAHFxcZQrVw6DwUBcXFw+W9xbb72Viy66CLCI89q1awssnW/RogV1nfhJu7LF3bZtG88++yxZWVmcPn2aTp06ub0OV3bH3loJu9rP/oYBEBERYbNf7tu3L927d+fkyZNkZWXRpk0bAB544IF8lsz22FsLu2p+0qpVK/r370/Pnj1t+/sCpbA985eJUFwaY7RV1zhLyWScsAhE2cgIr9I6OaY8lyWN4VQtE0pcve8XFqmu9AGOASONlVn20GxqxF7CVduOBJdVcREJK5H3tudkYdBa88ADDzg1zTIYDDZr2oiICJsJV0RERD5bXE9WuJDfAtjx/M7279+/P8nJyTRp0oT58+ezZs0aj9fi7DjeWgkXxXLY1Tnd4Y218Jw5c/jtt9/4+uuviY+PJzU11XYTLQqOAg4UaLDhbmL24ujy/J3pn64+gn+oeGgHAz4YyaNY5jheApYMnk1WFYtHvafFUqFEWKVrXJVDFaevY4cOHfjkk09sFryZmZk4Wi944vvvvyczM5OcnBySk5NtTT+8wZUt7qlTp7j00ksxmUxetZ9zZXfsrZWwq/0cMZvNthz7hx9+yE033UR0dDRVqlSxte/74IMPbKN6b3CMZffu3Vx//fU8//zzVKtWjb///tvNpz3jTWWNu4nZKlFluahCydbTC0XDcOo4973SnbUfjOQJYCGQvmwZX0/+wSbwVtwtlgolwmok72opcnH6Ol577bVMmjSJjh07YjabMRgMvPXWW17b+QLcdNNN3H///ezatYvevXsXyuWwc+fOpKamkpCQQNmyZbn99tt56aWXeOGFF7j++uupXbs2cXFxLr3drbz++uv07t2b119/nR49etje79OnD3fccYfNdteZlfBtt93G6NGjne7nSIUKFdi+fTvNmzcnOjra5s++YMEChgwZQnZ2NldccQXvv/++1z+DO+64g7vvvpsvvviCmTNnMmPGDHbu3InWmg4dOtCkmOVs3lTWuErJWCdmo8pZ2syd97JCRihZlOkct8x/nHGZGVwOfAlMbzeAH1dZ5rIOrfva6eeKkwUIFnxmUOYLfGE1XKSuLX5k/vz5bNy4kTfffDNgMZQkFStWtDUhCRXS//nXqYArFJdVteTlnaV0IpQitorlKdFxG8DhA3sYtOwf/wYvuEWbzSR8PonxuzbQGNgATIy7he23Dwdg35QugGtDstgYI+vGtC/BiItGUBuU+ZqSagEmhA8XR5d3KtIabZtcrRJVluxz58k8Y0JfmJ01a83BzBzbayG4uPLHBTyzfintgV1A35r1+bnPy6gIi9V0rF0a1x9ZgGAh7EQ+2Ojfvz/9+/cPdBglRqiN4gFbBY0zwbbPzZ/INhXYLgIffNTYtooRX79KL+Ao8ERUNF8NnEOesZLNDsJRwEOiEXgRCQmR96bJtSAUhypRZV1WyOTmmV0ueHKF1lpuACVMxYw/GbQwiaGACXgBWPTQbM5Ws3j9WCtkI5WiR/OCT/zhmgUIepEvX748x48f56KLLhKhF/yKu8lVb+rgrWitOZ/9L/uzfOuuKDjHcPIIPd8ZzKi881QE3gNm3TORrCua2/axv93mac2nmzJIqF01LEXdkaAX+csuu4yDBw/alscLgr/Izj1PVrYpX418hLI0tTiZbSLPy4G5RrM/y8TM3074J1ABAGU6S6d5j/Fc1j/EYvEyf7n9IDKuu9PjZ4u7SDKUCHqRNxgMTleCCoI/sFZnZWTlEKkUeVoTYzRw6tx58qQ8MijQZjMtPpnI+L2baITF43xgk0782fmxQh0nHMojvSHoRV4Q/IW7clv7SgtfN7UQik69Ve/y7O/JtAH+AnrHXsu63pNtFTOFoTiLJEMJEXmhVOLYMDsjK4fhi1NtXYSCvW1eaePiLSsY+e0b9AQOA8MqVOWbgbMwl6/osYGKIQLKREaGZXmkN4jIC2GFs3RLrMMoPTklgxGLU53WvsioPbiodGAbgz8awxDgHDAB+Oax9zgSdbHXxzCZ4d4WsaxOPxp25ZHeICIvhA2Oo3Ors6O92dTG/ZksXH8gYDEK3lH25GF6vj2I0dqMEXgHmNNzEifrxtO3Za1C/w5Xpx8NiZWr/kBEXggbnFlNW8kx5TFh2XZOykg9qInIzaHze4/y7L9HqAl8Bky7dSj/NLPYD1QoG8mnmwrfH6K0TLI6Q0ReCBs8dVeSVEzwos1mWi55jgn707gGWAs82PR2dnTM334y25SHszVpVosCV38DniZZg83zypeIyAthgzUHL4QW9VfO5dlNy2gNpAP31WrMr/dOQkUUdEJ39es9lJXDjHvjSVqa5rRXbrsG1V2e39kkfLh4yYOIvBBGiMCHFpekfkfS8jfpAfwDPFqpGt8OnI25rNFjxYwj0UaDXfnrFnIcGnm7W+Hqj45ywURYNQ0RSjexpaTuOdSpvD+NUVO7snb5m3QExilF28Hv8PUj8zGXLdrv0Op4ktg0lqoVyhXY7q4BiD86ygUTMpIXwgZndrFC8FAu8x/ue2cQSUA5YA4wt9dL/FurcbGPfSL7v/mWwop2zRij01x+uCyWkpG8EBZYJ85yTHlEipFdUBFx9gxdZj3Aj+8MYiLwLdC60zCmj/7KJwJvJTnFUnVT2DagSZ3qYzTkXzEbToulZCQvBA3eVjg4W/Bkj/MBUJgAACAASURBVOTmgwNtNnPjx08z4e9t1Ad+AvoldOOvDoP9cr6JX263/V1YbYWtuBPtcPaSBxF5IUjwpsIhOSWDCcu25yuFFEEPThqsmMVzKd/QCvgD6Fknnt/ued5pxYyvOJFtsqVtNP/5xzuueHZGuHrJg4i8UMK4Gq17qnBwvAkIwUnNzV8x6vs5JAKHgCHRF7NiwFuYy5YvdMVMcZHbvwURecHv2KdX7B+j7UfrribFMrJyXDZZFoKHmL2bGbpkHA8B2cCzEZEsGTyX3GjvPWb8RbjVvRcWEXnBrziOwB1HV9bRuqsKB/C8klUIHOWOH6TPu0MYCZQFZgFze0/l9OUNAxxZfsKp7r2wiMgLfsWdn4wV62pFSceEDhFnT9P13aE8e+YENYCPgVdve4IjjW8NWEyOk62OWJ8Ww9nCwBki8oJf8WZBSc0Yo9MKBxnBBx/anMdNi0Yz4VA69YDVQO8W3dnVbkCgQ7NNsrrzrwl3CwNnSJ284Fc8LShR/Jd3B1g3pj17p3QJmxrlcKLht2+weNqdLDqUzjmgxxXX0T9pWVAIPFgEft2Y9rx2b7zLund3E/zhit9H8kqpfcApIA84r7VO8Pc5heAhqVN9lw06wPkkLJDveyGwxG5IZvTqd+kGHAQGV6nJ9w++gTaUfMWMK+zr4N3VvY9YnOr08+FiYeCMkkrXtNNaHyuhcwlBRGLTWDbuz2TR+gP5hN5Z/jTHlMfEL7cTVbaM5OaDgJjdG3n0kwkMAE4DYyMNfDJ4LqbKrh0dA8Xk7nH50i2u6t4La2EQDvl7yckLfmdSYhwJtat6lW+3X9AiBAbj0f30nfcoI7AIxBvAe32nczq2QYAjc06s3ZyOJ5z5G7laDRsu+fuSEHkNrFBKaeBtrfVc+41KqcHAYIBatWqVQDhCILCOrKwjIyH4iMz5l27vDOHpnH+pDnwIvNblSY40Ct62eYX1mCmMhUG4WBCXhMi30lofUkrVAL5XSqVrrX+ybrwg+nMBEhISZJFaGCOrVoMTbc6jzcIkxv/zF1cCPwAvtuzJnjb9Ah2aW2KMBiZ0a1howfXWwiBcLIj9LvJa60MXvh5RSn0OtMDiVSSUMrypmRdKlkZfz2Dcth9oAaQB3a9qyaa7nvarx4yvOHfe0hjEX3nzcLEg9qvIK6UqABFa61MXvu8IPO/PcwrBS6iNgMKZWus/YfSP8+kCHAAGVr2MH/q/jjaUC5qKGStGQ0SBTk/w30T9WZPZL3nzwuTvgxl/j+QvBj5XFn/vMsCHWuvv/HxOIUiRBU6Bp+pf6xn2+ST6A/8Co8qU4/OH38FUsWqAI3NNbp7rLK6zSXpf5c3DxYLYryKvtd4DNPHnOYTgw9Xjs7uaeWnC7V+ijuzj/veHMRzLCshXgff7zeDMpfUCHJl7FJDnpDG3J3z11BgOFsRSQin4FE9lZ8NdLEbJ09rtknShaJQ5c5Ju7zzMM+dOcxHwAfDGHaM4eu3NgQ7NKzzJe4zRkK+/gJVQy5v7k+CfXRFCCk/Lxl0121ZAuwbViQi2hHCokmei7fuP8/2bfXj13Gk2Ae1u7MVzo78KGYH3hAImdGsY1q37fIGM5AWf4q7sLDklg8wz55xu18DC9Qf8GFnpofGX0xj/x480B1KAxPqtSE0cG+iwioQ7Z0lN+OTN/YmIvOBTXE2uxkQZLqRxClZJCL6h9q+LGfvTB3QG9gEPVq/N6gdeg0hDgCMrOhrXKRnrU2E45M39iaRrBJ+S1Ml553utkRp5P3HRjrVMnNqV1T99QAsgqayRW4d9wOoBb4W0wINFyCUlUzxkJC/4FFePz64mXIWiE3V4N/3nP8HjF15PB+b3f53si68MZFg+xT71IimZoiEiL3iFt6sKne0n+BbDmRMkzh3M2NwcYoD/A2YmjuV4/VaBDs2nxBgNtr8xSckUHRF5wS3JKRlMWLY9X07UsSzSXaPuEYtTg24FZciSZ6LdguFMOLqf2sC3wJTW97P/xnsDHZnPsVbOCMVHRF5wiTtDMfuySHeNurWT94TCE588mfE71tEU2AQ8ds3NbOk2KtBh+Q37yhmheIjICy7xZChmGdFvkYoZP1L350WM/eUjOgJ7gQdqXMGP/V4J+QlVTyh5/PMZIvKCS7xZGi4C7x+q//ETT3z5Mn2B48CT5SqybNDbnK8QHejQSgRxuPAdIvKCS8RQrOSp8M9OBvzfCIYBZmAy8MGDb5Jdo05gAxNCFhH5UoirShnH99s1qM6nmzKkvr0EMJzOpPvbgxhz/hzRwHzgzbueJfPqlgGOLDBUiQrvdFRJIiJfynBlILZxf2Y+Qc/IyuHTTRn0aB7LV2n/OF1xKBQfZTpHhwVPMO74QWoBXwFT2/Tn75Z3Bzq0gGGIVIy/QyprfIWIfCnDlYHYR7/9XcDqN8eUx1dp/9g68Ai+Q5vNNP/8JSbsWk9jYAPwSKP2bOvyZKBDCyixstDJ54jIlzJcTaa68nKXEbzvueLH/+OZ9UvoAOwG7r/0an7q8zIqUv4d140J3qbhoYr8VZUyXE2mStMO/1N9+2pGfPUKvYGjwHBjZb4cNIc8Y2VZMIZrG2qheIjIlzKc9a1UQMsrqrD5wEmZZPUDFTPSeWjhUzwKnAcmAQsfms3ZapcHOLLgwRCpxALDT4jIlzISm8aycX8mi9YfsK1E1cDmAyfp0TxWPN19iOHfo9wzdzCj80xUBOYBs+4ex4krWwQ6tKCjQtkykof3EyLypZDV6UcLWA3kmPJYnX6UCmUjOZMro/nioExn6fj+Yzx34h8uA5YBU9sNJKNFYqBDC1pOytyP3xCRL4W4mnyVhU/FQ5vNXPfp84zfs5E4YD0wqHFH/rztcU8fLfVIT1b/ISJfCpGVrL7nqjXv88xvn9IO2An0qdmAtX2moiIiPX201CMNQPyLiHwpJKlTfUYsThV3SB9QY+tKRn7zGvcCR4DHK1Thq4GzMZevKBUzXjK5e5zk4/2IiHwpJLFprHRqKiYV/97O4A9HMxTIBSYCHw6cw7mLLgtwZMGH0RBJeUMEJ7Kd92kVgfcvIvKllFhJ2RSJsicP03PuYEaZ84gC3gXm9HyerLrNAh1aUBJr1x3MsXRX0jQlg4h8KUX6rhaOiNyzdJz3KONOHqYm8Dkw7ZaHOdT8jkCHFvSMWJxKzRgjzWpFs37PCfK0JlIpejSXln4lgYh8KSWxaSwTv9zu9BFa+A9tNtNyyTjG70/lWuAXYEDT20jv+GigQwsJrE+LGVk5+Z4c87Tm000ZJNSuKkLvZ0TkQwRvG2kXhi6NL5XFT26o98O7PLcxmZuBHUCvyxvxy30voSIiAh1aWGBtISki719E5EMAV/bA4LwPpjc3hOSUDD7dlOH/4EOQS1K/46nlb3I38D/g0UrV+HbAW5jLV5CKGR/jTfcxoXiIyIcAruyBnY2CvLkhJKdkMHJJmhiSOVD5wBYe/uhpBgPngPEoPh40l3NVLw10aGGLLILyPyLyIYCr0Y6z9z3dEKw3ARH4/zAc3UevecNIAozAXGDOfS/yb+0mAY4sNHDlYKrA7VoMqa4pGUTkQwBXK1SdjYI83RCc3QRKK+rsaTq+fh8TgUuAT4DpHR/lf01vC3BkoYPREEmP5rEF2kRa31+dfjRfO0n719IcpGTwu8grpToDrwORwLta6yn+Pme44cwe2NUoyNUNIUIp6o75Wla5AtqcR7MZPZl8/hwNgJ+BXrWbsPe+FwMdWshRrkwECbWrklC7qs8LAwTfoLQfH9uVUpHAX8CtwEHgd6CX1voPZ/snJCTojRs3+i2eUMbb6hrHnLyQnys+HMMLf2+jFfAnMLaskZQnFkvFTDEwGiLFmiDAKKU2aa0TnG3z90i+BbBLa73nQiAfA3cCTkVecE1i0/wLR5JTMmg1ZZVT0S9viLCJvKe8aGmh+up5jN3wGd2Bf4BBwMonFqOlYqbYSClkcONvkY8F/rZ7fRC43n4HpdRgYDBArVq1/BxOeJCckkHS0jRMZot8Z2Tl8OTiVJ7+bAvZpvxNt0u7wFf8Yw2PfjmdQUA28Cyw+KFZmKrJ35ovkVLI4MXfIu9skJRPd7TWc7EUNJCQkFDaNckrJizbbhN4K2YoIPClmbL/202fBU8wEigHzAbeTnya0/VvDHBk4YmUQgYv/hb5g4B9I8vLgEN+PmfYkyVddFwSkf0vnWf2ZjxwMbAEmHrjfRxv3TfAkYUHra6sWqAXsJRCBjf+FvnfgXpKqbpABnAf0NvP5xRKIdqcR8vZD/HC6WNcDawBetZtzv6eEwMcWfhQJcrAL7sziYkyUK5MBCdzTFJJEwL4VeS11ueVUsOA5VhKKOdprbf785ylgSpRBjEWs+Pa72byXNpybgC2Ad2MlUkbtlAqZnyM9W/uRLYJoyGSGffGi7iHAH6vk9dafwN84+/zlCbG39GQpE/SMOWV7imM2I3LGPXDXO7E8pj4cMylrBgwE20oLxUzfkYqakIHWfEagiQ2jWXj/kwWrT9QKqtnYvZs4pGl43kIOA08HVmGpYPfwVS5eqBDK1VIRU1oICIfoqxOP1rqBL788b/p++5QnsTyhzsTeLfvNE7HXhPgyMKX2Bgj2bnnnaYHrRU1/rDBFnyHiHwIYf/PVJoEPjLnFHe8O4Rnsk9SHfgImNHlSY40ah/o0MKSWAehdraK2lpRU1gbbKHkEZEPEUqjXYE253HzwlGM/2cHVwGrgBdb3sPuNg8EOrSwJTbGyLox+W+eVrF2NlpvNWWV1zbYQmAQkQ8RSpt7ZMNvXmPc1pVcD2wBelzVgo13PSsVM37EXb27va2G9YlyxOJUl0+Ukq8PHkTkQwRnzpLhyGW/fcboNfO4A4sfxqCql7Gy/+toQzmpmPEjCrwyGfP2iVJWwAYPIvIhgqvGDOFC1Z2/8ehnL9AfOAWMKVOWTwe/g6nSRQGOrHSg8S6H7s0TpayADS5E5EOEcBV449F93D9vGCOACOA14P37X+F0TRGJYMRdGkaBVNcEISLyIUK4jeTLnDlJt3eH8PTZU1QDFgKv35HE0WvbBDq0UkmEl7kwV01pnE3YCsGBiHwQ4a7eOGwEPs/EzQuTmPC/XVwBfA+8dMN97L1ZDMQCidnLP6/CdCkTggMR+SDBU71xrIsRVCjR+MvpjPtjDQlAKnDX1TeyOXEsSsmUaqCJ9XKi1F05pRCciMgHCc4mtOzrjds1qM7C9QcCFF3xqLV+KWN+XMDtwH5gQLXarOr/GkQapGImCCjsSNyxS5kQ3IjIBwmuJrQOZeWQnJLBp5sySjii4lN1x688nvwi/YCTwKiyRj4f9DamilUDHZpwAcfVrUL4ISIfJLia0KoZYwy5hVBRR/bywPuP8cSF168A8/u9xplLrwpkWIId0ny79CDLB4OEpE71MRoi871nfYwOldWDZc6c5O4ZPVn7/mMkAYuBm+4czVujvxKBDzAxRgOxMUYUltG7CHzpQUbyQYK7Ca1py3cE96Rrnol2C0Yw/ug+6gDfAZNb92X/jfcFODABLPXrXZtcyqTEuECHIgQAEfkgwtWElrOytWChyRdTGZ/+M82AzcATDVqTdufoQIcl2KGBTzdlkFC7qozeSyEi8kGE2zp5sznA0eWnzrqPGLt2EZ2AvUD/Glewpt8rEGkIdGiCE8QZsvQiIh8kOKuTT/okjbGfbSHHFDwCXy39Zx7/Yip9gRPAyHIV+GLQXM5XiA50aIIHQmVuR/AtIvJBgrMKGlOeDpo+rhX+2cWD/zecYVge/18G/u/BmWTXqBvgyARvEWfI0omIfBCQnJIRtBOrhtOZ3DV3MGNNZ4kGFgAzE58hs/4NgQ6t1GOIVF4PAsR6oPQiIh9grGmaYEOZztF+wQjGHT9AbeBrYGqbBzjQ8p5Ah1aqiVQKs9a2OZuJX2532n+1QtlIYqLKivWAICJfUthPqkYbDSgFWdkmIoLMXVKbzTT7YjLj//qVeOB3YFjDdmztOjLQoQmAWWv2TumS772kT9LyjegNkYoX75I6eMGCiHwJ4DipmpXz38jLncAbIhQmb+0BfUDdnz7gmV8XcwuwB7j/0nr83OdlqZgJIhzz6mIYJnhCRL4EKIotgdVTpCQWQtXYvoYnvppOH+AYMMJYmWWD5pBnrOzX8wqFw1VeXQzDBHeIyJcAhS1ds/4zW/9xn1yS6rXfd2GoeGgHAz4YyaNAHvAisHDAW+RUr+37kwnFIlIpsSIQioSIfAngynzMGY7/zE9/tsXnAm84dZy73x7I6DwTlYD3gbd6jOPEVS18eyLBZ5i1LiDwySkZ+SZeY4wGJnRrKDcCIR8i8n4mOSWD7NzzXu+fpzXTlu+wvc724UIoZTrLre8/wXMnMrgcWAa83HYAB6/v7rNzCP7BMRefnJJRYMI1K8dE0tI0wLum3ELpQETejzj7R/SGjKwchi9OxVcNk7TZzHWfvcD43b8TB/wGPBx3C3/cPtw3JxD8irNc/LTlO5z+XZnMWuwLhHyIyPuRiV9uL9aKVV9UVl65Zj7P/PYJ7YGdQN+aDfi5z1RURKSnjwpBgrNcvLt5HrEvEOwRkfcjzhaplBQ1tv7AiG9m0As4AjwRFcNXA2eTZ6wkLfdCiCpRBqejcnfzPGJfINgjIh9mVPr7DwZ+OIqhgAl4AVg0cDZnL7o8wJEJhcUQqRh/R0On25I61XeaCjREKLEvEPLhN5FXSk0ABgFHL7z1tNb6G3+dLxiJMRryLXzyJ4aTR+g5dzCjzeeJAt4DZt0zkZNXNC+R8wtFJ8ZoIHV8R7dW045Y35fqGsETSvtpSf0FkT+ttZ7u7WcSEhL0xo0b/RJPIEhOySBpaZpfV61G5J6l47xhPHfyf8QCycC0DoPJSOjmt3MKRUdhcfG0Ir1WBV+glNqktU5wtk3SNX7E2ZLz7NzzPsnVa7OZFp9MYMLezTQEfgUGxnfmz07Din1soehEGSLclr1qLKuZxYJAKCn8LfLDlFL9gI3ASK31CT+fL+hwXHJed8zXxT5mvR/e5dmNybQBdgC9L2/EuvtelIqZABNliKBKhXJku6luiY0xsm5M+xKMSijtFEvklVIrgUucbHoGmI1l3k9f+PoKMMDJMQYDgwFq1apVnHBCgsKsfnXk4rTljPxuJj2Bw8CwilX55qHZmMtXkIqZIMBk1h5/t9m550lOyZDRu1Bi+C0nn+8kStUBvtJaN3K3X7jl5J3h6EjpDZUObOXhj8byMHAOmA58NOhtzlUVoQg2Ir2wjpY8vOBrApKTV0pdqrX+58LLu4Bt/jpXKOGYp4+JMrjM0ZfN+h/3zh3MKG3GCMwF3u45iZN140suYKFQ5GmN0RDp9iYuTbWFksSfOfmXlVLxWNI1+4CH/XiukMI+T5+cksHIpWnk2VXgROTm0Pm9R3ju36NcCnwKTL91KP806+L8gELQYG8RfSgrB1djelmVKpQUfhN5rfX9/jp2ODFt+Q6bwGuzmZaLn2XCgS1cA6wFHmzWhfRbhwY0RsE77C2irTfxVlNWOc3Ty6pUoaSQEsoAYxWA+t/P4bnNX3ETkA7cW7sx63tOQkVEBDQ+wT3WuvdYF+WQSZ3qF5iDkabaQkkiIh9gLkn5hqQVs+gB/AM8Urk63z00C3NZo1TMhAAz7o0nsWmsbbXqiMWp+erfpT2fEGhKpLrGW0pDdY2V7atW8WOHDgwGcoCXVQSLB88lN8ZZRaoQjCj+K4mVlaxCIJEVr37End+I47Z2Darz07pU2k/vx1NYjH3mAG/3msypWnGBvAyhCGj+S7c5DpWkgkYIFkTki4FjzXtGVg5jP9tq226/7Z//HSPrxXtYcjqTS4AlwMwuj1Przvs4tTszANEL/kYqaIRgQES+GExbvqNAPbR1BGf9XpvzaPXR04w/uJ36wI/A/QmJ7OwwEADz8Rz6tqzFR7/9TZ7WRCpFmUjFufO+a/snBAapoBGCARH5YuBqpGZ9v8HytxiX+i03AtuBnnWa8ds9E/JVzGRk5bBw/QFijAa6NrmU1elHi2x7IAQPUkEjBAsi8sXAlQ9Ngz++4+Ev3yQRyACGRF/CigFvYi5b3mXFTFaOiYXrD/gz3FKPdTIUKJS1hCFSce91l7M6/Wi++ZVPN2XkO4anckpBCAQi8sXAsQY6es8mHlk6noeAbOCZiDIsGTwXU3SNgMYpWNrojb8jf0ONact3ePfUpCGhdlUmJeafHE+oXVVKI4WgR0S+GFj/oWcsXEmHV/szEjAAbwG3rV3LdVF1WHNBBCK8MK4S/MdZB493aw27N4ZxJrN2WinjaCMtCMGIiHwxOHPiBFtva87Sw4epAXwMxL7/Pk/07w/A1ZCvnLKw7pP2xBgNnDp7Xm4URcRVSaPjYiXxmhHCDVkzXwTMeXnMuPFGDlWtynOHD7Md+GTUKO7TmtYXBN6RxKaxTO4eR2wRKi6MhkgmdGvIKz2byCrYYuBMqB3XMlSJMjj9rFTKCKGKjOQLyYKHHqL+vHmMALYCs7t2ZcgXX3jlMWP/eB8/cYVXTb4jlcq3cnLj/kyZoC0i0cb8Au5snYMhQmGIVJjy/hvTS6WMEMrISN5Lvpsxg2VK8cC8eVwGTLn6aupnZzP0yy+LZCI2oVtDjAbP7foilGbE4lRaTVlFckoGkxLjKBsp4/micOZCVyYrztY5mMyaCmXLEBtj8Q6KjTGKPYEQ0oh3jQfSvv2WDbffzgDgNPBa+fI8tmsXVWOL/0+fnJLBhGXbvRrRW6lXowI7j5wp9rlLK/Y9VuuO+dppDl4Be6eId78QOrjzrpGRvCtOnoRnnqHe7bfTD3gdOLp+PeNzcnwi8GBJ30zo1pCIQgzMReDzYzRE8tq98bx2b7xtviNSuf6B2uflXeXZJf8uhBMi8o7k5sIbb8CVV8JLL3GiTRu2LF7Mk1pz1fXX+/x005bvwBw8D1MhRZUog21xk3XyNDbGyCs9m7ic4LYX8KRO9QukzCT/LoQbMvFqxWyGpUvh6adhzx7o0AGmTiW2eXOcjdvduU8WBinNKzz2C5tcmcT1aB5bYEWqo4CL17tQGhCRB1i9GkaNgo0boXFj+O476NgRXDz2u3OfLKxAuLJGEPKzz0WO3JVJ3Or0o0zuHudRwGVBkxDulG6R37oVxoyBb76Byy+HBQugTx+IdF/14s59srCC4aw9nNUDxRABJjGjdIs7kzgRcEEorSJ/8CCMGwfz50N0NLz8Mjz2GJQv79XHPblPFgZPKYM6Y74u9DFLigigJO5B7iZSXT0JyeSpIFgoXSKflQVTp8Jrr1ly8E8+acnBV61aqMMURli8yd27G3FWiTJwItv7EsuSRDv2vPMTva6/3OU2aZQtCO4pHSJ/7hzMmgWTJsGJE9C3L7zwAtSuXaTDeSssnnL31htARlYOkRcMzBxtas8W0eumJPD3EotIpeh1/eUF3B/tkclTQXBPeIu82QwffwzPPAP79lkmU6dOhfj4Yh3WW2Hx1DnK/gZgNR5zvBHklNKkvP2iJU9I7l0QXBO+K15/+MFSMbN5s0XUX34Zbr3VN8f2EncrKj1V1UQqhVnrksiGBASjIZIezWNd+vDIqlNB8J7SteI1LQ06d4ZbboHjx2HhQti0qcQFHtyvqPQ0SZsXxgIP/5U5erNoSRCEohM+In/gADzwADRtChs2wCuvQHq6pSSyCAZivsDdisriiFhkhCLKULxrclOwUmIcysqRVaeC4GfCIyf/00+WfDtAUpKl9r1KlcDGhOfcfWGbiFjTPNZjNH1+RZErb4qTpTMaIjifZy52DX/NGKNMnAqCnwkPkW/RAoYOhREjoFatQEeTD1eTgvbiZl9dE+mmTaCjAGYFqLSyMJPBCouP+5nc8y492mXiVBD8R3iIfPnyMGNGQEPwphzSEWfi5q5NoGPlTUwRa+iNhkjKGyKcftZatrj497/ziXJRsK+Q8ZXXjyAIhSM8RD7AOAqzq3JIb3Ac4TtiX4J5+uz5QscaYzQwoVtDoGC6yGiItPm9FFfgnZmBiagLQskTPhOvAcRZPbwVe1H2lsSmsawb095lP9dDWTkWIXbiURxjNBBjdN6nFODcebPtHNaes44dkIpiz1AlyiDdlAQhCJGRvBu8TTF4EsWi2gm7qqWPUMpljf3JHBMz7o13mfLJMeUxfHEq05bvIKlTfacLjgrrjGk0RNqsfwVBCC6KNZJXSt2jlNqulDIrpRIcto1VSu1SSu1QSnUqXpgljzUFk5GVg+a/1It9j1Arnsohi1ou6ay8ECzpIFejfGvFirWZhivcXY+n8sW+LWvJqF0QQoTipmu2Ad2Bn+zfVEpdC9wHNAQ6A7OUUp67VgcRniwJ7HElxlC8mm+rWDtzYdRQQOgdK1ZcLTSy4up6EpvGUiXKeconNsbIpMQ41o1pz94pXVg3pr0IvCAEMcUSea31n1prZwnnO4GPtdbntNZ7gV1Ai+Kcq6QpjJ2wfX4b/rPG9cUoN7FpLGYXJZX6wjlcjajd3XysuLrO8Xc0LPYipeSUDFpNWUXdMV/Tasoqp08NgiD4F3/l5GOB9XavD154rwBKqcHAYIBaQVTjXlifcn9Wj7iKxZOJlzWeiV9ud1lqGe1ikra4i5R82T1LEISi41HklVIrgUucbHpGa/2Fq485ec/pcFRrPReYCxaDMk/xlBQl6VPuaYK3uLGcdbN46UzueZJTMlwu2CqqIPuye5YgCEXHo8hrrW8pwnEPAvadHi4DDhXhOAGjpJbbezPiLU4s7so7AUx52i/C68vuWYIgFB1/pWuWAR8qpV4FagL1gA1+OpffKIkFPN6OeIsaizei6g/hlbZ8ghAcFLeE8i6l1EHgBuBrpdRyAK31dmAJ8AfwHfCo1jp4WxwFkMKOeAs7memNqPpDeMVdUhCCg+JW13yutb5Ma11Oa32x1rqT3bYXtdZXaq3ra62/LX6o4Yk7z3lHClO7EiqMNQAAA3pJREFUb8VThY2/hNfdilpBEEoOWfEaYAozqVqUyUzHfH5MlAGtLStj/W0UJn41ghB4ROQDTGEmVYs6mSliKwilFxH5IMBbEZbJTEEQCou4UIYQMpkpCEJhkZF8CCGt8gRBKCwi8iGG5NcFQSgMkq4RBEEIY0TkBUEQwhgReUEQhDBGRF4QBCGMEZEXBEEIY5R20XUoECiljgL7Ax1HEagGHAt0ECWMXHPpoLRdc6heb22tdXVnG4JK5EMVpdRGrXWC5z3DB7nm0kFpu+ZwvF5J1wiCIIQxIvKCIAhhjIi8b5gb6AACgFxz6aC0XXPYXa/k5AVBEMIYGckLgiCEMSLygiAIYYyIvI9RSj2llNJKqWqBjsXfKKWmKaXSlVJblFKfK6ViAh2TP1BKdVZK7VBK7VJKjQl0PP5GKXW5Umq1UupPpdR2pdQTgY6ppFBKRSqlUpRSXwU6Fl8hIu9DlFKXA7cCBwIdSwnxPdBIa90Y+AsYG+B4fI5SKhJ4C7gNuBbopZS6NrBR+Z3zwEit9TVAS+DRUnDNVp4A/gx0EL5ERN63zABGAaViNltrvUJrff7Cy/XAZYGMx0+0AHZprfdorXOBj4E7AxyTX9Fa/6O13nzh+1NYRC/smxgopS4DugDvBjoWXyIi7yOUUt2ADK11WqBjCRADgG8DHYQfiAX+tnt9kFIgeFaUUnWApsBvgY2kRHgNyyDNHOhAfIl0hioESqmVwCVONj0DPA10LNmI/I+7a9Zaf3Fhn2ewPOIvKsnYSgjl5L1S8aSmlKoIfAoM11r/G+h4/IlSqitwRGu9SSnVNtDx+BIR+UKgtb7F2ftKqTigLpCmlAJL2mKzUqqF1vp/JRiiz3F1zVaUUg8AXYEOOjwXXRwELrd7fRlwKECxlBhKKQMWgV+ktf4s0PGUAK2Abkqp24HyQGWl1EKtdd8Ax1VsZDGUH1BK7QMStNah6GbnNUqpzsCrQBut9dFAx+MPlFJlsEwqdwAygN+B3lrr7QENzI8oy0hlAZCptR4e6HhKmgsj+ae01l0DHYsvkJy8UBzeBCoB3yulUpVScwIdkK+5MLE8DFiOZQJySTgL/AVaAfcD7S/8XlMvjHCFEERG8oIgCGGMjOQFQRDCGBF5QRCEMEZEXhAEIYwRkRcEQQhjROQFQRDCGBF5QRCEMEZEXhAEIYz5fxsgl4ZwrmBRAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.scatter(x, y_noisy, label='empirical data points')\n",
    "plt.plot(x, y, color='black', label='true relationship')\n",
    "plt.plot(inputs, outputs, color='red', label='predicted relationship (cpu)')\n",
    "plt.legend()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<a id=\"cuml\"></a>\n",
    "## cuML\n",
    "\n",
    "The mathematical operations underlying many machine learning algorithms are often matrix multiplications. These types of operations are highly parallelizable and can be greatly accelerated using a GPU. cuML makes it easy to build machine learning models in an accelerated fashion while still using an interface nearly identical to Scikit-Learn. The below shows how to accomplish the same Linear Regression model but on a GPU.\n",
    "\n",
    "First, let's convert our data from a NumPy representation to a cuDF representation."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "                    x                   y\n",
      "0 -0.2602409219007265 -1.6491150814315487\n",
      "1  -1.592967782370221  -2.911614740812639\n",
      "2  -1.849945489878636 -1.1900847094729925\n",
      "3 -0.2674443819927667  1.2062749267758428\n",
      "4   0.863977525145632   5.101137517220167\n"
     ]
    }
   ],
   "source": [
    "# create a cuDF DataFrame\n",
    "df = cudf.DataFrame({'x': x, 'y': y_noisy})\n",
    "print(df.head())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next, we'll load the GPU accelerated `LinearRegression` class from cuML, instantiate it, and fit it to our data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "cuML Version: 0.8.0a+1003.g7f9e070c\n"
     ]
    }
   ],
   "source": [
    "import cuml; print('cuML Version:', cuml.__version__)\n",
    "from cuml.linear_model import LinearRegression as LinearRegression_GPU\n",
    "\n",
    "\n",
    "# instantiate and fit model\n",
    "linear_regression_gpu = LinearRegression_GPU()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "CPU times: user 2.77 ms, sys: 8.56 ms, total: 11.3 ms\n",
      "Wall time: 10.3 ms\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "<cuml.linear_model.linear_regression.LinearRegression at 0x7f8dc8060fd0>"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "%%time\n",
    "\n",
    "linear_regression_gpu.fit(df[['x']], df['y'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can use this model to predict values for new data points, a step often called \"inference\" or \"scoring\". All model fitting and predicting steps are GPU accelerated."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "# create new data and perform inference\n",
    "new_data_df = cudf.DataFrame({'inputs': inputs})\n",
    "outputs_gpu = linear_regression_gpu.predict(new_data_df[['inputs']])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Lastly, we can overlay our predicted relationship using our GPU accelerated Linear Regression model (green line) over our empirical data points (light blue circles), the true relationship (blue line), and the predicted relationship from a model built on the CPU (red line). We see that our GPU accelerated model's estimate of the true relationship (green line) is identical to the CPU based model's estimate of the true relationship (red line)!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.legend.Legend at 0x7f8db81ec978>"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAD4CAYAAAAJmJb0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO2dd3xUVfr/32eSgSS0hA6hBBBCSYVQQ0eKghBBl8WKfhXLuru6iqJrwbaisOLaxfLDVXQRwYggnSC9JCShaACBAEkQAkkoSSBD5vz+mMwwSaYlmcmUnPfrxYvM3Dv3Pjflc899nud8jpBSolAoFArfROPuABQKhULhOpTIKxQKhQ+jRF6hUCh8GCXyCoVC4cMokVcoFAofxt/dAZjTvHlzGRYW5u4wFAqFwqtISUk5J6VsYWmbR4l8WFgYycnJ7g5DoVAovAohxAlr21S6RqFQKHwYJfIKhULhwyiRVygUCh/Go3LyltDpdGRlZXHlyhV3h6JQVImAgADatWuHVqt1dyiKOozHi3xWVhaNGjUiLCwMIYS7w1EoHEJKyfnz58nKyqJTp07uDkdRh/F4kb9y5YoSeIXXIYSgWbNm5ObmujsUryUxNZu5aw6RU1BM2+BAZo4NJyE21N1heR0eL/KAEniFV6J+b6tPYmo2zy7bT7GuFIDsgmKeXbYfQAl9FVGFV4VC4XHMXXPIJPBGinWlzF1zyE0ReS9K5O1QUFDAhx9+6O4wyMzMJCIiwu4+33zzjel1cnIyf/vb35waR1hYGOfOnav0/vLly5kzZ45Tz6Wou+QUFFfpfYV1lMjbwZbIl5aWWny/uly7dq1Gn68o8nFxcbz77rs1DcshJk6cyKxZs2rlXArfJjE1G42VVFfb4MBajsb7USJvh1mzZnH06FFiYmKYOXMmmzZtYsSIEdxxxx1ERkZWGmHPmzeP2bNnA3D06FHGjRtHnz59GDJkCBkZGZWOP3v2bGbMmMGYMWO45557KC0tZebMmfTt25eoqCg++eSTSp/JzMxkyJAh9O7dm969e7N9+3ZTrFu2bCEmJob58+ezadMmJkyYAEBeXh4JCQlERUUxYMAA9u3bZzr//fffz/Dhw+ncubPpplBYWMj48eOJjo4mIiKCxYsXm87/3nvv0bt3byIjI03XtHDhQh577DEApk+fzsMPP8yQIUPo1q0bK1asqOmPQeFDJKZmEz9nI51mrSR+zkYSU7PLbXt22X5KLaxYF6j1Y+bY8NoM1SfwisKrkccff5y0tDSnHjMmJoZ33nnH6vY5c+Zw4MAB03k3bdrE7t27OXDgAJ06dSIzM9PqZ2fMmMHHH39M165d2bVrF48++igbN26stF9KSgpbt24lMDCQBQsW0KRJE/bs2cPVq1eJj49nzJgx5Yp4LVu2ZN26dQQEBHDkyBGmTZtGcnIyc+bMYd68eSZR3bRpk+kzL730ErGxsSQmJrJx40buuece0zVlZGSQlJTEpUuXCA8P55FHHmH16tW0bduWlStXAnDhwgXTsZo3b87evXv58MMPmTdvHp999lmla8rMzOSXX37h6NGjjBgxgt9//52AgACr3ytF3cBeQdVSLh7ATwjemBypiq7VwKtE3lPo16+f3d7ny5cvs337dm6//XbTe1evXrW478SJEwkMNDyGrl27ln379vH9998DBnE9cuQI3bp1M+2v0+l47LHHSEtLw8/Pj8OHD9uNeevWrSxduhSAkSNHcv78eZNwjx8/nvr161O/fn1atmzJmTNniIyM5KmnnuKZZ55hwoQJDBkyxHSsyZMnA9CnTx+WLVtm8Xx/+tOf0Gg0dO3alc6dO5ORkUFMTIzdOBXeRVXbHK0VVJ/8Lh2wnnPXS6kEvpp4lcjbGnHXJg0aNDB97e/vj16vN702zszV6/UEBwc79ORhfjwpJe+99x5jx44tt4/5E8P8+fNp1aoV6enp6PV6h0bIlhZsNz4d1K9f3/Sen58f165do1u3bqSkpPDzzz/z7LPPMmbMGF588cVy+xv3tUTF9kHVTuh7VKfN0ZqIl0rJs8v20yRQS0GxrtJ2lYuvPionb4dGjRpx6dIlq9tbtWrF2bNnOX/+PFevXjWlSho3bkynTp1YsmQJYBDZ9PR0u+cbO3YsH330ETqd4Rf98OHDFBYWltvnwoULtGnTBo1Gw1dffWUqANuKdejQoSxatAgwpHGaN29O48aNrcaRk5NDUFAQd911F0899RR79+61G7s5S5YsQa/Xc/ToUY4dO0Z4uMql+hrVaXO0JdbFulKEMOTezVG5+JrhVSN5d9CsWTPi4+OJiIjgpptuYvz48eW2a7VaXnzxRfr370+nTp3o3r27aduiRYt45JFHeO2119DpdPz5z38mOjra5vkeeOABMjMz6d27N1JKWrRoQWJiYrl9Hn30UaZMmcKSJUsYMWKE6UkgKioKf39/oqOjmT59OrGxsabPzJ49m/vuu4+oqCiCgoL48ssvbcaxf/9+Zs6ciUajQavV8tFHHzn0/TISHh7OsGHDOHPmDB9//LHKx/sg2TbaHK2lcWaODS83+q9IQZGO+VNj1ExXJyIsPca7i7i4OFlx0ZDffvuNHj16uCkiRXWYPn06EyZM4LbbbnN3KG7HV39/E1OzeWJxGpbUIzhQy9Vr+nJCHqj1MxVOE1OzefK7dIsdNKHBgWybNdKFkfsmQogUKWWcpW0qXaNQKKrM3DWHLAq8AITAZhonITaUf/8pWqVlagmVrlE4nYULF7o7BIWLsVZAlRhSLvY+Y0y/qLSM61Eir1Aoqkzb4ECLOfnQssKqpW0Vi64JsaGVRF05TzofJfIKhQ9iLpZNArUIYRhhO0s4LRVQBTCiewviOjattM2RVIyznSfVDcOAEnmFwseoKJbmfefOsuxNiA0l+UQei3aeNOXmJbA0JZu4jk15Y3KkQwJrLsQaISoVY425fEdiNT9WcJCWy1euodNLp163N6JEXqHwMaxZAxgxn2FaE8FLysitVHw1ivK2WSOtivrs5QctTniy1G0DBoHu9eJqikpKkRgsDqb1b89rCZHljmt+Y8u3UBeoyg3Dl1DdNW6gYcOGgGHCkb02w3feeYeioqIqHd/cmKy6GGO0RkV3TkeupaoMHz6cii21UD0L5eLiYoYNG+ZU59CSkhKGDh1aY/dQZ+OIHa9xhqm5OZizzlPxfaMhWdislTy+OM2iwNujsEzgwRD71ztP8nziftN2ezc2ezG7GyklutKqf18cQYm8k6iOeLRt29bkUWON6oi8I9RUmCqKvCPX4iyqY6H8xRdfMHnyZPz8/Ozv7CD16tVj1KhR5Rw6XYEt10ZLOGoBYBzZVvX4xv2tzbBpEqgtt++zy/ZbnThVE77ddcr0taPi7Yn2CHvTVzP6uXY898aNLjm+Enk7ZGZm0r17d+69916ioqK47bbbTKIbFhbGK6+8wuDBg1myZIlVa+Hjx48zcOBA+vbtywsvvFDu2Eab4tLSUp566ikiIyOJiorivffe49133yUnJ4cRI0YwYsQIwGBgNnDgQHr37s3tt9/O5cuXAVi9ejXdu3dn8ODBVk3DFi5cyO23384tt9zCmDFjAJg7d67J1vill16q9JnLly8zatQok7Xwjz/+CFS2YDa/litXrnDfffcRGRlJbGwsSUlJpvNPnjyZcePG0bVrV55++mnTtU+fPp2IiAgiIyOZP3++6fxLliyhX79+dOvWjS1btgDln1Rmz57N3XffzciRI+natSuffvqpxWtftGgRkyZNMr1+6623iIyMJDo62uSDP3z4cB5//HEGDRpEREQEu3fvNp1j3rx5ps9GRESYvIQSEhJMdhGuwFwkJYbUxcwl6cS+staqKM8cG16pB90axly1+fFtjfAdEW2jTZFx0pMjI2xz/Bz0OTJP7zgi3p7Wh3/8ZDp3vtiLPok3kabPobPO9tNzdfGunPzjj4OTrYaJiQE7xmeHDh3i888/Jz4+nvvvv58PP/yQp556CoCAgAC2bt0KwKhRoyxaC//973/nkUce4Z577uGDDz6weI4FCxZw/PhxUlNT8ff3Jy8vj6ZNm/L222+TlJRE8+bNOXfuHK+99hrr16+nQYMGvPnmm7z99ts8/fTTPPjgg2zcuJEbbriBqVOnWr2WHTt2sG/fPpo2bcratWs5cuQIu3fvRkrJxIkT2bx5M0OHDjXtHxAQwA8//EDjxo05d+4cAwYMYOLEiZUsmM0N1IzXuH//fjIyMhgzZozJKTMtLY3U1FTq169PeHg4f/3rXzl79izZ2dkcOHAAMDwlGLl27Rq7d+/m559/5uWXX2b9+vWVrmnfvn3s3LmTwsJCYmNjGT9+PG3btjVtLykp4dixY4SFhQGwatUqEhMT2bVrF0FBQeTl5Zn2LSwsZPv27WzevJn777/fFJM1IiIi2LNnj819aoKlNIROL005Z0sFxYo96E0CtVy8okNvZehtyxXSPH9ta6aqOQVFOp5P3F+uKOsoxpmx1vL25pjfDCx1+2j9BA3q+XOh2HldRc7gfMFp/vXxHbxftAk/PTxX0Iun/7KIJuG2LU+qi3eJvJto37498fHxANx11128++67JpE3Cqota+Ft27aZbH7vvvtunnnmmUrnWL9+PQ8//DD+/oYfSdOmTSvts3PnTn799VdTLCUlJQwcOJCMjAw6depE165dTTEuWLDA4rWMHj3adOy1a9eydu1ak8fN5cuXOXLkSDmRl1Ly3HPPsXnzZjQaDdnZ2Zw5c8bm92vr1q389a9/BaB79+507NjRJPKjRo2iSZMmAPTs2ZMTJ07Qq1cvjh07xl//+lfGjx9vesqA8rbG1rz7J02aRGBgIIGBgYwYMYLdu3eTkJBg2n7u3DmCg4NNr9evX899991HUFAQUP57PW3aNMBg6Hbx4sVyNxxL+Pn5Ua9ePS5dukSjRo1s7lsdHElDWCooVuxBr1iYtEeplDy+OI3kE3m8lhBZJdGWwNc7Tzp0nopM6WPwlHckbz+tf3vT194wuar4aiHvffYg/8pZzCV/PffltuXluz4nNH6cS8/rXSLvJqthW7a5RnMwe9bC9qx2pZQO7TN69Gi+/fbbcu+npaU5bOVb0db42Wef5aGHHrK6/6JFi8jNzSUlJQWtVktYWJjJTtlWnNawZGscEhJCeno6a9as4YMPPuC7777jiy++KLd/TWyNAwMDy8Vs63tt6VjW7KSNXL161WUGbNYmHVXE3s3AUsujI3y982S1BbuqBGk1LE3JduhG1KCeX7nuGrA8ucoVVLX/vlRfyteLn+P59PlkBeoYn9eIOWPmEpEw43puy4WonLwDnDx5kh07dgDw7bffMnjw4Er72LIWjo+P53//+x+A1fztmDFj+Pjjj01CZkwhmNsHDxgwgG3btvH7778DUFRUxOHDh+nevTvHjx/n6NGjphgdYezYsXzxxRemvH52djZnz54tt8+FCxdo2bIlWq2WpKQkTpw4USmuipjbGh8+fJiTJ0/atBo+d+4cer2eKVOm8Oqrr1bZ1vjHH3/kypUrnD9/nk2bNtG3b99y20NCQigtLTWJ85gxY/jiiy9MtRXzdI2xiLp161aaNGlCkyZNCAsLM8W0d+9ejh8/btr//PnztGjRAq32erHRmTiaX3ckJ22p5dGTKNbpHX7SKCpx7vrKlrBUkLZUI7FWw5BSsnrdR/SeFcL0w2/R+jIktXmWFR8UEHHrQ7Ui8OAkkRdCfCGEOCuEOGD2XlMhxDohxJGy/0OccS530KNHD7788kuioqLIy8vjkUcesbjfokWL+Pzzz4mOjqZXr16mIuV//vMfPvjgA/r27VtuGT1zHnjgATp06EBUVBTR0dGmBblnzJjBTTfdxIgRI2jRogULFy5k2rRpprVaMzIyCAgIYMGCBYwfP57BgwfTsWNHh65rzJgx3HHHHQwcOJDIyEhuu+22SsJ95513kpycTFxcHIsWLTJZKZtbMM+cObPcZx599FFKS0uJjIxk6tSpLFy4sNwIviLZ2dkMHz6cmJgYpk+fzhtvvOFQ/Eb69evH+PHjGTBgAC+88EK5fLz5tRprJ+PGjWPixInExcURExNTrqgaEhLCoEGDePjhh/n8888BmDJlCnl5ecTExPDRRx+VW6UrKSmJm2++uUrxVoWE2FDemBxJaHAgAsNotyKOFhRd0eHiTKpyA9II4XA3UHWwJuazlx90yEN/b+rPjH42lJu2P8plXSH/q38nu+YWMHzGv0BTu2Nrp1gNCyGGApeB/0opI8reewvIk1LOEULMAkKklJWT0WZ4otVwZmYmEyZMsFuAU7iH2bNn07BhQ1ONxBqpqam8/fbbfPXVV1b3GT58OPPmzSMuzqJjq0UmT57MG2+8YfVJxZm/v5by6gK4c0CHSqkLS599fLGTmxY8BK2fYO5t0U5N1cTP2Vilm6IAjs8Zz/ET6Tz/2R184/8rzYrgRb+RPPzEIuq1aO202Cye34bVsFNy8lLKzUKIsApvTwKGl339JbAJsCnyCoWriI2NZcSIEZSWljqtV76kpISEhASnr3plLedrqdNGcr1fPCkj12qe2NZqTd6OrtRQJJ675pDTiq1VnTTVtkEhT/5rOO9f+QWNhOfyI3j6L9/QpKvtm29t4LRFQ8pEfoXZSL5AShlstj1fSlkpZSOEmAHMAOjQoUMfY87XiLtH8gpFTajq76+l0bqxrdDaIh2WMH4GDALv6akaZ2G+OElFqlIwtTaSDwnScsWsdiBlMWH575HefAuXtJLpuaG8fPfntBs4ttJnXYlHLxoipVwgpYyTUsa1aNHC3eEoFG7F2rqpT36XXm4mqT2KdaW8/NNBl802dRYCg3A6C2trzD6fuJ8nFqc5POnLUsE7UOvHS7f04o3JkbRtUo+2+V+g009jc+hmBuc3Ir33Aj7/MKvWBd4erhT5M0KINgBl/5+1s79CUeexliYolZLCkmtoNY53ZOQX6ao827S2kYCzVyC15J3ztYXWUVuLjlcseIcGB/LG5EgmxbQl7adnyc0ex462ywgt1pDU9jlWfJBPxKQHnXshTsKVffLLgXuBOWX//+jCcykUNSa/qIQzF65QUqqnnp+GVk0CCAmqV6sx2OqL15VKQoK0XCy+ZnfWqTdRHcMyW1RsJ529/KDVfW3l3iv23X/11VtEvPMMv3aGjv6wQDeJ/5v3LZoAz/PDMcdZLZTfAjuAcCFElhDi/zCI+2ghxBFgdNlrhcIjyS8qITu/mJJSw6SnklI92fnF5BeV1FoMianZ5BVetblPfpGOaf3bUzsd1t6HeTupsc/d1k3EkfkFW7etpPcUwT3HnuFMa5iwCtbdu4cHX0v0eIEHJ4m8lHKalLKNlFIrpWwnpfxcSnleSjlKStm17P88+0eqGyirYceoTavhMxeuoK8wOtZLyZkLtmf35ubmMm7cOPKLSsg4fZF9WQVknL5Y5ZtDYmo2M5ekU6zT2933m2p4wtQFBIYUzOOL0+j5wir+8V2a3XqErfkFxzJ/ZdAt9Ri1egK/9YDxm+H7vl/x005J11jH22zdjdsLr76Cshr2bqth4wje0feNtGjRgqYtWrFibVK5p4BTeUVk51v+uVmaSTl3zSHTKkb2sH8bqJuYf/eKdHqrhmz2uFR4gZuntCbmk17s6q1j2D74OOgZVmyQDL/tLqfEWpsokbeDshr2PathvV7Po48+Sq9evZgwYQI333wzSauWA3DTwCjm/+sl7pgwijsmjOL0SYOFwfTp08vdxMyfdOJHjmPFsu8qnfN8YQmnC4odmhbvyR0wvox54bW09Br3PdCXbi8EsyrqDNEn4aVT41n7o+Te57w32+xVBmWPr36ctD+cO2svpnUM74xTVsN1yWp42bJlZGZmsn//fs6ePUuPHj340533oCnzEmnYsBHfrNjAiqX/4+1X/8nYVT8DcPnqNTJOX6SkVI+Uhjx+SFA9ukVE886br1r8fpdKyol5fX+NxRZJPwvrmypcT05BMVJKZs++i6/PfsOx9hCZDQkHu/H+igP4uciTqDZRI3kHqGg1bBR1sGw1HBMTw0MPPcTp06cBg9Ww0cL27rvvtniOqloNx8TE8OWXX3LixIlyVsNCCO66y/ojpTWr4d69e5ORkcGRI0fK7W+0Go6KiuLGG2902GrYeJ3WrIYDAgJMVsOdO3c2WQ2vXr2axo0bm45VFavh5s2bm6yGzaloNbx161Zuv/12NBoNrVu3Nj0lGf2ibpp0G/4awQPT7yZ59y4ASq7pybtcYkrJSDAVZlu3akXumT9sfk/AIObWioBK4N1D6z+WETldwyuabyitD+N/DKSo/lfsG/ERPx3wja5vrxrJ2xtxuwplNex7VsPm6EoNAl5qTOIKwfUvDccq0RssY42f1+lKTIXZJvWhvhd0WSiuU+/sDhpdeJ2dUdCsCMavgpPR73FgTCfA8mIs3ooayTuAshr2LavhwYMHs3TpUvR6PWfOnGHL5l+QZmW7NT8tQy8lC7/6hoEDBwLQOrQDv+433MCT1vzMNZ1hRF5Sqic3K5MePXtWKWaFe9BcPkHL9ImcCH2dX3vAzVtg2cBvyU/YwMWQTuX2tTVZyptQIu8AymrYt6yGp0yZQrt27YiIiOChhx4iIqYPDRtdTxHpSkq485Yb+e+nH5mKwH++ezopO7dzx4RR7E9LJjDI8ERUz09DUlISN918M/5VmI2qqF2k7iItU+7gfMO/kNJPz9B9EJ08lZXrJUMT/mx1UlRVjco8EacZlDkDZTWsqCrVtRq+fPkyDRs25Pz588T2iWPhstU0b9mKmwZG8c3KJEKaNkMgaNc0kJCgeqbJUua99BohCA0JZPyYUcz95GsalS1raM6Zk8d4cPlp5160wmH0eh0t0/9Bdq/j/NEE4g+D5sQATg5+HoDMOeMB64ZkocGBbJs1slZjrg4utxpWKDydilbDEyZMoKCggJKSEp559p+0bNW60mQoiSQ73/CHHxJUj6Kr18gr1JlSO3op2XfkFH+a/rBFgVe4D4kk5OBbXGq3hZRBho6ZiM3tOTzyfUR7w1yJULPZrpYWAnd0MRZPR43kFQoM7ZBZecXlcvNGjD42FUfyjqBG8rVPg2P/QxPwNQe6QId86LklkEODP0MfeP1GbMmSuKprt3oSaiSvUNghJKgep/Isz1AtKdVbtD1QeBb1zmyj8cU3SImCpkVw8yr4PeI9fhvdCYHB9kACfkIwpU/lRb9rayHw2kaJvEJRRj0/jUUbA2vvKzwDzcVjtDj2d9L6SjQSxm2BrObPcnB4vGkf89tzqZQsTckmrmNTnxT1iqjuGoWijFZNAkyzXo1ohKBVkwD8VOeMxyF1F2mdPI28xn8jub9kyD7olTKN3+JWcCks3uZnfaU90hHUSF6hKMPoHW/0lBcI9FKSU1CMXg3kPQa9XkfrvU9wKiqTXUNg0GEQJwfxe/xzVTqOL7RHOoIaybuRQYMGVWubLWbPns28efPs7ldVK2FX8eKLL1r0ozFn06ZNbN++3enntmQPHBJUzzSiNxZhS/XSYkFWUbtIJE33vU69S7eyZ0gmzS/DyBUdyQr9kawqCjw45iXvCyiRdyOWhMtoWewKUasKtSXyr7zyCjfeeKPNfVwh8pYWCTmVV8TBnAuGkbsqsnoUDX//miY5t5Dafwe6+jBmeQMuB33D0VEfIDSW7aONaDVYXK/VF9ojHcHnRN6SV3dN+frrr+nXr5/JeMwoxA0bNuSZZ56hT58+3HjjjezevZvhw4fTuXNnli83WNcuXLiQSZMmMW7cOMLDw3n55ZdNxzWOpjdt2sSIESO44447iIyMLLcN4K233iIyMpLo6GhmzZoFwKeffkrfvn2Jjo5mypQpdj3nrdkdO2olbG2/ijRs2JAnn3yS3r17M2rUKHJzcwGDv86AAQOIiori1ltvJT8/Hyhv4RsWFsZLL71kOkdGRgaZmZl8/PHHzJ8/n5iYGLZs2cKSJUuIiIggOjq6nGMmlB+d78+6YHERj/yiEk7lFVkU8lK9vO5ho3A79XM20+LIBA5G/o/sNjBuFbTQf8Kh0YuRAY3tHwDQ6WFKn9BK67XWhaIr+JjIW/PqronQ//bbbyxevJht27aRlpaGn5+fyZelsLCQ4cOHk5KSQqNGjXj++edZt24dP/zwAy+++KLpGLt372bRokWkpaWxZMkSi6sd7d69m9dff51ff/213PurVq0iMTGRXbt2kZ6ebvJgnzx5Mnv27CE9PZ0ePXrw+eef27wOo93xnj17aN26tel9o5Xw3r17SUpK4sknn0RKyZw5c+jSpQtpaWnMnTvX6n4VKSwspHfv3uzdu5dhw4aZbmr33HMPb775Jvv27SMyMrLczc6c5s2bs3fvXh555BHmzZtHWFgYDz/8ME888QRpaWkMGTKEV155hTVr1pCenm66mULl0bkxxWK+lF92fpHVVkmF5+B38Sit0yaQ2fEtDvSAsVugzcnn+W34CsaNHFjl4yVl5LJt1kiOzxnPtlkj64zAg4+J/Nw1hyx6ddekir5hwwZSUlLo27cvMTExbNiwgWPHjgFQr149xo0bB0BkZCTDhg1Dq9USGRlZzhZ39OjRNGvWjMDAQCZPnlzOqthIv3796NSpU6X3169fz3333UdQUBBw3YL4wIEDDBkyhMjISBYtWsTBg9YXKwbrdseOWgk7up9GozHZLxttmS9cuEBBQQHDhg0D4N5772Xz5s0W43TEWjg+Pp7p06fz6aeflluRy1Yvu7GAer6w9tZsVVSDkgu02TOV843/zu4BMHgf9Ei7m4y4FVzuOIAG9fxYmlL1QVtdKbJawqe6a1xhMiSl5N5777VomqXVak1WtBqNxmTCpdFoytni2rPChfIWwBXPb2n/6dOnk5iYSHR0NAsXLmTTpk12r8XScRy1Eq6O5bC1c9rCEWvhjz/+mF27drFy5UpiYmJIS0ujWbNmdnvZVRrGc9HrdbRJ+Rsno0+xcygMOAwiazBHB84qt1+RrhRL93GjRYG1FbbsFVm9ebarPXxqJG/tB1mTKvqoUaP4/vvvTRa8eXl5JrtdR1m3bh15eXkUFxeTmJhoWoDEEcaMGcMXX3xhyi95kZQAACAASURBVLkbLYgvXbpEmzZt0Ol0Vu2LzbFmd+yolbC1/Sqi1+tNOfZvvvmGwYMH06RJE0JCQkzL93311VemUb0jVIzl6NGj9O/fn1deeYXmzZtz6tQpAASql93bkEiapb9G/Uu3snvoKZpfhuE/dyYndDk5FQQesCjwYBjIzRwbjtbKfIYR3VtYjcEVaV5PwqdEfubYcKdX0Xv27Mlrr73GmDFjiIqKYvTo0aYVnxxl8ODB3H333cTExDBlyhTi4hxf6X3cuHFMnDiRuLg4YmJiTO2Rr776Kv3792f06NEm+19bWLM7dtRK2Np+FWnQoAEHDx6kT58+bNy40VSb+PLLL5k5cyZRUVGkpaWVq1nY45ZbbuGHH34wFV5nzpxJZGQkERERDB06lOjoaADV5uhlNDr8FcE5t7B3wE5K6sPoFY243HAxx0e8i9BUTZqaBGpJiA1l7u3RBGorf3ZpSrZV0XZFmteT8DmDMk977Fq4cCHJycm8//77bouhNmnYsKFpEZLaxrj+qiehDMoqE5CzicaF89hT5jETtxmORX2ILqRDtY8ZEqQl9UXD4vRVtQ3uNGulxeGBAI6XWRF7OnXKoMxXTYYU9qmuU6SidvDP/50WmY+ztx9oJIzeAtmtX+TQsH41PnZ+0fW1c6tam2sbHGjxpuArk6V8Kl3jiUyfPr3OjOIBt43i84tKTN01KjfvYVzJp+3uP3Eu5HF2DYT4fdB933QOx62gsF3NBd6IMR1T1dqcK9K8noRXjOQdWeRa4f0YhbqkVG/ycDf6yVjbTyAq5eI9JTcvZd22Q9DrdbRNfoyTMdnsGGbomOH0cI71s72KV3V5+aeDzF1ziOyCYpOtsBFbom188vekNK8z8XiRDwgI4Pz58zRr1kwJvQ9TcXk94wQmuG4cll9UQk5BcblWSE8VUSkl14oucqJAZ39nH6R56stc6LyHXcOgVzZ029GF40PmI9q7LnmQX6QzpW0k1/3jQx0QbV9O83q8yLdr146srCzT9HiFd1NUco2Lxdco1Uv8NILGgf4E1fPnjwtXuGahjz33lKB1kwCKSq5RUKTDW1rdJZITBTre25Xv7lBqlSaHFiIafU/KIMOqTKNWNuHo4E+5PCyo1pNoXvKr4nI8XuS1Wq3FmaAK78HY8WTtMfqNyZE8kXjM6h9lqJXCmMJzCMraQKPi+eyOMXTMjF4tOBbzEb+PbOfu0Ex974DPjtZt4fEtlArvxjjRpGIfsjn2ZisqPBf/vEO0OvEkyWUdM4N3Qlab2RSHOj4XpLaw1kLpC9SpFkqFZ2FpoklFcgqKmT81xu7NQOFBXMmjXfpDHOhfzLG2MCwV8sX9HImb7LaQKj4lVsTYQulpc2lcjRJ5hUtxxDeobXCgxQ4HNbL3PPT6EkJ3/4WTsafZNhz6HwbOjOJ43BPuDs1UZLXlX1PxybIupHJUn7zCpdibUCIw/KHFz9kIYLKD9ZUeZV+iZcqLBFyazM4RpwkugiFrwjkdupw/PEDg4Xo65p2pMVb73n3dwsASLh/JCyEygUtAKXDNWt5I4ZvMHBvOE4vTrD5GG983H1EB5b5WuJcmv36GCE5kz2BDx8zIVSEcjf+EwqG13zFjDfM+eFt9708sTrP4eV+2Iq6tdM0IKeW5WjqXwoNIiA0l+UQei3aeLCf0lvKnxbpSXv7pIEH1/FVu3gMIOrmexlfeYVcfQ8fMqDUaMmM+5ujwtu4OrRIVV3qy1vdeVQsDX8jfq3SNwuW8lhDJ/Kkx5ZZfszayzy/SqVy8m6l3/lfapU7gaJd3SO8Jo7ZAyz9e5fehy7nW2PMEPtSspmOPqlgY+IoFcW2M5CWwVgghgU+klAvMNwohZgAzADp0qL4LncKzMY6sjCMjhechis/RLv0h9g+4yu+hMDQNCsSD/B43yd2hWaWqHjNVsTCwlb/3ptF8bYh8vJQyRwjRElgnhMiQUprWfisT/QVg6JOvhXgUbsKRnnlF7aPX6+iw6xGO9/6DrSOg32GQZ8eQ2edv7g7NJsGBWmZP7FVlwXXUwsAVK825A5eLvJQyp+z/s0KIH4B+gOUFPhU+jSM984rapVXy8+R3TWPbSOiZA5339ODkoDdd6jHjLK5eM6wd4Kq8ua9YELtU5IUQDQCNlPJS2ddjgFdceU6F5+JtIyBfpumBT5BNf2L3EGifD8NXNeV4/CcUDg70mI4ZI4FaDcW6yovBGAv1V3R6l/S9zxwbXunJ0xstiF09km8F/FDmHukPfCOlXO3icyo8FDXByf00PL6aRrr32dnX0DEzco0fJ3ov4PjwVu4OzSolpdazuOaLhRhxVt7cVyyIXSryUspjQLQrz6HwPKw9PtvqmfcTglIP8lHyNeqfO0jLU8+wu8xjZuQWyG73L44OjXJ3aDYRUM5a2lGc9dToCxbEytZA4VTsTRt/3MpklFIpldukC9AU5tJu30PsG1jCkXYwOA0u+D3C0TjvWLvUnrwHB2opKK48mve2vLkr8fzqisKrsDdtPNTKH58ARnRvgcbTEsJeiiy9Svtt93HV/z62jCwh/BT02XETJ7uv4EJX7xB4ewhg9sRePr10nzNQI3mFU7HVdpaYmk1e4VWL2yXw9c6TLoys7tB697Pkh+9n643QIwfCUnqRPehNaO/uyKqOLWdJie/kzV2JEnmFU7FWXA0O0palcSp3SSicQ7P9HyGbrWTXMEPHzNA1zcmM/4SiQfXdHVq1kVhPyRifCn0hb+5KVLpG4VSsTRuXEtUj7yIaHltJ28MT2NtvJZltYMQ6LfVKv+DE0IUIP+8VeDAIuUrJ1Aw1klc4FWuPz9YKrorqUz93H62ynmNXWcfM8C2Q0/ENjg2OdHdoTsM89aJSMtVDibzCIRydVWhpP4Vz8Ss8Q/v0h0gbdI0j7SE+DQrq/YXjcTe5OzSnEhyoNf2OqZRM9VEir7BJYmo2s5cfLJcTrdgWaW2h7uyCYp5YnOZxMyi9FVl6lY47ZnCs73l+GQV9D4M+bwKnoh92d2hOx9g5o6g5SuQVVrFlKGbeFmm+T8VOCGnhPUXVabvzac73+JUtow0dMx1SI8kZ8IZXdsw4gnnnjKJmKJFXWMWeoZhhRL9Pdcy4kOZp7yFbrmHHCEPHzJC1LTkx6COKBnh3QdUeQj3+OQ0l8gqrODI1XAm8a2j8+0801n/C9gEQXAzD1tfjVMwCTg5pXifSX8rhwnkokVdYRRmK1T4BZ9JonfM8O8o6ZoZuhdMd3yIzvqe7Q1N4KUrk6yDWOmUqvj+iewuWpmSr/vZawO/yH3RIf4jU+FIOd4RBaVBQ/69kxo11d2huISRI6+4QfAYl8nUMawZiySfyygl6dkExS1OymdInlBXppy3OOFQ4AV0xHXfO4Gi/fDbdCHGHobRgElmRD7o7Mreh9RO8dIvqrHEWSuTrGNYMxL7ddaqS1W+xrpQV6adNK/AonIfU62m382nO98xg8xjongP902P4o99rPtsx4wihaqKT01EiX8ewVky15uWuRvDOp8Xe/yBbr2P7KEPHTPz61pwc8CHF/eq5OzS3s23WSHeH4HMoka9jWCumqkU7XE/jIz/SWH7K9kGGjpmhG+pzKnYBWfHNlIkU1m2oFTVD/W7VMSwZiAlgQOeQSu8rnEPg6b102juBQz0/JaUnDN4KzfL+zYlBS9EHNnN3eB6B1k8oCwwXoUS+jpEQG8qUPqHleq0lsPfkBab0UXlQZ+J3MYcuWydyus2L/DII4vZD52N/52TcCnTNlKCZ06Cev8rDuwiVrqmDJGXkVrIaKNaVkpSRS4N6fhSWqJbJGqErotOOBznS/wIbR0OfI3Dtwq1k9/o/d0fmsVxQtR+XoUS+DmKt+KomPtUMqdfTYfuTnIs4wqaxho6Zvvv7cDbuZWjn7ug8G7Umq+tQIl8HUTNZnU+rlPmUttnA1tHQLh8GbWjLqQEfUBynJvXYQy0A4lpUTr4OMnNseJ3wP6kNmhxaRvvDE9gTv4FjbWDIxkD8+YrsQQvQaJTAO8IbkyNVPt6FKJGvgyTEhir73xoSmL2HLnsnkBHxBck9IX4rNM1/m5MDlyADQtwdnkcRqPWzalMQGhyoBN7FqHRNHSVUpWyqhf+FLML2P0LyYMmhLjAgDQoaPMmpuBHuDs0jCTVbHazi2gQqTVM7qJF8HUX9cVUNUVJE51/u4ErAw2wYLemSBZGpt5HTfQVF7ZXA2+KJxWnMXXOI3h2a4FdmFO8nBFP6qCX9agMl8nWUhNhQ5fTnAFKvp+PmvxNQ9CeSxl2kYRHEbevH+XYruNBjurvD83iyC4qRZf9vO5pnmlVdKiVLU7JJTM12b4B1ACXyXkJiajbxczbSadZK4udsdMofx/ioNk6IzHdpvWceLf+YyOaxR7kcAAOT2lHUJJHc3i+6OzSfwHwJSYXrUDl5L8CaPTBYXgfTml98xX2WpqhRlCVCfv2Oxv7/ZcsQg8dMfFIDsmM+IWdAsBoVORlHVh9T1Awl8l6ANXvguWsOWRRvezeExNRsnvwuXRmSVaBB1i7anH2Vrf1BSBi4Dc50/g9ZA7q4OzSfRU2Ccj1K5L0Aa6MdS+/buyEYbwJK4K/jl3uEzhlPkDwEfrsB+qdBfqOZZPcZ5u7QvAJrDqYCbLbqqu6a2kE9fXoB1kY7lt63d0OwdBOos1y5SOdVE7ja+AnWj4FOWRCRPpXT3VdwJVQJvCMEav2Y1r99JQfTQK0fdw7oQGhwIAJDK+VdFV6rSVC1g8tH8kKIccB/AD/gMynlHFef09eYOTbc4R5ja5YFGiHoNGulmgQFSH0p7RJvI3+kjqTJ0C0Hond2Jm/Yu+4Ozeuo768hrmNT4jo2tVsHUrgHIV342C6E8AMOA6OBLGAPME1K+aul/ePi4mRycrLL4vFmHCmmGvereENQXKfF6r8joo6ypwuE5kPrX7ScHb0EjUZlLqtLoNZPjcrdjBAiRUoZZ2mbq3+z+wG/SymPlQXyP2ASYFHkFdZJiC0/ccTYUmlJ9AO0GpPI28uL1hUa73ifps1WsznB0DEzYDWc7v8158aqjpmaYq0JQOEZuFrkQ4FTZq+zgP7mOwghZgAzADp06ODicHyDxNRsZi5JR6c3yHd2QTH/WJzGc8v2UaQrv+h2XRf4wEOraVf4PluGGjpmBmyDP8LmcXpYd3eH5lOoVkjPxdUib8nssJzuSCkXAAvAkK5xcTw+wezlB00Cb0QPlQS+LuN/JoPOh58ieQgcCoB+aZBf7zFy+oxzd2g+iWqF9FxcLfJZQHuz1+2AHBef0+cpUKvoWEUUFdDll7vIGAnrxkDMEWiVO44/Yh9zd2g+QXyXpuw9eUEZjXkRrhb5PUBXIUQnIBv4M3CHi8+pqINIfSlhP99N7uCLbJgMXXMgcldX8ofOV6syOYmQIC3bj+YRHKSlvr+GC8U61UnjBbhU5KWU14QQjwFrMLRQfiGlPOjKc9YFQoK05Bep0byRdjteR9dhB5tvN3TM9F4eSO6obykZqjpmnInxdy6/SEeg1o/5U2OUuHsBLm8skFL+LKXsJqXsIqV83dXnqwu8dEsvtH5qbadm+xfR6fAEto/YweG2MOCXRvj7fct51RLpcpS5mPeg/hK8kITYUJJP5LFo58k62T3TMHMrbfPmsLnMY6b/dsGZru9zul9Hd4dWp1AdNd6BEnkvJSkjt84JfL3zx+hy8G/sGgq/doe+aVAQ8k9O9x7o7tB8ltDgQIpKrllMDxo7ahydqKdwD0rkvQjzP6a6JPCa4gt02/4ABwcXs3osRB+Btlfv5kz3qe4OzScJrSDUlmZRGztqqmqDrah9lMh7CXXRrkDqS+ma9Bdy+mSxZoKhYybm8FDyI592d2g+S2hwINtmjSz3nlGsLY3W4+dsdNgGW+EelMh7CXXNPbL9tlfRhe1iwwRomw99t3bhj9h/UxKifmVdha1+d3NbDeMT5ROL06w+Uap8veeg/mK8BEvOkr5I87T/0jjoOzaNgibF0H9zE/6I+YSzfRoqjxkXIsAhkzFHnyjVDFjPQYm8l2BtYQZfofGxTbQtmMemso6Zfts1nO32Pn/0VX5GtYHEsRy6I0+UagasZ6FE3kvwVYGvn3uELr89wc6hcCAA+qTBxZDn+aP3AHeHprCArTSMANVd44EokfcSfG0k71eYT7cdD3JgyBVWjYWoI9BGdy+53W93d2h1Eo2Dc+usLUpjqWCr8AxUmtODMHrEd5q1kvg5G0lMzTZt8xmBL9Vxw4YZBOruZvUtVwgogug9w7nQbgWFnZTAuwu9g79eM8eGW1zqT6VnPBc1kvcQ7PUbh1oZQXkTHbe8REnnFFPHTJ9tXTkbOw9diJ/9DytcSqiDhVJb7ZQKz0SJvIdgqaBl3m88onsLvt550k3R1YyWqf+PRg2WkjTa0DHTb0sIf8R+wrneQepR0gOo6ki84iplCs9GibyHYK2glVNQTGJqNktTsi1u92Qa/76RdhffZmPZWmB9d/iT2+09zsS1t7iajKL2qTi7VeF7KJH3EKwVtNoGB3rdRKiAs4fpkvEPdph1zFxo+hJnYvu6OzRFGWrx7bqDelr2EGwVtLxl9qBfYT49106msMk/+HkstM+C7hn3kdt9BSUtlcC7k+BALaHBgQgMo3cl8HUHNZL3EGwVtOauOeTZRddSHeEbHia73xlWTYIbciDy2Cgu9HjC3ZEpMPSvT4huw2sJke4OReEGhPSg1ry4uDiZnJzs7jA8Dk82Jwv75QWu3pDKzhsMHTOtM8I5FzMXoR4SPQqVnvFthBApUso4S9vUSN6DsOXLXarXuzm68rRO+ZSGjX5k41hDx0zc1qacjfmY8zFBqqjqgShnyLqLGsl7CJZG61o/gb9GUKzzHIEPPrKO0Ev/YUNZx0zUHi253T5ABrV1b2AKuwjg+Jzx7g5D4QLUSN4LsNRBoyuV6Eo94yYc+Mdv3HB4JtuGwr4AiE2Hi81f5mxMH3eHpnAQ5QxZN1Ei7wEkpmZ7bGHV//I5wrc/yP5hOlaMhYgj0Kr0Ac6HJ7g7tDqP1k84PAhQ1gN1F1UdczPGNI2nIXRX6b7mfgKuTefnSTrqFUFE2hgutVtBcUcl8O7CTwhTG+Tc26IJCdJa3K9BPT/VMqkA1Ei+1jAvqjYJ1CIEFBTp0HiYu6TU6+nyy/Nc6baPNQnQJh9id/bkfPQcroWoMYG70UtZKa8+8/v0ciN6rZ/g9VuVqCsMKJGvBSoWVQuKdaZttgReqxHoHLUHdAJtdn9MoyYrWH8TNL4CfbY3JzfqQ/KiVceMp1Axr64MwxT2UCJfC1THlsDoKVIbE6GaZqymbdH7rI83vO6zqx654R9wLraNEncPwlpeXRmGKWyhRL4WqKotgfGP2fiH+4/v0hz2+64KQacPcMORWWwdCmkBEJMOF1u8Sm50rPNPpqgRfkKovLqiWiiRrwWsmY9ZouIf83PL9jld4P0v5dJz2wOkjSjlp7KOmZY8RF74Lc49kcJp6KWsJPCJqdm8/NNB8osM6b/gQC2zJ/ZSNwJFOZTIu5jE1GyKSq45vH+plMxdc8j0usiJE6GE7go91z/EyYHn+elW6JIDvTLHcSn8MaedQ+EaKubiE1OzKxVcC4p1zFySDji2KLeibqBE3oVY+kN0hOyCYh5fnIZwUkJc6vV0TXqW4u4H+XmyoWMmencE+ZH/Uh0zXoClXPzcNYcs/l7p9FLZFyjKoUTehbz808EazVh1Rmdl6M4PaRDyM+vHGzpmeu9oSW7UhxREBqiiqpdgKRdvq87jLdbUitpBibwLMeZK3UHTX1cSeuUj1g0xvO69uz7nur3P+Zg2agacFxESpLU4KrdV51H2BQpz1N+7j9EgK52YzRPI7PkRKwdBjwMQmvs656KWQkAbd4enqAJaP8FLt/SyuG3m2HC0fpWfxbQaoewLFOVwmcgLIWYLIbKFEGll/2521bk8leBAy1POXYH/hT+IXjWRyy3+yY9joW0W3HD0EfLDV1AaHF1rcSiqTnCglsw543lnakw5K4K5t0Vbza0nxIZWsjUIDtQy93brn1HUTVxmNSyEmA1cllLOc/QzvmY1nJiazcwl6S6dtaopuULP9Q9yYlA+B9pA59MQcO5mCm941GXnVFQfAZj/NqjFPBTOQFkNuwlLU86LSq45JVcv9Xq6bXia4p4ZrJxi6JiJSo6ioNdrFAarLJy7CNJqbLa9SgyjdGVBoKgtXC3yjwkh7gGSgSellPkuPp/HUXHKeadZK2t8zPbb3qVB87WsvcXQMRO7szXnIt/nQi/VMeNOgrQaQhrUp8hGd0tocCDbZo2sxagUdZ0aDfmEEOuFEAcs/JsEfAR0AWKA08C/rRxjhhAiWQiRnJubW5NwvIKadD40O/AjUXsnsGf4WjZFQMyeAJpc/Zy86M/QaAKcGKWiOuj00u7M5qKSaySmZtdSRApFLS3/J4QIA1ZIKSNs7edrOXlLVGdR7ganUgk//gKbhkJBAESlw8XWb6BvEunCSBXVwc8B62iVh1c4G1s5eVd215j3690KHHDVubyJhNhQ3pgcaeqisLboA4C2IIeYVbdwueUL/DAW2mTBDcf+QkH4CiXwHkqplARq/WzuY1xUW6GoDVyZk39LCBGDodaUCTzkwnN5FeZ5+sTUbJ5ckk6pWQeOpqSYiHUPcDz+Aj9Ohk6noUfWLRR2Ud9CT8fcIjqnoBhrY3o1K1VRW7hM5KWUd7vq2L7E3DWHTAIv9Xq6r3uSoogj/HQbtC6AyJQYLvR8hSLVMePxmFtEG2/i8XM2WszTq1mpitpCKYebMQpAh63z6XlkIusmHWF/KETvboNWu5SLPV9DqB+Tx2LsZrK2jurMseGV0jdqUW1FbaL65N1M833LaKf7gtUjDK+jkwPJ6/YhBZEtlLR7AfOnxpAQG2paw/eJxWnl+t/V8nwKd1Mr3TWOUhe6a4ykrl/Nay/eRFJZx0xkOlxs8yaysWWvEoXnIbhuFKZmsirciZrx6kKMIzhLo7SK20Z0b8Ev2/bSYs109o6EE2Ohx+/Q1P9vXAgf4+YrUVQVyfV0W8WhkrGDRom8wt0oka8BFXveswuKeXbZftN2822n/zjH/q8SKBhyia1TIOw0RBxKIGrE02w7mueW+BWuRXXQKDwBJfI1YO6aQ5UmNZn3QBfrSpH6Unqu+QeXo46y/E/QqgAi9vbhYo+XuBSsIfN8MXcN6MC3u05RKiV+QuDvJ7h6zXnL/incg+qgUXgCSuRrgLWRmvH9jpvnEdRqE2tuhUZXIHp3KHm93uVSj/qmrozsgmK+3nmS4EAtE6LbkJSR6/Ci3wrPRXXQKDwFJfI1wNrqPF0O/Ujjy5/y8yjD66jkBuR3+4CCyOZWO2YKinV8vfOk64JVmIqhQJWsJbR+gql925OUkVuuvrI0JbvcMYzF11DVQaPwIJTI14CZY8PLiUWj47vofupVNpR1zESkCy62fYuCiB7KHdLNhARpeemWXuWEd+6aQ449NUmI69iU1xLKW0nEdWyqWiMVHo8S+Rpg/IN+++u1tFx3P8kj4fue0P13WHjn/0M/abRJBDQOGFcpXMeVCh7vxh52RwzjdHppsVOmoo20QuGJKJGvAYX5+SQ+EEFB3wK2lHXMvN3yKZ54aa5pH/N2yqq6T5oTHKjl0pVr6kZRTay1NFacrKS8ZhS+hhL5aqAvLeXpm3uzu+U+tkw0dMz87ep45n+0HI2wnHU3F5OqFlYDtX7MnmiYJPXE4jSrQqSwjSWhrjiXwdrKXapTRuGtKJGvIm/PmMr63O9YPdDQMTPlVE++/jCZgHr2RcD88T7m5bUUFNtfBtBPiHIzJ5NP5KkCbTVpUmFhdUvzHLQagdZPoCu9fitVnTIKb0bZozjI92//iymjBc+1+o6NETDktxAOP5PN958ddEjgKzJ7Yi+7vuMAGiF5YnEa8XM2kpiazWsJkdTzU2Xc6lBYYVUmS/McdHpJg3r+Jr9/a8ZjCoW3oEbydtizcjlzX5/EhiGQHw9RBzT89/UdRIX3q9FxjaIxe/lBmyN6Y70wu6CYxxen8UHSEUpKVcKmOuhKyxdQreXZLxTrSHtJ2UwofAM1kreCviCfRS9M4vZ1k1gyFlrkwI/xX5H2fWmNBd5IQmwosyf2QlOFgfmRs4VOObevEKj1452pMbwzNYbQsry5n7D+DTUXdmt5dpV/V/gSSuQrUlLChn//hbh/tuAu/+U08m/I+x3+ScZXkltG3+X0081dcwi9GphXi5AgrWlyk7F4GhocyL//FG0S/IqYC7jyelfUBVS6xohez76v5vHMttmsDi2mY+MAvo6ZzbSJz1nsmLHlPlkVVGte1TGf2GTNJG5Kn9BKM1IrCrjyelfUBZTIA6dWLeaFpX/hv+3OE9zCj393fJBH7/wPAVrLo0Fb7pNVFQhr1giK8mTOGW/xfWsmcUkZubwxOdKugKsJTQpfp06LfEHKNt5YcA//aXkM2sDM4JuZ9cCXhDRsbvNzttwnqyoYFa0R4LoHilZzvfCqsIwtkzgl4ApFHRX5qyeO8sHb03g9YA/5beBu/z68OuMbOrTs5tDn7blPVgV7KYOwWSurfMzaQgPUxj3IViHV2pOQKp4qFAbqlMjr8/P437zp/LNoBZlNJWOvhfHmtC+JDh9apeNURVgcyd3bGnGGBGktzsD0BGTFNe9cxLT+7a1us/QkpIqnCsV16obIX73Khncf5+nMT9nbspRY2ZRPx7zLjQPvrNbhHBUWe7l74w0gu6AYvzIDs4o2tVeq6XVTG7jaRsdPCKb1b1/J/dEcVTxVKGzj2yKv17Pvy7d4ZvsrrG5X1jHT+xWmTZhl1WPGERwVFnsrR5nfAIzGYxVvBMV1NCkfGhzItlkjHdpX5d4VCusI6UGuhnFxcTI5Odkpxzr187e8sPQx/ts+j2CdH893+z9Dx4x/QoCnwQAAB9JJREFUgFOO7widZq20mM0Q2O+q8RMCvZQ+a0YWqPVjSp9Qqz48AjhupaNGoVCURwiRIqWMs7TN50byBXu28Man9/CfVpnQFp4KGc+zD3xJSINmtR6Lrdy9vSKtr1sKG9scQ1XhVKFwKT4z4/Xq8SO8/be+dPl+KHPbZjI1II7Dfz3EW39f4RaBB9szKmsiYn4aQZC2Zj86Gw0rtUZOQbGadapQuBifGMmnrPyU29bPILMZjC3txJt3fEl0tyHuDstu7r6qi4gY0zzGY8S+srbanTc1eVAI1Gq4VqqvcQ9/2+BAVThVKFyMT4h82KCbCdvRjgVj5jC6mh0zrsJaUbDiIiLG7ho/G8sEVhTAAje1VlalGCww+LgXllyz6tGuCqcKhevwCZFvFhJK0mun3BqDI+2QFbEkbraWCazYeRNczR76QK0fAVqNxc8a2xYX7zlVTpSrg3mHjLO8fhQKRdXwCZF3NxWF2Vo7pCPYWybQvAXz8pVrVY41OFBrWkrQUq+/0e+lpgJvyQxMibpCUfv4TOHVnVjqhzdiLsqOkhAbyrZZI7FWG80pKDYIsQWP4uBALcEVlrkz5+o1vekcb0yOtLgCUnXsGUKCtGo1JYXCA1EjeRs4mmKwJ4rVtRO21oKpEcJqj/2FYh3zp8ZYTfkU60p5fHEac9ccYubYcIsTjqrqjBmo9TNZ/yoUCs+iRiN5IcTtQoiDQgi9ECKuwrZnhRC/CyEOCSHG1izM2seYgskuKEZyPfVivkaoEXvtkNVtl7TUXgiGdJC1Ub6xY8W4mIY1bF2PvfbFuwZ0UKN2hcJLqGm65gAwGdhs/qYQoifwZ6AXMA74UAhhf9VqD8KeJYE51sQYatbzbRRrSy6MEioJfcWOFWurIxmxdj0JsaGEBFlO+YQGB/JaQiTbZo3k+JzxbJs1Ugm8QuHB1EjkpZS/SSktJZwnAf+TUl6VUh4HfgecszBqLVEVO2Hz/DZct8Z1xig3ITYUvZWWSll2Dmsjals3HyPWrvOlW3rVeJJSYmo28XM20mnWSuLnbLT41KBQKFyLq3LyocBOs9dZZe9VQggxA5gB0KFDBxeFU3Wq6lPuyu4Ra7HYM/EyxvPyTwettlo2sVKkrekkJWeunqVQKKqPXZEXQqwHWlvY9E8p5Y/WPmbhPYvDUSnlAmABGAzK7MVTW9SmT7m9Am9NY7liY/JSYck1ElOzrU7Yqq4gO3P1LIVCUX3siryU8sZqHDcLMF/poR2QU43juI3amm7vyIi3JrHYau8E0JVKlwivM1fPUigU1cdV6ZrlwDdCiLeBtkBXYLeLzuUyamMCj6Mj3urG4oioukJ41bJ8CoVnUNMWyluFEFnAQGClEGINgJTyIPAd8CuwGviLlNJzlzhyI1Ud8Va1mOmIqLpCeJW7pELhGdS0u+YHKWU7KWV9KWUrKeVYs22vSym7SCnDpZSrah6qb2JNYK2tF+to774Rex02rhJeWzNqFQpF7aFmvLqZqhRVq1PMrJjPDw7SIqVhZqyrjcKUX41C4X6UyLuZqhRVq1vMVGKrUNRdlMh7AI6KsCpmKhSKqqJcKL0IVcxUKBRVRY3kvQi1VJ5CoagqSuS9DJVfVygUVUGlaxQKhcKHUSKvUCgUPowSeYVCofBhlMgrFAqFD6NEXqFQKHwYIa2sOuQOhBC5wAl3x1ENmgPn3B1ELaOuuW5Q167ZW6+3o5SyhaUNHiXy3ooQIllKGWd/T99BXXPdoK5dsy9er0rXKBQKhQ+jRF6hUCh8GCXyzmGBuwNwA+qa6wZ17Zp97npVTl6hUCh8GDWSVygUCh9GibxCoVD4MErknYwQ4ikhhBRCNHd3LK5GCDFXCJEhhNgnhPhBCBHs7phcgRBinBDikBDidyHELHfH42qEEO2FEElCiN+EEAeFEH93d0y1hRDCTwiRKoRY4e5YnIUSeScihGgPjAZOujuWWmIdECGljAIOA8+6OR6nI4TwAz4AbgJ6AtOEED3dG5XLuQY8KaXsAQwA/lIHrtnI34Hf3B2EM1Ei71zmA08DdaKaLaVcK6W8VvZyJ9DOnfG4iH7A71LKY1LKEuB/wCQ3x+RSpJSnpZR7y76+hEH0fH4RAyFEO2A88Jm7Y3EmSuSdhBBiIpAtpUx3dyxu4n5glbuDcAGhwCmz11nUAcEzIoQIA2KBXe6NpFZ4B8MgTe/uQJyJWhmqCggh1gOtLWz6J/AcMKZ2I3I9tq5ZSvlj2T7/xPCIv6g2Y6slhIX36sSTmhCiIbAUeFxKedHd8bgSIcQE4KyUMkUIMdzd8TgTJfJVQEp5o6X3hRCRQCcgXQgBhrTFXiFEPynlH7UYotOxds1GhBD3AhOAUdI3J11kAe3NXrcDctwUS60hhNBiEPhFUspl7o6nFogHJgohbgYCgMZCiK+llHe5Oa4aoyZDuQAhRCYQJ6X0Rjc7hxFCjAPeBoZJKXPdHY8rEEL4YygqjwKygT3AHVLKg24NzIUIw0jlSyBPSvm4u+OpbcpG8k9JKSe4OxZnoHLyiprwPtAIWCeESBNCfOzugJxNWWH5MWANhgLkd74s8GXEA3cDI8t+rmllI1yFF6JG8gqFQuHDqJG8QqFQ+DBK5BUKhcKHUSKvUCgUPowSeYVCofBhlMgrFAqFD6NEXqFQKHwYJfIKhULhw/x/M6Z++coP12EAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.scatter(x, y_noisy, label='empirical data points')\n",
    "plt.plot(x, y, color='black', label='true relationship')\n",
    "plt.plot(inputs, outputs, color='red', label='predicted relationship (cpu)')\n",
    "plt.plot(inputs, outputs_gpu.to_array(), color='green', label='predicted relationship (gpu)')\n",
    "plt.legend()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<a id=\"dask\"></a>\n",
    "## Dask\n",
    "\n",
    "Dask is a library the allows facillitates distributed computing. Written in Python, it allows one to compose complex workflows using basic Python primitives like integers or strings as well as large data structures like those found in NumPy, Pandas, and cuDF. In the following examples and notebooks, we'll show how to use Dask with cuDF to accelerate common ETL tasks and train machine learning models like Linear Regression and XGBoost.\n",
    "\n",
    "To learn more about Dask, check out the documentation here: http://docs.dask.org/en/latest/\n",
    "\n",
    "#### Client/Workers\n",
    "\n",
    "Dask operates by creating a cluster composed of a \"client\" and multiple \"workers\". The client is responsible for scheduling work; the workers are responsible for actually executing that work. \n",
    "\n",
    "Typically, we set the number of workers to be equal to the number of computing resources we have available to us. For CPU based workflows, this might be the number of cores or threads on that particlular machine. For example, we might set `n_workers = 8` if we have 8 CPU cores or threads on our machine that can each operate in parallel. This allows us to take advantage of all of our computing resources and enjoy the most benefits from parallelization.\n",
    "\n",
    "To get started, we'll create a local cluster of workers and client to interact with that cluster."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Dask Version: 1.2.2\n"
     ]
    }
   ],
   "source": [
    "import dask; print('Dask Version:', dask.__version__)\n",
    "from dask.distributed import Client, LocalCluster\n",
    "\n",
    "\n",
    "# create a local cluster with 4 workers\n",
    "n_workers = 4\n",
    "cluster = LocalCluster(n_workers=n_workers)\n",
    "client = Client(cluster)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's inspect the `client` object to view our current Dask status. We should see the IP Address for our Scheduler as well as the the number of workers in our Cluster. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table style=\"border: 2px solid white;\">\n",
       "<tr>\n",
       "<td style=\"vertical-align: top; border: 0px solid white\">\n",
       "<h3>Client</h3>\n",
       "<ul>\n",
       "  <li><b>Scheduler: </b>tcp://127.0.0.1:36379\n",
       "  <li><b>Dashboard: </b><a href='http://127.0.0.1:8787/status' target='_blank'>http://127.0.0.1:8787/status</a>\n",
       "</ul>\n",
       "</td>\n",
       "<td style=\"vertical-align: top; border: 0px solid white\">\n",
       "<h3>Cluster</h3>\n",
       "<ul>\n",
       "  <li><b>Workers: </b>4</li>\n",
       "  <li><b>Cores: </b>40</li>\n",
       "  <li><b>Memory: </b>540.95 GB</li>\n",
       "</ul>\n",
       "</td>\n",
       "</tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Client: scheduler='tcp://127.0.0.1:36379' processes=4 cores=40>"
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# show current Dask status\n",
    "client"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "You can also see the status and more information at the Dashboard, found at `http://<ip_address>/status`. You can ignore this for now, we'll dive into this in subsequent tutorials.\n",
    "\n",
    "With our client and cluster of workers setup, it's time to execute our first distributed program. We'll define a function called `sleep_1` that sleeps for 1 second and returns the string \"Success!\". Executed in serial four times, this function should take around 4 seconds to execute."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "import time\n",
    "\n",
    "\n",
    "def sleep_1():\n",
    "    time.sleep(1)\n",
    "    return 'Success!'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "CPU times: user 179 ms, sys: 34.9 ms, total: 214 ms\n",
      "Wall time: 4 s\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "\n",
    "for _ in range(n_workers):\n",
    "    sleep_1()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "As expected, our workflow takes about 4 seconds to run. Now let's execute this same workflow in distributed fashion using Dask."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [],
   "source": [
    "from dask.delayed import delayed"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "['Success!', 'Success!', 'Success!', 'Success!']\n",
      "CPU times: user 64.5 ms, sys: 15.3 ms, total: 79.8 ms\n",
      "Wall time: 1.02 s\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "\n",
    "# define delayed execution graph\n",
    "sleep_operations = [delayed(sleep_1)() for _ in range(n_workers)]\n",
    "\n",
    "# use client to perform computations using execution graph\n",
    "sleep_futures = client.compute(sleep_operations, optimize_graph=False, fifo_timeout=\"0ms\")\n",
    "\n",
    "# collect and print results\n",
    "sleep_results = client.gather(sleep_futures)\n",
    "print(sleep_results)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Using Dask, we see that this whole workflow takes a little over a second - each worker is truly executing in parallel!"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<a id=\"daskcudf\"></a>\n",
    "## Dask cuDF\n",
    "\n",
    "In the previous example, we saw how we can use Dask with very basic objects to compose a graph that can be executed in a distributed fashion. However, we aren't limited to basic data types though. \n",
    "\n",
    "We can use Dask with objects such as Pandas DataFrames, NumPy arrays, and cuDF DataFrames to compose more complex workflows. With larger amounts of data and embarrasingly parallel algorithms, Dask allows us to scale ETL and Machine Learning workflows to Gigabytes or Terabytes of data. In the below example, we show how we can process 100 million rows by combining cuDF with Dask.\n",
    "\n",
    "Before we start working with cuDF DataFrames with Dask, we need to setup a Local CUDA Cluster and Client to work with our GPUs. This is very similar to how we setup a Local Cluster and Client in vanilla Dask."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Dask Version: 1.2.2\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/conda/envs/rapids/lib/python3.7/site-packages/distributed/bokeh/core.py:74: UserWarning: \n",
      "Port 8787 is already in use. \n",
      "Perhaps you already have a cluster running?\n",
      "Hosting the diagnostics dashboard on a random port instead.\n",
      "  warnings.warn(\"\\n\" + msg)\n"
     ]
    }
   ],
   "source": [
    "import dask; print('Dask Version:', dask.__version__)\n",
    "from dask.distributed import Client\n",
    "# import dask_cuda; print('Dask CUDA Version:', dask_cuda.__version__)\n",
    "from dask_cuda import LocalCUDACluster\n",
    "\n",
    "\n",
    "# create a local CUDA cluster\n",
    "cluster = LocalCUDACluster()\n",
    "client = Client(cluster)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's inspect our `client` object:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table style=\"border: 2px solid white;\">\n",
       "<tr>\n",
       "<td style=\"vertical-align: top; border: 0px solid white\">\n",
       "<h3>Client</h3>\n",
       "<ul>\n",
       "  <li><b>Scheduler: </b>tcp://127.0.0.1:44992\n",
       "  <li><b>Dashboard: </b><a href='http://127.0.0.1:45828/status' target='_blank'>http://127.0.0.1:45828/status</a>\n",
       "</ul>\n",
       "</td>\n",
       "<td style=\"vertical-align: top; border: 0px solid white\">\n",
       "<h3>Cluster</h3>\n",
       "<ul>\n",
       "  <li><b>Workers: </b>8</li>\n",
       "  <li><b>Cores: </b>8</li>\n",
       "  <li><b>Memory: </b>540.95 GB</li>\n",
       "</ul>\n",
       "</td>\n",
       "</tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Client: scheduler='tcp://127.0.0.1:44992' processes=8 cores=8>"
      ]
     },
     "execution_count": 27,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "client"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "As before, you can also see the status of the Client along with information on the Scheduler and Dashboard.\n",
    "\n",
    "With our client and workers setup, let's create our first distributed cuDF DataFrame using Dask. We'll instantiate our cuDF DataFrame in the same manner as the previous sections but instead we'll use significantly more data. Lastly, we'll pass the cuDF DataFrame to `dask_cudf.from_cudf` and create an object of type `dask_cudf.core.DataFrame`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Dask cuDF Version: 0.7.2+0.g3ebd286.dirty\n"
     ]
    }
   ],
   "source": [
    "import dask_cudf; print('Dask cuDF Version:', dask_cudf.__version__)\n",
    "\n",
    "\n",
    "# identify number of workers\n",
    "workers = client.has_what().keys()\n",
    "n_workers = len(workers)\n",
    "\n",
    "# create a cuDF DataFrame with two columns named \"key\" and \"value\"\n",
    "df = cudf.DataFrame()\n",
    "n_rows = 100000000  # let's process 100 million rows in a distributed parallel fashion\n",
    "df['key'] = np.random.binomial(1, 0.2, size=(n_rows))\n",
    "df['value'] = np.random.normal(size=(n_rows))\n",
    "\n",
    "# create a distributed cuDF DataFrame using Dask\n",
    "distributed_df = dask_cudf.from_cudf(df, npartitions=n_workers)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "---------------\n",
      "Type of our Dask cuDF DataFrame: <class 'dask_cudf.core.DataFrame'>\n",
      "---------------\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/conda/envs/rapids/lib/python3.7/site-packages/distributed/worker.py:3101: UserWarning: Large object of size 300.00 MB detected in task graph: \n",
      "  (<cudf.DataFrame ncols=2 nrows=12500000 >, 5)\n",
      "Consider scattering large objects ahead of time\n",
      "with client.scatter to reduce scheduler burden and \n",
      "keep data on workers\n",
      "\n",
      "    future = client.submit(func, big_data)    # bad\n",
      "\n",
      "    big_future = client.scatter(big_data)     # good\n",
      "    future = client.submit(func, big_future)  # good\n",
      "  % (format_bytes(len(b)), s)\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "   key                value\n",
      "0    0   0.3313946586367762\n",
      "1    0  -0.4015253450495195\n",
      "2    0  -1.4236005817346722\n",
      "3    0  -0.8342524790363371\n",
      "4    0  -0.7593642316512079\n"
     ]
    }
   ],
   "source": [
    "# inspect our distributed cuDF DataFrame using Dask\n",
    "print('-' * 15)\n",
    "print('Type of our Dask cuDF DataFrame:', type(distributed_df))\n",
    "print('-' * 15)\n",
    "print(distributed_df.head())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The above output shows the first several rows of our distributed cuDF DataFrame.\n",
    "\n",
    "With our Dask cuDF DataFrame defined, we can now perform the same `sum` operation as we did with our cuDF DataFrame. The key difference is that this operation is now distributed - meaning we can perform this operation using multiple GPUs or even multiple nodes, each of which may have multiple GPUs. This allows us to scale to larger and larger amounts of data!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2007.5109850552367\n"
     ]
    }
   ],
   "source": [
    "aggregation = distributed_df['value'].sum()\n",
    "print(aggregation.compute())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<a id=\"daskcuml\"></a>\n",
    "## Dask cuML (coming soon!)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<a id=\"conclusion\"></a>\n",
    "## Conclusion\n",
    "\n",
    "In this notebook, we showed at a high level what each of the packages in the RAPIDS are as well as what they do.\n",
    "\n",
    "To learn more about RAPIDS, be sure to check out: \n",
    "\n",
    "* [Open Source Website](http://rapids.ai)\n",
    "* [GitHub](https://github.com/rapidsai/)\n",
    "* [Press Release](https://nvidianews.nvidia.com/news/nvidia-introduces-rapids-open-source-gpu-acceleration-platform-for-large-scale-data-analytics-and-machine-learning)\n",
    "* [NVIDIA Blog](https://blogs.nvidia.com/blog/2018/10/10/rapids-data-science-open-source-community/)\n",
    "* [Developer Blog](https://devblogs.nvidia.com/gpu-accelerated-analytics-rapids/)\n",
    "* [NVIDIA Data Science Webpage](https://www.nvidia.com/en-us/deep-learning-ai/solutions/data-science/)\n"
   ]
  },
  {
   "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.7.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}