CNN_Autoencoder.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Convolutional Autoencoder\n",
    "\n",
    "Sticking with the MNIST dataset, let's improve our autoencoder's performance using convolutional layers. We'll build a convolutional autoencoder to compress the MNIST dataset. \n",
    "\n",
    ">The encoder portion will be made of convolutional and pooling layers and the decoder will be made of **transpose convolutional layers** that learn to \"upsample\" a compressed representation.\n",
    "\n",
    "<img src='notebook_ims/autoencoder_1.png' />\n",
    "\n",
    "### Compressed Representation\n",
    "\n",
    "A compressed representation can be great for saving and sharing any kind of data in a way that is more efficient than storing raw data. In practice, the compressed representation often holds key information about an input image and we can use it for denoising images or other kinds of reconstruction and transformation!\n",
    "\n",
    "<img src='notebook_ims/denoising.png' width=60%/>\n",
    "\n",
    "Let's get started by importing our libraries and getting the dataset."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import torch\n",
    "import numpy as np\n",
    "from torchvision import datasets\n",
    "import torchvision.transforms as transforms\n",
    "\n",
    "# convert data to torch.FloatTensor\n",
    "transform = transforms.ToTensor()\n",
    "\n",
    "# load the training and test datasets\n",
    "train_data = datasets.MNIST(root='data', train=True,\n",
    "                                   download=True, transform=transform)\n",
    "test_data = datasets.MNIST(root='data', train=False,\n",
    "                                  download=True, transform=transform)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Create training and test dataloaders\n",
    "\n",
    "num_workers = 0\n",
    "# how many samples per batch to load\n",
    "batch_size = 20\n",
    "\n",
    "# prepare data loaders\n",
    "train_loader = torch.utils.data.DataLoader(train_data, batch_size=batch_size, num_workers=num_workers)\n",
    "test_loader = torch.utils.data.DataLoader(test_data, batch_size=batch_size, num_workers=num_workers)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Visualize the Data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.image.AxesImage at 0x24064848548>"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAATEAAAEvCAYAAAAtufaDAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAPk0lEQVR4nO3db4xVdX7H8c+nqA9EFMhWJKyW1RgsGjs2iI2aqjGsf6LRUbdZEjc0GvGBJJhsSA1PVh9gSFW2IRoDG3HR7LJu4lrRNFUjKG1siAOiItRqDOuCE4gigvgvMN8+mGMy4Aznx7135swX3q+E3Ht/8+V3v8fDfDzn3N+ccUQIALL6q6YbAIB2EGIAUiPEAKRGiAFIjRADkBohBiC1E0byzWyzngNAqz6NiL8+fLCtIzHb19p+3/aHtu9rZy4AqPHnwQZbDjHbYyQ9Juk6SdMlzbY9vdX5AKAV7RyJzZT0YUR8FBHfSfqDpJs60xYAlGknxKZI+suA19urMQAYMe1c2PcgYz+4cG97rqS5bbwPAAypnRDbLunMAa9/LOmTw4siYrmk5RKfTgLovHZOJ9+UdK7tn9g+SdLPJa3uTFsAUKblI7GIOGB7nqSXJI2RtCIi3utYZwBQwCN5PzFOJwG0YUNEzDh8kB87ApAaIQYgNUIMQGqEGIDUCDEAqRFiAFIjxACkRogBSI0QA5AaIQYgNUIMQGqEGIDUCDEAqRFiAFIjxACkRogBSI0QA5AaIQYgNUIMQGqEGIDUCDEAqRFiAFIjxACkRogBSI0QA5AaIQYgNUIMQGqEGIDUCDEAqRFiAFIjxACkRogBSI0QA5AaIQYgNUIMQGqEGIDUCDEAqRFiAFIjxACkdkLTDSC3MWPG1NacdtppI9DJoebNm1dUd/LJJxfVTZs2rajunnvuqa15+OGHi+aaPXt2Ud0333xTW7N48eKiuR544IGiutGkrRCzvU3SPkkHJR2IiBmdaAoASnXiSOyqiPi0A/MAwFHjmhiA1NoNsZD0su0NtucOVmB7ru0e2z1tvhcA/EC7p5OXRcQntk+X9Irt/42IdQMLImK5pOWSZDvafD8AOERbR2IR8Un1uEvSc5JmdqIpACjVcojZHmt73PfPJf1U0uZONQYAJdo5nZwk6Tnb38/z+4j4z450BQCFWg6xiPhI0t91sBcM4ayzzqqtOemkk4rmuvTSS4vqLr/88qK68ePH19bceuutRXONZtu3by+qW7p0aW1Nd3d30Vz79u0rqnv77bdra15//fWiuTJiiQWA1AgxAKkRYgBSI8QApEaIAUiNEAOQGiEGIDVCDEBqhBiA1BwxcjeW4C4Wh+rq6iqqW7NmTW1NE7eAPhb09fUV1d1xxx1FdV9++WU77Ryit7e3qO7zzz+vrXn//ffbbWc02DDY3aM5EgOQGiEGIDVCDEBqhBiA1AgxAKkRYgBSI8QApEaIAUiNEAOQWru/dxJt+Pjjj4vqPvvss9qaY2HF/vr164vq9uzZU1tz1VVXFc313XffFdU9/fTTRXUYeRyJAUiNEAOQGiEGIDVCDEBqhBiA1AgxAKkRYgBSI8QApMZi1wbt3r27qG7BggW1NTfccEPRXG+99VZR3dKlS4vqSmzatKmobtasWUV1+/fvr605//zzi+aaP39+UR1GL47EAKRGiAFIjRADkBohBiA1QgxAaoQYgNQIMQCpEWIAUiPEAKTmiBi5N7NH7s2OM6eeempR3b59+4rqli1bVlR355131tbcfvvtRXOtWrWqqA7HrQ0RMePwwdojMdsrbO+yvXnA2ETbr9j+oHqc0OluAaBEyenkbyVde9jYfZJejYhzJb1avQaAEVcbYhGxTtLhP6l8k6SV1fOVkm7ucF8AUKTVC/uTIqJXkqrH0zvXEgCUG/Zb8dieK2nucL8PgONTq0diO21PlqTqcddQhRGxPCJmDPapAgC0q9UQWy1pTvV8jqTnO9MOABydkiUWqyT9j6RptrfbvlPSYkmzbH8gaVb1GgBGXO01sYiYPcSXru5wLwBw1LjH/jFi7969HZ3viy++6Nhcd911V1HdM888U1TX19fXTjs4xvCzkwBSI8QApEaIAUiNEAOQGiEGIDVCDEBqhBiA1AgxAKkRYgBS4x77GNTYsWOL6l544YXamiuuuKJoruuuu66o7uWXXy6qwzGntXvsA8BoRogBSI0QA5AaIQYgNUIMQGqEGIDUCDEAqRFiAFJjsSvacs4559TWbNy4sWiuPXv2FNWtXbu2tqanp6dorscee6yobiS/TzAkFrsCOPYQYgBSI8QApEaIAUiNEAOQGiEGIDVCDEBqhBiA1AgxAKmxYh/Drru7u6juySefLKobN25cO+0cYuHChUV1Tz31VFFdb29vO+3gyFixD+DYQ4gBSI0QA5AaIQYgNUIMQGqEGIDUCDEAqRFiAFIjxACkxop9jBoXXHBBUd2SJUtqa66++up22znEsmXLiuoWLVpUW7Njx4522zletbZi3/YK27tsbx4wdr/tHbY3VX+u73S3AFCi5HTyt5KuHWT81xHRVf35j862BQBlakMsItZJ2j0CvQDAUWvnwv482+9Up5sThiqyPdd2j+2yXwQIAEeh1RB7XNI5krok9Up6ZKjCiFgeETMGuyAHAO1qKcQiYmdEHIyIPkm/kTSzs20BQJmWQsz25AEvuyVtHqoWAIbTCXUFtldJulLSj2xvl/QrSVfa7pIUkrZJunsYewSAIbHYFemMHz++tubGG28smqv0lti2i+rWrFlTWzNr1qyiufAD3J4awLGHEAOQGiEGIDVCDEBqhBiA1AgxAKkRYgBSI8QApEaIAUiNFfs4rn377bdFdSecUPsTepKkAwcO1NZcc801RXO99tprRXXHEVbsAzj2EGIAUiPEAKRGiAFIjRADkBohBiA1QgxAaoQYgNQIMQCplS1DBkbAhRdeWFR322231dZcfPHFRXOVrsQvtWXLltqadevWdfQ9j3cciQFIjRADkBohBiA1QgxAaoQYgNQIMQCpEWIAUiPEAKRGiAFIjRX7aMu0adNqa+bNm1c01y233FJUd8YZZxTVddLBgweL6np7e2tr+vr62m0HA3AkBiA1QgxAaoQYgNQIMQCpEWIAUiPEAKRGiAFIjRADkBqLXY8zpQtFZ8+eXVRXspB16tSpRXM1oaenp6hu0aJFRXWrV69upx20oPZIzPaZttfa3mr7Pdvzq/GJtl+x/UH1OGH42wWAQ5WcTh6Q9MuI+FtJ/yDpHtvTJd0n6dWIOFfSq9VrABhRtSEWEb0RsbF6vk/SVklTJN0kaWVVtlLSzcPVJAAM5agu7NueKukiSeslTYqIXqk/6CSd3unmAKBO8YV926dIelbSvRGx13bp35sraW5r7QHAkRUdidk+Uf0B9ruI+FM1vNP25OrrkyXtGuzvRsTyiJgRETM60TAADFTy6aQlPSFpa0QsGfCl1ZLmVM/nSHq+8+0BwJGVnE5eJukXkt61vakaWyhpsaQ/2r5T0seSfjY8LQLA0GpDLCL+W9JQF8Cu7mw7AHB0WLGfwKRJk2prpk+fXjTXo48+WlR33nnnFdU1Yf369bU1Dz30UNFczz9fdhWEW0qPXvzsJIDUCDEAqRFiAFIjxACkRogBSI0QA5AaIQYgNUIMQGqEGIDUWLE/DCZOnFhUt2zZsqK6rq6u2pqzzz67aK4mvPHGG0V1jzzySFHdSy+9VFvz9ddfF82F/DgSA5AaIQYgNUIMQGqEGIDUCDEAqRFiAFIjxACkRogBSI3FrpVLLrmkqG7BggW1NTNnziyaa8qUKUV1Tfjqq6+K6pYuXVpb8+CDDxbNtX///qI6YCCOxACkRogBSI0QA5AaIQYgNUIMQGqEGIDUCDEAqRFiAFIjxACkxor9Snd3d0frOmnLli21NS+++GLRXAcOHCiqK71V9J49e4rqgOHCkRiA1AgxAKkRYgBSI8QApEaIAUiNEAOQGiEGIDVCDEBqhBiA1BwRI/dm9si9GYBjzYaImHH4YO2RmO0zba+1vdX2e7bnV+P3295he1P15/rh6BoAjqTkZycPSPplRGy0PU7SBtuvVF/7dUQ8PHztAcCR1YZYRPRK6q2e77O9VdLo/V1jAI4rR3Vh3/ZUSRdJWl8NzbP9ju0Vtid0uDcAqFUcYrZPkfSspHsjYq+kxyWdI6lL/Udqg967xfZc2z22ezrQLwAcoujTSdsnSnpR0ksRsWSQr0+V9GJEXFAzD59OAmhVy59OWtITkrYODDDbkweUdUva3IkuAeBolHw6eZmkX0h61/amamyhpNm2uySFpG2S7h6WDgHgCFjsCiCL1k4nAWA0I8QApEaIAUiNEAOQGiEGIDVCDEBqhBiA1AgxAKkRYgBSI8QApEaIAUiNEAOQGiEGIDVCDEBqhBiA1AgxAKkRYgBSI8QApEaIAUit5BeFdNKnkv582NiPqvGssvcv5d+G7P1L+bdhJPr/m8EGR/QXhQzagN0z2M3/s8jev5R/G7L3L+Xfhib753QSQGqEGIDURkOILW+6gTZl71/Kvw3Z+5fyb0Nj/Td+TQwA2jEajsQAoGWNhZjta22/b/tD2/c11Uc7bG+z/a7tTbZ7mu6nhO0VtnfZ3jxgbKLtV2x/UD1OaLLHIxmi//tt76j2wybb1zfZ45HYPtP2Wttbbb9ne341nmkfDLUNjeyHRk4nbY+R9H+SZknaLulNSbMjYsuIN9MG29skzYiINOt7bP+jpC8lPRURF1Rj/yppd0Qsrv6HMiEi/qXJPocyRP/3S/oyIh5usrcStidLmhwRG22Pk7RB0s2S/ll59sFQ2/BPamA/NHUkNlPShxHxUUR8J+kPkm5qqJfjSkSsk7T7sOGbJK2snq9U/z/IUWmI/tOIiN6I2Fg93ydpq6QpyrUPhtqGRjQVYlMk/WXA6+1q8D9CG0LSy7Y32J7bdDNtmBQRvVL/P1BJpzfcTyvm2X6nOt0ctadiA9meKukiSeuVdB8ctg1SA/uhqRDzIGMZPya9LCL+XtJ1ku6pTnUw8h6XdI6kLkm9kh5ptp16tk+R9KykeyNib9P9tGKQbWhkPzQVYtslnTng9Y8lfdJQLy2LiE+qx12SnlP/aXJGO6vrHN9f79jVcD9HJSJ2RsTBiOiT9BuN8v1g+0T1f/P/LiL+VA2n2geDbUNT+6GpEHtT0rm2f2L7JEk/l7S6oV5aYntsdVFTtsdK+qmkzUf+W6PWaklzqudzJD3fYC9H7ftv/kq3RvF+sG1JT0jaGhFLBnwpzT4Yahua2g+NLXatPn79N0ljJK2IiEWNNNIi22er/+hL6r8byO8zbIPtVZKuVP9dB3ZK+pWkf5f0R0lnSfpY0s8iYlRePB+i/yvVfwoTkrZJuvv760ujje3LJf2XpHcl9VXDC9V/TSnLPhhqG2argf3Ain0AqbFiH0BqhBiA1AgxAKkRYgBSI8QApEaIAUiNEAOQGiEGILX/BwIYAbUIKiJFAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 360x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline\n",
    "    \n",
    "# obtain one batch of training images\n",
    "dataiter = iter(train_loader)\n",
    "images, labels = dataiter.next()\n",
    "images = images.numpy()\n",
    "\n",
    "# get one image from the batch\n",
    "img = np.squeeze(images[0])\n",
    "\n",
    "fig = plt.figure(figsize = (5,5)) \n",
    "ax = fig.add_subplot(111)\n",
    "ax.imshow(img, cmap='gray')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---\n",
    "## Convolutional  Autoencoder\n",
    "\n",
    "#### Encoder\n",
    "The encoder part of the network will be a typical convolutional pyramid. Each convolutional layer will be followed by a max-pooling layer to reduce the dimensions of the layers. \n",
    "\n",
    "#### Decoder\n",
    "\n",
    "The decoder though might be something new to you. The decoder needs to convert from a narrow representation to a wide, reconstructed image. For example, the representation could be a 7x7x4 max-pool layer. This is the output of the encoder, but also the input to the decoder. We want to get a 28x28x1 image out from the decoder so we need to work our way back up from the compressed representation. A schematic of the network is shown below.\n",
    "\n",
    "<img src='notebook_ims/conv_enc_1.png' width=640px>\n",
    "\n",
    "Here our final encoder layer has size 7x7x4 = 196. The original images have size 28x28 = 784, so the encoded vector is 25% the size of the original image. These are just suggested sizes for each of the layers. Feel free to change the depths and sizes, in fact, you're encouraged to add additional layers to make this representation even smaller! Remember our goal here is to find a small representation of the input data.\n",
    "\n",
    "### Transpose Convolutions, Decoder\n",
    "\n",
    "This decoder uses **transposed convolutional** layers to increase the width and height of the input layers. They work almost exactly the same as convolutional layers, but in reverse. A stride in the input layer results in a larger stride in the transposed convolution layer. For example, if you have a 3x3 kernel, a 3x3 patch in the input layer will be reduced to one unit in a convolutional layer. Comparatively, one unit in the input layer will be expanded to a 3x3 path in a transposed convolution layer. PyTorch provides us with an easy way to create the layers, [`nn.ConvTranspose2d`](https://pytorch.org/docs/stable/nn.html#convtranspose2d). \n",
    "\n",
    "It is important to note that transpose convolution layers can lead to artifacts in the final images, such as checkerboard patterns. This is due to overlap in the kernels which can be avoided by setting the stride and kernel size equal. In [this Distill article](http://distill.pub/2016/deconv-checkerboard/) from Augustus Odena, *et al*, the authors show that these checkerboard artifacts can be avoided by resizing the layers using nearest neighbor or bilinear interpolation (upsampling) followed by a convolutional layer. \n",
    "\n",
    "> We'll show this approach in another notebook, so you can experiment with it and see the difference.\n",
    "\n",
    "\n",
    "#### TODO: Build the network shown above. \n",
    "> Build the encoder out of a series of convolutional and pooling layers. \n",
    "> When building the decoder, recall that transpose convolutional layers can upsample an input by a factor of 2 using a stride and kernel_size of 2. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "ConvAutoencoder(\n",
      "  (conv1): Conv2d(1, 16, kernel_size=(3, 3), stride=(1, 1), padding=(1, 1))\n",
      "  (conv2): Conv2d(16, 4, kernel_size=(3, 3), stride=(1, 1), padding=(1, 1))\n",
      "  (pool): MaxPool2d(kernel_size=2, stride=2, padding=0, dilation=1, ceil_mode=False)\n",
      "  (t_conv1): ConvTranspose2d(4, 16, kernel_size=(2, 2), stride=(2, 2))\n",
      "  (t_conv2): ConvTranspose2d(16, 1, kernel_size=(2, 2), stride=(2, 2))\n",
      ")\n"
     ]
    }
   ],
   "source": [
    "import torch.nn as nn\n",
    "import torch.nn.functional as F\n",
    "\n",
    "# define the NN architecture\n",
    "class ConvAutoencoder(nn.Module):\n",
    "    def __init__(self):\n",
    "        super(ConvAutoencoder, self).__init__()\n",
    "        ## encoder layers ##\n",
    "        # (increase depth from 1 to 16)\n",
    "        self.conv1 = nn.Conv2d(1, 16, 3, padding=1)\n",
    "        self.conv2 = nn.Conv2d(16, 4, 3, padding=1)\n",
    "        # (decrease size by half of length and width)\n",
    "        self.pool = nn.MaxPool2d(2, 2)\n",
    "        \n",
    "        ## decoder layers ##\n",
    "        ## a kernel of 2 and a stride of 2 will increase the spatial dims by 2\n",
    "        self.t_conv1 = nn.ConvTranspose2d(4, 16, 2, stride=2)\n",
    "        self.t_conv2 = nn.ConvTranspose2d(16, 1, 2, stride=2)\n",
    "\n",
    "    def forward(self, x):\n",
    "        ## encode ##\n",
    "        x = F.relu(self.conv1(x))\n",
    "        x = self.pool(x)\n",
    "        x = F.relu(self.conv2(x))\n",
    "        x = self.pool(x)\n",
    "        \n",
    "        ## decode ##\n",
    "        ## apply ReLu to all hidden layers *except for the output layer\n",
    "        ## apply a sigmoid to the output layer\n",
    "        x = F.relu(self.t_conv1(x))\n",
    "        x = F.sigmoid(self.t_conv2(x))\n",
    "                \n",
    "        return x\n",
    "\n",
    "# initialize the NN\n",
    "model = ConvAutoencoder()\n",
    "print(model)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---\n",
    "## Training\n",
    "\n",
    "Here I'll write a bit of code to train the network. I'm not too interested in validation here, so I'll just monitor the training loss and the test loss afterwards. \n",
    "\n",
    "We are not concerned with labels in this case, just images, which we can get from the `train_loader`. Because we're comparing pixel values in input and output images, it will be best to use a loss that is meant for a regression task. Regression is all about comparing quantities rather than probabilistic values. So, in this case, I'll use `MSELoss`. And compare output images and input images as follows:\n",
    "```\n",
    "loss = criterion(outputs, images)\n",
    "```\n",
    "\n",
    "Otherwise, this is pretty straightfoward training with PyTorch. Since this is a convlutional autoencoder, our images _do not_ need to be flattened before being passed in an input to our model."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "# specify loss function\n",
    "criterion = nn.MSELoss()\n",
    "\n",
    "# specify loss function\n",
    "optimizer = torch.optim.Adam(model.parameters(), lr=0.001)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "C:\\Users\\di.wei.tan\\AppData\\Local\\Continuum\\anaconda3\\envs\\deep-learning\\lib\\site-packages\\torch\\nn\\functional.py:1351: UserWarning: nn.functional.sigmoid is deprecated. Use torch.sigmoid instead.\n",
      "  warnings.warn(\"nn.functional.sigmoid is deprecated. Use torch.sigmoid instead.\")\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Epoch: 1 \tTraining Loss: 0.469362\n",
      "Epoch: 2 \tTraining Loss: 0.218012\n",
      "Epoch: 3 \tTraining Loss: 0.198438\n",
      "Epoch: 4 \tTraining Loss: 0.190665\n",
      "Epoch: 5 \tTraining Loss: 0.186450\n",
      "Epoch: 6 \tTraining Loss: 0.183504\n",
      "Epoch: 7 \tTraining Loss: 0.181051\n",
      "Epoch: 8 \tTraining Loss: 0.178758\n",
      "Epoch: 9 \tTraining Loss: 0.176872\n",
      "Epoch: 10 \tTraining Loss: 0.175204\n"
     ]
    }
   ],
   "source": [
    "# number of epochs to train the model\n",
    "n_epochs = 10\n",
    "\n",
    "for epoch in range(1, n_epochs+1):\n",
    "    # monitor training loss\n",
    "    train_loss = 0.0\n",
    "    \n",
    "    ###################\n",
    "    # train the model #\n",
    "    ###################\n",
    "    for data in train_loader:\n",
    "        # _ stands in for labels, here\n",
    "        # no need to flatten images\n",
    "        images, _ = data\n",
    "        # clear the gradients of all optimized variables\n",
    "        optimizer.zero_grad()\n",
    "        # forward pass: compute predicted outputs by passing inputs to the model\n",
    "        outputs = model(images)\n",
    "        # calculate the loss\n",
    "        loss = criterion(outputs, images)\n",
    "        # backward pass: compute gradient of the loss with respect to model parameters\n",
    "        loss.backward()\n",
    "        # perform a single optimization step (parameter update)\n",
    "        optimizer.step()\n",
    "        # update running training loss\n",
    "        train_loss += loss.item()*images.size(0)\n",
    "            \n",
    "    # print avg training statistics \n",
    "    train_loss = train_loss/len(train_loader)\n",
    "    print('Epoch: {} \\tTraining Loss: {:.6f}'.format(\n",
    "        epoch, \n",
    "        train_loss\n",
    "        ))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Checking out the results\n",
    "\n",
    "Below I've plotted some of the test images along with their reconstructions. These look a little rough around the edges, likely due to the checkerboard effect we mentioned above that tends to happen with transpose layers."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABXEAAADrCAYAAAAv1NW3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3dd8AU1fX4/4uKCNKbogKCBATBBmIvKMYKQaUYiTEKllgw2CPYsEUxlqgRjR2xgYjGggV7QQMKirSAAiK9SJMmPt8/fr/czzlH5j4zw+4+s8/zfv11Lmdn5+7Ofe7MDnvPViopKXEAAAAAAAAAgGzaqqw7AAAAAAAAAACIxk1cAAAAAAAAAMgwbuICAAAAAAAAQIZxExcAAAAAAAAAMoybuAAAAAAAAACQYdzEBQAAAAAAAIAM2ybJgytVqlSSr44gsSUlJSUNyroTcTBusqOkpKRSWfchDsZMpjDXIA3GDdJg3CANxg3SYNwgDcYNEuMzOFKInGv4Jm7xml3WHQBQITDXIA3GDdJg3CANxg3SYNwgDcYNgEKInGu4iQsAAAAAAAAAGcZNXAAAAAAAAADIMG7iAgAAAAAAAECGcRMXAAAAAAAAADKMm7gAAAAAAAAAkGHcxAUAAAAAAACADOMmLgAAAAAAAABkGDdxAQAAAAAAACDDuIkLAAAAAAAAABm2TVl3ACikyy67TLWrVq3q4z333FPlunfvHvk8DzzwgGp/+umnPh46dOiWdBEAAAAAAABQ+CYuAAAAAAAAAGQYN3EBAAAAAAAAIMMop4By77nnnvNxqESC9csvv0Tmzj33XNXu3Lmzj99//32VmzNnTux9ouJo2bKlj6dOnapyF198sY/vvffegvUJhbH99tur9uDBg31s55bx48erdo8ePXw8e/bsPPQOAAAAKC516tRR7SZNmsTazl5P9+/f38eTJk1SuenTp/t44sSJSbsI5ATfxAUAAAAAAACADOMmLgAAAAAAAABkGDdxAQAAAAAAACDDqImLckfWwHUufh1cW5f0jTfe8HHz5s1VrkuXLqq92267+bh3794qd+utt8baPyqWffbZx8e2/vLcuXML3R0UUKNGjVT77LPP9rEdC+3bt1ftE0880cf3339/HnqHsrTvvvuq9siRI32866675n3/v/3tb1V7ypQpPv7+++/zvn9ki7zWefnll1Xuwgsv9PGQIUNUbtOmTfntGFJr2LChj59//nmV++STT3z80EMPqdysWbPy2i+rVq1aqn3YYYf5ePTo0Sq3cePGgvQJQNk74YQTVLtr164+PuKII1SuRYsWsZ5T1rl1zrmmTZv6uEqVKpHbbb311rGeH8g1vokLAAAAAAAAABnGTVwAAAAAAAAAyDDKKaBc6NChg49POumkyMd98803qi2XYCxZskTlVq9e7eNtt91W5caOHavae+21l4/r1asXo8eo6Pbee28fr1mzRuVefPHFQncHedagQQMfP/HEE2XYE2TZMccco9qhZXz5YEsFnXXWWT4+9dRTC9oXFJ69fvnnP/8Z+dj77rvPx48++qjKrV27NrcdQ2p16tRRbXkdbEsWLFy40MeFLp/gnO7P+PHjVU6eQ22ZoRkzZuS3YwiqWbOmj20JubZt2/q4c+fOKkcZDPyPLEvonHMXXHCBj2XJMeecq1q1qmpXqlRpi/ffsmXLLX4OoJD4Ji4AAAAAAAAAZBg3cQEAAAAAAAAgw7iJCwAAAAAAAAAZVuY1cbt3767asu7JvHnzVG7dunU+HjZsmMotWLDAx9RGqngaNWrkY1sbR9b/svUG58+fH+v5L730UtVu06ZN5GNfffXVWM+JikXWBXPOuQsvvNDHQ4cOLXR3kGf9+vVT7W7duvm4Y8eOqZ/3sMMO8/FWW+n/h504caKPP/jgg9T7QGFts83/XYodf/zxZdiTX9ehvOSSS3y8/fbbq5yt5Y3iJ+cX55zbZZddIh/7zDPP+Fhen6Ps1a9f38fPPfecytWtW9fHtubxRRddlN+OlWLgwIE+btasmcqde+65PuZzXtnq3bu3at98880+bty4ceR2snauc84tXbo0tx1D0bLnmosvvjjv+5w6daqP7W/moPi0aNHCx/Ic6Nyvfy/piCOO8PEvv/yickOGDPHxxx9/rHJZOvfwTVwAAAAAAAAAyDBu4gIAAAAAAABAhpV5OYXbb79dtXfddddY28llNc45t2rVKh+XxVfi586d62P7msaNG1fo7lQ4//73v30sv07vnB4by5YtS/X8p556qmpXrlw51fOg4tp9991VWy5NtssdUfzuuusu1bbLddI6+eSTNxs759zs2bN93KtXL5Wzy+SRHZ06dfLxgQceqHL2eiLf6tSpo9qydFC1atVUjnIKxa9KlSqqPWDAgNjbyjJAJSUlOesTtty+++7rY7ls1Bo0aFABehNtjz32UG1ZuuzFF19UOa6TypZc7n733XerXL169Xwcmgvuvfde1ZZlxZxL/xkN2WGXscuyCHZp+ujRo328fv16lVuxYoWP7bWGLe305ptv+njSpEkq99lnn/n4yy+/VLm1a9dG7gPZJEsT2vlDfiay4zCJ/fff38c///yzyk2bNs3HH330kcrJsb5hw4bU+4+Lb+ICAAAAAAAAQIZxExcAAAAAAAAAMoybuAAAAAAAAACQYWVeE/fss89W7T333NPHU6ZMUbnWrVv7WNZ7ck7XfDrggANU7vvvv/dx48aNY/fN1sFYvHixjxs1ahS53Zw5c1SbmriFJetCbonLL7/cxy1btgw+VtbckTHwP1dccYVqy3HKHFE+vPbaaz7eaqvc/B/p0qVLVXv16tU+btq0qco1a9bMx59//rnKbb311jnpD7acrOnlnHPPPPOMj2fOnKlyt9xyS0H69D+/+93vCro/lK127dqpdvv27SMfa6+JX3/99bz0Cck1bNhQtU855ZTIx/bp08fH8nNNocg6uG+//Xbk42xNXPn7Fii8yy67zMd169ZN9Ry2Vv+xxx6r2jfffLOPbf3cQtSYRDqyRq2sT+ucc3vttZePTzrppMjnGDt2rGrL+zyzZs1SuSZNmqi2/F2iXP3+BMqOvBd4wQUXqJycQ2rWrBn5HD/88INqf/jhh6r93Xff+dh+Ppe/IdKxY0eVk3Pf8ccfr3ITJ0708ZAhQyL7lit8ExcAAAAAAAAAMoybuAAAAAAAAACQYWVeTmHMmDHBtjR69OjIXJ06dXy89957q5z8WvR+++0Xu2/r1q1T7enTp/vYlnqQX6+2yyFRHE488UTVHjRokI+33XZblVu0aJFq//Wvf/XxTz/9lIfeodjsuuuuqt2hQwfVlvPJmjVrCtEl5Njhhx+u2q1atfKxXdIVd4mXXYJjl6atWLHCx0ceeaTKDRgwIPJ5//znP/v4gQceiNUX5MfAgQNVWy5FtMtLZfmMfJHXL3ZMszSxfAstu7fsXITs+Pvf/67af/jDH3wsPwM559zw4cML0qcohx56qI932GEHlXv88cd9/NRTTxWqS9gMW67pzDPPjHzsV1995eOFCxeqXOfOnSO3q1WrlmrLkg3Dhg1TuQULFkR3FgVlPxM//fTTPpblE5zTJaFC5VMsW0JBsmUrUdwefPBB1ZZlN+rXrx+5nb1n+PXXX/v46quvVjl7T0866KCDVFt+Xnr00UdVTt5jtHPd/fff7+MXXnhB5fJRuohv4gIAAAAAAABAhnETFwAAAAAAAAAyjJu4AAAAAAAAAJBhZV4TN1eWL1/u43fffTfycaGau6WRtcNkDV7ndB2O5557LvU+UHZszVJb80eyx/j999/PS59QvGxtSSsf9XGQf7LW8bPPPqtyodpN0uzZs1Vb1k664YYbVC5UY9s+zznnnOPjBg0aqNztt9/u4+22207l7rvvPh9v3Lgxcn9Ir3v37j4+/vjjVW7GjBk+HjduXMH69D+ylrKtgfvee+/5+McffyxUl1Aghx12WDC/YcMGH4dqbqNslZSUqLb8O543b57KyWOaL1WrVvWxrU14/vnn+9j2+6yzzspvxxCb/X2ZGjVq+PjDDz9UOXm9a68vfv/73/vYjoXddttNtXfccUcfv/TSSyp33HHH+XjZsmXBviP3qlev7mP5OzDO6d+UWbJkicrdcccdPuY3YyouOy9cccUVPu7bt6/KVapUycf2s7L8TY/BgwerXNrfl6lXr55qb7311j6+/vrrVU7+PpetG15ofBMXAAAAAAAAADKMm7gAAAAAAAAAkGHlppxCPjRs2FC1//nPf/p4q630/e9Bgwb5mGUexWPUqFE+/u1vfxv5uCeffFK1Bw4cmLc+oXxo165dMC+Xt6N4bLPN/50245ZPcE6XXDn11FNVzi4/i8uWU7j11lt9fOedd6pctWrVfGzH3ssvv+zjmTNnpuoLwnr06OFjeSyc09cWhSBLgjjnXO/evX28adMmlbvpppt8TKmN8uGggw7abLw5cnnihAkT8tYn5M8JJ5yg2m+++aaPbYkUuVQ1CVs+6ogjjvDxAQccELndiBEjUu0P+VelShXVlqUv7rrrrsjt1q1bp9qPPfaYj+V50DnnmjdvHvk8dul9IcqAIFq3bt18fNVVV6ncnDlzfHzooYeq3IoVK/LbMRQFeU5wzrnLL7/cx7J8gnPO/fDDDz6WpUydc+7zzz9PtX9ZIsE55xo3buxje4/ntdde87EtnyrZfg8dOtTHhSg/xjdxAQAAAAAAACDDuIkLAAAAAAAAABnGTVwAAAAAAAAAyDBq4gZccMEFqt2gQQMfL1++XOWmTZtWkD5hyzRq1Ei1ZT04W/9J1qmUdQGdc2716tV56B2Knaz9duaZZ6rcl19+qdpvvfVWQfqEsjFu3DjVPuuss3yctgZuaWRtW1nn1Dnn9ttvv7zsE5tXq1Yt1Q7VhUxbhzKtc845R7VlbecpU6ao3LvvvluQPqFwkswFhR6bSOeee+5R7U6dOvl4p512UrnDDjvMx7amX9euXVPt3z6PrJ9qffvttz6++uqrU+0P+ff73/8+MmfrLMvfFwnp0KFD7P2PHTtWtfncVbZC9dPl55u5c+cWojsoMrYmrf39Bennn3/28f77769y3bt39/Huu+8e+Rxr165V7datW0e27WeyHXbYIfJ5pYULF6p2oX9Dgm/iAgAAAAAAAECGcRMXAAAAAAAAADKMcgrGwQcf7OOrrroq8nHdunVT7UmTJuWtT8idF154QbXr1asX+dinnnrKxzNnzsxbn1B+dO7c2cd169ZVudGjR6v2unXrCtIn5M9WW0X/P6hdAlQIckmr7Vuor9dff72PTz/99Jz3qyKy5Xl23nlnHz/zzDOF7o6y2267Rea4lin/Qkuaf/zxR9WmnEJxGD9+vGrvueeePt57771V7thjj/Xx5ZdfrnKLFy/28RNPPBF7/0OHDlXtiRMnRj72k08+8THX1tllz1Oy1IYtySKXNbdr107lTjrpJB/XqVNH5ex8I/Nnn322yskxNnny5GDfkXtyGbsl55TrrrtO5V566SUfT5gwIfcdQ1F45513VFuW6pKfnZ1zrkmTJj7+xz/+oXKhUj2yRIMt3xASKp/wyy+/qPaLL77o4379+qnc/PnzY+8zF/gmLgAAAAAAAABkGDdxAQAAAAAAACDDuIkLAAAAAAAAABlGTVzj+OOP93HlypVVbsyYMT7+9NNPC9YnbBlZx2nfffeNfNx7772n2rauD1Cavfbay8e2bs+IESMK3R3kwXnnnedjWyuprHXp0sXH++yzj8rJvtp+y5q4yI1Vq1aptqwFJ+tVOqfrZy9btiwv/WnYsKGPQ7XtPvroo7zsH2XnkEMOUe3TTjst8rErVqxQ7blz5+alT8iv5cuX+1jWHrTtK6+8Mif7a968uWrL+uy2DuZll12Wk30iv95++23VlnODrXsra9SGalba57zgggtU+5VXXvHxb37zG5WT9SfldRgKo0GDBj6215DyNwCuvfZalRs4cKCPhwwZonJjx471sayD6pxzM2bM8PE333wT7Nsee+zhY3t/hnNYNqxdu1a1Za3s2rVrq5z8TSr5W1XOObd06VIfz5kzR+XkOJSfx51zrmPHjgl7/P956KGHVPvqq6/2sa3pXWh8ExcAAAAAAAAAMoybuAAAAAAAAACQYdzEBQAAAAAAAIAMq/A1catWraraxx57rI83bNigcrJG6saNG/PbMaRWr1491Zb1S2ydY8nW7Vq9enVuO4ZyZ8cdd1TtQw891MfTpk1TuRdffLEgfUJ+ybqzZUHWJWvTpo3KybkuZPHixarN+Sz3bP2vmTNn+viUU05RuVdffdXHd955Z6r9tW3bVrVtjcpdd93Vx6GahVmr84wtZ6+Jttoq+vsbb731Vr67g3LI1sGUc4ytu2vPP8gmW5+9Z8+ePra/8VCrVq3I57n33nt9bMfCunXrVHvkyJE+lnUxnXPumGOO8fFuu+2mcvL8ivy44447fHzJJZfE3k6eb84//3yVs+1csPOL/L2bU089Nef7w5aztWXt334aTz75pGqHauLa37CQ4/vxxx9XuU2bNm1x33KFb+ICAAAAAAAAQIZxExcAAAAAAAAAMqzCl1O4/PLLVXufffbx8ejRo1Xuk08+KUifsGUuvfRS1d5vv/0iHztq1Cgfy3IZQBx/+tOfVLthw4Y+fv311wvcG1QEAwYM8PEFF1wQe7tZs2b5+IwzzlC5OXPmbHG/ECbPL5UqVVK5E044wcfPPPNMqudfsmSJatuSCfXr14/1PHbpGIpf9+7dI3N2GeODDz6Y7+6gHOjRo4dq//GPf1RtuTx16dKlBekT8uvtt9/2sZ1TTjvtNB/bOUWW2rDlE6wbb7zRx61bt1a5rl27bvY5nfv1NQ1yTy5xf+6551Tu6aef9vE22+hbS40bN/ZxqJRPrsiSY87psTpw4ECVu+mmm/LeHxTOFVdc4eMkpTPOO+881U57HV5ofBMXAAAAAAAAADKMm7gAAAAAAAAAkGHcxAUAAAAAAACADKtwNXFl7TnnnLvmmmtUe+XKlT4eNGhQQfqE3LrkkktiP/bCCy/08erVq/PRHZRjTZs2jcwtX768gD1BefXaa6+pdqtWrVI9z+TJk3380UcfbVGfkNzUqVN93LNnT5Xbe++9fdyiRYtUzz9ixIhg/oknnvBx7969Ix+3du3aVPtHtuyyyy4+lvUqrblz56r2uHHj8tYnlB/HHXdcMP/KK6/4+Isvvsh3d1Bgsj7u5tppyfOPrbsqa+J26tRJ5erWrevjZcuW5aQv0DZt2uRje55o2bJl5HZHHXWUjytXrqxy119/vY9Dv1+zJeRvELRv3z4v+0DZ6Nu3r2rLmse2NrP1zTff+HjkyJG57ViB8E1cAAAAAAAAAMgwbuICAAAAAAAAQIZViHIK9erV8/E//vEPldt6661VWy5dHTt2bH47hjInl+Bs3Lgx9fOsWLEi8nnk8pFatWpFPkft2rVVO25ZCLnExTnnrrzySh//9NNPsZ4D6Zx44omRuX//+98F7AkKRS7N2mqr6P8HDS03feihh1R7p512inys3ccvv/xSWhc3q0uXLqm2Q/5NmDBhs3Euffvtt7Ee17ZtW9WeNGlSPrqDPDvooIN8HJqnRo0aVYjuoJyx57c1a9ao9t///vdCdgfl0PPPP6/aspxCr169VE6WxqMUYraMGTMmMidLSdlyCj///LOPH3vsMZX717/+pdp/+ctffBwqH4Ti17FjRx/b80z16tUjt7MlM8877zwfr1+/Pke9Kyy+iQsAAAAAAAAAGcZNXAAAAAAAAADIMG7iAgAAAAAAAECGlcuauLbO7ejRo33crFkzlZs5c6ZqX3PNNfnrGDLnq6++ysnzDB8+3Mfz589XuR122MHHto5TPixYsMDHN998c973V9EccsghPt5xxx3LsCcoCw888ICPb7/99sjHvfLKK6odqmWbpM5t3McOGTIk9nOi/JO1nGVsUQO3fJC/BWEtWbLEx/fcc08huoNyQNYQlNe1zjm3aNEi1f7iiy8K0ieUX/ZaR15v/e53v1O56667zsfPPvusyk2fPj0PvUMuvPnmmz62n1e32eb/blGdffbZKteiRQvVPuKII2Ltb+7cuQl7iKyRv+9Ro0aNyMfZOu2yprZzzn388ce57VgZ4Ju4AAAAAAAAAJBh3MQFAAAAAAAAgAwrl+UUdtttN9Vu37595GMvueQS1bblFVB8XnvtNdW2y27yoUePHqm2+/nnn30cWib98ssvq/a4ceMiH/vhhx+m6gviOemkk3xsS7d8+eWXPv7ggw8K1icUzsiRI318+eWXq1yDBg3yvv/Fixf7eMqUKSp3zjnn+NiWdUHFVlJSstkY5dMxxxwTmZszZ46PV6xYUYjuoByQ5RTsHPLqq69GbmeXvNapU8fHciwCIRMmTPDxtddeq3KDBw/28S233KJyp59+uo/Xrl2bp94hDXkN+/zzz6tcz549I7fr1KlTZG7Tpk2qLeemq666KmkXUcbs+eOKK66Itd2wYcNU+7333stVlzKDb+ICAAAAAAAAQIZxExcAAAAAAAAAMoybuAAAAAAAAACQYeWmJm7Tpk19/Oabb0Y+ztYwfOWVV/LWJ5SNk08+WbVl/ZTKlSvHfp499tjDx7169Yq93aOPPqras2bNinzsCy+84OOpU6fG3gcKp1q1aqp9/PHHRz52xIgRPrZ1mVA+zJ4928ennnqqynXr1s3HF198cV72f/PNN/v4/vvvz8s+UP5st912kTnqBBY/e21jfxtCWrdunY83btyYtz6h4rDXO7179/Zx//79Ve6bb77x8RlnnJHfjqFcevLJJ1X73HPP9bH9DDho0CAff/XVV/ntGBKR1x5/+ctfVK569eo+7tChg8o1bNhQteXn7KFDh6rc9ddfv4W9RKHJYz958mSVC93HkX/fdjyVR3wTFwAAAAAAAAAyjJu4AAAAAAAAAJBh5aacwjnnnOPjJk2aRD7u/fffV+2SkpK89QnZcPvtt2/xc5x22mk56AmKkV1uunz5ch+//PLLKnfPPfcUpE/Ihg8++CCybcv6yHNUly5dVE6Oo4ceekjlKlWqpNp2aREQx5lnnunjH3/8UeVuvPHGQncHOfbLL7+o9rhx43zctm1blZsxY0ZB+oSKo2/fvqrdp08fHz/yyCMqx3yDLbV48WLV7ty5s49tCbsrr7zSx7LMB7Jl4cKFqi2vk08//XSVO+CAA1T7hhtu8PGiRYvy0DsU0pFHHunjXXbZReVC9+1k6R5ZNqq84pu4AAAAAAAAAJBh3MQFAAAAAAAAgAzjJi4AAAAAAAAAZFjR1sQ95JBDVPuiiy4qo54AKM9sTdyDDjqojHqCYjJ69OhgGyik//znPz6+8847Ve7dd98tdHeQY5s2bVLtAQMG+NjWkBs/fnxB+oTy5cILL/TxoEGDVM7Wh3/ggQd8LH9HwDnnNmzYkIfeoSKbM2eOj99++22V69q1q4/btGmjcvzGQHEYOnRosI3yRdZND9XAHTx4sGpXtGtZvokLAAAAAAAAABnGTVwAAAAAAAAAyLCiLadw6KGHqnb16tUjHztz5kwfr169Om99AgAAyJouXbqUdRdQQPPmzfPxWWedVYY9QXnx0Ucf+fjII48sw54A0bp3767aEydO9HGLFi1UjnIKQPbUrVvXx5UqVVK5RYsW+fjuu+8uWJ+yiG/iAgAAAAAAAECGcRMXAAAAAAAAADKMm7gAAAAAAAAAkGFFWxM3RNa/cc65o446ysfLli0rdHcAAAAAAECerFy5UrWbNWtWRj0BkMadd9652dg552688UYfz58/v2B9yiK+iQsAAAAAAAAAGcZNXAAAAAAAAADIsEolJSXxH1ypUvwHI9/Gl5SUdCjrTsTBuMmOkpKSSmXdhzgYM5nCXIM0GDdIg3GDNBg3SINxgzQYN0iMz+BIIXKu4Zu4AAAAAAAAAJBh3MQFAAAAAAAAgAzjJi4AAAAAAAAAZNg2CR+/xDk3Ox8dQWJNy7oDCTBusoExgzQYN0iDcYM0GDdIg3GDNBg3SINxg6QYM0gjctwk+mEzAAAAAAAAAEBhUU4BAAAAAAAAADKMm7gAAAAAAAAAkGHcxAUAAAAAAACADOMmLgAAAAAAAABkGDdxAQAAAAAAACDDuIkLAAAAAAAAABnGTVwAAAAAAAAAyDBu4gIAAAAAAABAhnETFwAAAAAAAAAyjJu4AAAAAAAAAJBh3MQFAAAAAAAAgAzjJi4AAAAAAAAAZBg3cQEAAAAAAAAgw7iJCwAAAAAAAAAZxk1cAAAAAAAAAMgwbuICAAAAAAAAQIZxExcAAAAAAAAAMmybJA+uVKlSSb46gsSWlJSUNCjrTsTBuMmOkpKSSmXdhzgYM5nCXIM0GDdIg3GDNBg3SINxgzQYN0iMz+BIIXKuSXQTN6RSpUqR7V9++aUotkuybQa2mx25YRHJ1XHcaqv/+1L5pk2bctS73LP9LimJnidzPW5C2xSTLM0Z+ZprQtsVcsz8/9sx1xRwu9KeJ3T8027HuImWpeMf9xhurs21TWEVy3zDeSpbmG+Yb9LI0rzBZ/DikaXjX4jr4rTSnt/iblNMimXMVIS5hnIKAAAAAAAAAJBh3MQFAAAAAAAAgAzLWTkF+zXxuF8bL5bt7GOzvF0xydVxzEWpgNByiajlEZsTeg2Mmy1XLHPGlsw1abdjzESrCMc/7XaMm2jFcvy5tsmWYplvOE9lS7EcR+abbCmWeYNxky3FdPxzsV1adn+hsjPldaz8T7GMmYow1/BNXAAAAAAAAADIMG7iAgAAAAAAAECGJS6nkOZX90JL0bO0XVnsc0v6Wkzy/WuNodIHoX1UrVrVx/Xq1VO5NWvWRD5n3bp1VXvu3Lk+3rBhg8ol+fVEqaKPm2Kca5K892n3udVWW0U+rqKPGecKO27S2pL9Ffr4M26y855ybZM95X3cMN/kR3kfN8XU12LCuMnGdsWmvI+bQp+n5P7sZ/zyMm7K+5gppr5G4Zu4AAAAAAAAAJBh3MQFAAAAAAAAgAzjJi4AAAAAAAAAZFjimrhbb721j2UdEFkfxDld38HmsrqdzYjijOQAACAASURBVNs6J6GaK/nYLlQ/I1RnNYvSjJsktWS32WabyJzcR4sWLVTukksu8XHr1q1VrmbNmj5u3Lixys2bN0+1u3bt6uNZs2ap3M8//+xj+5qi3hfb71yMm4owZspizogSGoc2HzqGoWMvx1Zp24XGWkWfa9Ie/1DtLbudbNv3W+4vdCzs8+Tq+HOOyu58I4XGDdc2ZSOr4yZX5zd5bbVp0yaVK22u2tz+StuOcZP74yj35Zw+jsw3jJusbGfzjJvcKW/jxpKvL3SessdUtpOc3+R5cf369ZHbFXN93PI2ZsrjXMM3cQEAAAAAAAAgw7iJCwAAAAAAAAAZlricgvxab+hr4qGlwnGXIYSWqiaRZFm+fGyS7eIusbVC24WWABSbqNcZWtoQElqqXqVKFZUbNmyYj48++miVk8snqlWrFtnnb7/9VuUWL16s2vfee6+PBw0apHJjx4799QvYjCTju6KMm/9JO2byvZ3NyfFU2v5q1Kjh46uvvjpy/7fddpvKrVy5MrLfof1VlDGTZq5J8prjjht7/ggtb+3Zs6eP+/fvr3KtWrVS7Z133tnHtpxC3PNQkvIdcXPFPm6kfIwbecw3btwY+Tj7nE2aNPHx/PnzVU7ON6WdP+xYiepb7dq1VW7NmjWb3Z9z8a+RGDdle56y7788bnYukstEbd7m5Di2+5DXU6tXr47sN/NN+uMY2i7ufGOPqSw59v3336vcunXrIp8nNMZCc5Nd8hqap+RrDI3TDRs2RPatmJc3O5f/65vQkt+4833oGG7J/kPiLmGvqNfFUtpxE/dvJ+29myTjJnTvxu5PftbaaaedVK558+aROdmfbbfdVuVk2UT7Gd+e74pVLu7bxP1btOeBuNe2ofNOadKeF9JuF/faJgm+iQsAAAAAAAAAGcZNXAAAAAAAAADIMG7iAgAAAAAAAECGJa6JG1ULIkldq5BQzaO0NTJC0tZHCtV8SfKcofodoRqGxSaq/2lfl32P5dgYMGCAyh1xxBE+tnVtZG3b//73vyr32Wef+djW6unYsaNqt2/f3sfHHHOMyoVq4oZef5KazFJ5GTe2/uL/pH1Nad/PuH+jpW1n2zfccIOPW7durXKyxvLy5csj97El/YnarpjHjHPp5ppQ/dpc1Huy7VNOOUXlHnvssdjPs8suu/h49uzZKpekZmHcXEUZN2nmmyTjJu77U7NmTdWWNdb79u0bud2W1HqU29r3oXHjxj62NTLXr18f+ZyMm+I8T1nPPPOMj21t1YkTJ/r44YcfVjlZSzlJjcSKMm7SnKfSXm/YusdNmzb1sT0XyRrcl156aeT+SzumaX8LRNZStr9TIX/v4qefflI5WYcyV7+nkkW5uL6R0l4XhK5v2rVrp3Lz5s3z8YoVK1TOzin5mDfj/h5Aku2KTa6vb3JVdzgueR3inHNnnHGGau+5554+PvTQQ1VOzhuhcVu9enWVk5/75fnMOefeeOMNH8vzoHPOrVq16tcvoAilmWuSfF4K1TGOO0625O+ycuXKPk7y+yJppT0nhvBNXAAAAAAAAADIMG7iAgAAAAAAAECGJS6nIIVKJqRdjiqXz9SuXVvlli1b5mO7jE9+pXqrraLvTZf2NeW0S2flPpN8LTzue1jsol5nrspZyOUS5557rsrJMTV16lSVO/vss30syyc451z9+vV93Lt3b5U77bTTVFuWaShtGb2Ui7Ij5XXcpBkz9m8/7fKfqOdIIlRixjnnzjnnHB/bcfn+++/72C6DkktAkiwlqQhjxrl04ya0xCpX703z5s19fO2116qcLB1kx7BdxnXYYYf5ePjw4SoXKqfAOSosF+NGSvLeyOXOf/7zn1WuR48ePj7vvPNUzi4By4W2bduqthxv99xzT+znYdwU53nKLr3v1q2bj+25SF6jP/jggyonr4ns0veQijJuouTqOEr9+vVT7f79+/vYnjM++OADH2/Jsui451Cb69Onj4/32Wcflfvuu+98LMt8OKdfRz7mxawo5HkqyfGXJToef/xxlXvvvfd8LMuIOffr8gqFvk6pKPNNLs5TUj7emwYNGqi2LE144403qtyOO+6o2vJzke23fO32HBb625DX5XZOOfDAA30cuu4uZmU518Q9f9hjHSrRIMeIc7p00KJFi1ROXrPYMSP3GSo5YuVjruGbuAAAAAAAAACQYdzEBQAAAAAAAIAM4yYuAAAAAAAAAGTYFtXElWxdCllXy9a6kLk99thD5WStQFkD1znnRo8e7eMpU6aonKwpaWtUyHoWtt5XqH6HrYMht7XPI+tZJKm5EqqREapHk6QWZhZEvU77GmXbvo+SrW3SrFkzHy9ZskTlpk+f7uMBAwao3NixYyP3t3jx4s323znnqlevrtpr166NfJ7tt9/ex6tXr1a50NgIHf8046YijJlQnZl8bBeqvWa3q1Wrlmr/+OOPPn722WdVLlRjKe58Esox14TnmlBNdzn32/dG5tq0aaNyTz31lI/t/CGfx55b7JzRq1cvH8+aNUvlxo8f7+OlS5eqnKzvZfsdGhuMm3R//3bchN5/eR1k66/L57HnvQ0bNmy2z6W1bd+qVq3qY1kD1znnTjnlFB+/+OKLKvftt99uti92H4ybwp+nQnUAJTkunfv18ZfPY/fxxRdf+NjOU3KOC9XWtBg3ubnekDWJbZ3tOnXq+Pjee+9Vudtvvz3U7cj9pa1Daue0iy++2Mfy9yycc+5vf/ubj+fPnx/5nOVZWV7fyLnC1it++umnfdy4cWOVk+0nn3xS5SZMmBB7/6G5MO68kaQeeUWfb0KvMXTPx85ToXOI/Pu3c5Gsxy5/92Zz5D5Dc4G9TpH7t9vJvxs7TmUNaDuH5eM3NcpCIa9tQvsOXdvI85y16667qnaXLl1U+/jjj/fxxIkTVe7zzz/38UsvvaRycqxt3LhR5ULvRdzfR7KCNYEjMwAAAAAAAACAMsdNXAAAAAAAAADIsMTlFKK+Xp2kLIFcVtynTx+Va926tY+32247lZPLDO3Xi1etWuVjuzysZs2akX2zy5blPu1Xn7/55hsfz507V+UWLVrk4/vuu0/l5NfwQ8uhKqLQMrtQzo6padOm+Xj33XdXOXkc7diQy0jXrFmjcg0bNvSxLNfgnHMzZsxQbbm0Qi4xdc659evX+7jQS76KeYlZVN+TjJm028l22vfQLrE59thjVVvONZMnT1Y5u0RDCi2NzcXxLuYxExI6/raEQagET2jZmDy3vfzyyyq38847R/bN7kOypReOPPLIzcbO6XEjy8E451yTJk18nGQpYHkdD1aa+caeT+S4CS0Bs8/ZvHnzzT6/c86NGDHCx3a8yeOYZNmefWxo3DRq1MjHdpnsf//739j7jFLs4yur56m472uNGjVU25ackux5SV732vOSLK+QjyWlxT5ucv2e2HNI586dfWyvX+WS9uuuu07l5DGWpcCc09eyW9J/eezs+JNLYO11kSypJ0uYORd/aa5VzMudpVydp+Q4sn/TLVu29PFrr72mcrJEh33OevXq+XjMmDEqN2rUKNXu16+fj5Oc7+T1duj6xr6m0LiRinm+qVSpkhoD8m88bcmIUDkL+9lHlmS67LLLVO6dd97ZbL+c+/W1iCTLJDrn3MqVK338xBNPqNykSZN8vHz5cpWT947seJPvjT3+8r0JfSarCEJzaNzrntC9MTt/tWjRwseyhJxz+p6OfU55n8g559q2bevjQw89VOUWLFjg46OPPlrl5GNfffVVlbv66qt9bEtMyb7Z81fa8xDfxAUAAAAAAACADOMmLgAAAAAAAABkGDdxAQAAAAAAACDDEtfElXUbZGxrp4RqqYRql8hao7NmzVK5gw8+2Me2jlK1atV83K5dO5WTNZZsTQ5b23b48OE+tjVSW7Vq5eP99ttP5eQ+R44cqXJffPGFixKq8VNeajU5Fz0ektRmjDveQu+brbkja6bY2kC/+93vfHzaaaepnK0/Jus8vfjii8F9SnFrDoXqD4W2S/L+Zk1UnaG0ryntdmn/Dm1dqD/84Q+qLceFrdUT6mvcsV4Rx4xz6eaa0Ptox6F8HvuePvDAAz7eaaedVC5UXyvk+++/V21Zf9uS+7C1dBs0aODjZcuWqZwci/a9iDsWi12ux43dTj72wAMPVLmePXv6eMOGDSr34IMP+ljWpCxt/0lUqVLFx2effbbKyXls/vz5Kpd2bJSn+aYsz1Oh4x93jrG1utu3bx/5PLbG3PPPP+9j+z6kHZsVZdxE9T/t5wBZk9Q552699VYf2+tVOd+H6jna3wzJFdnXgQMHqpysFXjzzTer3MyZM32cdl4udmnOU6FzkX2v5HiQNRydc+6mm27ysa2XLM9b2267rcr99NNPPrZ1Iu2ccumll/p42LBhKifvCdhxGxrHFf08VVJSoo5P6PjHZT8v165d28fyt42c0zVE3333XZW79957fWznsPvvv9/H9l6N3b98fUuXLlW50NgIvRfleR6JI805Ksk1sWR/A6tDhw4+PvPMM1VO/oaErZssf8tI3mt0TtdGdk7/xoPtm5yz7GcueU/xT3/6k8rJ+4bXXHONyslxmKuxxTdxAQAAAAAAACDDuIkLAAAAAAAAABmWuJyCJJdZJflqsFw+8eWXX6rcV1995WO7dHDUqFGb3bdzetlHrVq1VE5+7d5+ZVqWYXBOl3OoX7++yl199dU+ll+nthYuXKjaoa+Qx13iXOyixkra5YBp35vQMoAWLVqoXO/evX0sl5s659zkyZNV+6KLLvKxXA62JeK+N+V1CUiaMWOX2OR7O0v2eY899lC5ww47TLXnzZvnY7tUKB9/FxVhzDiXbtyEhLY79thjVXv//ff38Tbb6NOrPJ/ZUhty6c5nn32mcna5qTx/XnzxxSrXuXNnH9vzYNeuXX0sl0E759yKFStclLTn+fIiV+NGzjFDhw5VObnMql+/fir30Ucf+dgui07LziFy3NpyVHJp4pw5c1SO+aaw56kkZVji7r9Hjx4qZ5dQS3ZusuU1ovaR5PhX9HFjxT2O55xzjsrJ8m/2OSZOnOjjfC0Tl6/PLq+Xn6X69u2rcuPHj/fxSy+9pHKhZdFyf/bvpCKMm3xc39jPPfI8YbeT1xQdO3ZUOdnnevXqqdyRRx6p2nKs2PONLbEoyb+F0DixKsp8EyXta7TL3x9//HEf2/PbO++8s9nHOaevL+RnIueSLT9Pe52a9vOUfY1xtysmuZ5rLPm3/q9//UvlTjzxxMjtpk6d6uNbbrlF5WSZFzsP2GubF154wcfHHXecysl5qWXLlion3xc7R3br1s3HgwYNUjk5Zuz9zbT4Ji4AAAAAAAAAZBg3cQEAAAAAAAAgw7iJCwAAAAAAAAAZtkU1cSVbH0TWA9ywYYPKyToVtmZF2romsg7pqlWrYm9n+y3rQy1fvjzysbYezOzZs328bt06lZP1Dzdu3Khyodcr92frtiSp+ZMFUa/Tvv+yLWs/hp7PtpPU+JI1WWxNMVn32L7fssaPc7qWT9p6PHY7Wf8wVOMrNG6kfNU+y5c0YyZU0yq0Xei9STJG5Xg69NBDg8/z+eefR+5THvskc6TcLjS3xR0zdrtikO+5RtZR/8tf/qJyDRs2jHweeRxvu+02lXvooYd8vGzZMpWz5xPZn0svvVTlZG1uW1tV5uy57e23347MhWqGlaea7rkeN7YmsqyrZWsQP/vssz7+8MMPY/Y4XAfS9lv21V6/HH300ZHPM2HCBB+H5g07T8S9trEqwnyT9jyVZJ6S489ed9atW9fHJ5xwQvB5fvzxRx/feuutkY+1r0nu3x5TmQuN0/I8bqLY1yjP6fb4N2vWzMd//OMfI7ezc/p3333nY1uvVn5eC9VutrUAba3TJUuW+Lh79+4q94c//CFy/7IGfJK6gaG5KFRLnM9S+u9dHo9evXqp3IIFC3y8cuVKlbv22mt9vOeee6pc27ZtfWxrcO+4446qXadOncj9jxgxwserV69WubTXKRXlPBV1PRY6b9i/Dfk7C3fddZfKyWuIRx55ROWGDx/uYztupNA4LU3c6007F8Sts21/x0Ief/sc8j2Uca7qoJa1tHONPWfIv29bG1t+zpo+fbrKyc8y8ne0LFsDV54vndM1t+3nrDZt2vjY/r6NZK/z5TiR11nO6TkrV3Xb+SYuAAAAAAAAAGQYN3EBAAAAAAAAIMMSl1OIuwRZfr3cft1Ybhda5mC/si2fM7TMIbQ/+1V6+1j5dXe7zKNdu3Y+/v7771Xuqquu8vGUKVNUbkuWCJR3oa+U26/ey/cx9HX+0JIFu4z0lltu8bFdfiz78t5776mcXVaY7yVZab9qH/r7yrqovofGTJIlVXI7Ow/IsWbnjNDyRvmcTZs2jXxO55y74447fGyXcuRiqVba5yjmMROSZNzInF1GJZeNHnDAASonlwDJEj/OOde1a1cfz5w5U+Xmz5/v49KOm+zb0qVLVW7UqFE+7tatm8o1b97cx3fffbfKyaVwcvm0c+HSHhVR3LlYjgXnnBswYICPW7VqpXJyudhPP/2Uk76Eyhs0atRI5Y4//vjI55QlQ9asWRPcRxrFPt8U8jxlz0Xy/bfnsFA5BXltu9tuu6mcvbaSy+vlMnzbN0vu0z5OziMV9TwVt/+hcSNLFkybNk3l5N+4vZZ+/fXXI/cnz1v2uvfwww+P3O6LL76I3L889zinj/lLL72kcuPGjfNxkuve0DiK+7dXzJJc38q/P3t9I9v9+/ePzMnj5JxecmyvL+S58OOPP1a5p59+WrVlqakffvhB5ex1shS6J5CL0k7lddzYc4osmWCvL4899lgfn3zyySony3CceuqpKvfXv/7Vx/ZcJIVKwFlpS3nZa1j5+kP7s/2OW2ojdB7MuqhSI0n+FuT7ZD+TDB482Me2HE/os/RZZ53l45tuuknl5Ji15z37maxTp04+vvLKK1VOlsi0171HHXWUizJjxgwf2+MtnydXY4Fv4gIAAAAAAABAhnETFwAAAAAAAAAyjJu4AAAAAAAAAJBhiWviRtUMsfUd5ONs7RD52FBdCJuT7VBtviTPaZ+ndu3aPpY17JzT9Q8XL16sch988IGPbe2UJHXrorbLRe25shTV/9DrkrXYnMvN+9GiRQvV7tOnj4/tOF21apWPbb2U5cuXq3baekmh15H2NZaXcRP1N572NSV5r0NzVGju2WGHHXwcqhnlnHPffvtt5HOmnTMq+phxLt1cE3q/bS1jOS/UrFkzcjtZ59Y5595//30f2/ki7Xtux41sN2zYUOVkXTpbI0/WLJw8eXJwH1Kx1fgKiXttE/c1V69eXbU7duzoY3v85fVEknpuSfopz2+HHHKIysnzoj0PynGcjxqVxT7fFPI8lWSeCo2V9u3b+3j77bcP9uehhx7y8Zw5c2L3J+3xj1vrsNjHTdzPQfI4hmrs2fq1O++8s49/85vfRO7bks9pa2LKmoJ2nrC1vOU5xR5T+Tw33HCDyoVqZoZUlPkm19c39v2W403WhXTOuV122cXHdr6pW7euj88991yV69mzp4/3228/lZO/NeOcrlEqa7U7p8+ptnZ/6O8kpKKMmyj2PLFy5crIx3bo0MHH9vOMfK9s7XQpSa34JJ+ro+q3libuPGxVhHEj+572NYXu28nPUnL+KG0f8vOLvGfnnK7b3rdvX5Vr1qyZasvPb3Xq1FG5vfbay8e2Jq58TZMmTVK5M88808f2M2A+6mrzTVwAAAAAAAAAyDBu4gIAAAAAAABAhiUupxD1lfW0yyrLejvbbtOmjY+7dOmickuWLPHxPffco3Lya+FJhJYAlKelqmU5brbbbjsfDx48WOXk8hz7nDfeeKOPv/zyy+A+4i4BDG2X5PjLpWzlddykGTNpS7eEhN5fuyz9iiuu8HGNGjVUburUqaqddtkgYyYs33ONLJlhySVARx11lMrJ5Yd2nOaKXBotyyeUJvRYubwxtES72OV63MjlxM7p5Vrr169XuQULFmzx/ix7TqpSpYqPZRkh+9iZM2eqXNpjzrVNds5TdumrXN4s/76d+/XYlOUUkoyFtMe/oowbKe1SYHm+efjhh1VuzJgxPu7WrZvKyWWltmzZunXrfPzGG2+onFzCbkuTzZ49W7Xla7Lj5q9//auPv/rqK5WT11ShpddWRRk3+b6+kWPq3XffVTm5nL5x48Yq16RJEx/bcdO6dWsfy+XOzjnXqFGjyL7Jc5Zzzu26664+njBhgsqFrlMYN9FC5eJsqZ2WLVv62J5T5Dwyfvx4lZMlM8aOHatycrzZ4y0/I9nyHZY9j0XtI7Qd56loac9Rcruvv/5a5eRYsOUM5DWS/XuWc0///v1VTs5R8n6ec78+18nPPfZcI/e/aNEilZOf5Xv16qVytgSNlHashfBNXAAAAAAAAADIMG7iAgAAAAAAAECGcRMXAAAAAAAAADIsZzVxbY0v2Q7VI7Hbyee328n92foVoe1CbE3L3//+9z6eP3++yg0fPtzHDzzwgMqFaoWFalGG6oyE6iYmeY1ZkGbchGrg2OMvj6M9FrL+2xFHHBH5PJ999pnKPfLII5vt8+aExrusg2Jfr+xr6PgnqaEXNW4qwpgJ1ZnJ1VirXLmyj3faaSeVk3XBli9frnJ33323aieZF6P6k+8xU1o/syjXc42ttdW+fXsf2/fmo48+8vHChQsjnzNurS3bT+fC5xo55lavXq1ysr6ZrHvonHM777xz5P5D9b7T1skqJknGjRwrNidrZdlxI89fsoa7c78+VnH7WbduXdWWNXplLXjndB3UUaNGqZw8xnYsyNdrc/I1Mt8U/jwl6wvaenMNGjSI3M7WxF2xYkXsfUbl7OsNjZuKcp6KYl+jnBtsTT/5mu1xkvUHbS1C+b7a85u8vrF1L7/44gsfyzHk3K9rxffu3dvHr7/+usqNHDnSRZH9scdUvjcVddzk+vrGqlq1qo8//PBDlZPnlCOPPFLlDjnkEB8//vjjKidrU9arVy9yf87p6xv7+Vz+Tkn37t1VLlQTNe08VZ7GTRT73sjf8thxxx1VrmnTppHbyfdu7733Vjk5V8jfFnLOuY8//tjHU6ZMUblly5b52B4LWZ/bOX087LVvqD536O8k1+Om2MZM3LlG/u3Z33qRr/nbb79VuWuvvdbH9jdjZP3lPffcU+Xkb1B16tRJ5WS9XHv+sv2Wx9sea7mPoUOHqtyQIUN8LMdoaeR8FrpPZPsS+v0cvokLAAAAAAAAABnGTVwAAAAAAAAAyDBu4gIAAAAAAABAhm1RTdzQv4fqZ8iaJKGaI7aehawLYevfyOexdaNC/bS14c477zwf25qGsj6QrRsWt95ZWqH3qZglGTeyhkjo/bD1eLp27erjd999V+WWLl3q42nTpqncmjVrIvdhyX6Hajnbejyhujqh54y7XTFLM9fErWFkt7O1t+RxCo1DWWfUOV1P7tNPP1W5sWPHRvYtyTG0fZVyUaO0Is41IfYY9+jRw8d2bMh6TKHzUIjtl50z5HlRnq+cc+7iiy/2sa2JKs2ZM0e1n3/+eR/b9yluPb20729WRP3tJHldofpfw4YN8/Ell1yicpMnT458flnr8s0331S5k08+2cdy7nHu13W0QrVt5VitVq2ayoXqvYdqvIVqFEb1qxgV8jyVJCePTdu2bVWuWbNmPrbXOf/85z9VW/79J/mbjntcK+p5Ku7vhsj53/6Ny+1C9TvtdvLv3V5PyL9bW3dStleuXKly9lwkzZgxQ7XltqH6f6HrXiv0+uM+RzFLe/4NXd/KeqXO6Zr7thbl/vvv72P7GUzWZC7t+kLWojzmmGMi+2brVMoazFaxXYsUUujvRtbxd86522+/3cd33HGHytWsWdPH++23n8p17Ngxch99+vSJzMnatnYu+uabb1Rb/lbFTz/9pHKy1u4tt9yicuPHj4/cf0W/dxP386T8u7TXfaFa5XPnzvXx3/72N5WT84I9R8m2vYfXqlUrH9uxJT+fOefcUUcd5WNbR/nCCy/08YgRI1Qu9HsPoeMtz212u7j3tyy+iQsAAAAAAAAAGcZNXAAAAAAAAADIsMTlFOzSzv8JLQcKCW1nl6OGvpYd13bbbafaAwYMiHzs559/rtpffvmlj+37kLY/odcfWk5fbKL6H/f1l7adfOz111+vckcffbSP7dfUFy9e7OMbbrhB5eIuI95cf6L6lpbtd9yleMW8jKiQc40VtwTG4MGDVa5du3Y+HjVqlMotWrRItaNeX2lC26V9byr6XBNil+vYpcmSXFIWKiMU+ru029myQnJ+u+yyy1ROzll2zpDLhf785z+r3Pz5831cXueT0uRi3Mj3Ry49dc654cOH+/iAAw5QuQ4dOkQ+pyzLIMeXc3rZoC2DII+pc841adIkch+yPNR//vMflUtyHpQqyrVNIc9Tob+/UBmWxx9/XOVq1arlY1tO4a233lLttGVhOE+FRfU/tMwyScktyf4Nx72+CalTp45q9+rVS7WXLVvm4yeffFLl5JhLck5J8pkgarvyOm5yNd+E/t5nz57t44cffljl9txzTx+PGTNG5eQ1jS3t0bBhQ9WWZQw/+OADlevfv7+P9913X5WT10lJPp/H/fsqr+zxt8vKpffee8/Ho0ePVrndd9/dxzvvvLPKyWNsl9vL/dtzkSx3aK/Dmzdvrtq1a9febOycXn5vx428z5PkeFeE+SbNOSpU1iV0jRJ67+35S+7PltmQ88txxx2ncvb6We7/kUceUbkXXnjBx6HSZFahz1F8ExcAAAAAAAAAMoybuAAAAAAAAACQYdzEBQAAAAAAAIAMS1wTV9aCiFuj1tb4y0Vt2yTbyf3/8Y9/VLl+/fqptqwxN3DgQJVbu3ZtrH0kqasT9X6Wtl2xSTNukpDPc9BBB6mcPDa2tsr06dN9LOvvMiMFAQAAEk5JREFUJBWqURt6jbI+UGg7W8t5++239/GKFStULm0Nw6wp5FwTtw6cfZ5OnTqp3Lp163z8/fffq5w9TrJOU5Ixk3bOSDtHFZtczzWNGzdW7Ro1akQ+tkqVKj629Zfk/GLrgsmxsMMOO6jcAw88oNp77LGHj+14l/u356vrrrvOx+PGjVO5UA0m+X6G6l0Vu1yPGzsPz50718fnnHOOysnxIMeCc7oWnD3eDRo08LGtgStrCzqna4fZuofyusf2W/bH1gbj2ia718SyRuROO+2kcvLv3dZA/uyzz1Rbjs0k1xacp8LSjJu0uZC0n1fOPvtslbN1KF988UUfz5w5M9U+Sutr1HaMm/RCzynrpf7rX/+K7Feorq68RnZOn3vsPr/77rvI55F1vZ1zrlGjRrG2s9KOt/IiSZ1SeU1hf1dBnm/seaJ+/fo+btWqlcoNGzbMx/baR37OtdfTMuecPo722lc+rx2bst9c32i5nmtytZ3sV7169VTumGOO8bEdI/a66+uvv/bxtddeq3JxfwsgyWeifIwZvokLAAAAAAAAABnGTVwAAAAAAAAAyLCclVOwX1OOuyQqtF1o6Zb9CrP8urxcDuCcc23atPHx9ddfr3KVK1dW7ZEjR/p41qxZkfu0r0nu3+bkcjT7euXX9+129rFSaPlrFqUZN0mOv1wSsWbNmsjc8uXLVW748OGRzyn7YseJFfoqvBwbdrmIXBJkcwcccICP5RIB55ybNGmSj19//XWVW7RokY/lEtqKMGYKMdfI8WSXLMslzV999VXkczqnj7fdf+g1heYhuV2SJbvMNdHvv1xC6JwuwWKX8silgg899JDK9ejRw8f77befyl1xxRU+7tq1q8rZ0gtyjNlznRxTtjzM0KFDfbx+/XoXV2i5oXyf7ONs37Iu1+eo0DLFxYsXx+qHc/rcE3eJ1+b897//9fG8efNUTj6vXfpctWpVH9txE5qLJOab/J+n5HnJOX39YMeNvEZ47LHHVM4ud65WrZqP7XLT0HGMe77hPBX/+Ns5VW6Xq2Xh8nxj9ydL/dhyCnb/Tz31VOTzhITmFNk3e14MLQOv6OMmyfVt6Hwj//7tXJCWPY6hz8TysfYz2eGHH+5je10ul9fb9yLu9bRVbOMmij3+ofI5oeMvzxv2fZSfyefMmaNy/fv39/G5556rcvLcY/tiyyvI8nXymsU+zwknnKByI0aM8LEt7REaG2nGTbGNmXzPNbIdem/sHCFLyrVu3VrlTj/9dB/bcgrWGWec4eNVq1YFHyvFPUfZazI5n+ZqruGbuAAAAAAAAACQYdzEBQAAAAAAAIAM4yYuAAAAAAAAAGTYFtXEDf173FpNoe1sjYhQXQi5na2H0q9fPx/bOjpjxoxR7YsuusjHtj5KiKwPY/sZ970ISbtd1iU5/jJn68fK2iO2bp+si1KjRg2Vu++++3z897//XeVkHR3LHmNZDyi0na0XI+t4ValSReXk65U1nZzT4/iSSy5Ruccff9zHjzzyiI9tHcSsSzPXJHm+uNvZejy9evXycdOmTVVO1ji2QvWw7XiKW3/bilujMqQizjWWPB62Rqisq/7yyy+rXP369X188sknq5z827fzQKhuun3sjjvuGNlvWWv1mWeeUbkVK1ZEbpereorFLM18Y+cGOW5CNa5srSx5HghtF7omsn2x1zoPP/ywj6+55hqVk2Ps4IMPVrlvvvlms/0sje1PlGKfb7J6njrxxBMjc9ttt52Pv/322+D+5PWUPU+Frsk5T4XlYr6JO2+E/m7tc1avXt3Htk7gcccd52P7OcuOhbFjx0buM1QfWo4bW3dTPjZJfWap2MdNlLTnqVAte/s3HNouVN8y9Bk49Fj7uzR33323j+3nrMGDB/t4ypQpKvfpp5+6KPI12hrAofrcxUaOgdBxlK/TXkPEnW8suT87b7zzzjs+vuWWW1SucePGm+3X5tgauZI8ritXrlQ5+ZsXoVr1aZXH+SZXY0a2Q7/3YJ/z/PPP93GfPn1UTo5ze/6yvzci54kkxzp0/SLbaX/DIsmY4Zu4AAAAAAAAAJBh3MQFAAAAAAAAgAxLXE4hallOaFlVSNwSCaXlZPuGG25QOfkVarlUyDnnPvjgA9WWX79O8tX6uK8jyVe25WPTvr9ZEdX/tMffLjGWpS9kaQPn9LIr+7V8uVTQfvVdLvsobYmjXR4bxX69X5ZJsMtBfvzxRx/bpWlyOeR1112ncnKJQGgpZNYVcq4JsWOmU6dOkY+dPHmyjz/++GOVs2M2tBwpNPZDSyPTvsaKPteE2HlhwoQJPp49e7bKtWzZ0sd2CaOcQ+wyxdDyHDtnyOe1ywYvvfRSH3/55ZcqJ8dfkvNQ3Dm62MdNWV7bxH0fk1yD2HF75ZVX+tguafzhhx98/P7776ucnbfiYtwU9jxl55vddtstMidLN5VWTiFJCY2423Geys64sdutWbMm8rH77LOPj+U16OYccsghPrYliRYsWOBje70u35d8fAYr9nGT6+sb+56GSgLl4rNskmNqS2Z88sknPn700UdV7sEHH/TxwIEDVU6WQLPXU6HxVp5EjY/Q8Q8dq9D8HnofbZlKeYztPLHTTjv5uLTyTKG+ybJnF154YeT+uS7W0sw1ofIoaa8t2rRpo9qnnXaaj+24kPsYNWqUytn7fWll6X4f38QFAAAAAAAAgAzjJi4AAAAAAAAAZBg3cQEAAAAAAAAgwxLXxI2q3Req/ZCrujohsp7oRRddpHKyZsaYMWNU7s4778zJ/kM1DeO+N0m2KzZpxk1IqI5P9+7dVe7ggw/2sa0TKGtTyrqAzjl3xhln+LhVq1Yq16JFC9WePn26jxcvXqxyX3zxhY9tTZalS5f6eN68eSq3cuVKH9saKdWqVdvs4zb32GKVlbmmYcOGqn344Yf72NZ3kvW4V69eHXxe+fpKq7kctR1zza/leq6xdZtkTbXDDjtM5WRN2lNOOUXlGjVq5GNbm132zdYktH/fch9vvPGGyoXGXOj4h6Qdb8WmLOebXLyPpT1HnTp1fGxrDY4bN87HU6dOVTl5jrTbhTBuCnue2n777VW7cePGkY9dsmSJj5cvXx583rTnDc5TYVkZN6FraXmd6ZxzrVu39rH9rQB7LSR/i2TatGkq99xzz/lY/i6Ec/o1JqnHzHxTtuNGstezaY+NzcnfCZF1Tp1zbsiQIT6W9cCd079VMH78eJWT9wfseKvo4yYkV9vJa4qePXuq3L333utj+TneuV+f7+Tru+KKK1Tuscce87Gt8S+vbwrxOayYlOWYkcf3rrvuUjl7XpLkuaZPnz4qZ++NyL/9Yv0MzjdxAQAAAAAAACDDuIkLAAAAAAAAABmWs3IKdrmGbIeWd4e2s197l6pWraral112mY/tsh65dMyWWrBLV6tUqRK7r1Lo69Vxv7Jttwvtr9iWzKcZN6Hjb99H+dV7WdrAtu12oXEqS2/YpWNyeY7dR5Ilp2m/Qm/LQkQJjcusSzNmQksUkow1ucSmdu3aKifLbvznP/9ROblUq7TlGaF+h+YMOU5Dz5lkCR1zTfy5Rr4f8tzinHMDBgzw8SOPPKJy//jHP3wsl7Y759z777/v47lz56rcyJEjVXv+/Pk+Tvs3HZoHt912W5WTcw3jJv61TWhJaZL3Ju4cbvdXo0YN1ZZj9bvvvlO5t956y8f2/CVLf6xfvz5y/7ZvoX5X9HGTq/NUKPfee+/5uE2bNio3evRoH9tr3tC5MMn5Rj4P56lfK8txE9UP55zbbrvtfNy0aVOVq1u3buTz2LIIspyC/Zwl55gnnnhC5dIuY2W+yf24SfLZPer5S9tOzi/OhT/byP7YMjByn7vssovKybKJ3bp1U7l169ZF7q+YPz9ZWb2+kecfez19+umn+7h9+/Yqd/7556v2bbfd5uMpU6aoXNzXYV+T/BwmY+fin9/i9iOL8v1ZSt5XsaVM+vbt6+MDDzwwcn8LFy5UuRNPPDFWX5zTxzT0Gdy+XtnXsr7fxzdxAQAAAAAAACDDuIkLAAAAAAAAABnGTVwAAAAAAAAAyLBKSeq8VKpUqUTWr5F1GmztBylJPR75PLaehaw1cfXVV6vcVVdd5eNq1aqpnKy/dNddd6ncjTfeqNqyPk5p9TSkqPfFufB7I4W2s8+xadOm8SUlJR1id7AMpR03ltwuNG5s/Vp5HEPb2b7I2ky2bpOtjyP7ZmsKhurqyNoqoXp3obqcccfNpk2bXElJSfw3vAylHTOh9ym0XZLaX6E6PknmDDkWkswZodpTzDW5HzehOcPOC/Jv3+ZCtcbkdqH9OafPn3a8xT3+oXqCSf5OGDfp5pu024Wupey5xbZlzbghQ4aonKzDHKp7m2R+C413xk3uz1N2vpG/G/HKK6+oXL9+/Xw8bdo0lQvV9Lfnu1Bf4x7/EMbNr7ZT7bjb2eMm2fONrIles2ZNlWvSpImP27VrF/mczjn3/PPP+9jWHY07x+Xjuohxk93tbN4eR/na7ee8+vXr+/jiiy9WuS5duvh41apVKifrbYbOfZvBuImQ9p5P6LO7vZ6xzyM/ryc5T4XGW6hOqcRn8P8TOtfYYyh/y8jmnn32WR937txZ5eRxOffcc1VO1li39bVD9WuTHPssfQbnm7gAAAAAAAAAkGHcxAUAAAAAAACADNum9Ido8ivA8qvJScoyhJbLxN13jRo1VE4uAbPP+fbbb/v45ptvVrnVq1erdtyvz4f6FlpKkETa9ymLCjlubC4f739oibPdn9zWLkcNbRe337kab1mTZswkOS5Rjwv1wzm9RGNL3vu0c4Z8jcw1v5aLuUYKLQ8K/T2HtgtJ8v6nHTeh5ynEuTyLCjnfpN3OCo032x48eLCPFy5cGPlYu395TZR2frMYN7k/T9m5SF7bnnXWWSr33Xff+bi066XQNXHcY8d56tcKOW5C7Hby+mb58uUqt2LFCh9PmjRJ5exS6LVr1+a8b6Fc2hINxaYYzlNbMk5Defk53853y5Yt8/GDDz6oci+99JKP7XyW9jqt2JTlfJOLz+5JSjklOd/k4vwW2n8xzzf5/iz1008/+bhVq1Yq16xZMx/bYyT3L/+2nQt/PreSjEUpS5/B+SYuAAAAAAAAAGQYN3EBAAAAAAAAIMO4iQsAAAAAAAAAGZa4Jm6UJHVVQvUkQmRdjO+//17lZsyY4WNbT6JPnz4+XrlyZXAfsm+hejBJxK3HUl5rm4YUYtzk4n3ckmMTqiUTt46XlY9xWixyUfdzS/Yn26F6SltyHPLxGplrynbchPaXZB+hvsbdLm3NMMaNlov5uxDbOefcggULfJykDmDac00+6oYVk7Keb2Rb1sB1Lv3xz9VxY76JlqVxk6Qm7bp161Q7F5/zSttn3P0xbornPBUSGptr1qzxsbwf4Jxz06dPj9xu22239XFF+yzlXGGOf9y5YEve77TXxSGcpzYvV69P3quz1yiffvqpj1etWqVyX3/9tY9lnfak0p4zsvQZnG/iAgAAAAAAAECGcRMXAAAAAAAAADIscTmFqK8f2yUw8nG2vIFkt5Ptn3/+WeXkErCnn35a5T7//HMf//jjjyq3cOHCyP1bW2+9tY/t15tlzn6dWr5Gu518TUm+hh1aVhR6T7Mo1+MmtFwjyXuTdsmXXILjnHMbNmzI+Xayb3LsORd/vIX+nrIuzZix4m4XWl5qnz/ucj97zCpXrqzaGzdu3Ox2dh9p5wzmmvRzjdwuV0sD5f6TzG12HIXGqtyHHW/y759x82uFnG/SbpdknrLiLqEPzXdxzzXOxT9HWYybdMc/1K8k29n5Ju5jQ0uRmW9+LSvjJvT3nmTchPoW2uc22+iPoaF5IxfjzWLcZGu7uGPOjpu0n2/kdZJ9jtD9gC352ygLhRw3Sa6n0352t/uPu4zdfgbP9+cwqSLMNUnONfJY2HI8ffv29XH16tVVTu4j7uehze1fKtbP4HwTFwAAAAAAAAAyjJu4AAAAAAAAAJBh3MQFAAAAAAAAgAyrlKSmX6VKlRY752bnrztIoGlJSUmDsu5EHIybzGDMIA3GDdJg3CANxg3SYNwgDcYN0mDcICnGDNKIHDeJbuICAAAAAAAAAAqLcgoAAAAAAAAAkGHcxAUAAAAAAACADOMmLgAAAAAAAABkGDdxAQAAAAAAACDDuIkLAAAAAAAAABnGTVwAAAAAAAAAyDBu4gIAAAAAAABAhnETFwAAAAAAAAAyjJu4AAAAAAAAAJBh/w+ttlXVpP7DywAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 1800x288 with 20 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# obtain one batch of test images\n",
    "dataiter = iter(test_loader)\n",
    "images, labels = dataiter.next()\n",
    "\n",
    "# get sample outputs\n",
    "output = model(images)\n",
    "# prep images for display\n",
    "images = images.numpy()\n",
    "\n",
    "# output is resized into a batch of iages\n",
    "output = output.view(batch_size, 1, 28, 28)\n",
    "# use detach when it's an output that requires_grad\n",
    "output = output.detach().numpy()\n",
    "\n",
    "# plot the first ten input images and then reconstructed images\n",
    "fig, axes = plt.subplots(nrows=2, ncols=10, sharex=True, sharey=True, figsize=(25,4))\n",
    "\n",
    "# input images on top row, reconstructions on bottom\n",
    "for images, row in zip([images, output], axes):\n",
    "    for img, ax in zip(images, row):\n",
    "        ax.imshow(np.squeeze(img), cmap='gray')\n",
    "        ax.get_xaxis().set_visible(False)\n",
    "        ax.get_yaxis().set_visible(False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}