Video+3+-+Autograd.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Autograd\n",
    "\n",
    "PyTorch's *Autograd* feature is part of what make PyTorch flexible and fast for building machine learning projects. It allows for the rapid and easy computation of multiple partial derivatives (also referred to as *gradients)* over a complex computation. This operation is central to backpropagation-based neural network learning.\n",
    "\n",
    "The power of autograd comes from the fact that it traces your computation dynamically *at runtime,* meaning that if your model has decision branches, or loops whose lengths are not known until runtime, the computation will still be traced correctly, and you'll get correct gradients to drive learning. This, combined with the fact that your models are built in Python, offers far more flexibility than frameworks that rely on static analysis of a more rigidly-structured model for computing gradients.\n",
    "\n",
    "## What Do We Need Autograd For?\n",
    "\n",
    "A machine learning model is a *function*, with inputs and outputs. For this discussion, we'll treat the inputs a as an *i*-dimensional vector $\\vec x$, with elements $x_{i}$. We can then express the model, *M*, as a vector-valued function of the input: *$\\vec y$ = $\\vec M$($\\vec x$)*. (We treat the value of M's output as a vector because in general, a model may have any number of outputs.)\n",
    "\n",
    "Since we'll mostly be discussing autograd in the context of training, our output of interest will be the model's loss. The *loss function* *L($\\vec y$) = L($\\vec M$($\\vec x$))* is a single-valued scalar function of the model's output. This function expresses how far off our model's prediction was from a particular input's *ideal* output. *Note: After this point, we will often omit the vector sign where it should be contextually clear - e.g., $y$ instead of $\\vec y$.*\n",
    "\n",
    "In training a model, we want to minimize the loss. In the idealized case of a perfect model, that means adjusting its learning weights - that is, the adjustable parameters of the function - such that loss is zero for all inputs. In the real world, it means an iterative process of nudging the learning weights until we see that we get a tolerable loss for a wide variety of inputs.\n",
    "\n",
    "How do we decide how far and in which direction to nudge the weights? We want to *minimize* the loss, which means making its first derivative with respect to the input equal to 0: $\\frac{\\partial L}{\\partial x} = 0$.\n",
    "\n",
    "Recall, though, that the loss is not *directly* derived from the input, but a function of the model's output (which is a function of the input directly), $\\frac{\\partial L}{\\partial x}$ = $\\frac{\\partial {L({\\vec y})}}{\\partial x}$. By the chain rule of differential calculus, we have $\\frac{\\partial {L({\\vec y})}}{\\partial x}$ = $\\frac{\\partial L}{\\partial y}\\frac{\\partial y}{\\partial x}$ = $\\frac{\\partial L}{\\partial y}\\frac{\\partial M(x)}{\\partial x}$.\n",
    "\n",
    "$\\frac{\\partial M(x)}{\\partial x}$ is where things get complex. The partial derivatives of the model's outputs with respect to its inputs, if we were to expand the expression using the chain rule again, would involve many local partial derivatives over every multiplied learning weight, every activation function, and every other mathematical transformation in the model. The full expression for each such partial derivative is the sum of the products of the local gradient of *every possible path* through the computation graph that ends with the variable whose gradient we are trying to measure.\n",
    "\n",
    "In particular, the gradients over the learning weights are of interest to us - they tell us *what direction to change each weight* to get the loss function closer to zero.\n",
    "\n",
    "Since the number of such local derivatives (each corresponding to a separate path through the model's computation graph) will tend to go up exponentially with the depth of a neural network, so does the complexity in computing them. This is where autograd comes in: It tracks the history of every computation. Every computed tensor in your PyTorch model carries a history of its input tensors and the function used to create it. Combined with the fact that PyTorch functions meant to act on tensors each have a built-in implementation for computing their own derivatives, this greatly speeds the computation of the local derivatives needed for learning.\n",
    "\n",
    "## A Simple Example\n",
    "\n",
    "That was a lot of theory - but what does it look like to use autograd in practice?\n",
    "\n",
    "Let's start with a straightforward example. First, we'll do some imports to let us graph our results:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import torch\n",
    "\n",
    "import matplotlib.pyplot as plt\n",
    "import matplotlib.ticker as ticker\n",
    "import math"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next, we'll create an input tensor full of evenly spaced values on the interval $[0, 2{\\pi}]$, and specify `requires_grad=True`. (Like most functions that create tensors, `torch.linspace()` accepts an optional `requires_grad` option.) Setting this flag means that in every computation that follows, autograd will be accumulating the history of the computation in the output tensors of that computation."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "tensor([0.0000, 0.2618, 0.5236, 0.7854, 1.0472, 1.3090, 1.5708, 1.8326, 2.0944,\n",
      "        2.3562, 2.6180, 2.8798, 3.1416, 3.4034, 3.6652, 3.9270, 4.1888, 4.4506,\n",
      "        4.7124, 4.9742, 5.2360, 5.4978, 5.7596, 6.0214, 6.2832],\n",
      "       requires_grad=True)\n"
     ]
    }
   ],
   "source": [
    "a = torch.linspace(0., 2. * math.pi, steps=25, requires_grad=True)\n",
    "print(a)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next, we'll perform a computation, and plot its output in terms of its inputs:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x19bac017430>]"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjgAAAGdCAYAAAAfTAk2AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAABaSklEQVR4nO3de1iUZcIG8HtmgBkQGM4nOXoCPCuewNRMQy0rsxLLUMssd2tLzV1z2/222i3XdivbSkvTLCu1UsvKE5hn8ACKR0RUEEQOgjADCAzMvN8foxQhCMrwzOH+Xddc1+fLOy/3jPs1t8+8z/PIJEmSQERERGRF5KIDEBEREbU1FhwiIiKyOiw4REREZHVYcIiIiMjqsOAQERGR1WHBISIiIqvDgkNERERWhwWHiIiIrI6d6AAiGAwGXL58GS4uLpDJZKLjEBERUQtIkoTy8nIEBARALm9+jMYmC87ly5cRFBQkOgYRERHdhtzcXAQGBjZ7jk0WHBcXFwDGN8jV1VVwGiIiImoJrVaLoKCg+s/x5thkwbnxtZSrqysLDhERkYVpye0lvMmYiIiIrA4LDhEREVkdFhwiIiKyOiw4REREZHVYcIiIiMjqsOAQERGR1WHBISIiIqvDgkNERERWhwWHiIiIrI5JC86ePXvwwAMPICAgADKZDN9///0tn7N7925ERUVBpVKhU6dO+Pjjjxuds379enTv3h1KpRLdu3fHxo0bTZCeiIiILJVJC05lZSX69OmDDz/8sEXnZ2Vl4b777sOwYcNw9OhR/PWvf8WLL76I9evX15+TnJyMuLg4xMfH49ixY4iPj8ekSZNw8OBBU70MIiIisjAySZKkdvlFMhk2btyICRMmNHnO/PnzsWnTJqSnp9cfmzVrFo4dO4bk5GQAQFxcHLRaLbZs2VJ/ztixY+Hu7o41a9a0KItWq4VarYZGo+FeVERERBaiNZ/fZnUPTnJyMmJjYxscGzNmDFJSUlBbW9vsOUlJSU1et6amBlqttsGD6Pdq6vTYc/YK3k04i03HLkNTVSs6EhER3Saz2k28oKAAvr6+DY75+vqirq4OxcXF8Pf3b/KcgoKCJq+7cOFCvP766ybJTJat7JoOOzOKkHi6CLvPXkFFTV39z+zkMgzu5IHRkb4YHemLIA8ngUmJiKg1zKrgAI23QL/xDdpvj9/snOa2Tl+wYAHmzp1b/2etVougoKC2iEsWKLu4EonphUg4XYiUi6XQG379ltbbRYnoTp5Iz9cis6gC+8+VYP+5Erz+42lE+LkYy053X/TuqIZc3vT/5oiISCyzKjh+fn6NRmKKiopgZ2cHT0/PZs/5/ajObymVSiiVyrYPTBZBb5CQlluKhNNFSEwvxLmiigY/b6q43ChCiemFOJxdijMF5ThTUI4Pd56Dt4sSoyN9MDrSF0O7eEFlrxDx0oiIqAlmVXCio6Px448/Nji2fft2DBgwAPb29vXnJCQkYM6cOQ3OiYmJadesZN6u6eqwN7MYiacL8cuZIpRU6up/1tKvnkK9OuCZYZ3wzLBOKLumw66MK0hIL8TujCu4Ul6DNYdyseZQLlT2cgzr6o17I30xMsIH3i4s00REopm04FRUVODcuXP1f87KykJaWho8PDwQHByMBQsWIC8vD1988QUA44ypDz/8EHPnzsXMmTORnJyMFStWNJgd9dJLL2H48OFYtGgRHnroIfzwww9ITEzEvn37TPlSyAIUaquxI904SrPvXDF0dYb6n7mo7DAy3Aeju/tiRDdvqB3tW3VtNycHTOjXERP6dYSuzoCDWSVIPF2IxPQi5JVVIeG08SsvmQzoF+SG0d19cW+kL7r4ODf79SkREZmGSaeJ79q1CyNHjmx0fNq0aVi1ahWmT5+O7Oxs7Nq1q/5nu3fvxpw5c3Dq1CkEBARg/vz5mDVrVoPnf/fdd/jb3/6GCxcuoHPnznjzzTcxceLEFufiNHHrkldWhbnr0nAw62qD44Hujrj3etEYGOYBe0XbTxqUJAnp+eX1X2Udv6Rp8PNOXh3w1sReGNLJs81/NxGRrWnN53e7rYNjTlhwrEd6vhbTPzuEQm0NAKBvkBvu7W786qmbb/uPnhRoqrHjTCESTxdi//kS6OoMcFDI8W5cH4zvHdCuWYiIrA0Lzi2w4FiH5PMlePaLFJTX1KGrjzNWTBuIYE/zmcpdUVOHP397DFtOFkAmA/5+f3c8fVeY6FhERBbLYhf6I2qpH49dxrSVh1BeU4dBoR74blaMWZUbAHBW2uHDJ/pjanQIJAl446fTWLg5HQaDzf2bgoio3bHgkMVZsS8Lf1pzFDq9AeN6+uGLGYOgdmrdTcPtRSGX4fUHe+AvY8MBAJ/suYC536Q1uAGaiIjaHgsOWQyDQcJbm9Pxz59OAwCmRofgwyf6m/0aNDKZDH+8uwveeawP7OQyfJ92GU+vOozyam4FQURkKiw4ZBF0dQbM+SYNy/ZcAAD8ZWw4Xn+wBxQWtJrwI1GBWDF9IJwcFNh3rhhxnxxAkbZadCwiIqvEgkNmr7y6Fk+tOoQf0i7DTi7DO4/1wR/v7mKR68uM6OaNdc9Gw8vZAafztZi4NAnnr1Tc+olERNQqLDhk1oq01Zj0yQHsP1cCJwcFVkwfiEeiAkXHuiO9AtVY/4cYhHo64VJpFR5dmoQjOaWiYxERWRUWHDJb569U4OElSUjP18LL2QHrno3GiG7eomO1iRDPDvjuDzHoE6hG6bVaPLH8ABJOF4qORURkNVhwyCylXizFI0uTkFdWhVBPJ2z4w1D0ClSLjtWmvJyVWPPsEIwM90Z1rQHPrU7BmkM5omMREVkFFhwyOwmnCzHl0wMou1aLPte/zjG3NW7aipODHZZPHYBJAwJhkIAFG07gvYSzsMH1N4mI2hQLDpmVrw/m4LnVKaiuNWBkuDfWPDsEns7WvTu3nUKORY/0xov3dAEAvL8jE6+sP4E6PdfKISK6XSw4ZBYkScK7CWfx140nYJCASQMCsXzqADg5mHTDe7Mhk8kwNzYcbz7cE3IZsC4lF8+uTsU1XZ3oaEREFokFh4Sr0xvwyvoT+N+OTADAi/d0waJHesPOBLt/m7spg0Pw8ZNRUNrJ8cuZIjyx/CBKKmpExyIisji29wlCZuWarg7Prk7FupRcyGXAmw/3xNzYcItc46atxPbww9czB8PNyR5puWV49ONk5JRcEx2LiMiisOCQMCUVNXh8+UH8cqYISjs5Pn4yClMGh4iOZRaiQjzw3axodHRzRFZxJSYuTcLJPI3oWEREFoMFh4S4pqvD5GUHcCy3DG5O9vh65mDE9vATHcusdPFxwYY/xiDS3xXFFTWI+yQZmYXlomMREVkEFhwS4l8/pyOzqAI+Lkp8NysGUSEeoiOZJV9XFdY9NwSDQj1QqdPjpbXciZyIqCVYcKjdJZ4uxNcHjQvavRfXF118nAUnMm+uKnt8+EQ/eHQw7l/1bsJZ0ZGIiMweCw61qyvlNZi//jgAYOawMAzt4iU4kWXwcVVh4cReAIBP9pzHgQslghMREZk3FhxqN5IkYf764yip1CHCzwXzxoSLjmRRxvTww+SBQZAkYO66NGiqakVHIiIyWyw41G6+PJiDX84UwcFOjsWT+0JppxAdyeL8fXx3hHg64bKmGv/3w0nRcYiIzBYLDrWLc0UVePPn0wCA+WMjEOHnKjiRZeqgtMN7cX2hkMvwQ9pl/JCWJzoSEZFZYsEhk9PVGTB73VFU1xpwVxcvPBUTKjqSResf7I4/Xd+36m/fn0ReWZXgRERE5ocFh0zu/R1ncTJPCzcne/z3sT6Qy213leK28sLILugX7Iby6jrMXZcGvYG7jxMR/RYLDpnUoayrWLLrPADgrYd7wU+tEpzIOtgp5Fgc1xdODgoczLqKT/deEB2JiMissOCQyWirazFnXRokCXg0KhD39fIXHcmqhHh2wGsP9AAA/Hd7BrdyICL6DRYcMpnXNp1CXlkVgjwc8Y8HuouOY5UeGxCIMT18UauXMHtdGqpr9aIjERGZBRYcMomfjl/GhiN5kMuA9yb1hYvKXnQkqySTybBwYm94uyhxrqgC/95yRnQkIiKzwIJDbS5fU4VXNxrXaHl+ZBcMCOU+U6bk0cEB/32sDwBgVVI2dmUUCU5ERCQeCw61KYNBwrxvj0FTVYvegWq8OKqr6Eg2YUQ3b0y/Pv3+z98dx9VKndhARESCseBQm1q5Pwv7z5XA0V6BxXF9Ya/g/8TayyvjItDVxxlXymvwyvrjkCROHSci29Uunz5LlixBWFgYVCoVoqKisHfv3ibPnT59OmQyWaNHjx496s9ZtWrVTc+prq5uj5dDTUjP1+LtrRkAjFsKdPLmLuHtSWWvwOLJfWGvkGH76UJ8m3JJdCQiImFMXnDWrVuH2bNn49VXX8XRo0cxbNgwjBs3Djk5OTc9//3330d+fn79Izc3Fx4eHnjssccanOfq6trgvPz8fKhUXGNFlOpaPWavTYNOb8DoSB88PihIdCSb1CNAjXmxxk1MX/vxFLKLKwUnIiISw+QF591338WMGTPwzDPPIDIyEosXL0ZQUBCWLl160/PVajX8/PzqHykpKSgtLcVTTz3V4DyZTNbgPD8/P1O/FGrGf7ZlIKOwHF7ODvj3I70hk3G1YlGeGdYJQzp54JpOj9nr0lCnN4iORETU7kxacHQ6HVJTUxEbG9vgeGxsLJKSklp0jRUrVmD06NEICQlpcLyiogIhISEIDAzE+PHjcfTo0SavUVNTA61W2+BBbWdfZjFW7MsCALz9aG94OSsFJ7JtCrkM70zqCxeVHdJyy/DhznOiIxERtTuTFpzi4mLo9Xr4+vo2OO7r64uCgoJbPj8/Px9btmzBM8880+B4REQEVq1ahU2bNmHNmjVQqVQYOnQoMjMzb3qdhQsXQq1W1z+Cgvj1SVspu6bDy9+mAQCeHBKMeyJ8m38CtYuObo7414SeAIAPfjmHIzmlghMREbWvdrnJ+PdfV0iS1KKvMFatWgU3NzdMmDChwfEhQ4bgySefRJ8+fTBs2DB888036NatGz744IObXmfBggXQaDT1j9zc3Nt+LfQrSZLw140nUKitQSfvDnj1Pq5WbE4e6tsRD/UNgN4gYc66NFTW1ImORETUbkxacLy8vKBQKBqN1hQVFTUa1fk9SZKwcuVKxMfHw8HBodlz5XI5Bg4c2OQIjlKphKura4MH3bn1R/Kw+UQB7OQyvB/XD44OCtGR6HfeeKgnOro54mLJNfzzp9Oi4xARtRuTFhwHBwdERUUhISGhwfGEhATExMQ0+9zdu3fj3LlzmDFjxi1/jyRJSEtLg78/N3NsLzkl1/CPH4yrFc+5txt6BaoFJ6KbUTva451JfSCTAWsP52LryVt/NUxEZA1M/hXV3Llz8emnn2LlypVIT0/HnDlzkJOTg1mzZgEwfn00derURs9bsWIFBg8ejJ49ezb62euvv45t27bhwoULSEtLw4wZM5CWllZ/TTKtOr0Bc79JQ6VOj4Gh7pg1orPoSNSMIZ088dxw49/Rgg3HUaTlelFEZP3sTP0L4uLiUFJSgjfeeAP5+fno2bMnNm/eXD8rKj8/v9GaOBqNBuvXr8f7779/02uWlZXh2WefRUFBAdRqNfr164c9e/Zg0KBBpn45BODj3eeRcrEUzko7vDupLxRyTgk3d3Pv7YY9Z6/gdL4W8747js+fGsip/ERk1WSSDa7nrtVqoVarodFoeD9OKx3LLcMjS5NQZ5Dw7qQ+mNg/UHQkaqHMwnKM/2AfauoMeP3BHph2fe8qIiJL0ZrPb24URC2mN0iYv/446gwS7u/tj4f7dRQdiVqhq68L/npfJABg4ZZ0FGj4VRURWS8WHGqxjUfzcKagHK4qO/zroZ78isMCTY0OwcBQd1TXGrA48azoOEREJsOCQy1SXavHu9uNG2n+cWQXuHdofuo+mSeZTIZXxkUAAL5JycW5onLBiYiITIMFh1rki+RsXNZUw1+twnTeu2HRokI8MKaHLwwSsOj67u9ERNaGBYduSXOtFh/tPA/AuOaNyp4L+lm6P4+JgFwGJJwuREr2VdFxiIjaHAsO3dLS3eehqapFN19nPMJZU1ahi48z4gYa92T795YzsMHJlERk5VhwqFn5mip8tt+4U/hfxkRwzRsr8tKoblDZy5FysRSJ6UWi4xARtSkWHGrW4oRM1NQZMCjUA6MifUTHoTbkp1bh6aFhAIC3t55Bnd4gOBERUdthwaEmZRaW49tU487r88dFcFq4FZp1d2e4Odkjs6gC649cEh2HiKjNsOBQkxZtzYBBAsb08EVUiLvoOGQCrip7vDCyCwDgvYRMVOn0ghMREbUNFhy6qcPZV5GYXgiFXIY/j4kQHYdMKD46BB3dHFGgrcaqpGzRcYiI2gQLDjUiSRL+veUMAGDSgCB08XEWnIhMSWmnwMux3QAAS3adQ9k1neBERER3jgWHGtl+uhCpF0uhspdj9uiuouNQO5jQtyMi/V1RXl2Hj3aeEx2HiOiOseBQA3V6A97eahy9mXFXGHxdVYITUXuQy2WYPzYcAPB50kVcKr0mOBER0Z1hwaEGvku9hPNXKuHuZI/nRnQWHYfa0Yhu3oju5Amd3oD3EjJFxyEiuiMsOFSvSqfHe9d3mH7hnq5wVdkLTkTt6bcbcW44eglnCrSCExER3T4WHKq3cn8WCrU16OjmiCeHBIuOQwL0CXLD/b39IUnAous3mhMRWSIWHAIAlFbq8PEu44aa88Z0g9KOG2raqnmx4bCTy7Az4wqSz5eIjkNEdFtYcAgA8NHOcyivqUOkvyse6tNRdBwSKMyrAx4fZBzB+/dWbsRJRJaJBYdwqfQavki+CAB4ZVwE5NxQ0+b9aVQXODkocCy3DFtPFoiOQ0TUaiw4hHe3n4VOb0BMZ08M7+olOg6ZAR8XFZ4Z1gkA8Pa2DNRyI04isjAsODbu9GUtNqblAQDmj+WGmvSrZ4d3gmcHB2QVV2Ld4VzRcYiIWoUFx8a9ve0MJAm4v7c/+gS5iY5DZsRZaYcXRxlXsl6cmInKmjrBiYiIWo4Fx4YlnS/GrowrsJPL8OfYcNFxyAw9PigYwR5OKK6owcp9WaLjEBG1GAuOjZIkqX6dkycGByPUq4PgRGSOHOzkmDfGWH4/2XMBJRU1ghMREbUMC46N2nyiAMcuaeDkoMCf7uGGmtS08b380aujGhU1dfjgF27ESUSWgQXHBtXqDfjPNuPozcxhneDtohSciMyZXP7rFg5fHbyInBJuxElE5o8FxwatPZyL7JJr8HJ2wMzhnUTHIQswtIsXhnX1Qq1ewjsJGaLjEBHdEguOjamsqcP7icadol8c1RXOSjvBichSzB9rHMX5Ie0yTuZpBKchImoeC46N+XRvFoorahDi6YTJA7mhJrVcz45qPNQ3AACwaCs34iQi88aCY0OKK2qwbI9xQ80/jwmHgx3/+ql15sWGw14hw97MYuzNvCI6DhFRk9rlE27JkiUICwuDSqVCVFQU9u7d2+S5u3btgkwma/Q4c6bhvxjXr1+P7t27Q6lUonv37ti4caOpX4bF+/CXc6jU6dE7UI37evqLjkMWKMjDCU8OCQFgHMUxGLgRJxGZJ5MXnHXr1mH27Nl49dVXcfToUQwbNgzjxo1DTk5Os8/LyMhAfn5+/aNr11+nMicnJyMuLg7x8fE4duwY4uPjMWnSJBw8eNDUL8diXSypxFcHr2+oOZYbatLte2FkFzgr7XAyT4ufTuSLjkNEdFMySZJM+k+wwYMHo3///li6dGn9scjISEyYMAELFy5sdP6uXbswcuRIlJaWws3N7abXjIuLg1arxZYtW+qPjR07Fu7u7lizZs0tM2m1WqjVamg0Gri6urb+RVmgP605ih+PXcbwbt744ulBouOQhftgRybeSTiLYA8nJM4dwa87iahdtObz26T/VdLpdEhNTUVsbGyD47GxsUhKSmr2uf369YO/vz9GjRqFnTt3NvhZcnJyo2uOGTOmyWvW1NRAq9U2eNiSE5c0+PHYZchkxtEbojs1Y1gYvF2UyLl6DV9fHxkkIjInJi04xcXF0Ov18PX1bXDc19cXBQUFN32Ov78/li1bhvXr12PDhg0IDw/HqFGjsGfPnvpzCgoKWnXNhQsXQq1W1z+CgoLu8JVZlvd3nAUAPNQnAN0DbGPEikzLycEOs0cbvzb+aNd5VNfqBSciImqoXRZBkcka3u8hSVKjYzeEh4cjPPzXjR+jo6ORm5uL//73vxg+fPhtXXPBggWYO3du/Z+1Wq3NlJzMwnIkphdBJkP9ztBEbeGxqCB8+Ms55GuqsfFoHh4fxGUHiMh8mHQEx8vLCwqFotHISlFRUaMRmOYMGTIEmZmZ9X/28/Nr1TWVSiVcXV0bPGzFJ3suAADGdPdDJ29nwWnImjjYyTHjrjAAwPI9F6DnjCoiMiMmLTgODg6IiopCQkJCg+MJCQmIiYlp8XWOHj0Kf/9fpzVHR0c3uub27dtbdU1bkK+pwg9peQCA50ZwSwZqe5MHBcNVZYcLxZVIOH3zr4iJiEQw+VdUc+fORXx8PAYMGIDo6GgsW7YMOTk5mDVrFgDj10d5eXn44osvAACLFy9GaGgoevToAZ1Ohy+//BLr16/H+vXr66/50ksvYfjw4Vi0aBEeeugh/PDDD0hMTMS+fftM/XIsyoq9WajVSxgc5oF+we6i45AVclbaYWp0KD7ceQ5Ld1/AmB5+TX5VTETUnkxecOLi4lBSUoI33ngD+fn56NmzJzZv3oyQEONiYfn5+Q3WxNHpdJg3bx7y8vLg6OiIHj164Oeff8Z9991Xf05MTAzWrl2Lv/3tb/j73/+Ozp07Y926dRg8eLCpX47F0FyrxZpDxvd11t2dBachazZ9aCiW772AY7llOJh1FUM6eYqORERk+nVwzJEtrIPz0c5z+M+2DET4uWDLS8P4r2oyqb99fwJfHsjB3eHeWPUU11kiItMwm3VwSIzqWj0+258FwHjvDcsNmdrMYZ0glwG7Mq4gPd+21pkiIvPEgmOF1h+5hOIKHTq6OWJ87wDRccgGhHh2wLhexokAn+w+LzgNERELjtXRGyQsvz41fMZdYbBX8K+Y2ses4cZ7vX48no9LpdcEpyEiW8dPPyuz9WQBskuuwc3JHpMH2cZihmQeegWqMbSLJ/QGCZ/uzRIdh4hsHAuOFZEkCR9f/3pganQonBzaZaFqonqzRhhHcdYdzkVppU5wGiKyZSw4ViT5fAlO5GmgspdjWnSI6Dhkg+7q4oUeAa6oqtXji2RuwklE4rDgWJGl10dvJg0IgqezUnAaskUymax+FOfz5GxU6bgJJxGJwYJjJU5d1mBvZjEUchlmDuO2DCTOuJ5+CPJwxNVKHb5JyRUdh4hsFAuOlfhkt3Hm1P29/BHk4SQ4DdkyO4Ucz14v2cv3XkCd3iA4ERHZIhYcK5B79Rp+On4ZAPDscI7ekHiPDQiCZwcHXCqtws8n8kXHISIbxIJjBZbvvQCDBAzr6oWeHdWi4xBBZa/AtJhQAMDHuy/ABneEISLBWHAsXElFTf19Dn8YwU01yXxMjQ6Bk4MC6fla7MksFh2HiGwMC46F+zz5IqprDejVUY3oztzFmcyHm5MDJg8MBsDtG4io/bHgWLBrujp8kZwNwLjAGjfVJHMzY1gY7OQyJJ0vwfFLZaLjEJENYcGxYOsO56LsWi1CPJ0wtqef6DhEjXR0c8SDfYwbvn7MURwiakcsOBaqVm+o3+9n5rBOUMg5ekPm6bnr94ZtOVmArOJKwWmIyFaw4Fion45fRl5ZFbycHfBoVKDoOERNCvdzwT0RPpAk44w/IqL2wIJjgSRJql/Y76mhYVDZKwQnImrec9fXZ/ou9RKKyqsFpyEiW8CCY4F2nb2CMwXl6OCgwJODuakmmb9BYR7oF+wGXZ0Bq/Zni45DRDaABccCfbzLeLPm44OCoXayF5yG6NZ+uwnn6gMXUVFTJzgREVk7FhwLczSnFAezrsJeIcOMYWGi4xC12L2Rvujk3QHl1XVYczBHdBwisnIsOBbmxr03D/XtCH+1o+A0RC0nl8vq78VZsS8LujpuwklEpsOCY0HOX6nAttMFAH69aZPIkkzo1xG+rkoUaKvxfVqe6DhEZMVYcCzI8j0XIEnA6EgfdPV1ER2HqNWUdgo8PdT41eqyPRdgMHATTiIyDRYcC1GkrcaGI8Z/8c7ipppkwR4fHAwXpR3OFVVgx5ki0XGIyEqx4FiIlfuzodMbEBXijgGhHqLjEN02V5U9pgwxLm/A7RuIyFRYcCyAtroWXx24CICjN2Qdnh4aCgeFHKkXS5GSfVV0HCKyQiw4FmDNwRyU19Shi48zRkX4iI5DdMd8XFWY2L8jAI7iEJFpsOCYuZo6PVbsM26q+ezwTpBzU02yEs8O7wSZDEhML8LZwnLRcYjIyrDgmLnvj+ahqLwGfq4qTOjbUXQcojbTydsZY7r7Afh1fSciorbCgmPGDAYJn+wx/od/xl1hcLDjXxdZl+dGGNdz+iEtD/maKsFpiMiatMsn5pIlSxAWFgaVSoWoqCjs3bu3yXM3bNiAe++9F97e3nB1dUV0dDS2bdvW4JxVq1ZBJpM1elRXW9cuxQnphbhwpRIuKjtMHhQkOg5Rm+sX7I7BYR6oM0hYsTdLdBwisiImLzjr1q3D7Nmz8eqrr+Lo0aMYNmwYxo0bh5ycm+9Fs2fPHtx7773YvHkzUlNTMXLkSDzwwAM4evRog/NcXV2Rn5/f4KFSqUz9ctqNJEn1N1/GDwmBi4qbapJ1mnW3cWbgmkM50FyrFZyGiKyFyQvOu+++ixkzZuCZZ55BZGQkFi9ejKCgICxduvSm5y9evBh/+ctfMHDgQHTt2hVvvfUWunbtih9//LHBeTKZDH5+fg0e1uRwdimO5pTBwU6O6UNDRcchMpm7u3kjws8FlTo9vjx4UXQcIrISJi04Op0OqampiI2NbXA8NjYWSUlJLbqGwWBAeXk5PDwaLm5XUVGBkJAQBAYGYvz48Y1GeH6rpqYGWq22wcPcrbw+c2piv47wcbGekSmi35PJZHj2+t5qXyRno1bPTTiJ6M6ZtOAUFxdDr9fD19e3wXFfX18UFBS06BrvvPMOKisrMWnSpPpjERERWLVqFTZt2oQ1a9ZApVJh6NChyMzMvOk1Fi5cCLVaXf8ICjLv+1kulV7D9uubaj51fd8eImt2f29/eDkrUaitwZaTLftvAxFRc9rlJmOZrOHaLZIkNTp2M2vWrMFrr72GdevWwcfn1wXuhgwZgieffBJ9+vTBsGHD8M0336Bbt2744IMPbnqdBQsWQKPR1D9yc3Pv7AWZ2OrkizBIwNAungj346aaZP2Udgo8OSQYAPDZft5sTER3zqQFx8vLCwqFotFoTVFRUaNRnd9bt24dZsyYgW+++QajR49u9ly5XI6BAwc2OYKjVCrh6ura4GGurunqsOaQ8Qbsp2I4ekO2Y8rgENgrZDiaU4a03DLRcYjIwpm04Dg4OCAqKgoJCQkNjickJCAmJqbJ561ZswbTp0/H119/jfvvv/+Wv0eSJKSlpcHf3/+OM4u28WgetNV1CPZwwkhuy0A2xNtFiQd6BwAAVnEUh4jukMm/opo7dy4+/fRTrFy5Eunp6ZgzZw5ycnIwa9YsAMavj6ZOnVp//po1azB16lS88847GDJkCAoKClBQUACNRlN/zuuvv45t27bhwoULSEtLw4wZM5CWllZ/TUslSRJW7c8GAEyLCYWC2zKQjblxz9nPJ/JRpLWuda2IqH2ZvODExcVh8eLFeOONN9C3b1/s2bMHmzdvRkhICAAgPz+/wZo4n3zyCerq6vD888/D39+//vHSSy/Vn1NWVoZnn30WkZGRiI2NRV5eHvbs2YNBgwaZ+uWY1P5zJcgsqkAHBwUeGxAoOg5Ru+sVqEZUiDtq9RK+PHjztbKIiFpCJkmSJDpEe9NqtVCr1dBoNGZ1P86MVYex40wRpkWH4PWHeoqOQyTET8cv44Wvj8LL2QH7X7kHSjuF6EhEZCZa8/nNzY3MRHZxJX7JKAJg/HqKyFaN6eEHf7UKxRU6/HQsX3QcIrJQLDhm4vPkbEgScHe4Nzp5O4uOQySMvUKOJ4cYv8L+LCkLNjjITERtgAXHDJRX1+LblEsAuLAfEQA8MSgYSjs5TuZpkXKxVHQcIrJALDhmYH3qJVTU1KGzdwcM7+olOg6RcO4dHPBwv44AUD+zkIioNVhwBDMYJHyebNxgcHpMaItWeCayBTc2md16qgCXy6rEhiEii8OCI9ius0XIKq6Ei8oOE/tzajjRDRF+roju5Am9QcIXydxlnIhahwVHsM+uD7/HDQhCB6Wd2DBEZubGKM7awzmo0unFhiEii8KCI9C5onLszSyGXMap4UQ3MzrSF0Eejii7Vovv0/JExyEiC8KCI9CN0Rvjf8SdxIYhMkMKuQzTokMBGG825pRxImopFhxBNNdqseGI8V+kN4bhiaixxwYEwclBgYzCciSfLxEdh4gsBAuOIOtSclBVq0eEnwuiO3mKjkNkttSO9njk+g34KzllnIhaiAVHAL1BwudJnBpO1FI37lHbcaYQOSXXxIYhIovAgiNAwulC5JVVwd3JHhOuL2ZGRE3r4uOM4d28IUnGbU2IiG6FBUeAz/ZnAQAeHxQMlT13SiZqiaeu36v2zeFcVNbUiQ1DRGaPBaednb6sxcGsq1DIZfUbChLRrY3o6o1OXh1QXlOH9UcuiY5DRGaOBaedrUoyjt6M7emHADdHwWmILIdcLqu/F2fV/mwYDJwyTkRNY8FpRyUVNfg+7TIA4GlODSdqtUeiAuGitMOF4krsybwiOg4RmTEWnHa09nAudHUG9OqoRv9gd9FxiCyOs9IOjw0IAvDrQplERDfDgtNOavUGrL6+YeBTQzk1nOh2TYsJgUwG7D57BeevVIiOQ0RmigWnnWw9WYACbTW8nJW4v7e/6DhEFivEswNGRfgAAD5PyhYbhojMFgtOO7kxNXzK4GAo7Tg1nOhOPDU0DADwXeolaKpqBachInPEgtMOjuWW4UhOGewVMkwZEiw6DpHFi+nsiW6+zrim0+PblFzRcYjIDLHgtINV14fRx/cOgI+LSmwYIisgk8kwPcY4ivN5cjb0nDJORL/DgmNiRdpq/HTcODV8+vU1PIjozj3cryPUjvbIvVqFX84UiY5DRGaGBcfEvjyYg1q9hP7BbugT5CY6DpHVcHRQYPKgG1PGswSnISJzw4JjQjV1enx98MbU8DDBaYisz9ToUCjkMiSdL8GZAq3oOERkRlhwTOinY/kortDBz1WFsT39RMchsjod3RwxpocvAE4ZJ6KGWHBMRJKk+puL46NDYK/gW01kCjduNt5wJA+llTrBaYjIXPBT10RSL5biRJ4GDnZyPD6IU8OJTGVgqDt6BLiips6ANYdzRMchIjPBgmMiN/bJmdA3AB4dHMSGIbJixinjoQCA1ckXUac3iA1ERGaBBccELpdVYeupAgC8uZioPTzQJwCeHRyQr6nGtlOFouMQkRlol4KzZMkShIWFQaVSISoqCnv37m32/N27dyMqKgoqlQqdOnXCxx9/3Oic9evXo3v37lAqlejevTs2btxoqvittvrARegNEoZ08kCkv6voOERWT2WvwJTBxq+CVyVxyjgRtUPBWbduHWbPno1XX30VR48exbBhwzBu3Djk5Nz8u/KsrCzcd999GDZsGI4ePYq//vWvePHFF7F+/fr6c5KTkxEXF4f4+HgcO3YM8fHxmDRpEg4ePGjql3NL1bV6rDlkfG03bn4kItObMiQEdnIZDmeX4mSeRnQcIhJMJkmSSdc4Hzx4MPr374+lS5fWH4uMjMSECROwcOHCRufPnz8fmzZtQnp6ev2xWbNm4dixY0hOTgYAxMXFQavVYsuWLfXnjB07Fu7u7lizZs0tM2m1WqjVamg0Gri6tu0Iy9pDOXhlwwkEujti959HQiGXten1iahpL609ih/SLmNi/454d1Jf0XGIqI215vPbpCM4Op0OqampiI2NbXA8NjYWSUlJN31OcnJyo/PHjBmDlJQU1NbWNntOU9esqamBVqtt8DAFSZLqby6edn0BMiJqPzduNv7pWD6ulNeIDUNko2r1BsxYdRgbj15CrcCb/k1acIqLi6HX6+Hr69vguK+vLwoKCm76nIKCgpueX1dXh+Li4mbPaeqaCxcuhFqtrn8EBQXd7ktqVvKFEmQUlsPRXoFJA0zzO4ioaf2C3dE3yA06vQFfH+SUcSIRNp/Ix44zRXhr8xmY9jui5rXLTcYyWcORDEmSGh271fm/P96aay5YsAAajab+kZub26r8LdWzoxp/uz8Sf7y7M9RO9ib5HUTUvKeGhgIAvjx4Ebo6Thknam83FrmdMjgYDnbiJmvbmfLiXl5eUCgUjUZWioqKGo3A3ODn53fT8+3s7ODp6dnsOU1dU6lUQqlU3u7LaDFXlT2eGdbJ5L+HiJo2rqc/3nRJR1F5DTafyMeEfh1FRyKyGWm5ZTiaUwYHhRxTBocIzWLSauXg4ICoqCgkJCQ0OJ6QkICYmJibPic6OrrR+du3b8eAAQNgb2/f7DlNXZOIbIeDnRzxQ4z/Yf1sfxZMPI+CiH5j1X7jMg3j+/jD28X0AwvNMfnY0dy5c/Hpp59i5cqVSE9Px5w5c5CTk4NZs2YBMH59NHXq1PrzZ82ahYsXL2Lu3LlIT0/HypUrsWLFCsybN6/+nJdeegnbt2/HokWLcObMGSxatAiJiYmYPXu2qV8OEVmAxwcHw0Ehx7FLGhzNLRMdh8gmFGmr8fOJfADAU2awTIrJC05cXBwWL16MN954A3379sWePXuwefNmhIQY/4WVn5/fYE2csLAwbN68Gbt27ULfvn3xz3/+E//73//wyCOP1J8TExODtWvX4rPPPkPv3r2xatUqrFu3DoMHDzb1yyEiC+DlrMSDfQMA/LptChGZ1pcHLqJWL2FAiDt6BapFxzH9OjjmyJTr4BCReTiZp8H4D/bBTi7Dvvn3wE+tEh2JyGrV1OkRs/AXlFTq8NET/XF/b3+T/B6zWQeHiEiUnh3VGBTqgTqDhC8PXBQdh8iq/XgsHyWVOvirVYjtcfMJP+2NBYeIrNaNKeNfH8pBda1ebBgiK2Vc5NZ4c3F8dAjsFeZRLcwjBRGRCdzb3Rcd3RxxtVKHTccui45DZJVSLpbi1GUtlHZyPD4wWHSceiw4RGS17BRyxEffmDKezSnjRCZwY/Tm4X4d4d7BQXCaX7HgEJFVmzwwCCp7OdLztTiUdVV0HCKrkldWhW2nCgEA069/JWwuWHCIyKq5OTng4X6BADhlnKitrU6+CL1BQnQnT0T4mdesZBYcIrJ6N2423n66ALlXr4kNQ2QlqnR6rDlkXMfuKTMbvQFYcIjIBnTzdcHQLp4wSOCUcaI28n1aHjRVtQjycMSoSPOYGv5bLDhEZBNuLB2/5lAOrunqBKchsmy/nRo+LToUCrlMcKLGWHCIyCbcE+GDEE8naKvrsPFonug4RBYt+XwJzhZWwMlBgccGBImOc1MsOERkE+RyGaZGhwIAVnHKONEdWXn9hv1H+gdC7WgvNkwTWHCIyGY8NiAQHRwUyCyqwL5zxaLjEFmkiyWV2HHGODV8Wkyo2DDNYMEhIpvhqrKvH05fxSnjRLfli+SLkCRgeDdvdPFxFh2nSSw4RGRTpl5f2fiXjCJkF1cKTkNkWSpq6vDN4VwA5jk1/LdYcIjIpnTydsbIcG9IErAqKVt0HCKLsuHIJZTX1KGTVweM6OotOk6zWHCIyOZMH2qcMv5d6iWUV9cKTkNkGQwGqf6r3WkxoZCb4dTw32LBISKbM7yrFzp7d0BFTR2+S70kOg6RRdideQUXiivhorTDI1GBouPcEgsOEdkcmUxWP4rzeVI2DAZOGSe6lRujN5MGBsFZaSc2TAuw4BCRTZrYryNcVHbILrmGXWeLRMchMmvniiqw++wVyGTGlYstAQsOEdmkDko7TB5onDLOXcaJmvf59RvyR0X4ItjTSWyYFmLBISKbNTU6FHIZsDezGJmF5aLjEJklTVUt1h8x3qtm7lPDf4sFh4hsVpCHE0Zf3wWZU8aJbu7blFxc0+kR7uuCmM6eouO0GAsOEdm0p67fbLzhSB401zhlnOi39AYJnydnAwCmDw2FTGbeU8N/iwWHiGzakE4eiPBzQVWtHutSckTHITIrO9ILkXu1Cm5O9pjQt6PoOK3CgkNENk0mk9XfV/B50kXU6Q1iAxGZkRs34E8eGAxHB4XYMK3EgkNENu+hvh3h7mSPvLIqJKZzyjgRAJwp0CL5QgkUchnir+/hZklYcIjI5qnsFXh8UDAA4LP9WYLTEJmHGwv7jenhi45ujmLD3AYWHCIiAPHRIVDIZTiYdRWnL2tFxyESqrRSh41H8wD8eiO+pWHBISIC4K92xNiefgCAVUkcxSHbtuZwDmrqDOgR4IoBIe6i49wWFhwiouuevn6z8fdpl1FSUSM2DJEgtXoDVidfBGAcvbGkqeG/xYJDRHRd/2B39Oqohq7OgLWHc0XHIRJi+6lC5Guq4dnBAeN7+4uOc9tMWnBKS0sRHx8PtVoNtVqN+Ph4lJWVNXl+bW0t5s+fj169eqFDhw4ICAjA1KlTcfny5Qbn3X333ZDJZA0ekydPNuVLISIb8Nsp46uTL6KWU8bJBt240X7K4GCo7C1ravhvmbTgPPHEE0hLS8PWrVuxdetWpKWlIT4+vsnzr127hiNHjuDvf/87jhw5gg0bNuDs2bN48MEHG507c+ZM5Ofn1z8++eQTU74UIrIR9/f2h5ezEgXaamw9WSA6DlG7OnFJg5SLpbCTy/DkEMubGv5bdqa6cHp6OrZu3YoDBw5g8ODBAIDly5cjOjoaGRkZCA8Pb/QctVqNhISEBsc++OADDBo0CDk5OQgODq4/7uTkBD8/P1PFJyIbpbRTYMrgYLy/IxOf7c/CA30CREciajefXb/B/v7e/vBxVQlOc2dMNoKTnJwMtVpdX24AYMiQIVCr1UhKSmrxdTQaDWQyGdzc3Boc/+qrr+Dl5YUePXpg3rx5KC9veifgmpoaaLXaBg8ioqZMGRIMe4UMR3LKcCy3THQconZxpbwGPx3LB2C5U8N/y2QFp6CgAD4+Po2O+/j4oKCgZcO+1dXVeOWVV/DEE0/A1dW1/viUKVOwZs0a7Nq1C3//+9+xfv16TJw4scnrLFy4sP4+ILVajaCgoNa/ICKyGT4uKjzQ2zhyw13GyVZ8fTAHOr0B/YLd0DfITXScO9bqgvPaa681usH394+UlBQAuOnUMkmSWjTlrLa2FpMnT4bBYMCSJUsa/GzmzJkYPXo0evbsicmTJ+O7775DYmIijhw5ctNrLViwABqNpv6Rm8vZEUTUvOnXbzb+6fhlFGqrxYYhMrGaOj2+PGicGj49JlRsmDbS6ntwXnjhhVvOWAoNDcXx48dRWFjY6GdXrlyBr69vs8+vra3FpEmTkJWVhV9++aXB6M3N9O/fH/b29sjMzET//v0b/VypVEKpVDZ7DSKi3+od6IYBIe5IuViKz/Zn45VxEaIjEZnM90fzcKW8Bn6uKtzXy3Knhv9WqwuOl5cXvLy8bnledHQ0NBoNDh06hEGDBgEADh48CI1Gg5iYmCafd6PcZGZmYufOnfD09Lzl7zp16hRqa2vh728dfylEZB6eG9EZKV+k4KsDF/HHkZ3hqrIXHYmozRkMEj7ZcwEAMOOuMNgrrGOJPJO9isjISIwdOxYzZ87EgQMHcODAAcycORPjx49vMIMqIiICGzduBADU1dXh0UcfRUpKCr766ivo9XoUFBSgoKAAOp0OAHD+/Hm88cYbSElJQXZ2NjZv3ozHHnsM/fr1w9ChQ031cojIBo2K8EFXH2eU19Th64M5ouMQmcT204W4cKUSrio7PD44+NZPsBAmrWlfffUVevXqhdjYWMTGxqJ3795YvXp1g3MyMjKg0WgAAJcuXcKmTZtw6dIl9O3bF/7+/vWPGzOvHBwcsGPHDowZMwbh4eF48cUXERsbi8TERCgUlrsgERGZH7lchmeHdwIArNyXhZo6veBERG1LkiR8vPs8AOOGs85Kk60e0+5kkiRJokO0N61WC7VaDY1Gc8v7e4jItunqDBj+9k4UaKux6JFeiBtoPf/CJTp4oQRxyw7AwU6O/fPvgbeLed+v2prPb+v4oo2IyEQc7OSYcZdxTZBP9lyAwWBz/yYkK3Zj9ObRqECzLzetxYJDRHQLjw8OhqvKDheuVGL76cazQ4ks0ZkCLXZmXIFcBjw7rJPoOG2OBYeI6BaclXaIjzbuy/Px7vOwwW/2yQot222cOTWupz9CvToITtP2WHCIiFpgekwYHOzkSMstw6Gsq6LjEN2RvLIqbDp2GQDw3AjrG70BWHCIiFrE20WJR6MCAfx63wKRpfp07wXUGSTEdPZE70A30XFMggWHiKiFnh3WCTIZsDPjCs4UcNNeskyllTqsPWTcsmjWiM6C05gOCw4RUQuFenXAuJ5+AH69f4HI0qw+cBFVtXp093fFsK633pnAUrHgEBG1wo1/8W46dhl5ZVWC0xC1TpVOj1VJ2QCM9960ZPNrS8WCQ0TUCr0D3RDT2RN1Bgmf7uUoDlmWb1NzcbVShyAPR9xvJZtqNoUFh4iolW6M4qw9lIvSSp3gNEQtU6c3YPn1Uj5zWCfYWcmmmk2x7ldHRGQCw7p6obu/K6pq9Vh94KLoOEQtsvlkAXKvVsGjgwMeiwoSHcfkWHCIiFpJJpPVrx2yKikbVTpuwknmTZIkfLzLuLzBtOhQODpY/+bULDhERLfh/l7+CHR3xNVKHb5LzRUdh6hZ+84V43S+Fo72Cky9viq3tWPBISK6DXYKOWZe379n2d4LqNMbBCciatqNxSnjBgbBvYOD4DTtgwWHiOg2TRoQBI8ODsi9WoXNJwtExyG6qROXNNh/rgQKuQzPDAsTHafdsOAQEd0mRwcFpkWHAgA+3sVNOMk83Ri9ebBPAALdnQSnaT8sOEREd2BqdAgc7RU4na/FvnPFouMQNZBdXIktJ/MBWO+mmk1hwSEiugPuHRwQN9A45ZabcJK5Wb73AgwScHe4NyL8XEXHaVcsOEREd+iZYWFQyGXYf64EJy5pRMchAgBcKa/Bt6mXAFj3pppNYcEhIrpDge5OeLBPAADg4z0cxSHz8HlSNnR1BvQNcsPgMA/RcdodCw4RURu4cX/DlhP5uFhSKTgN2bqKmjp8kZwNwDh6Y82bajaFBYeIqA1E+Lni7nBvGCRg2R5uwklirT2UA211HTp5dcC93X1FxxGCBYeIqI3cuM/h29RLuFJeIzgN2SpdnQGf7s0CADw7vBMUctsbvQFYcIiI2szgMA/0DXKDrs6Az5OyRcchG7Xp2GUUaKvh7aLEw/07io4jDAsOEVEbkclkmHX9XpwvkrNRUVMnOBHZGoNBwifXlyt4emgYlHbWv6lmU1hwiIja0L3d/dDJqwO01XVYeyhHdByyMb+cKUJmUQVclHaYMiRYdByhWHCIiNqQQi7Ds8ONozgr9mVBV8dNOKn9fHJ9mYInhgTDVWUvOI1YLDhERG3s4f4d4e2iRL6mGpuOXRYdh2xE6sWrOJxdCgeFHE8PtZ1NNZvCgkNE1MaUdor6D5hPdp+HwcBNOMn0lu4yLk/wcL+O8HVVCU4jHgsOEZEJTBkSDBelHTKLKvDLmSLRccjKZRaWIzG9EDIZ8KyNbarZFBYcIiITcFXZ44nrN3l+wu0byMRuLC4Z290Xnb2dBacxDyYtOKWlpYiPj4darYZarUZ8fDzKysqafc706dMhk8kaPIYMGdLgnJqaGvzpT3+Cl5cXOnTogAcffBCXLl0y4SshImq9p4eGwUEhx+HsUqRevCo6DlmpfE0Vvk/LAwA8Z4ObajbFpAXniSeeQFpaGrZu3YqtW7ciLS0N8fHxt3ze2LFjkZ+fX//YvHlzg5/Pnj0bGzduxNq1a7Fv3z5UVFRg/Pjx0Ov1pnopRESt5uuqwsP9jAut3bg/gqitrdyXhVq9hEFhHugf7C46jtmwM9WF09PTsXXrVhw4cACDBw8GACxfvhzR0dHIyMhAeHh4k89VKpXw8/O76c80Gg1WrFiB1atXY/To0QCAL7/8EkFBQUhMTMSYMWPa/sUQEd2mZ0d0wjepuUhML8S5onJ08XERHYmsiKaqFl8fNK639AeO3jRgshGc5ORkqNXq+nIDAEOGDIFarUZSUlKzz921axd8fHzQrVs3zJw5E0VFv96gl5qaitraWsTGxtYfCwgIQM+ePZu8bk1NDbRabYMHEVF76OztjNjrmx1+spujONS2vjxwEZU6PcJ9XXB3uLfoOGbFZAWnoKAAPj4+jY77+PigoKCgyeeNGzcOX331FX755Re88847OHz4MO655x7U1NTUX9fBwQHu7g2H4Xx9fZu87sKFC+vvA1Kr1QgKCrqDV0ZE1Do37ovYeDQP2cWVgtOQtSivrsWKfcZNNZ8b0QkymW1uqtmUVhec1157rdFNwL9/pKSkAMBN32xJkpr9S4iLi8P999+Pnj174oEHHsCWLVtw9uxZ/Pzzz83mau66CxYsgEajqX/k5ua24hUTEd2Z/sHuuDvcG3UGCf/dniE6DlmJZXsu4GqlDp28O+DBPgGi45idVt+D88ILL2Dy5MnNnhMaGorjx4+jsLCw0c+uXLkCX1/fFv8+f39/hISEIDMzEwDg5+cHnU6H0tLSBqM4RUVFiImJuek1lEollEpli38nEVFbmz82ArvPXsFPx/Mxc1gZ+gS5iY5EFqxIW41P9xpHb/4yJgJ2Cq768nutfke8vLwQERHR7EOlUiE6OhoajQaHDh2qf+7Bgweh0WiaLCI3U1JSgtzcXPj7+wMAoqKiYG9vj4SEhPpz8vPzcfLkyVZdl4ioPUX6u9bPqPr3ljOQJK5uTLdv8Y5MVNXq0T/YDWN6tHzQwJaYrPJFRkZi7NixmDlzJg4cOIADBw5g5syZGD9+fIMZVBEREdi4cSMAoKKiAvPmzUNycjKys7Oxa9cuPPDAA/Dy8sLDDz8MAFCr1ZgxYwZefvll7NixA0ePHsWTTz6JXr161c+qIiIyR3Pv7QYHhRzJF0qwJ7NYdByyUOevVGDdYeOtFq+Mi+S9N00w6ZjWV199hV69eiE2NhaxsbHo3bs3Vq9e3eCcjIwMaDQaAIBCocCJEyfw0EMPoVu3bpg2bRq6deuG5ORkuLj8OrXyvffew4QJEzBp0iQMHToUTk5O+PHHH6FQKEz5coiI7kiguxOmRocAMI7icI8quh3/3ZYBvUHC6EgfDArzEB3HbMkkGxwn1Wq1UKvV0Gg0cHV1FR2HiGxIaaUOw9/eifKaOiyO64sJ17+2ImqJIzmlmLgkCXIZsOWl4Qj3s611lVrz+c27koiI2pF7BwfMuts4bfy/2zNQU8cV2KllJEnCv7ecAQA80j/Q5spNa7HgEBG1s6eHhsHXVYlLpVX48kCO6DhkIXZmFOFQ1lUo7eSYc2830XHMHgsOEVE7c3RQYM5o4wfUh79kQltdKzgRmTu9QcKiLcY1lKYPDUWAm6PgROaPBYeISIBHowLR2bsDSq/VYhm3cKBb2Hg0DxmF5XBV2eGPI7qIjmMRWHCIiASwU8jxl7ERAIBP911AobZacCIyV9W1erx7fQXs50d2gdrJXnAiy8CCQ0QkSGx3X0SFuKO61oDFiZmi45CZ+iI5G5c11QhQqzAtJlR0HIvBgkNEJIhMJsMr44yjON+k5OJcUYXgRGRuNNdq8dHO8wCAOfd2g8qe6721FAsOEZFAA0M9MDrSF3qDhP9u40ac1NDS3eehqapFN19nTOwfKDqORWHBISIS7C9jwyGXAVtPFSD1YqnoOGQmLpdV4bP9xg0154+NgELOLRlagwWHiEiwbr4ueDTK+K/zRdyIk65bnHgWNXUGDAr1wD0RPqLjWBwWHCIiMzDn3m5Q2slxKPsqfjlTJDoOCXa2sBzfpV4CALxyXwQ31LwNLDhERGbAX+2Ip4aGAQAWbT0DPTfitGlvb82AQQLG9vBD/2B30XEsEgsOEZGZ+MOIzlA72uNsYQU2HLkkOg4Jcjj7KhLTC6GQy/DnseGi41gsFhwiIjOhdrLH8yONG3G+m3AW1bXciNPWSJKEhZvTAQCTBgShs7ez4ESWiwWHiMiMTI0ORYBahXxNNT5PyhYdh9rZ9tOFOJJTBkd7BeaM7io6jkVjwSEiMiMqewXmxhq/lvho5zlornEjTltRpzfg7a1nAAAz7gqDj6tKcCLLxoJDRGRmHu7XEeG+LtBW12HJ7nOi41A7+S71Es5fqYS7kz2eHdFJdByLx4JDRGRmFHIZ5o8zjuJ8tj8bl8uqBCciU6vS6fFe4lkAwAv3dIWrihtq3ikWHCIiMzQy3AeDwjygqzPgvYSzouOQia3cn4VCbQ0C3R3x5JBg0XGsAgsOEZEZkslkWHB9I871Ry4ho6BccCIyldJKHT7eZdxQc15sOJR23FCzLbDgEBGZqX7B7hjX0w8GCfjPtjOi45CJfLTzHMpr6hDp74oH+wSIjmM1WHCIiMzYvDHhUMhlSEwvwqGsq6LjUBu7VHoNXyRfBAC8Mi4Ccm6o2WZYcIiIzFhnb2fEDQwCACzcks6NOK3Mu9vPQqc3IKazJ4Z39RIdx6qw4BARmbnZo7rC0V6Bozll2HaqUHQcaiOnL2uxMS0PgHH0hhtqti0WHCIiM+fjqsIzw4wbcb697Qzq9AbBiagtvL3tDCQJGN/bH70D3UTHsTosOEREFuDZ4Z3g7mSPC1cq8W0qN+K0dEnni7Er4wrs5DLMi+WGmqbAgkNEZAFcVPb40z3GvYneSziLa7o6wYnodkmShEVbjLPinhgcjFCvDoITWScWHCIiCzFlSDAC3R1RVF6DFXuzRMeh2/TT8Xwcu6RBBwdFfWmltseCQ0RkIZR2Cvx5jPHrjA92nsOZAq3gRNRaV8pr8NqmUwCAmcM7wdtFKTiR9WLBISKyIA/2CcCoCB/o6gyYvTYN1bV60ZGohSRJwvz1x1FSqUOEnwv+cHdn0ZGsGgsOEZEFkclkWPRob3g5O+BMQTn+uy1DdCRqoa8O5uCXM0VwsJNj8eS+3JLBxExacEpLSxEfHw+1Wg21Wo34+HiUlZU1+xyZTHbTx3/+85/6c+6+++5GP588ebIpXwoRkdnwclZi0SO9AQCf7svC/nPFghPRrZy/UoF//XwaADB/bAQi/FwFJ7J+Ji04TzzxBNLS0rB161Zs3boVaWlpiI+Pb/Y5+fn5DR4rV66ETCbDI4880uC8mTNnNjjvk08+MeVLISIyK6MifTFlsHHX6Ze/OYayazrBiagptfobXycacFcXLzwVEyo6kk2wM9WF09PTsXXrVhw4cACDBw8GACxfvhzR0dHIyMhAePjN5/37+fk1+PMPP/yAkSNHolOnTg2OOzk5NTqXiMiWvHp/JJLPl+BCcSVe3XgSHz7Rj6vhmqH3EzNxIk8DtaM9/vtYH+431U5MNoKTnJwMtVpdX24AYMiQIVCr1UhKSmrRNQoLC/Hzzz9jxowZjX721VdfwcvLCz169MC8efNQXl7e5HVqamqg1WobPIiILJ2Tgx0WT+4LO7kMP5/Ix8ajeaIj0e8czr6KJbvOAQAWTuwFP7VKcCLbYbKCU1BQAB8fn0bHfXx8UFBQ0KJrfP7553BxccHEiRMbHJ8yZQrWrFmDXbt24e9//zvWr1/f6JzfWrhwYf19QGq1GkFBQa17MUREZqp3oBvm3NsNAPB/P5xC7tVrghPRDdrqWsxZlwaDBDzSPxD39fIXHcmmtLrgvPbaa03eCHzjkZKSAgA3HSqVJKnFQ6grV67ElClToFI1bLwzZ87E6NGj0bNnT0yePBnfffcdEhMTceTIkZteZ8GCBdBoNPWP3NzcVr5qIiLzNWtEZwwIcUdFTR3mrEuD3sAdx83Ba5tO4VJpFYI8HPHag91Fx7E5rb4H54UXXrjljKXQ0FAcP34chYWNd729cuUKfH19b/l79u7di4yMDKxbt+6W5/bv3x/29vbIzMxE//79G/1cqVRCqeRiSkRknRRyGd6L64tx7+9FysVSfLz7PJ4f2UV0LJv20/HL2HAkD3IZ8N6kvnBR2YuOZHNaXXC8vLzg5eV1y/Oio6Oh0Whw6NAhDBo0CABw8OBBaDQaxMTE3PL5K1asQFRUFPr06XPLc0+dOoXa2lr4+3P4j4hsU5CHE15/sAde/vYY3ks4i2FdvbhDtSD5miq8uvEkAOD5kV0wINRDcCLbZLJ7cCIjIzF27FjMnDkTBw4cwIEDBzBz5kyMHz++wQyqiIgIbNy4scFztVotvv32WzzzzDONrnv+/Hm88cYbSElJQXZ2NjZv3ozHHnsM/fr1w9ChQ031coiIzN7E/h1xfy9/1BkkzF6bxg05BTAYJMz79hg0VbXoHajGi6O415QoJl0H56uvvkKvXr0QGxuL2NhY9O7dG6tXr25wTkZGBjQaTYNja9euhSRJePzxxxtd08HBATt27MCYMWMQHh6OF198EbGxsUhMTIRCwVUhich2yWQyvPlwT/i5qnChuBJvbU4XHcnmrNyfhf3nSuBor8DiuL6wV3DDAFFkkiTZ3N1oWq0WarUaGo0Grq5cTZKIrMv+c8WY8ulBAMDK6QNwT8St73ukO3emQIsHP9gPnd6ANx/uiSmDQ0RHsjqt+fxmtSQisjJDu3jhmbvCAAB/+e44iitqBCeyftW1ery0Jg06vQGjInzwxKBg0ZFsHgsOEZEVmjcmHBF+Liiu0GH+d8dhg4P17eo/2zKQUVgOL2cHLHq0N1eUNgMsOEREVkhlr8DiyX3hoJBjx5kifH0oR3Qkq7Uvsxgr9mUBABY90htezlyWxByw4BARWakIP1f8Zaxx1uo/fzqN81cqBCeyPmXXdHj52zQAwJNDgjEqkvc7mQsWHCIiK/b00DDc1cUL1bUGzFmXhlq9QXQkqyFJEv668QQKtTXo5N0Br97H1YrNCQsOEZEVk8tl+O9jfaB2tMfxSxr8b0em6EhWY8ORPGw+UQA7uQyL4/rC0YFLlZgTFhwiIivnp1Zh4cReAICPdp7D4eyrghNZvtyr1/CPTacAAHPu7cZVo80QCw4RkQ24r5c/HukfCIMEzFmXhvLqWtGRLFad3vh1X0VNHQaEuGPWiM6iI9FNsOAQEdmI1x7sjkB3R1wqrcJrm06LjmOxPt59HikXS+GstMN7cX2hkHNKuDliwSEishEuKnu8F9cXchmw/sgl/Hw8X3Qki3MstwyLE433Mb3+YA8EeTgJTkRNYcEhIrIhA0M98Me7uwAA/rrxBAo01YITWY5rujrMWZeGOoOE+3v5Y2L/jqIjUTNYcIiIbMxLo7uid6AamqpazPv2GAwGrnLcEm/+nI4LxZXwc1XhzYd7crViM8eCQ0RkY+wVcuO0ZnsF9p0rxsIt6Sw5t/DlgYv46qBxNeh3JvWBm5OD4ER0Kyw4REQ2qJO3M15/sAcAYPneLLz87THo6rgI4O9JkoR3tmfgb9+fBAD84e7OGNrFS3Aqagk70QGIiEiMSQODIJMBr2w4gY1H81BcUYOlT0bBWcmPBgCo1Rvw6sYT+CblEgBg9uiueGlUV8GpqKU4gkNEZMMeGxCEFdMGwMlBgb2ZxZi8LBlF5bzx+JquDs9+kYJvUi5BLgMWTuyF2aO78b4bC8KCQ0Rk4+4O98GamUPg2cEBJ/O0eGRpEi7Y8MacJRU1eHzZAezMuAKVvRzL4gfg8UHBomNRK7HgEBER+gS5YcMfYxDi6YTcq1V49ONkHM0pFR2r3V0sqcQjS5Nw7JIG7k72+HrmEIzuzh3CLRELDhERAQBCPDtg/R9i0DtQjauVOjy+/AB2pBeKjtVuTlzS4JGlScguuYZAd0d894cY9A92Fx2LbhMLDhER1fNyVmLNzCG4O9wb1bUGPLs6FWsP5YiOZXK7z15B3LJkFFfo0N3fFRv+EIPO3s6iY9EdYMEhIqIGOijtsHzqADwaFQi9QcIrG07g/cRMSJJ1rpWzPvUSZqw6jGs6Pe7q4oV1zw2Bj6tKdCy6Qyw4RETUiL1Cjv882hsvjDRu6/Be4ln8deNJ1OmtZ60cSZKwZNc5vPztMdQZJEzoG4CV0wfCRWUvOhq1ARYcIiK6KZlMhnljwvHPCT0hkwFrDuVg1pdHUKXTi452x/QGCf/YdApvb80AADw3vBPendQXDnb8WLQW/JskIqJmxQ8JwdIpUVDayZGYXogpnx5AaaVOdKzbVl2rxwtfH8EXyRchkwH/N747FtwXCbmca9xYExYcIiK6pbE9/fDVM4OhdrTHkZwyPPJxEnKvXhMdq9U012oxdcUhbDlZAAeFHB883g9P3xUmOhaZAAsOERG1yIBQD3w3KxoBahUuXKnExKVJOHVZIzpWi10uq8KjHyfhUPZVuKjs8PnTgzC+d4DoWGQiLDhERNRiXX1dsOGPQxHh54Ir5TWI++QA9p8rFh3rls4UaDFxSRIyiyrg56rCt7OiEd3ZU3QsMiEWHCIiahU/tQrfzIrGkE4eqKipw/TPDuGHtDzRsZp04EIJHvs4GQXaanT1ccaGP8Ygws9VdCwyMRYcIiJqNVeV/fWvePxRq5fw0to0LN11Hro685lGrjdI+P5oHqauOITy6joMDHXHt7OiEeDmKDoatQOZZK0rNzVDq9VCrVZDo9HA1ZUtnojodhkMEt7cnI4V+7IAAM5KO4wI98a9kb64O9wbbk4O7Zrnmq4OezOLkXi6EL+cKULJ9dleY3v4YfHkvlDZK9o1D7Wt1nx+s+Cw4BAR3bHPk7Lx4c5zuFJeU39MIZdhYKg7Rkf64t7uvgjx7GCS312orcaO9CIkphdi37niBqNIrio7TBkSgnmx4VBwGrjFM5uC8+abb+Lnn39GWloaHBwcUFZWdsvnSJKE119/HcuWLUNpaSkGDx6Mjz76CD169Kg/p6amBvPmzcOaNWtQVVWFUaNGYcmSJQgMDGxRLhYcIqK2ZzBIOJ6nQeLpQiSmF+JMQXmDn3f1ccbo7r4YHemLvkFut104JEnCmYLy+t9z7FLDmVxBHo64N9IPo7v7YGCoB+wVvBvDWphNwfnHP/4BNzc3XLp0CStWrGhRwVm0aBHefPNNrFq1Ct26dcO//vUv7NmzBxkZGXBxcQEA/OEPf8CPP/6IVatWwdPTEy+//DKuXr2K1NRUKBS3Hn5kwSEiMr3cq9eQmG4sIQcvXEWd4dePGy9nB9wT4YPRkb64q6sXnBzsmr2Wrs6AQ1lXkZheiITThcgrq6r/mUwG9A1yqx8p6urjDJmMozXWyGwKzg2rVq3C7Nmzb1lwJElCQEAAZs+ejfnz5wMwjtb4+vpi0aJFeO6556DRaODt7Y3Vq1cjLi4OAHD58mUEBQVh8+bNGDNmzC3zsOAQEbUvTVUtdp+9gsTThdiZUYTy6rr6nynt5LirixdGd/fFqAif+o0uNddqsetsERJOF2J3xhWU1/z6HJW9HHd18ca93X0wMsIHPi7cHNMWtObzu/nK3M6ysrJQUFCA2NjY+mNKpRIjRoxAUlISnnvuOaSmpqK2trbBOQEBAejZsyeSkpJuWnBqampQU/Pr98Jarda0L4SIiBpQO9rjwT4BeLBPAGr1BhzOuoqE66M7uVersONMEXacKQIA9Alyg6O9HIezS6FvMOqjxOhI46jP0C5ecHTgDcPUNLMqOAUFBQAAX1/fBsd9fX1x8eLF+nMcHBzg7u7e6Jwbz/+9hQsX4vXXXzdBYiIiai17hRwxXbwQ08UL/ze+O84WVtR/9ZSWW4ZjuWX154b7umB0d2Op6RPoxv2iqMVaXXBee+21W5aFw4cPY8CAAbcd6vffnUqSdMvvU5s7Z8GCBZg7d279n7VaLYKCgm47HxERtQ2ZTIZwPxeE+7ng+ZFdUFRejV1nrqCmTo8R3XwQ7OkkOiJZqFYXnBdeeAGTJ09u9pzQ0NDbCuPn5wfAOErj7+9ff7yoqKh+VMfPzw86nQ6lpaUNRnGKiooQExNz0+sqlUoolcrbykRERO3Hx0WFSQP5D1C6c60uOF5eXvDy8jJFFoSFhcHPzw8JCQno168fAECn02H37t1YtGgRACAqKgr29vZISEjApEmTAAD5+fk4efIk3n77bZPkIiIiIsti0ntwcnJycPXqVeTk5ECv1yMtLQ0A0KVLFzg7OwMAIiIisHDhQjz88MOQyWSYPXs23nrrLXTt2hVdu3bFW2+9BScnJzzxxBMAALVajRkzZuDll1+Gp6cnPDw8MG/ePPTq1QujR4825cshIiIiC2HSgvN///d/+Pzzz+v/fGNUZufOnbj77rsBABkZGdBofl2k6S9/+Quqqqrwxz/+sX6hv+3bt9evgQMA7733Huzs7DBp0qT6hf5WrVrVojVwiIiIyPpxqwaug0NERGQRWvP5zfWriYiIyOqw4BAREZHVYcEhIiIiq8OCQ0RERFaHBYeIiIisDgsOERERWR0WHCIiIrI6LDhERERkdVhwiIiIyOqYdKsGc3Vj8WatVis4CREREbXUjc/tlmzCYJMFp7y8HAAQFBQkOAkRERG1Vnl5OdRqdbPn2OReVAaDAZcvX4aLiwtkMlmbXlur1SIoKAi5ubnc5+om+P40je9N8/j+NI/vT/P4/jTNkt4bSZJQXl6OgIAAyOXN32VjkyM4crkcgYGBJv0drq6uZv8/FJH4/jSN703z+P40j+9P8/j+NM1S3ptbjdzcwJuMiYiIyOqw4BAREZHVYcFpY0qlEv/4xz+gVCpFRzFLfH+axvemeXx/msf3p3l8f5pmre+NTd5kTERERNaNIzhERERkdVhwiIiIyOqw4BAREZHVYcEhIiIiq8OC04aWLFmCsLAwqFQqREVFYe/evaIjmY09e/bggQceQEBAAGQyGb7//nvRkczGwoULMXDgQLi4uMDHxwcTJkxARkaG6FhmY+nSpejdu3f9ImTR0dHYsmWL6FhmaeHChZDJZJg9e7boKGbhtddeg0wma/Dw8/MTHcus5OXl4cknn4SnpyecnJzQt29fpKamio7VJlhw2si6deswe/ZsvPrqqzh69CiGDRuGcePGIScnR3Q0s1BZWYk+ffrgww8/FB3F7OzevRvPP/88Dhw4gISEBNTV1SE2NhaVlZWio5mFwMBA/Pvf/0ZKSgpSUlJwzz334KGHHsKpU6dERzMrhw8fxrJly9C7d2/RUcxKjx49kJ+fX/84ceKE6Ehmo7S0FEOHDoW9vT22bNmC06dP45133oGbm5voaG2C08TbyODBg9G/f38sXbq0/lhkZCQmTJiAhQsXCkxmfmQyGTZu3IgJEyaIjmKWrly5Ah8fH+zevRvDhw8XHccseXh44D//+Q9mzJghOopZqKioQP/+/bFkyRL861//Qt++fbF48WLRsYR77bXX8P333yMtLU10FLP0yiuvYP/+/Vb7bQNHcNqATqdDamoqYmNjGxyPjY1FUlKSoFRkqTQaDQDjhzg1pNfrsXbtWlRWViI6Olp0HLPx/PPP4/7778fo0aNFRzE7mZmZCAgIQFhYGCZPnowLFy6IjmQ2Nm3ahAEDBuCxxx6Dj48P+vXrh+XLl4uO1WZYcNpAcXEx9Ho9fH19Gxz39fVFQUGBoFRkiSRJwty5c3HXXXehZ8+eouOYjRMnTsDZ2RlKpRKzZs3Cxo0b0b17d9GxzMLatWuRmprKkeKbGDx4ML744gts27YNy5cvR0FBAWJiYlBSUiI6mlm4cOECli5diq5du2Lbtm2YNWsWXnzxRXzxxReio7UJm9xN3FRkMlmDP0uS1OgYUXNeeOEFHD9+HPv27RMdxayEh4cjLS0NZWVlWL9+PaZNm4bdu3fbfMnJzc3FSy+9hO3bt0OlUomOY3bGjRtX/3/36tUL0dHR6Ny5Mz7//HPMnTtXYDLzYDAYMGDAALz11lsAgH79+uHUqVNYunQppk6dKjjdneMIThvw8vKCQqFoNFpTVFTUaFSHqCl/+tOfsGnTJuzcuROBgYGi45gVBwcHdOnSBQMGDMDChQvRp08fvP/++6JjCZeamoqioiJERUXBzs4OdnZ22L17N/73v//Bzs4Oer1edESz0qFDB/Tq1QuZmZmio5gFf3//Rv9IiIyMtJrJMSw4bcDBwQFRUVFISEhocDwhIQExMTGCUpGlkCQJL7zwAjZs2IBffvkFYWFhoiOZPUmSUFNTIzqGcKNGjcKJEyeQlpZW/xgwYACmTJmCtLQ0KBQK0RHNSk1NDdLT0+Hv7y86ilkYOnRooyUpzp49i5CQEEGJ2ha/omojc+fORXx8PAYMGIDo6GgsW7YMOTk5mDVrluhoZqGiogLnzp2r/3NWVhbS0tLg4eGB4OBggcnEe/755/H111/jhx9+gIuLS/1IoFqthqOjo+B04v31r3/FuHHjEBQUhPLycqxduxa7du3C1q1bRUcTzsXFpdG9Wh06dICnpyfv4QIwb948PPDAAwgODkZRURH+9a9/QavVYtq0aaKjmYU5c+YgJiYGb731FiZNmoRDhw5h2bJlWLZsmehobUOiNvPRRx9JISEhkoODg9S/f39p9+7doiOZjZ07d0oAGj2mTZsmOppwN3tfAEifffaZ6Ghm4emnn67//ytvb29p1KhR0vbt20XHMlsjRoyQXnrpJdExzEJcXJzk7+8v2dvbSwEBAdLEiROlU6dOiY5lVn788UepZ8+eklKplCIiIqRly5aJjtRmuA4OERERWR3eg0NERERWhwWHiIiIrA4LDhEREVkdFhwiIiKyOiw4REREZHVYcIiIiMjqsOAQERGR1WHBISIiIqvDgkNERERWhwWHiIiIrA4LDhEREVkdFhwiIiKyOv8P2SR1rxmevpgAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "b = torch.sin(a)\n",
    "plt.plot(a.detach(), b.detach())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's have a closer look at the tensor `b`. When we print it, we see an indicator that it is tracking its computation history:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "tensor([ 0.0000e+00,  2.5882e-01,  5.0000e-01,  7.0711e-01,  8.6603e-01,\n",
      "         9.6593e-01,  1.0000e+00,  9.6593e-01,  8.6603e-01,  7.0711e-01,\n",
      "         5.0000e-01,  2.5882e-01, -8.7423e-08, -2.5882e-01, -5.0000e-01,\n",
      "        -7.0711e-01, -8.6603e-01, -9.6593e-01, -1.0000e+00, -9.6593e-01,\n",
      "        -8.6603e-01, -7.0711e-01, -5.0000e-01, -2.5882e-01,  1.7485e-07],\n",
      "       grad_fn=<SinBackward0>)\n"
     ]
    }
   ],
   "source": [
    "print(b)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This `grad_fn` gives us a hint that when we execute the backpropagation step and compute gradients, we'll need to compute the derivative of $sin(x)$ for all this tensor's inputs.\n",
    "\n",
    "Let's perform some more computations:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "tensor([ 0.0000e+00,  5.1764e-01,  1.0000e+00,  1.4142e+00,  1.7321e+00,\n",
      "         1.9319e+00,  2.0000e+00,  1.9319e+00,  1.7321e+00,  1.4142e+00,\n",
      "         1.0000e+00,  5.1764e-01, -1.7485e-07, -5.1764e-01, -1.0000e+00,\n",
      "        -1.4142e+00, -1.7321e+00, -1.9319e+00, -2.0000e+00, -1.9319e+00,\n",
      "        -1.7321e+00, -1.4142e+00, -1.0000e+00, -5.1764e-01,  3.4969e-07],\n",
      "       grad_fn=<MulBackward0>)\n",
      "tensor([ 1.0000e+00,  1.5176e+00,  2.0000e+00,  2.4142e+00,  2.7321e+00,\n",
      "         2.9319e+00,  3.0000e+00,  2.9319e+00,  2.7321e+00,  2.4142e+00,\n",
      "         2.0000e+00,  1.5176e+00,  1.0000e+00,  4.8236e-01, -3.5763e-07,\n",
      "        -4.1421e-01, -7.3205e-01, -9.3185e-01, -1.0000e+00, -9.3185e-01,\n",
      "        -7.3205e-01, -4.1421e-01,  4.7684e-07,  4.8236e-01,  1.0000e+00],\n",
      "       grad_fn=<AddBackward0>)\n"
     ]
    }
   ],
   "source": [
    "c = 2 * b\n",
    "print(c)\n",
    "\n",
    "d = c + 1\n",
    "print(d)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Finally, let's compute a single-element output. When you call `.backward()` on a tensor with no arguments, it expects the calling tensor to contain only a single element, as is the case when computing a loss function."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "tensor(25., grad_fn=<SumBackward0>)\n"
     ]
    }
   ],
   "source": [
    "out = d.sum()\n",
    "print(out)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Each `grad_fn` stored with our tensors allows you to walk the computation all the way back to its inputs with its `next_functions` property. We can see below that drilling down on this property on `d` shows us the gradient functions for all the prior tensors. Note that `a.grad_fn` is reported as `None`, indicating that this was an input to the function with no history of its own."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "d:\n",
      "<AddBackward0 object at 0x0000019BAD20E1F0>\n",
      "((<MulBackward0 object at 0x0000019BAD20E2E0>, 0), (None, 0))\n",
      "((<SinBackward0 object at 0x0000019BAD20E1F0>, 0), (None, 0))\n",
      "((<AccumulateGrad object at 0x0000019BAD20E2E0>, 0),)\n",
      "()\n",
      "\n",
      "c:\n",
      "<MulBackward0 object at 0x0000019BAD20E340>\n",
      "\n",
      "b:\n",
      "<SinBackward0 object at 0x0000019BAD20E340>\n",
      "\n",
      "a:\n",
      "None\n"
     ]
    }
   ],
   "source": [
    "print('d:')\n",
    "print(d.grad_fn)\n",
    "print(d.grad_fn.next_functions)\n",
    "print(d.grad_fn.next_functions[0][0].next_functions)\n",
    "print(d.grad_fn.next_functions[0][0].next_functions[0][0].next_functions)\n",
    "print(d.grad_fn.next_functions[0][0].next_functions[0][0].next_functions[0][0].next_functions)\n",
    "print('\\nc:')\n",
    "print(c.grad_fn)\n",
    "print('\\nb:')\n",
    "print(b.grad_fn)\n",
    "print('\\na:')\n",
    "print(a.grad_fn)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "With all this machinery in place, how do we get derivatives out? You call the `backward()` method on the output, and check the input's `grad` property to inspect the gradients:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "tensor([ 2.0000e+00,  1.9319e+00,  1.7321e+00,  1.4142e+00,  1.0000e+00,\n",
      "         5.1764e-01, -8.7423e-08, -5.1764e-01, -1.0000e+00, -1.4142e+00,\n",
      "        -1.7321e+00, -1.9319e+00, -2.0000e+00, -1.9319e+00, -1.7321e+00,\n",
      "        -1.4142e+00, -1.0000e+00, -5.1764e-01,  2.3850e-08,  5.1764e-01,\n",
      "         1.0000e+00,  1.4142e+00,  1.7321e+00,  1.9319e+00,  2.0000e+00])\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x7fa0004bbfd0>]"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD4CAYAAADvsV2wAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO3dd3yV5f3/8dcng4SVAFkEws5gzxRU0CIgIKJUq1brqG0t5asiVmu1am31W7usHU5EcIF1fFVELYoLRYZoQJQNIawww0rIzkk+vz9y9EdpAoRzkuuMz/PxyCMn59zkfsfx5sp17vu6RFUxxhgT+iJcBzDGGNM0rPCNMSZMWOEbY0yYsMI3xpgwYYVvjDFhIsp1gBNJTEzUrl27uo5hjDFBY8WKFQdUNamu1wK68Lt27UpOTo7rGMYYEzREZHt9r9mUjjHGhAkrfGOMCRNW+MYYEyas8I0xJkxY4RtjTJjwufBFpJOILBSR9SKyVkSm1XGMiMjDIpIrIl+LyGBfz2uMMaZh/HFZpge4TVVXikhrYIWIvK+q64455nwgw/sxDHjC+9kYY0wT8XmEr6p7VHWl9/FRYD3Q8bjDJgHPa63PgDYikurruevzyIebefOr3RwormisUxhjjN/tLypn7pf5TP9kS6N8f7/eeCUiXYFBwPLjXuoI7Dzm63zvc3vq+B6TgckAnTt3bnCG8qpqnl6ylcOlVQD0So1jRHoCw9MTGdqtHS2aBfS9ZsaYMFJc4WF53kEW5x5gSe4BNu0rBqB9XCw/O7s7kRHi1/P5rf1EpBXwGnCLqhYd/3Idf6TOnVdUdQYwAyA7O7vBu7PERkeSc895rNlV+O0/xOeWbuepT7cSHSkM6tyWEemJDE9PZEBaPFGR9r61MaZpVFXXsGrnERZvru2mVTuP4KlRYqIiGNqtHd8fnMbw9ER6p8YR4eeyBxB/7HglItHA28ACVf1bHa8/CXysqi96v94IjFTV/xrhHys7O1v9sbRCWWU1OdsPffsXwNrdRahC65gohnVPYER6AmdnJtEjqZXP5zLGmGNt3neURd6CX553kJLKaiIE+qW1+Xb2YXDntsRGR/rlfCKyQlWz63rN5xG+iAgwC1hfV9l7vQncJCIvUftmbeHJyt6fmjeL5OyMJM7OqF1P6HBJJcuO+TXqg/X7AJh8Tnd+NS7LRv3GGJ9VeKr5/dvrmf1Z7dI23RNbcol3BH9m9wTiW0Q3eSZ/TOkMB64BVovIKu9zdwGdAVR1OjAfmADkAqXAj/1w3tPWtmUzJvRLZUK/2veNdx4q5clFW5ixKI+vdh7h0R8OJql1jMuIxpggtvtIGTe8sJJVO49w/Yhu/HhENzq2ae46ln+mdBqLv6Z0TtXrK/O5a+5q4ptH8/hVgxnSpV2TndsYExqW5B5g6otfUlFVzV8vG8D5/RrtgsQ6nWhKx+YujnHJ4DTm3jCc2OhIfvDkZzyzZCuB/BeiMSZwqCqPf5zLNbOW065lM+bdNKLJy/5krPCP0ys1jjdvGsHIrGTue2sd015aRUmFx3UsY0wAKyqvYvLsFfzl3Y1M6JfKvBuHk54ceBeB2EXpdYhvHs2Ma4bwxCdbeOi9jWzYW8T0q4fQ3a7iMcYcZ8PeIqbMXkH+4TLundibHw/vSu21LIHHRvj1iIgQbjw3ned/MowDxZVc9OgS3l2z13UsY0wAeePLXXzvsSWUVlbz4uQz+MmIbgFb9mCFf1IjMhJ5e+oIeiS3YsqcFfzxnfV4qmtcxzLGOFTpqeHeeWu45eVV9E9rw9s3j+A7XQP/Ig8r/FPQoU1zXvn5GVw1rDNPfpLHNbM+t3V6jAlTewrL+MGMZTy/bDvXj+jGC9cPI7l1rOtYp8QK/xTFREXywMX9eOiyAazccZiJDy9mxfbDrmMZY5rQ0i0HmPjwYjbtPcpjPxzMPRN7Ex1EN2oGT9IA8f0htZduNouK4OqZy1m/5/hlg4wxoejLHYe57ukvaNuyGfNuGs4F/QPrkstTYYV/Gnp3iOPVKWcS1zyKybNzOFxS6TqSMaYR7S8qZ8qcFaTEx/B/Pz+T9OTWriOdFiv805QcF8v0q4ewr7CCm15caW/kGhOiKjzVTJmzgqIyDzOuyaZty2auI502K3wfDOrclgcu7suS3IP88Z0NruMYY/xMVfntvLWs3HGEhy4fQK/UONeRfGI3XvnosuxOrN1dxKzFW+mdGsf3h6S5jmSM8ZM5n23npS92ctO56d8uthjMbITvB3df0Iszuyfw67mr+Tr/iOs4xhg/WJ53kPveWsfonsncel6m6zh+YYXvB9GRETx21WCSWsXw89krKDhq1+gbE8x2eZc37pzQgr9fMbBRdp9ywQrfT9q1bMaMa4dwuLSS/5mzgkqPvYlrTDAqq6zm57NzqPTU8NS12cTFNv1GJY3FCt+P+nSI58FLB5Cz/TD3vbXWdRxjTAOpKne+/jVrdxfxzysHhty2p34pfBF5WkT2i8iael4fKSKFIrLK+3GvP84biC4c0IEp3+3BC8t38K/lO1zHMcY0wMxPtzJv1W5+OTaLUT1TXMfxO3+N8J8Fxp/kmE9VdaD3434/nTcg3T4ui5FZSfz2zTXkbDvkOo4x5hQs2lTAH99Zz4R+7blhZA/XcRqFXwpfVRcB1mxekRHCP68YRFrbFkyZs5I9hWWuIxljTmD7wRKmvvglmSmtefDSAQG9xLEvmnIO/0wR+UpE3hGRPvUdJCKTRSRHRHIKCgqaMJ5/fbOJSlmlhymzV1BeVe06kjGmDsUVHn72fA4i8NS12bSMCd3bk5qq8FcCXVR1APAI8EZ9B6rqDFXNVtXspKSkJorXODJSWvP3Hwzkq/xC7p67xvbHNSbA1NQot72yii0FJTz2w8F0atfCdaRG1SSFr6pFqlrsfTwfiBaRxKY4t2tj+7TnF2MyeW1lPs8s2eY6jjHmGI8uzGXB2n3cNaEXw9NDv5KapPBFpL14J8VEZKj3vAeb4tyBYOqodMb1SeGB+etZuuWA6zjGGODD9fv42/ubuGRQR34yvKvrOE3CX5dlvggsA7JEJF9EfioiU0RkiveQS4E1IvIV8DBwhYbR/EZEhPDQ5QPpktCC2//va8oqbT7fGJcKy6q447XV9OkQxx8u6Reyb9Iezy/vTqjqlSd5/VHgUX+cK1i1ioniT5f05/Inl/HPDzdz5/k9XUcyJmw9uGADh0oqePbH3yE2OtJ1nCZjd9o2oaHd2nF5dhozP81j496jruMYE5a+3HGYF5bv4LqzutG3Y7zrOE3KCr+J3Xl+L1rHRnH33NXU1ITNrJYxAcFTXcNdc9eQ0jqWW8eGxgqYDWGF38TatWzGXRN6kbP9MK/k7HQdx5iw8uzSbazfU8TvLupNqxC+3r4+VvgOXDokjaHd2vHHdzZwsNiWUjamKew+Usbf3t/EqJ7JjOvT3nUcJ6zwHRAR/nBxX0orPTwwf73rOMaEhd+9uZYaVe67qE/YXJVzPCt8R9KTW/Pzc3rw+spddm2+MY3s/XX7eG/dPqaNzgz5u2lPxArfoZtGpdO5XQvueWMNFR67Nt+YxlBa6eF3b64lM6UV15/dzXUcp6zwHYqNjuT+SX3IKyjhyU/yXMcxJiT944PN7DpSxh8u7kd0ZHhXXnj/9AFgZFYyF/RP5dGFuWw7UOI6jjEhZf2eImYt3soV3+lEdtd2ruM4Z4UfAH47sTcxkRH8Zp6tqGmMv9TUKHfNXU1882juGG93toMVfkBIjovl9vFZfLr5AG99vcd1HGNCwktf7OTLHUe4e0Iv2rZs5jpOQLDCDxBXDetC/7R47n9rHYVlVa7jGBPUCo5W8Kd31nNG93ZcMrij6zgBwwo/QERGCH+4uB+HSip4cMEG13GMCWoP/HsdZVXV/P574bMS5qmwwg8gfTvG86OzuvLC8h18ueOw6zjGBKUluQd4Y9Vu/ue7PUhPbuU6TkCxwg8wt43NIqV1LHfPXYOnusZ1HGOCSnlVNfe8sYYuCS244dx013ECjhV+gGkVE8XvLurNuj1FPLt0m+s4xgSV6Z9sYeuBEv53Ut+wWuf+VPlrx6unRWS/iKyp53URkYdFJFdEvhaRwf44b6ga16c9o3om87f3N7H7SJnrOMYEhbyCYh5fuIULB3TgnMwk13ECkr9G+M8C40/w+vlAhvdjMvCEn84bkkSE+y7qQ40qv3tzres4xgQ8VeWeN9YQEx3Bbyb2ch0nYPml8FV1EXDoBIdMAp7XWp8BbUQk1R/nDlWd2rVg2uhM3lu3j4837ncdx5iANn/1XpZuOcivxmWR3DrWdZyA1VRz+B2BY3f7yPc+919EZLKI5IhITkFBQZOEC1Q/HdGNzu1a8Od3N9ruWMbUo6q6hr++t5GslNb8cFgX13ECWlMVfl0XwtbZYKo6Q1WzVTU7KSm85+GaRUVw29hM1u8p4q2vd7uOY0xAeiVnJ1sPlHD7uCwiI+ya+xNpqsLPBzod83UaYA12Ci7s34Ge7Vvz0HubqPTYZZrGHKusspp/frCZIV3aMrpXsus4Aa+pCv9N4Frv1TpnAIWqaovGnIKICOGO8T3ZcaiUl20PXGP+w7NLt7H/aAV3jO9pd9SeAn9dlvkisAzIEpF8EfmpiEwRkSneQ+YDeUAu8BRwgz/OGy5GZiUxtGs7Hv5wM6WVHtdxjAkIhaVVPPFxLqN6JjO0my19fCr8sm27ql55ktcVuNEf5wpHIsId52fx/SeW8cySbdxodxAaw/RFWzha4eH2cVmuowQNu9M2SAzp0o4xvVKY/vEWDpdUuo5jjFP7isp5ZslWJg3oQK/UONdxgoYVfhC5fVwWxZUepn+yxXUUY5z654eb8VQrt55no/uGsMIPIlntW3PxoI48u3QbewptyQUTnrYeKOHlL3byw2Gd6ZzQwnWcoGKFH2R+MSaTGlUe/nCz6yjGOPHQexuJiYpg6qgM11GCjhV+kOnUrgVXDevCKzn5bCkodh3HmCa1Zlchb3+9h5+O6EZS6xjXcYKOFX4QumlUOrFRETz03kbXUYxpUn9ZsJE2LaL52TndXUcJSlb4QSixVQzXn92d+av38tXOI67jGNMklm45wKJNBdw4Mp242GjXcYKSFX6Quv7sbrRr2YwHF9go34Q+VeXP724kNT6Wa860BdJOlxV+kGodG82N56azOPcAizcfcB3HmEa1YO0+vtp5hF+MybSdrHxghR/ErhrWmY5tmvPndzdQezOzMaHH413+uEdSSy4ZXOeq6uYUWeEHsdjoSG4Zk8HqXYW8s2av6zjGNIrXv9xF7v5ibh+XRVSkVZYv7J9ekLtkcBoZya3464KNeKpt+WQTWsqrqvnH+5sY0KkN4/q0dx0n6FnhB7nICOH2cVnkHSjh1RX5ruMY41dzPtvO7sJy7hifZcsf+4EVfgg4r3cKgzq34R8fbKa8qtp1HGP8oqi8iscW5nJ2RiJn9Uh0HSckWOGHAJHaTVL2FpXz3NJtruMY4xczF+VxuLSKX43r6TpKyPDXBijjRWSjiOSKyJ11vH6diBSIyCrvx/X+OK/5/87onsB3M5N4/OMtFJZVuY5jjE8KjlYwc/FWLuifSr+0eNdxQobPhS8ikcBjwPlAb+BKEeldx6Evq+pA78dMX89r/tvt47IoLKtixiJbPtkEt8cW5lLhqeG28zJdRwkp/hjhDwVyVTVPVSuBl4BJfvi+poH6doznwgEdeHrxNg4UV7iOY8xp2XWkjBeWb+fy7E50T2rlOk5I8UfhdwSO3V073/vc8b4vIl+LyKsi0qm+byYik0UkR0RyCgoK/BAvvEwbnUG5p5qnPs1zHcWY0/L4wlygdpFA41/+KPy6rpU6/rbPt4Cuqtof+AB4rr5vpqozVDVbVbOTkpL8EC+8pCe34sL+HZi9bDuHbCtEE2R2HynjlZydXJbdiY5tmruOE3L8Ufj5wLEj9jRg97EHqOpBVf1mjuEpYIgfzmvqcfPodMqqbJRvgs8323feMLKH4yShyR+F/wWQISLdRKQZcAXw5rEHiEjqMV9eBKz3w3lNPdKTW3NBv1SeX7rNNjw3QWNvYTkvfb6TS4ekkdbWti5sDD4Xvqp6gJuABdQW+SuqulZE7heRi7yH3Swia0XkK+Bm4Dpfz2tO7ObRGZRWVTNzsY3yTXCY/skWalS5YaTN3TeWKH98E1WdD8w/7rl7j3n8a+DX/jiXOTWZKa2Z0DeV55Zu52dnd6dNi2auIxlTr31F5fzr8x18f3AandrZ6L6x2J22IWzq6HSKKzzMWrzVdRRjTmj6J1uorlFuPNdG943JCj+E9Wwfx/l92/Pskm0UltrdtyYw7S8q51/Ld3DJoI50TrDRfWOywg9xU0dlcLTCw6wlNso3genJRXl4atSuu28CVvghrneHOMb2TuGZJVttjR0TcAqOVvDC8u1MGtiBLgktXccJeVb4YeDm0RkcLffw7JJtrqMY8x9mLNpCpaeGqaMyXEcJC1b4YaBvx3jG9Eph1uI8isptlG8Cw4HiCmZ/tp1JAzvSLdFG903BCj9MTBudQVG5h+dslG8CxFOL8qj01NjcfROywg8T/dLiGd0zmZmLt3LURvnGsYPFFTy/bDsXDuhAD1sRs8lY4YeRaWMyKCyr4vll211HMWFu5uKtlHuqmWqj+yZlhR9G+qe14dysJJ76NI/iCo/rOCZMHS6p5Pml25jYvwPpya1dxwkrVvhhZtqYTI6UVvH8sm2uo5gwNXNxHqVV1dxso/smZ4UfZgZ2asN3M5N4alEeJTbKN03sSGklzy3dzoR+qWSk2Oi+qVnhh6FpYzI4XFrFnM9sLt80rVmLt1Jc4eFmu+7eCSv8MDS4c1vOzkhkxqI8SittlG+aRmFpFc8u2caEfu3Jam+jexes8MPUtNEZHCyp5IXPdriOYsLErCVbOVrhsbtqHbLCD1PZXdsxPD2BJxdtoayy2nUcE+IKy6p4ZslWxvVJoVdqnOs4YcsvhS8i40Vko4jkisiddbweIyIve19fLiJd/XFe45tpozM5UFzJC8ttLt80rmeWbOVouYebR9vo3iWfC19EIoHHgPOB3sCVItL7uMN+ChxW1XTg78CffT2v8d3Qbu04s3sCTy7Ko7zKRvmmcRSVV/H04q2c1zuFPh3iXccJa/4Y4Q8FclU1T1UrgZeASccdMwl4zvv4VWC0iIgfzm18NG1MBgVHK/jXcpvLN43j2SXbKCr3MM1G9875o/A7AjuP+Trf+1ydx3g3PS8EEur6ZiIyWURyRCSnoKDAD/HMiZzRPYFh3dox/ZMtNso3fne0vIpZi7cyplcyfTva6N41fxR+XSN1PY1jap9UnaGq2aqanZSU5HM4c3LTxmSw/2gFL31uo3zjX88t3UZhWRXTRme6jmLwT+HnA52O+ToN2F3fMSISBcQDh/xwbuMHZ3ZPYGjXdjxho3zjR8UVHmYu3sqonsn0S7PRfSDwR+F/AWSISDcRaQZcAbx53DFvAj/yPr4U+EhV6xzhm6YnIkwbk8G+ogpeydl58j9gzCl4ftk2jpRW2dx9APG58L1z8jcBC4D1wCuqulZE7heRi7yHzQISRCQXuBX4r0s3jVtn9Uggu0tbnvh4CxUeG+Ub35RUeHhqUR4js5IY0KmN6zjGyy/X4avqfFXNVNUeqvqA97l7VfVN7+NyVb1MVdNVdaiq5vnjvMZ/vhnl7yks5/9y8l3HMUFu9mfbOVxaZdfdBxi709Z8a0R6IoM6t+GJj2s3ljbmdJRW1o7uz85IZHDntq7jmGNY4ZtviQjTRmew60gZr66wUb45PXM+287BkkpuGWOj+0BjhW/+w3cza+dcH1uYa6N802BlldXMWJTHiPREhnRp5zqOOY4VvvkPIsIt3lH+6yttlG8a5oXl2zlQXMk0G90HJCt8819GZiXRPy2eRxfmUlVto3xzasoqq5n+SR5n9UjgO11tdB+IrPDNf/lmLj//cBlzV+5yHccEiX99voMDxRV23X0As8I3dRrVM5l+HW2Ub05NeVU10z/Zwhnd2zGse53LZJkAYIVv6iQi3Dw6gx2HSnnjSxvlmxN76fMdFBytsDVzApwVvqnXmF7J9OkQx6MLc/HYKN/Uo7yqmic+2VK7v0IPG90HMit8U69vRvnbD5Yyb9Xx6+EZU+uVnJ3sK6rgFpu7D3hW+OaExvau3YPURvmmLhWeap74eAvf6drWRvdBwArfnFDtFTvpbD1Qwltf2yjf/KdXcvLZU1jOtNGZ2CZ2gc8K35zU2N7t6dm+NY98lEt1ja1qbWpVeKp5YmEuQ7q0ZXi6je6DgRW+OamIiNq5/LyCEt62Ub7xenVFPrsLy5k2OsNG90HCCt+ckvF92pOZ0oqHP9xso3xDpaeGxxduYWCnNpydkeg6jjlFVvjmlERECFNHZbCloIR/r97jOo5x7LWV+ew6Usa0MTa6DyY+Fb6ItBOR90Vks/dznYtfi0i1iKzyfhy//aEJEhP6pZKe3IpHPtxMjY3yw1ZVdQ2PLcxlQFo8IzOTXMcxDeDrCP9O4ENVzQA+pP6tC8tUdaD346J6jjEBLjJCmDoqnc37i3lnzV7XcYwjc1fuIv+wje6Dka+FPwl4zvv4OeB7Pn4/E+Am9u9Aj6SWPGyj/LBUVV3Dowtz6Z8Wz7lZya7jmAbytfBTVHUPgPdzff8FxIpIjoh8JiIn/EtBRCZ7j80pKCjwMZ7xt0jvXP7GfUeZv8bm8sPN3JW72HGolJtH2eg+GJ208EXkAxFZU8fHpAacp7OqZgM/BP4hIj3qO1BVZ6hqtqpmJyXZ/GAgunBABzKSW/HQe5tsJc0wUl5Vzd8/2MSATm0Y3ctG98HopIWvqmNUtW8dH/OAfSKSCuD9vL+e77Hb+zkP+BgY5LefwDS5yAjh9nFZbD1QYnvfhpE5n21nT2E5d4zPstF9kPJ1SudN4Efexz8C5h1/gIi0FZEY7+NEYDiwzsfzGsfO653C4M5t+McHmyivqnYdxzSyovIqHluYy9kZiZzVw667D1a+Fv6fgPNEZDNwnvdrRCRbRGZ6j+kF5IjIV8BC4E+qaoUf5ESEO8b3ZF9RBc8t3eY6jmlkMxflcbi0ijvG93Qdxfggypc/rKoHgdF1PJ8DXO99vBTo58t5TGAa1j2BkVlJPP7xFq4Y2pn45tGuI5lGUHC0gpmLtzKxfyp9O8a7jmN8YHfaGp/cPi6LwrIqnvxki+soppE8+tFmKjw13DY2y3UU4yMrfOOTPh3iuWhAB55espX9ReWu4xg/23GwlH99voMffKcT3RJbuo5jfGSFb3x263mZeKqVhz/a7DqK8bO/f7CJyAhhmu1mFRKs8I3Puia25MqhnXnp851sO1DiOo7xk/V7inhj1S5+PLwbKXGxruMYP7DCN34xdVQ60ZER/O39Ta6jGD95cMFGWsdEMeWceu+TNEHGCt/4RXJcLD8Z0ZU3v9rN2t2FruMYH32+9RAfbdjP/4xMJ76FXX0VKqzwjd9MPqcH8c2j+cu7G11HMT5QVf7y7gZS4mK47qyuruMYP7LCN34T3zyaG8/twSebCli25aDrOOY0fbRhPznbDzNtdCbNm0W6jmP8yArf+NW1Z3alfVwsf1mwAVVbPjnYVNcof3l3I90SW3JZdprrOMbPrPCNX8VGR3LLmAy+3HGE99ftcx3HNNC8VbvYuO8ot43NJDrS6iHU2L9R43eXDkmje2JLHlyw0TY8DyIVnmr+9v4m+naMY0LfVNdxTCOwwjd+FxUZwS/HZbF5fzFzv9zlOo45RS8u30H+4TJ+Na4nERG2/HEossI3jeL8vu3pnxbP39+35ZODQXGFh0c+yuWsHgmcnWHLH4cqK3zTKL5ZPnnXkTJeWL7DdRxzErM+3crBkkp+Nb6nbW4SwqzwTaMZnp7IiPREHluYy9HyKtdxTD0OFlfw1Kd5jO/TnoGd2riOYxqRT4UvIpeJyFoRqRGR7BMcN15ENopIrojc6cs5TXC5fVwWh0oqmfnpVtdRTD0e/3gLpZUefjku03UU08h8HeGvAS4BFtV3gIhEAo8B5wO9gStFpLeP5zVBYkCnNkzo156Zn+ZxoLjCdRxznF1Hypi9bDuXDkkjPbm16zimkflU+Kq6XlVPdh/9UCBXVfNUtRJ4CZjky3lNcLltbBblnhoeW5jrOoo5zj/e3wQCt4yx0X04aIo5/I7AzmO+zvc+VycRmSwiOSKSU1BQ0OjhTOPrkdSKy7PTmPPZdnL3F7uOY7xW5xfy2sp8rj2jCx3aNHcdxzSBkxa+iHwgImvq+DjVUXpdb/nXezeOqs5Q1WxVzU5KSjrFU5hAd9vYLFo0i+LuuattyYUAUF2j3P3GahJaxXDzGNvcJFycdBNzVR3j4znygU7HfJ0G7Pbxe5ogk9gqhjvP78mvX1/Nayt3cekQW6fFpdnLtvF1fiGPXDmIuFhb/jhcNMWUzhdAhoh0E5FmwBXAm01wXhNgfpDdiSFd2vKH+es5XFLpOk7Y2ldUzl/f28TZGYlM7G9LKIQTXy/LvFhE8oEzgX+LyALv8x1EZD6AqnqAm4AFwHrgFVVd61tsE4wiIoQHLu5LYVkVf3png+s4Yev+t9ZRWV3D77/X126yCjO+XqUzV1XTVDVGVVNUdZz3+d2qOuGY4+araqaq9lDVB3wNbYJXz/ZxXD+iGy/n7OSLbYdcxwk7Czfu59+r9zD13HS6JLR0Hcc0MbvT1jS5aWMy6NimOXfPXU2lp8Z1nLBRVlnNvfPW0COpJZO/2911HOOAFb5pci2aRXH/pD5s2lfMzMV5ruOEjUc+2szOQ2U8cHE/YqJsJ6twZIVvnBjdK4VxfVJ4+MPN7DxU6jpOyNu07ygzFuXx/cFpnNE9wXUc44gVvnHmdxf1IVKEe+etsWvzG1FNjXLP3DW0io3irgk9XccxDlnhG2dS45vzi/MyWbixgHfW7HUdJ2S9uiKfz7cd4tfn9yShVYzrOMYhK3zj1HVndaV3ahz3vbXWllBuBIdKKvnDO+v5Tte2XDak08n/gAlpVvjGqajICP5wST/2H63gofc2uY4Tcv4wfz3F5R4euLifbVtorPCNewM7teHqYV14fuR30l8AAAyaSURBVNk2VucXuo4TMj7LO8irK/L52TndyUyxpY+NFb4JELePzyKhVQx3v7Ga6hp7A9dXlZ4a7nljDWltm3PzKFsczdSywjcBIS42mt9M7M3X+YXMXrbNdZygN2PRFnL3F/O/k/rSvJldc29qWeGbgHFh/1TOzkjkr+9tYl9Rues4QWv7wRIe+SiXCf3ac27PZNdxTACxwjcBQ0T430l9qayu4f631rmOE5RUld/MW0t0ZAT3TuzjOo4JMFb4JqB0TWzJ1HPT+ffqPSzcuN91nKDz9td7WLSpgNvGZtI+PtZ1HBNgrPBNwJn83e70SGrJvfPWUFZZ7TpO0Cgqr+L+t9fRr2M8157Z1XUcE4Cs8E3AiYmK5Pff68fOQ2X8/QO7Nv9U/XH+Bg4WV/DAxX2JtGvuTR2s8E1AOrNHAlcN68yMRXnMX73HdZyA90rOTl78fAc/O7s7/dPauI5jApSvO15dJiJrRaRGRLJPcNw2EVktIqtEJMeXc5rwce+FvRncuQ23vfIV6/cUuY4TsFbuOMw9c9cwIj2R28dluY5jApivI/w1wCXAolM49lxVHaiq9f7FYMyxYqIimX71EOKaRzF5do7tg1uH/UXlTJm9gpT4GB65chBRkfZLu6mfr1scrlfVjf4KY8zxkuNimX71EPYVVnDTiyvxVNsOWd+o8FTz8zkrOFruYcY12bRt2cx1JBPgmmo4oMB7IrJCRCaf6EARmSwiOSKSU1BQ0ETxTCAb1Lktv7+4L0tyD/JH2/wcqL3e/t431vLljiM8dPkAeqXGuY5kgkDUyQ4QkQ+A9nW8dLeqzjvF8wxX1d0ikgy8LyIbVLXOaSBVnQHMAMjOzrZFVQwAl2d3Yt3uImYt3kqfDnFcMjjNdSSn5ny2nZdzdnLTuelM6JfqOo4JEictfFUd4+tJVHW39/N+EZkLDOXU5v2N+dbdF/Ri496j3Pn6atKTW4Xt1SjL8w5y31vrGNUzmVvPy3QdxwSRRp/SEZGWItL6m8fAWGrf7DWmQaIjI3j0h4NIahXDz2evoOBohetITW7XkTJueGElnRNa8I8rBtoa96ZBfL0s82IRyQfOBP4tIgu8z3cQkfnew1KAxSLyFfA58G9VfdeX85rwldAqhhnXDuFwaSU3vLCCSk/4vIlbVlnN5OdzqPTU8NS12cTFRruOZIKMr1fpzFXVNFWNUdUUVR3nfX63qk7wPs5T1QHejz6q+oA/gpvw1adDPA9eOoAvth3mvrfWuo7TJFSVO1//mnV7ivjnlQPpkdTKdSQThE46h29MILpwQAfW7i5i+idb6NMhnh8O6+w6UqN66tM85q3azS/HZjKqZ4rrOCZI2V0aJmjdPi6L72Ym8ds315Cz7ZDrOI1m0aYC/vTOBib0a8+N56a7jmOCmBW+CVqREcLDVwyiY5vmTJmzkj2FZa4j+d32gyVMffFLMlNa8+ClAxCxN2nN6bPCN0EtvkU0T12bTVmlhymzV1BeFTrLKRdXePjZ8zmIwIxrsmkZYzOwxjdW+CboZaS05u8/GMhX+YXcPXcNqsF/v15NjXLbK6vI3V/Mo1cOpnNCC9eRTAiwwjchYWyf9twyJoPXVuZzzxtrqPAE70i/rLKaX776FQvW7uOuCb0YkZHoOpIJEfY7ogkZN4/KoKyqmic/yWPt7iIev2owHdo0dx2rQbYfLGHKnJVs2FvELWMy+OmIbq4jmRBiI3wTMiIihF+f34vpVw8md38xEx9ZzJLcA65jnbIP1+9j4iOL2X2kjKev+w63jMm0N2mNX1nhm5Azvm8q824aTkLLZlwzazmPLcylpiZw5/Wra5S/LtjIT5/LoUtCC96eOoJzs5JdxzIhyArfhKQeSa1448bhXNC/Aw8u2Mjk2SsoLKtyHeu/HCqp5LpnPufRhblcnp3Gq1POolM7e4PWNA4rfBOyWsZE8fAVA/nthb35eON+Jj26mA17A2erxK92HuHCRxazfOsh/nRJP/5y6QBioyNdxzIhzArfhDQR4cfDu/HS5DMorazme48t4Y0vdznNpKr8a/kOLpu+DIBXp5zJFUNDe2kIExis8E1YyO7ajrdvHkH/tDbc8vIq7p23xslKm+VV1dz+6tfcNXc1Z/RI4O2pI8J2XX/T9OyyTBM2klvH8sL1w/jLuxt46tOtrN5VyONXDSY1vmku3dxxsJQpc1awbk8RN4/OYNroDCJtPXvThGyEb8JKdGQEd1/Qm8evGsymvUeZ+PBiPtlU0Kh359bUKO+t3cvERz4l/3ApT1+Xza3nZVrZmyZnI3wTlib0SyUzpTVT5qzgR09/TmKrGEakJ3BWeiIj0hN9vmFr56FSluQeYHHuAZZuOcihkkp6p8Yx/eohtkyCccanwheRB4ELgUpgC/BjVT1Sx3HjgX8CkcBMVf2TL+c1xh/Sk1sx78bh/Hv1nm/L+Y1VuwHontiS4emJDE9P5MweCcQ3P/HuUodLKlmWd5DFuQdYknuA7QdLAUiJi2FkVhIj0hOZ0C/VrsIxTokvv8qKyFjgI1X1iMifAVT1juOOiQQ2AecB+cAXwJWquu5k3z87O1tzcnJOO58xDaGqbNx3lMWba0t7+dZDlFZWEyHQL60NI9ITGJ6eyJAubVGFnG2Hvy34NbsLUYVWMVGc0T2BEekJjMhIpEdSK7tb1jQpEVmhqtl1vuavuUsRuRi4VFWvOu75M4HffbP9oYj8GkBV/3iy72mFb1yq9NSwaueRb0t91c4jVNcosdER1Gjt69GRwqDObRnh/W1gQFo8UZH21phx50SF7885/J8AL9fxfEdg5zFf5wPD6vsmIjIZmAzQubNdm2zcaRYVwdBu7RjarR23npfJ0fIqPt96iCW5B4kQGJ6RyNCu7WydehM0Tvpfqoh8ALSv46W7VXWe95i7AQ/wQl3foo7n6v21QlVnADOgdoR/snzGNJXWsdGM7pXC6F62p6wJTictfFUdc6LXReRHwERgtNY9P5QPdDrm6zRgd0NCGmOM8Z1Pk43eq2/uAC5S1dJ6DvsCyBCRbiLSDLgCeNOX8xpjjGk4X99dehRoDbwvIqtEZDqAiHQQkfkAquoBbgIWAOuBV1R1rY/nNcYY00A+vdukqun1PL8bmHDM1/OB+b6cyxhjjG/s+jFjjAkTVvjGGBMmrPCNMSZMWOEbY0yY8NvSCo1BRAqA7af5xxOBA36M09SCPT8E/88Q7Pkh+H8Gy99wXVQ1qa4XArrwfSEiOfWtJxEMgj0/BP/PEOz5Ifh/BsvvXzalY4wxYcIK3xhjwkQoF/4M1wF8FOz5Ifh/hmDPD8H/M1h+PwrZOXxjjDH/KZRH+MYYY45hhW+MMWEi5ApfRMaLyEYRyRWRO13naSgReVpE9ovIGtdZToeIdBKRhSKyXkTWisg015kaSkRiReRzEfnK+zPc5zrT6RCRSBH5UkTedp3ldIjINhFZ7V2JN+j2OhWRNiLyqohs8P7/cKbzTKE0h+/LhumBQkTOAYqB51W1r+s8DSUiqUCqqq4UkdbACuB7QfbvQICWqlosItHAYmCaqn7mOFqDiMitQDYQp6oTXedpKBHZBmSralDeeCUizwGfqupM714gLVT1iMtMoTbCHwrkqmqeqlYCLwGTHGdqEFVdBBxyneN0qeoeVV3pfXyU2j0QOrpN1TBaq9j7ZbT3I6hGRiKSBlwAzHSdJRyJSBxwDjALQFUrXZc9hF7h17VhelCVTSgRka7AIGC52yQN550OWQXsB95X1WD7Gf4B/AqocR3EBwq8JyIrRGSy6zAN1B0oAJ7xTqvNFJGWrkOFWuE3aMN003hEpBXwGnCLqha5ztNQqlqtqgOp3YN5qIgEzfSaiEwE9qvqCtdZfDRcVQcD5wM3eqc7g0UUMBh4QlUHASWA8/cUQ63wbcP0AOCd934NeEFVX3edxxfeX8M/BsY7jtIQw4GLvHPgLwGjRGSO20gN5905D1XdD8yldso2WOQD+cf8ZvgqtX8BOBVqhW8bpjvmfcNzFrBeVf/mOs/pEJEkEWnjfdwcGANscJvq1Knqr1U1TVW7Uvv/wEeqerXjWA0iIi29b/rjnQoZCwTNlWuquhfYKSJZ3qdGA84vXPBpT9tAo6oeEflmw/RI4Olg2zBdRF4ERgKJIpIP/FZVZ7lN1SDDgWuA1d45cIC7vPsaB4tU4DnvVV8RwCuqGpSXNgaxFGBu7fiBKOBfqvqu20gNNhV4wTv4zAN+7DhPaF2WaYwxpn6hNqVjjDGmHlb4xhgTJqzwjTEmTFjhG2NMmLDCN8aYMGGFb4wxYcIK3xhjwsT/A4h8Ue4ZlLYFAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "out.backward()\n",
    "print(a.grad)\n",
    "plt.plot(a.detach(), a.grad.detach())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Recall the computation steps we took to get here:\n",
    "\n",
    "```\n",
    "a = torch.linspace(0., 2. * math.pi, steps=25, requires_grad=True)\n",
    "b = torch.sin(a)\n",
    "c = 2 * b\n",
    "d = c + 1\n",
    "out = d.sum()\n",
    "```\n",
    "\n",
    "Adding a constant, as we did to compute `d`, does not change the derivative. That leaves $c = 2 * b = 2 * sin(a)$, the derivative of which should be $2 * cos(a)$. Looking at the graph above, that's just what we see.\n",
    "\n",
    "Be aware than only *leaf nodes* of the computation have their gradients computed. If you tried, for example, `print(c.grad)` you'd get back `None`. In this simple example, only the input is a leaf node, so only it has gradients computed.\n",
    "\n",
    "## Autograd in Training\n",
    "\n",
    "We've had a brief look at how autograd works, but how does it look when it's used for its intended purpose? Let's define a small model and examine how it changes after a single training batch. First, define a few constants, our model, and some stand-ins for inputs and outputs:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "BATCH_SIZE = 16\n",
    "DIM_IN = 1000\n",
    "HIDDEN_SIZE = 100\n",
    "DIM_OUT = 10\n",
    "\n",
    "class TinyModel(torch.nn.Module):\n",
    "\n",
    "    def __init__(self):\n",
    "        super(TinyModel, self).__init__()\n",
    "        \n",
    "        self.layer1 = torch.nn.Linear(1000, 100)\n",
    "        self.relu = torch.nn.ReLU()\n",
    "        self.layer2 = torch.nn.Linear(100, 10)\n",
    "    \n",
    "    def forward(self, x):\n",
    "        x = self.layer1(x)\n",
    "        x = self.relu(x)\n",
    "        x = self.layer2(x)\n",
    "        return x\n",
    "    \n",
    "some_input = torch.randn(BATCH_SIZE, DIM_IN, requires_grad=False)\n",
    "ideal_output = torch.randn(BATCH_SIZE, DIM_OUT, requires_grad=False)\n",
    "\n",
    "model = TinyModel()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "One thing you might notice is that we never specify `requires_grad=True` for the model's layers. Within a subclass of `torch.nn.module`, it's assumed that we want to track gradients on the layers' weights for learning.\n",
    "\n",
    "If we look at the layers of the model, we can examine the values of the weights, and verify that no gradients have been computed yet:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "tensor([ 0.0244,  0.0843,  0.0378, -0.0212, -0.0114,  0.0351, -0.0460,  0.0561,\n",
      "         0.0040,  0.0179], grad_fn=<SliceBackward>)\n",
      "None\n"
     ]
    }
   ],
   "source": [
    "print(model.layer2.weight[0][0:10]) # just a small slice\n",
    "print(model.layer2.weight.grad)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's see how this changes when we run through one training batch. For a loss function, we'll just use the square of the Euclidean distance between our `prediction` and the `ideal_output`, and we'll use a basic stochastic gradient descent optimizer."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "tensor(157.9841, grad_fn=<SumBackward0>)\n"
     ]
    }
   ],
   "source": [
    "optimizer = torch.optim.SGD(model.parameters(), lr=0.001)\n",
    "\n",
    "prediction = model(some_input)\n",
    "\n",
    "loss = (ideal_output - prediction).pow(2).sum()\n",
    "print(loss)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now, let's call `loss.backward()` and see what happens:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "tensor([ 0.0244,  0.0843,  0.0378, -0.0212, -0.0114,  0.0351, -0.0460,  0.0561,\n",
      "         0.0040,  0.0179], grad_fn=<SliceBackward>)\n",
      "tensor([ 2.8054,  1.4428,  4.5022,  4.4852, 10.5636,  8.2099,  7.4645,  1.9516,\n",
      "         1.9949,  3.5033])\n"
     ]
    }
   ],
   "source": [
    "loss.backward()\n",
    "print(model.layer2.weight[0][0:10])\n",
    "print(model.layer2.weight.grad[0][0:10])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can see that the gradients have been computed for each learning weight, but the weights remain unchanged, because we haven't run the optimizer yet. The optimizer is responsible for updating model weights based on the computed gradients."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "tensor([ 0.0216,  0.0829,  0.0332, -0.0257, -0.0220,  0.0269, -0.0535,  0.0541,\n",
      "         0.0020,  0.0144], grad_fn=<SliceBackward>)\n",
      "tensor([ 2.8054,  1.4428,  4.5022,  4.4852, 10.5636,  8.2099,  7.4645,  1.9516,\n",
      "         1.9949,  3.5033])\n"
     ]
    }
   ],
   "source": [
    "optimizer.step()\n",
    "print(model.layer2.weight[0][0:10])\n",
    "print(model.layer2.weight.grad[0][0:10])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "You should see that `layer2`'s weights have changed.\n",
    "\n",
    "One important thing about the process: After calling `optimizer.step()`, you need to call `optimizer.zero_grad()`, or else every time you run `loss.backward()`, the gradients on the learning weights will accumulate:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "tensor([ 2.8054,  1.4428,  4.5022,  4.4852, 10.5636,  8.2099,  7.4645,  1.9516,\n",
      "         1.9949,  3.5033])\n",
      "tensor([ 4.6037, -1.8794,  6.1310, 20.7784, 36.1072, 11.4850, 43.7955,  1.5995,\n",
      "         3.7064, -0.5014])\n",
      "tensor([0., 0., 0., 0., 0., 0., 0., 0., 0., 0.])\n"
     ]
    }
   ],
   "source": [
    "print(model.layer2.weight.grad[0][0:10])\n",
    "\n",
    "for i in range(0, 5):\n",
    "    prediction = model(some_input)\n",
    "    loss = (ideal_output - prediction).pow(2).sum()\n",
    "    loss.backward()\n",
    "    \n",
    "print(model.layer2.weight.grad[0][0:10])\n",
    "\n",
    "optimizer.zero_grad()\n",
    "\n",
    "print(model.layer2.weight.grad[0][0:10])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "After running the cell above, you should see that after running `loss.backward()` multiple times, the magnitudes of most of the gradients will be much larger. Failing to zero the gradients before running your next training batch will cause the gradients to blow up in this manner, causing incorrect and unpredictable learning results.\n",
    "\n",
    "## Turning Autograd Off and On\n",
    "\n",
    "There are situations where you will need fine-grained control over whether autograd is enabled. There are multiple ways to do this, depending on the situation.\n",
    "\n",
    "The simplest is to change the `requires_grad` flag on a tensor directly:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "tensor([[1., 1., 1.],\n",
      "        [1., 1., 1.]], requires_grad=True)\n",
      "tensor([[2., 2., 2.],\n",
      "        [2., 2., 2.]], grad_fn=<MulBackward0>)\n",
      "tensor([[2., 2., 2.],\n",
      "        [2., 2., 2.]])\n"
     ]
    }
   ],
   "source": [
    "a = torch.ones(2, 3, requires_grad=True)\n",
    "print(a)\n",
    "\n",
    "b1 = 2 * a\n",
    "print(b1)\n",
    "\n",
    "a.requires_grad = False\n",
    "b2 = 2 * a\n",
    "print(b2)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In the cell above, we see that `b1` has a `grad_fn` (i.e., a traced computation history), which is what we expect, since it was derived from a tensor, `a`, that had autograd turned on. When we turn off autograd explicitly with `a.requires_grad = False`, computation history is no longer tracked, as we see when we compute `b2`.\n",
    "\n",
    "If you only need autograd turned off temporarily, a better way is to use the `torch.no_grad()`:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "tensor([[5., 5., 5.],\n",
      "        [5., 5., 5.]], grad_fn=<AddBackward0>)\n",
      "tensor([[5., 5., 5.],\n",
      "        [5., 5., 5.]])\n",
      "tensor([[6., 6., 6.],\n",
      "        [6., 6., 6.]], grad_fn=<MulBackward0>)\n"
     ]
    }
   ],
   "source": [
    "a = torch.ones(2, 3, requires_grad=True) * 2\n",
    "b = torch.ones(2, 3, requires_grad=True) * 3\n",
    "\n",
    "c1 = a + b\n",
    "print(c1)\n",
    "\n",
    "with torch.no_grad():\n",
    "    c2 = a + b\n",
    "\n",
    "print(c2)\n",
    "\n",
    "c3 = a * b\n",
    "print(c3)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "`torch.no_grad()` can also be used as a function or method dectorator:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "tensor([[5., 5., 5.],\n",
      "        [5., 5., 5.]], grad_fn=<AddBackward0>)\n",
      "tensor([[5., 5., 5.],\n",
      "        [5., 5., 5.]])\n"
     ]
    }
   ],
   "source": [
    "def add_tensors1(x, y):\n",
    "    return x + y\n",
    "\n",
    "@torch.no_grad()\n",
    "def add_tensors2(x, y):\n",
    "    return x + y\n",
    "\n",
    "\n",
    "a = torch.ones(2, 3, requires_grad=True) * 2\n",
    "b = torch.ones(2, 3, requires_grad=True) * 3\n",
    "\n",
    "c1 = add_tensors1(a, b)\n",
    "print(c1)\n",
    "\n",
    "c2 = add_tensors2(a, b)\n",
    "print(c2)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "There's a corresponding context manager, `torch.enable_grad()`, for turning autograd on when it isn't already. It may also be used as a decorator.\n",
    "\n",
    "Finally, you may have a tensor that requires gradient tracking, but you want a copy that does not. For this we have the `Tensor` object's `detach()` method - it creates a copy of the tensor that is *detached* from the computation history:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "tensor([0.5443, 0.0012, 0.8133, 0.8857, 0.7629], requires_grad=True)\n",
      "tensor([0.5443, 0.0012, 0.8133, 0.8857, 0.7629])\n"
     ]
    }
   ],
   "source": [
    "x = torch.rand(5, requires_grad=True)\n",
    "y = x.detach()\n",
    "\n",
    "print(x)\n",
    "print(y)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We did this above when we wanted to graph some of our tensors. This is because `matplotlib` expects a NumPy array as input, and the implicit conversion from a PyTorch tensor to a NumPy array is not enabled for tensors with requires_grad=True. Making a detached copy lets us move forward.\n",
    "\n",
    "### Autograd and In-place Operations\n",
    "\n",
    "In every example in this notebook so far, we've used variables to capture the intermediate values of a computation. Autograd needs these intermediate values to perform gradient computations. *For this reason, you must be careful about using in-place operations when using autograd.* Doing so can destroy information you need to compute derivatives in the `backward()` call. PyTorch will even stop you if you attempt an in-place operation on leaf variable that requires autograd, as shown below.\n",
    "\n",
    "*Note: The following code cell throws a runtime error. This is expected.*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "ename": "RuntimeError",
     "evalue": "a leaf Variable that requires grad is being used in an in-place operation.",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mRuntimeError\u001b[0m                              Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-20-cffbbcd7d649>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m\u001b[0m\n\u001b[1;32m      1\u001b[0m \u001b[0ma\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mtorch\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mlinspace\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m0.\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;36m2.\u001b[0m \u001b[0;34m*\u001b[0m \u001b[0mmath\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpi\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0msteps\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;36m25\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mrequires_grad\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mTrue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 2\u001b[0;31m \u001b[0mtorch\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msin_\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0ma\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[0;31mRuntimeError\u001b[0m: a leaf Variable that requires grad is being used in an in-place operation."
     ]
    }
   ],
   "source": [
    "a = torch.linspace(0., 2. * math.pi, steps=25, requires_grad=True)\n",
    "torch.sin_(a)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Autograd Profiler\n",
    "\n",
    "Autograd tracks every step of your computation in detail. Such a computation history, combined with timing information, would make a handy profiler - and autograd has that feature baked in. Here's a quick example usage:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "-------  ---------------  ---------------  ---------------  ---------------  ---------------  ---------------  \n",
      "Name     Self CPU total %  Self CPU total   CPU total %      CPU total        CPU time avg     Number of Calls  \n",
      "-------  ---------------  ---------------  ---------------  ---------------  ---------------  ---------------  \n",
      "mul      51.18%           7.985ms          51.18%           7.985ms          7.985us          1000             \n",
      "div      48.82%           7.616ms          48.82%           7.616ms          7.616us          1000             \n",
      "-------  ---------------  ---------------  ---------------  ---------------  ---------------  ---------------  \n",
      "Self CPU time total: 15.601ms\n",
      "\n"
     ]
    }
   ],
   "source": [
    "device = torch.device('cpu')\n",
    "run_on_gpu = False\n",
    "if torch.cuda.is_available():\n",
    "    device = torch.device('cuda')\n",
    "    run_on_gpu = True\n",
    "    \n",
    "x = torch.randn(2, 3, requires_grad=True)\n",
    "y = torch.rand(2, 3, requires_grad=True)\n",
    "z = torch.ones(2, 3, requires_grad=True)\n",
    "\n",
    "with torch.autograd.profiler.profile(use_cuda=run_on_gpu) as prf:\n",
    "    for _ in range(1000):\n",
    "        z = (z / x) * y\n",
    "        \n",
    "print(prf.key_averages().table(sort_by='self_cpu_time_total'))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The profiler can also label individual sub-blocks of code, break out the data by input tensor shape, and export data as a Chrome tracing tools file. For full details of the API, see the [documentation](https://pytorch.org/docs/stable/autograd.html#profiler)\n",
    "\n",
    "## Advanced Topic: More Autograd Detail and the High-Level API\n",
    "\n",
    "If you have a function with an n-dimensional input and m-dimensional output, $\\vec{y}=f(\\vec{x})$, the complete gradient is a matrix of the derivative of every output with respect to every input, called the *Jacobian:*\n",
    "\n",
    "\\begin{align}J=\\left(\\begin{array}{ccc}\n",
    "   \\frac{\\partial y_{1}}{\\partial x_{1}} & \\cdots & \\frac{\\partial y_{1}}{\\partial x_{n}}\\\\\n",
    "   \\vdots & \\ddots & \\vdots\\\\\n",
    "   \\frac{\\partial y_{m}}{\\partial x_{1}} & \\cdots & \\frac{\\partial y_{m}}{\\partial x_{n}}\n",
    "   \\end{array}\\right)\\end{align}\n",
    "\n",
    "If you have a second function, $l=g\\left(\\vec{y}\\right)$ that takes m-dimensional input (that is, the same dimensionality as the output above), and returns a scalar output, you can express its gradients with respect to $\\vec{y}$ as a column vector, $v=\\left(\\begin{array}{ccc}\\frac{\\partial l}{\\partial y_{1}} & \\cdots & \\frac{\\partial l}{\\partial y_{m}}\\end{array}\\right)^{T}$ - which is really just a one-column Jacobian.\n",
    "\n",
    "More concretely, imagine the first function as your PyTorch model (with potentially many inputs and many outputs) and the second function as a loss function (with the model's output as input, and the loss value as the scalar output).\n",
    "\n",
    "If we multiply the first function's Jacobian by the gradient of the second function, and apply the chain rule, we get:\n",
    "\n",
    "\\begin{align}J^{T}\\cdot v=\\left(\\begin{array}{ccc}\n",
    "   \\frac{\\partial y_{1}}{\\partial x_{1}} & \\cdots & \\frac{\\partial y_{m}}{\\partial x_{1}}\\\\\n",
    "   \\vdots & \\ddots & \\vdots\\\\\n",
    "   \\frac{\\partial y_{1}}{\\partial x_{n}} & \\cdots & \\frac{\\partial y_{m}}{\\partial x_{n}}\n",
    "   \\end{array}\\right)\\left(\\begin{array}{c}\n",
    "   \\frac{\\partial l}{\\partial y_{1}}\\\\\n",
    "   \\vdots\\\\\n",
    "   \\frac{\\partial l}{\\partial y_{m}}\n",
    "   \\end{array}\\right)=\\left(\\begin{array}{c}\n",
    "   \\frac{\\partial l}{\\partial x_{1}}\\\\\n",
    "   \\vdots\\\\\n",
    "   \\frac{\\partial l}{\\partial x_{n}}\n",
    "   \\end{array}\\right)\\end{align}\n",
    "\n",
    "*Note: You could also use the equivalent operation $v^{T}\\cdot J$, and get back a row vector.*\n",
    "\n",
    "The resulting column vector is the *gradient of the second function with respect to the inputs of the first* - or in the case of our model and loss function, the gradient of the loss with respect to the model inputs.\n",
    "\n",
    "**`torch.autograd` is an engine for computing these products.** This is how we accumulate the gradients over the learning weights during the backward pass.\n",
    "\n",
    "For this reason, the `backward()` call can *also* take an optional vector input. This vector represents a set of gradients over the tensor, which are multiplied by the Jacobian of the autograd-traced tensor that precedes it. Let's try a specific example with a small vector:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "tensor([-663.2701, 1493.5994,  789.5186], grad_fn=<MulBackward0>)\n"
     ]
    }
   ],
   "source": [
    "x = torch.randn(3, requires_grad=True)\n",
    "\n",
    "y = x * 2\n",
    "while y.data.norm() < 1000:\n",
    "    y = y * 2\n",
    "\n",
    "print(y)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "If we tried to call `y.backward()` now, we'd get a runtime error and a message that gradients can only be *implicitly* computed for scalar outputs. For a multi-dimensional output, autograd expects us to provide gradients for those three outputs that it can multiply into the Jacobian:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "tensor([1.0240e+02, 1.0240e+03, 1.0240e-01])\n"
     ]
    }
   ],
   "source": [
    "v = torch.tensor([0.1, 1.0, 0.0001], dtype=torch.float) # stand-in for gradients\n",
    "y.backward(v)\n",
    "\n",
    "print(x.grad)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "(Note that the output gradients are all related to powers of two - which we'd expect from a repeated doubling operation.)\n",
    "\n",
    "### The High-Level API\n",
    "\n",
    "There is an API on autograd that gives you direct access to important differential matrix and vector operations. In particular, it allows you to calculate the Jacobian and the *Hessian* matrices of a particular function for particular inputs. (The Hessian is like the Jacobian, but expresses all partial *second* derivatives.) It also provides methods for taking vector products with these matrices.\n",
    "\n",
    "Let's take the Jacobian of a simple function, evaluated for a 2 single-element inputs:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(tensor([0.7327]), tensor([0.8072]))\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "(tensor([[4.1616]]), tensor([[3.]]))"
      ]
     },
     "execution_count": 24,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def exp_adder(x, y):\n",
    "    return 2 * x.exp() + 3 * y\n",
    "\n",
    "inputs = (torch.rand(1), torch.rand(1)) # arguments for the function\n",
    "print(inputs)\n",
    "torch.autograd.functional.jacobian(exp_adder, inputs)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "If you look closely, the first output should equal $2e^x$ (since the derivative of $e^x$ is $e^x$), and the second value should be 3.\n",
    "\n",
    "You can, of course, do this with higher-order tensors:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(tensor([0.3723, 0.7282, 0.0756]), tensor([0.9602, 0.6986, 0.2260]))\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "(tensor([[2.9022, 0.0000, 0.0000],\n",
       "         [0.0000, 4.1428, 0.0000],\n",
       "         [0.0000, 0.0000, 2.1570]]),\n",
       " tensor([[3., 0., 0.],\n",
       "         [0., 3., 0.],\n",
       "         [0., 0., 3.]]))"
      ]
     },
     "execution_count": 25,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "inputs = (torch.rand(3), torch.rand(3)) # arguments for the function\n",
    "print(inputs)\n",
    "torch.autograd.functional.jacobian(exp_adder, inputs)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The `torch.autograd.functional.hessian()` method works identically (assuming your function is twice differentiable), but returns a matrix of all second derivatives.\n",
    "\n",
    "There is also a function to directly compute the vector-Jacobian product, if you provide the vector:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(tensor([ -52.9634, -894.0098, 1167.3909]),\n",
       " tensor([5.1200e+01, 5.1200e+02, 5.1200e-02]))"
      ]
     },
     "execution_count": 26,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def do_some_doubling(x):\n",
    "    y = x * 2\n",
    "    while y.data.norm() < 1000:\n",
    "        y = y * 2\n",
    "    return y\n",
    "\n",
    "inputs = torch.randn(3)\n",
    "my_gradients = torch.tensor([0.1, 1.0, 0.0001])\n",
    "torch.autograd.functional.vjp(do_some_doubling, inputs, v=my_gradients)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The `torch.autograd.functional.jvp()` method performs the same matrix multiplication as `vjp()` with the operands reversed. The `vhp()` and `hvp()` methods do the same for a vector-Hessian product.\n",
    "\n",
    "For more information, including preformance notes on the [docs for the functional API](https://pytorch.org/docs/stable/autograd.html#functional-higher-level-api)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Advanced Topic: Higher-Order Derivatives\n",
    "+ computing your own chain-rule derivatives with the high-level API"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.8.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}