{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\n",
    "Custom Workloads with Dask Delayed\n",
    "==================================\n",
    "\n",
    "<img src=\"http://dask.readthedocs.io/en/latest/_images/dask_horizontal.svg\" \n",
    "     width=\"30%\" \n",
    "     align=right\n",
    "     alt=\"Dask logo\">\n",
    "\n",
    "*Because not all problems are dataframes*\n",
    "\n",
    "This notebook shows using [dask.delayed](http://dask.pydata.org/en/latest/delayed.html) to parallelize generic Python code.  \n",
    "\n",
    "Dask.delayed is a simple and powerful way to parallelize existing code.  It allows users to delay function calls into a task graph with dependencies.  Dask.delayed doesn't provide any fancy parallel algorithms like Dask.dataframe, but it does give the user complete control over what they want to build.\n",
    "\n",
    "Systems like Dask.dataframe are built with Dask.delayed.  If you have a problem that is paralellizable, but isn't as simple as just a big array or a big dataframe, then dask.delayed may be the right choice for you."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Start Dask Client for Dashboard\n",
    "\n",
    "Starting the Dask Client is optional.  It will provide a dashboard which \n",
    "is useful to gain insight on the computation.  \n",
    "\n",
    "The link to the dashboard will become visible when you create the client below.  We recommend having it open on one side of your screen while using your notebook on the other side.  This can take some effort to arrange your windows, but seeing them both at the same is very useful when learning."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "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 style=\"text-align: left;\">Client</h3>\n",
       "<ul style=\"text-align: left; list-style: none; margin: 0; padding: 0;\">\n",
       "  <li><b>Scheduler: </b>tcp://127.0.0.1:44986</li>\n",
       "  <li><b>Dashboard: </b><a href='/user/eglez/proxy/8787/status' target='_blank'>/user/eglez/proxy/8787/status</a></li>\n",
       "</ul>\n",
       "</td>\n",
       "<td style=\"vertical-align: top; border: 0px solid white\">\n",
       "<h3 style=\"text-align: left;\">Cluster</h3>\n",
       "<ul style=\"text-align: left; list-style:none; margin: 0; padding: 0;\">\n",
       "  <li><b>Workers: </b>1</li>\n",
       "  <li><b>Cores: </b>4</li>\n",
       "  <li><b>Memory: </b>34.36 GB</li>\n",
       "</ul>\n",
       "</td>\n",
       "</tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Client: 'tcp://127.0.0.1:44986' processes=1 threads=4, memory=34.36 GB>"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from dask.distributed import Client, progress\n",
    "client = Client(threads_per_worker=4, n_workers=1)\n",
    "client"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Create simple functions\n",
    "\n",
    "These functions do simple operations like add two numbers together, but they sleep for a random amount of time to simulate real work."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import time\n",
    "import random\n",
    "\n",
    "def inc(x):\n",
    "    time.sleep(random.random())\n",
    "    return x + 1\n",
    "\n",
    "def dec(x):\n",
    "    time.sleep(random.random())\n",
    "    return x - 1\n",
    "    \n",
    "def add(x, y):\n",
    "    time.sleep(random.random())\n",
    "    return x + y "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can run them like normal Python functions below"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "CPU times: user 66.4 ms, sys: 13.6 ms, total: 80 ms\n",
      "Wall time: 1.91 s\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "3"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "%%time\n",
    "x = inc(1)\n",
    "y = dec(2)\n",
    "z = add(x, y)\n",
    "z"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "These ran one after the other, in sequence.  Note though that the first two lines `inc(1)` and `dec(2)` don't depend on each other, we *could* have called them in parallel had we been clever."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Annotate functions with Dask Delayed to make them lazy\n",
    "\n",
    "We can call `dask.delayed` on our funtions to make them lazy.  Rather than compute their results immediately, they record what we want to compute as a task into a graph that we'll run later on parallel hardware."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "import dask\n",
    "inc = dask.delayed(inc)\n",
    "dec = dask.delayed(dec)\n",
    "add = dask.delayed(add)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Calling these lazy functions is now almost free.  We're just constructing a graph"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Delayed('add-3e8e3f73-0449-4351-a58d-5dca1bdaf28d')"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "x = inc(1)\n",
    "y = dec(2)\n",
    "z = add(x, y)\n",
    "z"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Visualize computation\n",
    "\n",
    "You will need graphviz installed for this to work"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAb4AAACtCAYAAADPuBpoAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deViN+f8/8OdpX04rbYSynFJEKrRoJ0aISWUJw2gMk2Tsw8duZMbSmLGv2VUTIoxIpTL2vYXRZKuklMppv39/zK++Y0oq9zn3qV6P6+q6cs7d+/28zxiv7vt+LzyGYRgQQgghrUOlFNcJCCGEEHGiwkcIIaRVocJHCCGkVWkWhS8/Px+lpaVcxyCEENICyHAdAACeP3+OhIQEpKSkIDU1FampqcjMzERRURGKiopqjpOVlQWfz0fbtm1haGgIIyMjGBsbo1evXujXrx9kZCTidAghhEgwHhejOsvKyhAVFYWzZ88iJiYGjx8/hqysLDp37gxjY2MIBAJ06NABysrKUFFRgZqaGoRCIYqLi1FUVITXr1/j6dOnNYUyLy8PfD4fAwYMgLOzM7y8vNCxY0dxnxYhhBDJVynWwnfjxg3s378fR44cwdu3b9GvXz84OzvD2dkZ1tbWUFRUbFK7jx8/RkxMDC5duoTo6Gi8ffsWDg4OmDBhAkaPHg1lZWWWz4QQQkgzJZ7Cd+XKFQQFBeH06dMwNjaGt7c3Jk6cCENDQ9b7qqysRExMDEJCQhAeHg4lJSXMmDEDAQEB0NDQYL0/QgghzYpoC9/Vq1cRGBiIq1evwtXVFYsWLYKTk5OouqslJycHwcHB+PXXX8Hj8TB37lx8//33kJeXF1sGQgghEkU0E9hzc3Ph5+cHW1tbKCsr4+rVq7hw4YJYix4AaGlpYdWqVcjIyMDs2bOxZs0amJmZ4cKFC2LNQQghRHKwXvgiIyNhbGyMM2fO4PDhw4iOjka/fv3Y7qZR1NTUsGTJEiQnJ6NHjx4YNGgQJk2a9MGIUUIIIa0Da4WvrKwMs2fPxogRIzBixAikpKTA29ubreZZ0aFDB4SHhyMyMhJRUVGwsrLCvXv3uI5FCCFEjFgpfDk5OXBwcMDOnTsREhKCXbt2QUVFhY2mRcLd3R23b9+GtrY2+vfvj6NHj3IdiRBCiJh89ozvjIwMuLm5oby8HNevX4exsTEbuUSuffv2uHjxIubMmYNx48YhJycH/v7+XMcihBAiYp9V+JKTkzFw4EC0bdsWMTEx0NPTYyuXWMjIyGDTpk3Q09NDQEAAcnJysGLFCq5jEUIIEaEmF76MjAwMHDgQhoaGOH36NNTU1NjMJVbz58+HlpYWpk6dCmVlZcyfP5/rSIQQQkSkSYXvzZs3GDx4MDQ1NXHq1KlmXfSqTZ48GaWlpZgxYwbatm2LKVOmcB2JEEKICDS68JWVlWH48OEoKytDTExMi1oN5dtvv8WLFy8wbdo0GBgYwMXFhetIhBBCWNbolVu+//577NixAzdu3ICRkZGocnGGYRh4e3sjPj4et2/fhq6uLteRCCGEsKdxS5ZFRUXB3d0d+/fvh6+vryiDcaqoqAiWlpbQ09NDdHQ0pKWluY5ECCGEHQ1fsiw3NxcTJ07E5MmTW3TRAwA+n4/Dhw8jKSkJmzdv5joOIYQQFjX4is/Pzw+RkZFISUlpEYNZGmLp0qXYuHEjkpOT0b59e67jsCYsLAznz5/nOoZEMzQ0xKJFi7iOQQhhX8NudV6/fh39+/fH4cOHJW4ZMlEqLS2FmZkZzM3NW9TqLnPnzsWOHTvQt29frqNIpMePH0NLSwvXr1/nOgohhH0NK3w2NjZQUlJCdHS0OEJJlDNnzsDd3R0JCQmwsbHhOg4r5s6di8uXL9M/7B9Bnw8hLdqnn/FdvHgRSUlJrXZFk6FDh8LOzg6rV6/mOgohhBAWfLLwrV69Gq6uri3maqcpFi5ciKioKNy4cYPrKIQQQj5TvYXv5s2biImJafUP+b/44gtYWFhg/fr1XEchhBDymeotfHv37kX37t3FvnO6JJo2bRpOnDiB/Px8rqMQQgj5DB8tfGVlZTh27BgmTJggzjwSy8vLC1JSUggNDeU6CiGEkM/w0cJ39uxZ5OXlYdy4ceLMI7FUVVUxfPhwHDhwgOsohBBCPsNHC19UVBT69euHDh06iDOPRPP09ERiYiLd7iSEkGbso4Xv0qVLtDvBfzg5OYFhGMTFxXEdhRBCSBPVWfieP3+OJ0+eiHRQy9GjR8Hj8cDj8aCgoCCyftikqakJMzMzxMTEcB2FEEJIE9VZ+BISEiArKwtra2uRdezj4wOGYZrdVaWjoyOuXLnCdQxCCCFNVGfhS0lJQZcuXaCoqCjuPBLP1NQUKSkpaOQ2hoQQQiTERwtfS9xklg1GRkYoKirCq1evuI5CCCGkCeosfGlpaRAIBOLO0ixU/0KQkpLCcRJCCCFNUWfhy8zMhL6+PqsdpaSkwMPDA2pqalBWVsaAAQPqfVaWk5ODmTNnwsDAAHJyctDS0sKoUaNw586dWsfm5uZi9uzZ6NKlC+Tl5aGvrw9XV1fs27cPQqGQ1fPQ1taGnJwcMjMzWW2XEEKIeNRZ+AoLC6GqqspaJ0+ePIG1tTVu3LiBsLAwZGdnY8uWLVi5ciX++uuvWsdnZmbCysoKx48fx5YtW5CXl4fLly8jLy8P1tbWSEpKqjk2KysLVlZWOHLkCIKDg/HmzRvcvHkTjo6O+Oqrr7B9+3bWzqOaiooKCgsLWW+3KZYvX44///yT6xiEENJs1Cp8VVVVeP/+Pfh8PmudLFq0CPn5+QgODsbAgQPB5/PRs2dP7N27t84rp4ULFyIjIwMbNmzAF198AT6fD1NTUxw9ehQMw8Df3/+DY9PT0xEcHAx3d3eoqKhAR0cHixcvxuDBg1k7h3+TpMK3bds29O/fH6ampti6dSvevXvHdSRCCJFotQpfcXExGIZhtfCdO3cOAODm5vbB6+3atavzWeKJEycgJSUFd3f3D17X1dWFqakpbt68iRcvXgAAIiIiAABDhgyp1c7Zs2cxa9YsVs7h3ySp8BUXFwMAkpOT4e/vD21tbUyZMoU2USWEkI+Q+e8LcnJyAP5ZpJoNpaWlKCwshIKCQp3FVFtbG2lpaR8cX1BQAABQU1P7aLuPHz+GlpYWCgoKoKCgABUVFVbyNkRpaSlyc3Pxxx9/QFFREUpKSlBXV4eioiIUFRWhoaEhtizVzzAZhkFlZSUqKytx4MAB7NmzB926dcOUKVPg5+cn1kyEECLJahU+eXl5yMnJsXZFIy8vX3OFVFRUVKv45eXl1TpeXV0dRUVFEAqFkJGpFfEDampqKCgoQGFhodiKX2FhIW7evImtW7fWe5yCggI0NDSgqKjYqO8/9l6bNm0gLy9f035ZWRkqKipq9VteXg7gn2erP/zwA5YsWQJvb2/MmjULFhYW7H4YhBDSzNRZVfh8PoqKiljrZMiQITh+/DjOnTsHT0/PmtffvHmD1NTUWsePGjUKe/bsQUJCAhwcHD54LygoCL/99huePn0KGRkZjBw5Evv27UNUVBS8vb0/ONbc3ByOjo7YuHEja+cCAO/evcO0adNw4cIFCIVCFBcXo6CgAEKhEO/fv0d+fj6EQiGEQiHevn1b831+fj7ev38PoVCIzMzMmu8LCgpQXFwMoVD4yWd08vLyNVeY1VfnH/Pvq8CjR4/i4MGDMDMzg66uLqqqqtj8SAghpNmos/BpaGggNzeXtU7WrFmD6OhozJo1C2pqarC2tsazZ88we/Zs8Pn8mlub1X788UfExsZi8uTJ+PXXX2FjY4PKykqEhoZixYoV2Lt3b82VYPWxgYGB4PP5sLe3R0FBAdasWYPMzEwEBgaydh4AUFJSgvfv30NdXR2qqqqsjn6tVn21W1hYiMLCQgiFQhQVFeHdu3cfFNrs7Ow6f3GoS/WV4b1793Dv3j3IysrizJkzGDp0KOv5CSFEojF1GDRoEDNp0qS63mqy1NRUxsPDg1FVVWUUFRUZKysr5vTp04yLiwsDgAHATJkypeb43NxcZvbs2Uznzp0ZWVlZRktLixk0aBBz4cKFWm2/efOGmTVrFmNoaMjIysoyenp6jI+PD5OWlsbqOTAMw9y/f58BwNy/f5/1thsrOTm55rOr70tWVpYBwGhpaTFTp05lRo4cyfTp04fr+BJrzpw5jKWlJdcxCCGiUcFjmNqLTs6cORM3b95EQkKCmMpv8xEeHg4vLy8UFxdzvqvEjRs3YGVlVet1aWlp8Hg8VFRUoFu3bhg1ahTc3d1ha2sLHo+HuXPn4vLlyzTy8yPo8yGkRaus81ankZERDh06JO4wzUJaWho6derEedED/m8qAwDIyMigsrIScnJycHZ2hoeHB9zd3dGuXTsOExJCiOSps/CZm5sjLy8Pjx8/Rrdu3cSdSaJdvXoV5ubmXMcAgJoBSFpaWhg5ciSGDRsGFxcX2lWDEELqUWfh69u3L1RUVBATE0OF718qKysRFxeHlStXch0FAGBiYoIbN26gT58+4PF4XMchBMA/o4mfPXuG9PR0vHnzBsXFxSguLkZhYSE0NDSgrKwMZWVlaGlpQSAQQEtLi+vIpJWps/DJyMjA1tYWMTEx8PPzE3cmiXXz5k3k5+fD2dmZ6ygAAENDQxgaGnIdg7RyGRkZiImJwaVLl3D//n2kpaXh/fv3Ne8rKipCWVkZqqqqePv2LYqLiz9YIENTUxMCgQBWVlZwcnKCg4MDNDU1uTgV0kp8dHa4i4sL1q5di/LycsjKyoozk8Q6d+4c9PT00L17d66jEMKpBw8eICQkBOHh4Xj69CkUFRVhbW0NNzc3+Pv7w9jYGF26dIG2tnaddyPKy8uRlZWF1NRUpKWlISUlBYmJidiyZQsYhoG5uTl8fHwwduxYek5NWPfRwufl5YX58+fj7NmzGD58uDgzSazDhw/Dx8eHbiuSVkkoFGLv3r3YtWsXbt++DQMDA4wZMwaDBg1C//79GzXgS1ZWFh06dECHDh3g6upa8/rbt28RGxuLs2fPYs2aNViwYAEGDhyIadOmYfjw4fT/HmFFndsSAUDHjh3h4OCAAwcOiDOPxEpKSkJqaip8fX25jkKIWBUWFmLdunUwNDTEnDlzYG5ujtjYWDx9+hRr1qyBo6Mja6OcNTQ04OHhge3bt+PVq1c4evQopKWlMXLkSPTq1QtHjhxBZWUlK32R1uujhQ8AJkyYgMjISLx580ZceSTWvn37YGpqKjEjOgkRtaqqKuzcuROGhoZYtWoVJk2ahPT0dOzevRv29vYiv/pSUFCAp6cnTp8+jXv37qFnz57w9fWFmZkZYmJiRNo3adnqLXxeXl5QVVXFpk2bxJVHImVmZiIkJAQzZszgOgohYnH79m3Y2tpi+vTpmDBhAjIyMrB27Vro6OhwkqdHjx44dOgQkpOT0blzZ7i4uGD8+PHIysriJA9p3uotfEpKSpg1axY2b96M/Px8cWWSOD///DPU1dUxadIkrqMQIlIMwyA4OBj9+vWDjIwMbt68iQ0bNkjMtlbdunVDZGQkTp48icTERBgbGyMsLIzrWKSZqbfwAcCMGTMgJSWFzZs3iyOPxHn9+jW2b9+OefPm0cTwVqSkpKTVXU3k5ORgyJAhmDdvHoKCghAXFwczMzOuY9Vp2LBhePDgAb788kt4eXlhzpw5NdtxEfIpnyx8ampqmDt3LoKCgvD8+XNxZJIo8+bNg6amJr755huuoxAxevXqFfT19TFw4EDs3bu3xd/xePjwIfr06YPU1FTEx8cjMDBQ4kdQKikpYffu3di/fz+2b98OV1fXFv/fibDjk4UPAObMmYMOHTqwvsWPpEtISEBISAg2bdoEJSUlruMQMTIwMEBERAT09PTg7+8PXV1dDBs2DCEhIazuVSkJEhMTYW9vD0NDQ9y6dQt9+/blOlKj+Pr64urVq0hPT4e9vT1evXrFdSQi4RpU+OTk5PDLL78gPDwcp0+fFnUmiVBaWopvv/0WgwcPxqhRo7iOQ8RMSkqqptC9fPkSO3bsAAB8/fXX0NbWxrBhwxAaGvrBCiTN0fnz5zFw4EAMGDAA58+fl5hneY1lamqKhIQEVFRUwNbWFn/99RfXkYgEa1DhA4CBAwdi4sSJmDx5Ml6+fCnKTBLh+++/x7Nnz/Dbb79xHYVwTE1NrWZqT1ZWFrZt24aSkhJ4e3tDV1e35r3qzX6bi6SkJIwaNQqenp4ICwtr9s+wO3TogPj4eLRp0wZubm6t7hktabg69+P7mOLiYlhZWaFt27a4dOlSzS7oLU1YWBhGjx6NQ4cOYezYsVzHYR3tN1e/hn4+z58/x++//47Q0FAkJCSgffv2+PLLLzF69OiavQ8l1cOHD2Fvbw9bW1v8/vvvLer/5Tdv3sDOzg5ycnKIi4uDuro615FEgmEYREdHcx1D4vXq1Qva2tr/fqmyzh3Y63P37l1GUVGRCQgIYHNHXIlx7949Rk1NjZk+fTrXUUSGdhivX1M+n0ePHjFLly5lBAIBA4Dp1KkTM3/+fCY5OVlEKZvu9evXTIcOHRh7e3vm/fv3XMcRiadPnzJ6enrMoEGDmMrKSq7jiER5eTkDgL4+8RUeHv7fj66i0b/mmZmZYc+ePRg3bhz09PQwf/78xjYhsTIyMjBkyBCYm5tjw4YNXMchzUj37t2xbNkyLFu2DA8fPkRoaChCQkIQFBQEExMTjB49Gr6+vujSpQunOauqqjBhwgTIyMjgxIkTzf725scYGhri1KlTsLOzw+rVq7FkyRKuI4lMSEgI3N3duY4hcSorKz+65VWT7m/4+PggJycHAQEB0NLSwuTJkz8roCTIycnB4MGD0bZtW5w4cQLy8vJcRyLNlKmpKUxNTfG///0PiYmJCA0NxdatW7F8+XJYWFjA19cXPj4+TVoF5dGjRzAxMWlytrVr1+LixYuIjY1ttgNZGsrS0hLr1q1DYGAgbGxs4OLiwnUkkVBWVm7x/y2bor5n7g0e3PJf/v7+WLx4MaZOnYpffvmlqc1IhIyMDAwYMADl5eU4e/Ys1NTUuI5EWgApKSnY2dkhODgYr169woULF2BiYoIlS5agffv2sLOzw44dO/Du3bsGtVdRUQF7e3t89913TZqsfevWLSxduhTr1q2DtbV1o3++OZo5cyY8PDwwadIkFBYWch2HSIgmFz4AWLFiBTZs2IBZs2ZhwYIFYBo+TkZiPHz4EHZ2dpCVlUVsbCz09PS4jkRaIGlpabi6uiIkJASvX79GREQE2rVrB39/f+jo6NRMnSguLv5oG9HR0cjNzcW2bdvg4ODQqFGLVVVV+O6779C/f38EBASwcUrNxo4dO1BSUoLly5dzHYVIiM8qfAAQEBCAPXv2YP369fD29m7wb6+SIDQ0FDY2NujatSuuXLmC9u3bcx2JtAIKCgoYNmwYjh8/juzsbGzfvh0AMGXKFGhra8PLywuRkZG1ruoOHToEWVlZVFZW4saNGzAzM8PVq1cb1OfOnTtx/fp1/PrrrxI92lQU2rRpg9WrVyM4OBh3797lOg6RAI2azlCfixcvYvz48eDz+Th27Bj69OnDRrMiUVpaiu+//x6//fYbpk+fjg0bNrSqZ3pz587Fvn374OHhwXUUiXTt2jXIycmJfbpHbm4uwsPDERISgsTERGhoaGDo0KGYMGECbG1toaWl9cEVobS0NHg8Hn777Tf4+fl9tN13796hc+fOmDhxItavXy+OU5E4VVVVsLGxgbKyMi5evMh1HFZUVFRAVlYW4eHhtMhGHer5fCpZm7zj4uKC27dvY/z48bCxscGCBQuwYMEC1jaoZMuVK1cwffp0PHv2DMePH8fo0aO5jiR2AoEAvXv3xt9//811FImkra0NIyMjsffbpk0b+Pn5wc/PD0+fPsWRI0dw5MgRHDhwANra2nj//v0Hx1dvyDpt2jTExcVh586ddY7S3LJlC8rLy7F48WKxnIckkpKSwtq1a+Hk5IQrV67Azs6O60iES2zPLamsrGQ2btzIqKqqMl26dGHOnDnDdhdNkpWVxUycOJHh8XjMkCFDmKdPn3IdiZAGuXfvHtOrVy9GRkbmo3OVZGRkmF69ejEZGRkf/KxQKGT09PSYBQsWcJRestjZ2TFffPEF1zFYUT2Pr455aoSp9/Op+OxnfP8lJSWFWbNmISUlBTY2Nhg6dChsbW0RGRnJyeCX169fY9myZRAIBLhw4QL27duHqKgoGBoaij0LIU3RqVMnJCcn1zs8u6KiAo8ePULv3r0/2J189+7dKCgoaHULzH/MDz/8gLNnz+LOnTtcRyEcYr3wVdPT00NISAgSEhKgrq6OESNGwMrKCnv27BHLAJikpCR888036NSpE3bu3Illy5bh8ePHmDBhgsj7JoRNv//+e4PWAS0vL8e7d+/g6uqKoKAgAMCuXbswduzY/y7Z1GoNHjwYxsbG2L17N9dRCIdEVviq2djY4MyZM7h58ya6deuGGTNmQFdXF2PGjEF4eDhyc3NZ6aeyshLXr1/HihUrIBAIYGNjg4SEBGzYsAFPnz5FYGAgbS1EmqVDhw41+NjKykpUVVVhwYIF8PDwwJ07d+iXvf/w9fXFkSNHmv3OGqTpWBvV2VD5+fkIDQ3FgQMHkJiYCIZh0KtXLzg4OMDU1BRGRkYwMjKq9zfU0tJSPH78GKmpqUhLS8PVq1cRFxeH/Px86OnpwdvbG76+vhI9spSQhsjJyYGenl7NQJbGkpeXR1paGjp27Mhysubr2bNnMDQ0REREBIYPH851nCajUZ31E8uozoZSV1fH1KlTMXXqVOTn5yMuLg6XLl1CfHw8du3aVbPJp5ycHPh8PtTV1aGiogKhUIjCwkIUFRXVrMAgJSWFTp06oXfv3li5ciWcnJxgamoq7lMiRGT4fD7S0tKa9LOOjo4YPnx4i12Ps6k6duwIOzs7hIeHN+vCR5qO071I1NXVMXz48A/+8r148QKpqanIyspCUVER8vPz8e7dOygpKYHP54PP56NNmzbo3LkzBAKBxE2XIIRNioqK6Ny5c6N/LiMjA8+fP8fo0aM/ulBvazZw4EBs3bqV6xiEIxK3CZe+vj709fW5jkFIs3bx4kUoKSmhf//+XEcBABw9ehRjxowB8M/t15KSElaPbyxnZ2csWbIEaWlpEAgErLZNJJ/IB7cQQsTv8uXLsLGxkZgViXx8fMAwTIN3SGjs8Y3Vt29fqKiofDD1g7QeVPgIaYHu3LkDKysrrmNILBkZGfTu3ZvW7mylqPAR0sJUVVXhyZMnnCy71pwYGRkhNTWV6xgAgIKCAq4jtCpU+AhpYZ49ewahUEiF7xMkqfBt3boVBgYGWLx4MR4+fMh1nBaPCh8hLcyTJ08AAF27dmWtzYqKChw7dgwDBw6Erq4uFBUV0bNnTwQHB6OqqqrW8SkpKfDw8ICamhqUlZUxYMAAXLly5aPtN/Z4NggEArx69arW4t9cef78OdatW4cePXqge/fuCAoKQkZGBtexWiQqfIS0MG/fvgWPx4OmpiZrbZ47dw4+Pj5wdnZGcnIynj9/Dj8/P8yePRvz58//4NgnT57A2toaN27cQFhYGLKzs7FlyxasXLkSf/31V622G3s8WzQ1NcEwjMTcZpSVla3ZgzElJQVLliyBgYEBevXqheDg4EZtPEw+QazLZRNCRG7Pnj2MsrIyq21GRkYyjo6OtV4fP348IysryxQUFNS8Nnr0aAYAExYW9sGxL1++ZOTl5Rl5efkPXm/s8Wy5c+cOA4BJTU0VSfuN8eOPPzLy8vJ17rzB4/EYaWlpRkpKiunfvz+zfft2pqCggHZn+IT6dmeQuHl8hJDPU1hYCBUVFVbbdHd3h7u7e63Xe/XqhYMHD+Lhw4ewtrYG8M/VIQC4ubl9cGy7du0gEAhqrUTT2OPZUv0ZXbx4Ec+ePRNJHw1VfXu6LgzD1CxZd+3aNVy7dg3+/v4YMmQIANRcJZKGo8JHSAtTWlrK+opGBQUFWL9+PSIiIvDixQvk5+d/8H71c7LS0lIUFhZCQUEBfD6/Vjva2tofFLLGHs+m6qXcpk+fLpL2G6shcy6rn6eWlZXh5MmTAICffvoJdnZ2aN++vUjztSRU+AhpYZSUlFBcXMxqm8OGDUN8fDyCg4MxZswYtG3bFjweD5s2bUJgYGDNXpvy8vJQUVGpWVf3v8UsLy/vgz839ng2Va8LfPnyZZiZmYmsn4bYtGlTzVZS9ZGTk0NZWRm6du2KMWPGYOXKlVi0aBEVvUaiwkdIC8Pn82v+UWdDZWUlEhISoKuri5kzZ37wnlAorHX8kCFDcPz4cZw7dw6enp41r79586bO6QONPZ4t1Yvd6+vrQ0NDQ2T9NER9C4lXD3rR1dWFl5cXJk2aBHNzc1RUVGDlypViTNly0KhOQlqY6t1MGrJ5bUNIS0vD0dERWVlZ+Omnn/DmzRsIhULExMRg27ZttY5fs2YNNDU1MWvWLFy4cAFFRUV49OgRxo8fX+ftzMYez5bqwsf281A2yMnJAQDU1NTw1VdfIT4+Hq9evUJwcDDMzc05Ttf8UeEjpIWp3o0hOzubtTaPHTuGb775Bps3b0a7du1gaGiIkJAQjB07FsA/ux1YWloCALp06YKkpCRYWVnB09MT2tramDRpEvz9/dGzZ0+UlpaCx+Ph66+/btLxbMnOzoa0tDTU1dVZbbepysrKwOPxoKysjHHjxiE6Ohq5ubnYvn077OzswOPxuI7YYtCtTkJamG7dugEAHj9+zNqzn7Zt29Z5dQcAP/74Y63XBAIBIiIiar0+dOjQOtto7PFsSE1NhYGBQc3VFZeUlJTg5eWFMWPGYPDgwRKzuHhLRVd8hLQwurq6UFdXl5jluCRVWlqaxCzrNnPmTBw9ehQjRoygoicGVPgIaYEEAgFSUlK4jiHRUlNTJabwEfGiwkdIC2RpaYnExESuY0is9+/f4+7duzXPJUnrQs/4CJFg169fx99//w11dXXIy8tDWVkZKioqkJeXh6qqKhQVFeucrO7k5ITt27cjPz9fYgZvSJK4uDiUl5fDycmJ6yiEA1T4CJFgd+7cgRofnJsAABd8SURBVJ+f3yeP4/P5kJOTg6KiIlRUVCAjI4OqqipMnDgRR48erXeeWGt06dIlmJiYQE9Pj+sohANU+Fqhx48f4++//+Y6hkRTU1ND3759uY4Bd3d38Hi8mpVRPqauCesaGhpo27YtFb06nD17Fi4uLlzHIByhwtcK7dixAz///DPXMSSapaUlrl+/znUM6OnpoXfv3rhz584nix/wzyofPB4P69atg1AoxE8//YQtW7bQSMF/uXPnDh48eIAdO3ZwHYVwRGIKn1AoxMOHD5GamoqUlBSkp6ejuLgYRUVFyM/Ph6KiYs3zDW1tbRgbG0MgEMDExAT6+vpcx292zM3NcfHiRa5jSKSlS5ciKSmJ6xg1Ro0ahQcPHnxyFX4ZGRkYGhoiLCwMPXv2xMuXL7F48WJERUVh5MiRYkor+Q4cOIAuXbqgf//+XEchHOG08F29ehXnzp1DTEwM/vzzT5SWlkJOTg5du3ZF165doaysDG1tbaipqUEoFNYUwjt37uDYsWN48+YNAMDQ0BBOTk5wdnaGu7s71NTUuDytZkFaWprz9QkllaRcHVVUVCAuLg6pqan1Fj0pKSlUVVXBx8cH27dvh5KSEgCgffv2cHJywr59+6jw/X9lZWU4fPgwpk2bRiuhtGJiL3wZGRnYv38/Dhw4gCdPnsDAwABOTk74+uuvYW1tDQMDA8jINCxWXl4e7t27h8uXL+PSpUs4ePAgpKSk4OHhAV9fXwwePBhSUjRjgzQf+fn5OHfuHE6dOoWzZ88iPz8fpqam0NLSQk5OTq3jZWVlIScnhz179sDLy6vW+99++y1Gjx6N5ORkdO/eXRynINEOHDiAvLw8TJ48mesohEvi2g33wYMHjK+vLyMjI8Noamoyfn5+THx8PKt95OfnM/v372dcXV0ZHo/HdO7cmdm0aRNTUlLCaj/N3Zw5cxhLS0uuY0gscX8+6enpzPbt2xl3d3dGTk6OkZaWZmxtbZm1a9cyycnJDMMwjL+/PyMnJ/fBztzS0tKMhYUF8/Tp04+2XVVVxfTo0YPx9fUV1+lIrIqKCqZbt27M1KlTuY7CCtqBvX717cAu8sL3/Plzxtvbm5GSkmJMTU2ZgwcPMuXl5aLulklJSWEmTZrEyMrKMgYGBszx48dF3mdzQYWvfqL+fCorK5kbN24wS5cuZSwsLBgAjLKyMuPu7s7s37+fefv2ba2f+eOPPz4oeDwej5k5cyZTVlb2yf4OHDjAyMjIME+ePBHF6TQbBw8eZGRkZJjHjx9zHYUVVPjqV1/hE9l9wPLycqxfvx7du3fHzZs3ERoainv37mHcuHENvpX5OYyMjLB3716kpaXB0dER3t7eGDx4MB4/fizyvgn5L6FQiOjoaAQEBKBDhw6wtLTE/v37YWFhgVOnTiEvLw+RkZGYMGFCnRPOHRwcoKSkBB6Ph7Zt2yImJgbBwcGQlZX9ZN8+Pj7o2rUr5syZI4pTaxaKi4uxaNEijB8/Hl27duU6DuGYSArfs2fP4OjoiEWLFmHGjBm4f/8+Ro0axcnzNgMDA+zduxfXrl1Dbm4uzMzMEBwcLPYcpPXJyclBSEgIvLy8oK2tDTc3NyQkJGDq1Km4ceMG0tPTsX37dgwbNuyTOwTIycnhiy++wJAhQ/Dw4UM4ODg0OIeMjAw2b96MEydO4PTp0597Ws3SypUrkZ+fjzVr1nAdhUgA1i+9IiIiMGXKFOjr6+Pu3bswNjZmu4smsbS0RFJSEv73v/9h9uzZuH79OrZt2ybSjS5J6/Pw4UOcPn0akZGRSEpKgry8PGxtbbFq1SqMHj0a7dq1a3LbmzZtQrt27Zo0GtHV1RXe3t4ICAiAi4tLq5rU/ujRI2zcuBEbNmyglVoIAJav+H7++Wd8+eWX8PT0xJ9//ikxRa+ajIwM1qxZg6ioKFy4cAEODg6sbtZJWp/KykpcuXIFCxYsgJGREXr06IH169ejc+fOOHr0KF6/fo0LFy4gICDgs4oe8M/0hM8Zgr9hwwa8ffsWAQEBn5WjOSkpKcG4cePQu3dvTJs2jes4REKwcsXHMAyWL1+OFStWYN26dRL/LMHNzQ1//vkn3NzcYG1tjfPnz9ds3knIp+Tl5eHixYuIjIzEqVOnUFBQABMTE4wcORLu7u6wtbWVyDli7dq1Q0hICIYPH44BAwbA19eX60giFxgYiPT0dNy6dQvS0tJcxyESgpXCN2vWLGzduhUHDhzAuHHj2GhS5AwMDBAfH48vvvgCjo6OSEhIgIGBAdexiIR6+vQpIiMjcfr0acTGxoJhGPTr1w8LFy7EqFGjms0vTu7u7pg5cyZmzJgBS0vLFj237/Dhw9i+fTuOHz+Ozp07cx1HZO7duwcVFRWuY0icqqqqj7/5uUNGV61axUhLSzOhoaGf2xQn8vPzmd69ezMCgYB5/fo113HEgqYz1K/68/nvlANNTU1m9OjRzP79+5n8/HyuYzZZaWkpY2Njw3Ts2JF5/vw513FEIiYmhpGXl2cCAwO5jiIy1cP16av+r7qmM3zWFd++ffuwZMkSbNmyBZ6enp/TFGfU1NRw9uxZ2NraYujQoYiLi6tzfzPSujx9+hSWlpYQCAQYPnw4NmzYAFtb2xZxu0xOTg5RUVFwcHCAm5sb4uPjoampyXUs1ty/fx8jR46Eh4dHi16MXUZGBnl5eVzHkHh1DWBscuG7f/8+pk+fjgULFjT7h8a6uro4d+4crKysEBgYiK1bt3IdiXBMT08PSUlJEAgEXEcRCTU1NZw+fRq2trYYPHgwoqKi0LZtW65jfba7d+9i8ODBsLCwQEhISItfspDW222aJv2tKC4uhre3NywtLbFixQq2M3GiW7du2LVrF7Zt24bDhw9zHYdwTFFRscUWvWr6+vqIjo5GTk4OBgwYgIyMDK4jfZbLly/DwcEBpqamiIiI+OTcSNJ6NanwzZ8/Hzk5OThy5IhYVmERF09PT8yYMQMzZsxAVlYW13EIEblu3bohISEBcnJysLW1xe3bt7mO1CRHjhzB4MGD4ebmhjNnztBgD1KvRhe+mzdvYtu2bVi/fj3at28vikyc+umnn6ChoYF58+ZxHYUQsWjXrh1iY2NhbGwMGxsbbNmyhetIDVZSUoLp06dj7Nix+Pbbb3HkyBGJ2VaKSK5GFb6qqip8++23sLW1bbFzgBQVFbFhwwYcPHgQsbGxXMchRCzU1dVx/vx5zJ8/HzNnzoSXl5fED5xITk6GtbU1Dh8+jNDQUGzcuLHFP9Mj7GjU35Lw8HDcvHkTv/32m0RO0GWLh4cH3NzcsGDBAq6jECI20tLSWLZsGf744w8kJCTAyMgIe/bsAcMwXEf7QHFxMRYuXIjevXtDVlYWt27darajygk3GlX4goKC4OnpiR49eogqj8RYtmwZrl69ikuXLnEdhRCxcnZ2RnJyMsaNG4dvvvkGdnZ2iI+P5zoWKisrcejQIZiYmNQ8bklKSmrRk9OJaDS48EVFReHWrVtYuHChKPNIjH79+sHV1RWrV6/mOgohYqeqqopNmzbh+vXrkJWVhb29PRwcHHDhwgWxZykvL8fevXvRvXt3TJw4EU5OTkhJScF3333XIuZVEvFrcOHbunUrBg4ciN69e4syj0SZN28eLl26hNTUVK6jEMKJ3r174/Lly4iNjYWCggIGDRqEHj16YN26dXj58qVI+7537x7mzJmDTp06Ydq0abC3t0dKSgr27dsHHR0dkfZNWrYGFb7Xr1/j/PnzmDRpkojjSBYXFxd06NABBw4c4DpKnXbv3o01a9YgMTER5eXlXMchLZi9vT3Onz+PGzduwN7eHkFBQejYsSOcnZ3x448/4s8//0RlZeVn9fH+/Xv88ccfWLBgAczMzNCrVy9ERETAz88PT548wa5du2gTWcKKBk3CO3ToEBQVFTFixAhR55EoUlJSGDt2LA4ePIgVK1ZI3Iixly9fYunSpQAABQUF2NjYwNXVFY6OjrCysmpRcyyJZLCwsICFhQU2bdqEM2fOIDw8HJs3b8aiRYugqqoKCwsLCASCmq82bdqAz+eDz+dDRUUF7969Q0FBAYqKipCdnY20tDSkpqYiJSUFd+7cQVlZGbp37w4XFxds2bJFYne6IM1bg/5lPHXqFEaMGAElJSVR55E4Pj4+CAoKwp07d9CnTx+u43xAR0cHMjIyqKioQElJCWJiYnDlyhWUlZVBXl4e1tbWcHBwgJ2dHezt7WklC8IaOTk5jBw5EiNHjgTwz9SCmJgY3L17F6mpqTh58uQnF4GQkpKCgYEBBAIBbGxs4O/vD2dn58/et5CQT/lk4SspKcHVq1fx1VdfiSOPxOnVqxe0tbURExMjkYWvoqKi5s8Mw6CsrAwAUFpaitjYWCQkJKC8vByKioqws7ODq6srMjMzJW6IOmneunfvXmuLo+LiYuTl5aGoqAjFxcV49+4d1NTUoKqqCj6fD01NTZpsTjjxycKXmJiIkpISODg4sNJhUVERzM3NYWRkhNOnT7PSpijxeDw4ODggJiYG33//PddxPvCpB/wMw9Q8+xMKhYiOjq4ZlScnJ4eYmBg4OTmJPCdpnZSVlaGsrMx1DEJq+WThS0pKgqGhITp16sRKhwzDoKqqqv5NAiWMg4MDFi9eLNY+GYZBTk4OXr9+jezsbGRlZeH169fIyspCdnY2Xr9+3ahFhaWkpMAwDHR1dWFsbIyCggIqeoSQVumThS85ORkmJiasdaiiooK//vqLtfbEwcTEBPn5+cjOzv6sYdQMw+D169fIyclBVlYWsrKykJOTg8zMTGRnZyMnJwevXr2qOebftzFlZWWhra0NHR0d6OrqQkdHByYmJnj06FG9fVY/AzQyMsKCBQswZswYLFq0CJcvX27yeRBCSHP2ycKXmpoKe3t7cWSRWEZGRgD++SzqKnxv377Fq1ev8PbtW2RmZn70++fPn9eadqChoQE9PT20a9cOenp66NatW833GhoaNd/r6OjUOVl3y5YtEAqFtV6vLngWFhb44Ycf4O7uTqPjCCEEDSh8T548wZQpU1jp7MSJEzWjwIB/njspKCjUej09PR3z58/H+fPnIS0tDWtrawQHB6NLly4ftJebm4vVq1fj5MmTePHiBbS0tGBsbIzx48fD29sbioqKrOTW09ODkpISVq1aBV1dXWRnZyMzM7PmVuS/b9vKy8tDS0urplhpaWnBzMwM2tra0NbWRrt27aClpQUdHR20adPms7O1adMGL168qPmzrKwsqqqqMGLECMyfPx9WVlaf3QchhLQk9Ra+qqoqFBQUQEtLi5XOPDw8wDAMPDw8cPLkyY++PmvWLMyfPx979uxBUlIShg8fjjFjxuDatWs1P5OVlQUbGxsIhULs3LkTDg4OeP/+PXbu3ImvvvoK+fn5mDVrFiu5eTwelJWVkZKSAnl5eejr66Nv3741BUxPTw9aWlrQ1dUV+47Ienp6ePHiBaSlpSEvL4/p06cjICAA+vr6Ys1BCCHNRb2Fr6ioCAzDgM/niysPAODrr7+GtbU1AMDV1RVDhw5FWFgY3rx5g7Zt2wIAFi5ciPT0dBw7dgzu7u4A/nl+uHjxYiQkJLCeSUdHByNHjpS4Hef19fXx7NkzzJkzB35+flBVVeU6EiGESLR6C19hYSEAiH034//enuvQoQMA4NWrVzWFLyIiAgAwZMiQWj9/9uxZ1jOpqqrWfB6SZMWKFTAyMoKsrCzXUQghpFmot/BVT3IW96AINTW1D/5cveJI9bO00tJSFBQUQEFBQWxFWUpK6rPXIhSF1rBFFCGEsKnexSerb3EWFRWJJUxDycvLQ01NDSUlJWK7Cnv37h3dRiSEkBag3sJXfTUlibf4qkeBRkVF1XrP3NwcgYGBrPZXWFgo9lu+hBBC2Fdv4ZOWloaKigry8vLElafBfvzxRxgaGiIwMBBnzpxBYWEhXrx4genTpyMzM5P1wvf27Vuoq6uz2iYhhBDx++Q+O126dMHjx49Z6ezEiRPg8Xg1UxkUFRUxfvx4XL16tdbr1UuE8Xg8BAUFAfjnSq56BKeuri6uX78Ob29v+Pv7o02bNujbty/evn2L+Ph4dOzYkZXMwD/7Eebn59NeYIQQ0gJ8cgK7kZER0tLSWOmser5eXRr7OvDP5O2NGzdi48aNrOT7mOod2I2NjUXaDyGEENH75BWfsbExkpOTxZFFYqWkpIDP59M+YYQQ0gJ8svD17dsXaWlpyM7OFkceiRQfH49+/frRWpeEENICfLLw2dvbQ0ZGBrGxseLII5Fo3zpCCGk5Pln4+Hw+LC0tcfHiRXHkkTgpKSl48eIFnJ2duY5CCCGEBZ8sfADg7u6OiIiIWlvqtAahoaHQ0dGhXQ4IIaSFaFDhGz9+PHJzc0WyBqakO3ToEMaMGQMZmU8OgCWEENIMNOhf844dO8LBwQH79+/H8OHDRZ1JYiQmJiI1NRWHDx/mOgrr3r17hwsXLnAdQyJlZGRwHYEQIkINvozx8/PD+PHj8eTJk1YzkXv9+vWwtLREnz59uI7CurS0NAwaNIjrGBLL0tKS6wiEEBHhMfXNEP+XyspKmJiYwN7eHjt37hR1Ls4lJyejR48eCAsL+2B3+JZAKBSipKSE6xgSTUZGhtZmJaRlqmxw4QOA3bt3Y/r06UhJSYGhoaEog3HOx8cHDx8+xL1792j+HiGEtByNK3xlZWXo1asXBAJBzbqaLVFsbCycnJwQERGBESNGcB2HEEIIexpX+ID/KwonT57EsGHDRBWMMxUVFbCwsEC7du1a5ShWQghp4Rpf+ABg7NixSExMxO3bt6GhoSGKYJxZvnw5goKC8ODBA3Tu3JnrOIQQQtjVtMKXl5cHc3NzmJmZ4dSpUy3mGVhsbCxcXFwQHByMGTNmcB2HEEII+5pW+AAgLi4OLi4uWLduHeubvnIhKysLffr0gY2NDcLCwriOQwghRDSaXvgAYO3atVi8eDFCQ0Ob9ZD/wsJCODo6orCwENeuXaOd1gkhpOWq/Kx1uBYsWIDMzEz4+PjgzJkzcHV1ZSuY2JSVlcHT0xMvXrzAlStXqOgRQkgL91lXfMA/E9t9fHxw/vx5nDhxolntYlBcXAwvLy8kJCTg8uXL6N27N9eRCCGEiFZlgxapro+0tDQOHToEd3d3DBkyBMePH2cjmMjl5eXBzc0NV69exdmzZ6noEUJIK8HKlgNycnI4ePAgtLS0MHbsWKSlpWHRokWQkvrsuioSd+7cwejRo1FRUYGkpCQIBAKuIxFCCBET1iqTlJQUgoODsWnTJqxatQpubm7Izs5mq3nWbN26FdbW1tDX10diYiIVPUIIaWVYvyT77rvvkJCQgPT0dPTs2RN79+7FZz5GZMWTJ08wZMgQfPfdd5g3bx6io6Ohp6fHdSxCCCFiJpJ7kRYWFrh16xbGjBmDqVOnYsCAAbh+/boouvqkd+/e4X//+x969uyJly9fIi4uDsuXL4e0tDQneQghhHBLZA/hVFVVERwcjBs3boBhGPTt2xdDhgxBfHy8qLr8QF5eHpYuXQoDAwP88ssvWL16NW7dugVbW1ux9E8IIUQyffZ0hob6448/sHr1asTFxcHc3BwTJkzAmDFjoKOjw1oflZWViI6OxsGDBxEREQEFBQUEBATA39+f5ucRQggBPnfllqZITEzEnj17EBYWhuLiYlhbW8PZ2RnOzs7o27cvFBQUGtXe06dPcfnyZVy6dAnR0dHIzs6GjY0NfH19MX78ePD5fBGdCSGEkGZI/IWvmlAoRGRkJM6fP49Lly7h77//hpSUFAwMDCAQCNC5c2eoqKhAQ0MDfD4fQqEQRUVFKCoqwqtXr5CWloa0tDQUFhZCSUkJtra2cHZ2hqenJ7p27crFKRFCCJF83BW+/0pPT8etW7eQlpaG1NRUpKeno6ioCPn5+SgsLISioiL4fD5UVFSgo6MDgUAAgUAAExMTWFpaQl5enutTIIQQIvkkp/ARQgghYvD5S5YRQgghzQkVPkIIIa0KFT5CCCGtyv8Dz3TQrkuzKYIAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "z.visualize(rankdir='LR')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Run in parallel\n",
    "\n",
    "Call `.compute()` when you want your result as a normal Python object\n",
    "\n",
    "If you started `Client()` above then you may want to watch the status page during computation."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "3"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "z.compute()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Parallelize Normal Python code\n",
    "\n",
    "Now we use Dask in normal for-loopy Python code.  This generates graphs instead of doing computations directly, but still looks like the code we had before.  Dask is a convenient way to add parallelism to existing workflows."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "zs = []"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "CPU times: user 36.8 ms, sys: 16.3 ms, total: 53.1 ms\n",
      "Wall time: 37.6 ms\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "for i in range(256):\n",
    "    x = inc(i)\n",
    "    y = dec(x)\n",
    "    z = add(x, y)\n",
    "    zs.append(z)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "zs = dask.persist(*zs)  # trigger computation in the background"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To make this go faster, add additional workers.\n",
    "\n",
    "(although we're still only working on our local machine, this is more practical when using an actual cluster)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "client.cluster.scale(10)  # ask for ten 4-thread workers"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "By looking at the Dask dashboard we can see that Dask spreads this work around our cluster, managing load balancing, dependencies, etc.."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Custom computation: Tree summation\n",
    "\n",
    "As an example of a non-trivial algorithm, consider the classic tree reduction.  We accomplish this with a nested for loop and a bit of normal Python logic.\n",
    "\n",
    "```\n",
    "finish           total             single output\n",
    "    ^          /        \\\n",
    "    |        c1          c2        neighbors merge\n",
    "    |       /  \\        /  \\\n",
    "    |     b1    b2    b3    b4     neighbors merge\n",
    "    ^    / \\   / \\   / \\   / \\\n",
    "start   a1 a2 a3 a4 a5 a6 a7 a8    many inputs\n",
    "```"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "([65536],)"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "L = zs\n",
    "while len(L) > 1:\n",
    "    new_L = []\n",
    "    for i in range(0, len(L), 2):\n",
    "        lazy = add(L[i], L[i + 1])  # add neighbors\n",
    "        new_L.append(lazy)\n",
    "    L = new_L                       # swap old list for new\n",
    "\n",
    "dask.compute(L)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "If you're watching the [dashboard's status page](../proxy/8787/status) then you may want to note two things:\n",
    "\n",
    "1.  The red bars are for inter-worker communication.  They happen as different workers need to combine their intermediate values\n",
    "2.  There is lots of parallelism at the beginning but less towards the end as we reach the top of the tree where there is less work to do.\n",
    "\n",
    "Alternativley you may want to navigate to the [dashboard's graph page](../proxy/8787/graph) and then run the cell above again.  You will be able to see the task graph evolve during the computation."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Further Reading\n",
    "\n",
    "For a more in-depth introdution to delayed and lazy operation in Dask, see the [dask tutorial](https://github.com/dask/dask-tutorial), notebooks 01 and 01x."
   ]
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "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.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}