dmnist_demo.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%%sh\n",
    "pip install torch -q \n",
    "pip install ipywidgets -q"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import os\n",
    "import sys\n",
    "import data\n",
    "\n",
    "import torch\n",
    "import torch.nn as nn\n",
    "import torch.nn.functional as F\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "from collections import OrderedDict"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "def uint82bin(n, count=8):\n",
    "    \"\"\"returns the binary of integer n, count refers to amount of bits\"\"\"\n",
    "    return ''.join([str((n >> y) & 1) for y in range(count - 1, -1, -1)])\n",
    "\n",
    "def labelcolormap(N):\n",
    "    cmap = np.zeros((N, 3), dtype=np.uint8)\n",
    "    for i in range(N):\n",
    "        r, g, b = 0, 0, 0\n",
    "        id = i + 1  # let's give 0 a color\n",
    "        for j in range(7):\n",
    "            str_id = uint82bin(id)\n",
    "            r = r ^ (np.uint8(str_id[-1]) << (7 - j))\n",
    "            g = g ^ (np.uint8(str_id[-2]) << (7 - j))\n",
    "            b = b ^ (np.uint8(str_id[-3]) << (7 - j))\n",
    "            id = id >> 3\n",
    "        cmap[i, 0] = r\n",
    "        cmap[i, 1] = g\n",
    "        cmap[i, 2] = b\n",
    "\n",
    "    \n",
    "    cmap[0] = (255, 255, 255)\n",
    "\n",
    "    return cmap"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "IPython.OutputArea.prototype._should_scroll = function(lines) {\n",
       "    return false;\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%%javascript\n",
    "IPython.OutputArea.prototype._should_scroll = function(lines) {\n",
    "    return false;\n",
    "}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\n",
       "\n",
       "\n",
       "mpl.get_websocket_type = function() {\n",
       "    if (typeof(WebSocket) !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof(MozWebSocket) !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert('Your browser does not have WebSocket support. ' +\n",
       "              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
       "              'Firefox 4 and 5 are also supported but you ' +\n",
       "              'have to enable WebSockets in about:config.');\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n",
       "    this.id = figure_id;\n",
       "\n",
       "    this.ws = websocket;\n",
       "\n",
       "    this.supports_binary = (this.ws.binaryType != undefined);\n",
       "\n",
       "    if (!this.supports_binary) {\n",
       "        var warnings = document.getElementById(\"mpl-warnings\");\n",
       "        if (warnings) {\n",
       "            warnings.style.display = 'block';\n",
       "            warnings.textContent = (\n",
       "                \"This browser does not support binary websocket messages. \" +\n",
       "                    \"Performance may be slow.\");\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.imageObj = new Image();\n",
       "\n",
       "    this.context = undefined;\n",
       "    this.message = undefined;\n",
       "    this.canvas = undefined;\n",
       "    this.rubberband_canvas = undefined;\n",
       "    this.rubberband_context = undefined;\n",
       "    this.format_dropdown = undefined;\n",
       "\n",
       "    this.image_mode = 'full';\n",
       "\n",
       "    this.root = $('<div/>');\n",
       "    this._root_extra_style(this.root)\n",
       "    this.root.attr('style', 'display: inline-block');\n",
       "\n",
       "    $(parent_element).append(this.root);\n",
       "\n",
       "    this._init_header(this);\n",
       "    this._init_canvas(this);\n",
       "    this._init_toolbar(this);\n",
       "\n",
       "    var fig = this;\n",
       "\n",
       "    this.waiting = false;\n",
       "\n",
       "    this.ws.onopen =  function () {\n",
       "            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n",
       "            fig.send_message(\"send_image_mode\", {});\n",
       "            if (mpl.ratio != 1) {\n",
       "                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n",
       "            }\n",
       "            fig.send_message(\"refresh\", {});\n",
       "        }\n",
       "\n",
       "    this.imageObj.onload = function() {\n",
       "            if (fig.image_mode == 'full') {\n",
       "                // Full images could contain transparency (where diff images\n",
       "                // almost always do), so we need to clear the canvas so that\n",
       "                // there is no ghosting.\n",
       "                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "            }\n",
       "            fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "        };\n",
       "\n",
       "    this.imageObj.onunload = function() {\n",
       "        fig.ws.close();\n",
       "    }\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_header = function() {\n",
       "    var titlebar = $(\n",
       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
       "        'ui-helper-clearfix\"/>');\n",
       "    var titletext = $(\n",
       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
       "        'text-align: center; padding: 3px;\"/>');\n",
       "    titlebar.append(titletext)\n",
       "    this.root.append(titlebar);\n",
       "    this.header = titletext[0];\n",
       "}\n",
       "\n",
       "\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = $('<div/>');\n",
       "\n",
       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
       "\n",
       "    function canvas_keyboard_event(event) {\n",
       "        return fig.key_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
       "    this.canvas_div = canvas_div\n",
       "    this._canvas_extra_style(canvas_div)\n",
       "    this.root.append(canvas_div);\n",
       "\n",
       "    var canvas = $('<canvas/>');\n",
       "    canvas.addClass('mpl-canvas');\n",
       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
       "\n",
       "    this.canvas = canvas[0];\n",
       "    this.context = canvas[0].getContext(\"2d\");\n",
       "\n",
       "    var backingStore = this.context.backingStorePixelRatio ||\n",
       "\tthis.context.webkitBackingStorePixelRatio ||\n",
       "\tthis.context.mozBackingStorePixelRatio ||\n",
       "\tthis.context.msBackingStorePixelRatio ||\n",
       "\tthis.context.oBackingStorePixelRatio ||\n",
       "\tthis.context.backingStorePixelRatio || 1;\n",
       "\n",
       "    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband = $('<canvas/>');\n",
       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
       "\n",
       "    var pass_mouse_events = true;\n",
       "\n",
       "    canvas_div.resizable({\n",
       "        start: function(event, ui) {\n",
       "            pass_mouse_events = false;\n",
       "        },\n",
       "        resize: function(event, ui) {\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "        stop: function(event, ui) {\n",
       "            pass_mouse_events = true;\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "    });\n",
       "\n",
       "    function mouse_event_fn(event) {\n",
       "        if (pass_mouse_events)\n",
       "            return fig.mouse_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
       "\n",
       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
       "\n",
       "    canvas_div.on(\"wheel\", function (event) {\n",
       "        event = event.originalEvent;\n",
       "        event['data'] = 'scroll'\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        mouse_event_fn(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.append(canvas);\n",
       "    canvas_div.append(rubberband);\n",
       "\n",
       "    this.rubberband = rubberband;\n",
       "    this.rubberband_canvas = rubberband[0];\n",
       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
       "\n",
       "    this._resize_canvas = function(width, height) {\n",
       "        // Keep the size of the canvas, canvas container, and rubber band\n",
       "        // canvas in synch.\n",
       "        canvas_div.css('width', width)\n",
       "        canvas_div.css('height', height)\n",
       "\n",
       "        canvas.attr('width', width * mpl.ratio);\n",
       "        canvas.attr('height', height * mpl.ratio);\n",
       "        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n",
       "\n",
       "        rubberband.attr('width', width);\n",
       "        rubberband.attr('height', height);\n",
       "    }\n",
       "\n",
       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
       "    // upon first draw.\n",
       "    this._resize_canvas(600, 600);\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n",
       "        return false;\n",
       "    });\n",
       "\n",
       "    function set_focus () {\n",
       "        canvas.focus();\n",
       "        canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    window.setTimeout(set_focus, 100);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>');\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            // put a spacer in here.\n",
       "            continue;\n",
       "        }\n",
       "        var button = $('<button/>');\n",
       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
       "                        'ui-button-icon-only');\n",
       "        button.attr('role', 'button');\n",
       "        button.attr('aria-disabled', 'false');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "\n",
       "        var icon_img = $('<span/>');\n",
       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
       "        icon_img.addClass(image);\n",
       "        icon_img.addClass('ui-corner-all');\n",
       "\n",
       "        var tooltip_span = $('<span/>');\n",
       "        tooltip_span.addClass('ui-button-text');\n",
       "        tooltip_span.html(tooltip);\n",
       "\n",
       "        button.append(icon_img);\n",
       "        button.append(tooltip_span);\n",
       "\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    var fmt_picker_span = $('<span/>');\n",
       "\n",
       "    var fmt_picker = $('<select/>');\n",
       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
       "    fmt_picker_span.append(fmt_picker);\n",
       "    nav_element.append(fmt_picker_span);\n",
       "    this.format_dropdown = fmt_picker[0];\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = $(\n",
       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
       "        fmt_picker.append(option);\n",
       "    }\n",
       "\n",
       "    // Add hover states to the ui-buttons\n",
       "    $( \".ui-button\" ).hover(\n",
       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
       "    );\n",
       "\n",
       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n",
       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
       "    // which will in turn request a refresh of the image.\n",
       "    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_message = function(type, properties) {\n",
       "    properties['type'] = type;\n",
       "    properties['figure_id'] = this.id;\n",
       "    this.ws.send(JSON.stringify(properties));\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_draw_message = function() {\n",
       "    if (!this.waiting) {\n",
       "        this.waiting = true;\n",
       "        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n",
       "    }\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    var format_dropdown = fig.format_dropdown;\n",
       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
       "    fig.ondownload(fig, format);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_resize = function(fig, msg) {\n",
       "    var size = msg['size'];\n",
       "    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n",
       "        fig._resize_canvas(size[0], size[1]);\n",
       "        fig.send_message(\"refresh\", {});\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
       "    var x0 = msg['x0'] / mpl.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
       "    var x1 = msg['x1'] / mpl.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n",
       "    x0 = Math.floor(x0) + 0.5;\n",
       "    y0 = Math.floor(y0) + 0.5;\n",
       "    x1 = Math.floor(x1) + 0.5;\n",
       "    y1 = Math.floor(y1) + 0.5;\n",
       "    var min_x = Math.min(x0, x1);\n",
       "    var min_y = Math.min(y0, y1);\n",
       "    var width = Math.abs(x1 - x0);\n",
       "    var height = Math.abs(y1 - y0);\n",
       "\n",
       "    fig.rubberband_context.clearRect(\n",
       "        0, 0, fig.canvas.width / mpl.ratio, fig.canvas.height / mpl.ratio);\n",
       "\n",
       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_figure_label = function(fig, msg) {\n",
       "    // Updates the figure title.\n",
       "    fig.header.textContent = msg['label'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_cursor = function(fig, msg) {\n",
       "    var cursor = msg['cursor'];\n",
       "    switch(cursor)\n",
       "    {\n",
       "    case 0:\n",
       "        cursor = 'pointer';\n",
       "        break;\n",
       "    case 1:\n",
       "        cursor = 'default';\n",
       "        break;\n",
       "    case 2:\n",
       "        cursor = 'crosshair';\n",
       "        break;\n",
       "    case 3:\n",
       "        cursor = 'move';\n",
       "        break;\n",
       "    }\n",
       "    fig.rubberband_canvas.style.cursor = cursor;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_message = function(fig, msg) {\n",
       "    fig.message.textContent = msg['message'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_draw = function(fig, msg) {\n",
       "    // Request the server to send over a new figure.\n",
       "    fig.send_draw_message();\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_image_mode = function(fig, msg) {\n",
       "    fig.image_mode = msg['mode'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message(\"ack\", {});\n",
       "}\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function(fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            /* FIXME: We get \"Resource interpreted as Image but\n",
       "             * transferred with MIME type text/plain:\" errors on\n",
       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "             * to be part of the websocket stream */\n",
       "            evt.data.type = \"image/png\";\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src);\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                evt.data);\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n",
       "            fig.imageObj.src = evt.data;\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        var msg = JSON.parse(evt.data);\n",
       "        var msg_type = msg['type'];\n",
       "\n",
       "        // Call the  \"handle_{type}\" callback, which takes\n",
       "        // the figure and JSON message as its only arguments.\n",
       "        try {\n",
       "            var callback = fig[\"handle_\" + msg_type];\n",
       "        } catch (e) {\n",
       "            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        if (callback) {\n",
       "            try {\n",
       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
       "                callback(fig, msg);\n",
       "            } catch (e) {\n",
       "                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "}\n",
       "\n",
       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
       "mpl.findpos = function(e) {\n",
       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
       "    var targ;\n",
       "    if (!e)\n",
       "        e = window.event;\n",
       "    if (e.target)\n",
       "        targ = e.target;\n",
       "    else if (e.srcElement)\n",
       "        targ = e.srcElement;\n",
       "    if (targ.nodeType == 3) // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "\n",
       "    // jQuery normalizes the pageX and pageY\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    // offset() returns the position of the element relative to the document\n",
       "    var x = e.pageX - $(targ).offset().left;\n",
       "    var y = e.pageY - $(targ).offset().top;\n",
       "\n",
       "    return {\"x\": x, \"y\": y};\n",
       "};\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * http://stackoverflow.com/a/24161582/3208463\n",
       " */\n",
       "function simpleKeys (original) {\n",
       "  return Object.keys(original).reduce(function (obj, key) {\n",
       "    if (typeof original[key] !== 'object')\n",
       "        obj[key] = original[key]\n",
       "    return obj;\n",
       "  }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function(event, name) {\n",
       "    var canvas_pos = mpl.findpos(event)\n",
       "\n",
       "    if (name === 'button_press')\n",
       "    {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * mpl.ratio;\n",
       "    var y = canvas_pos.y * mpl.ratio;\n",
       "\n",
       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
       "                             step: event.step,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "\n",
       "    /* This prevents the web browser from automatically changing to\n",
       "     * the text insertion cursor when the button is pressed.  We want\n",
       "     * to control all of the cursor setting manually through the\n",
       "     * 'cursor' event from matplotlib */\n",
       "    event.preventDefault();\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    // Handle any extra behaviour associated with a key event\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.key_event = function(event, name) {\n",
       "\n",
       "    // Prevent repeat events\n",
       "    if (name == 'key_press')\n",
       "    {\n",
       "        if (event.which === this._key)\n",
       "            return;\n",
       "        else\n",
       "            this._key = event.which;\n",
       "    }\n",
       "    if (name == 'key_release')\n",
       "        this._key = null;\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which != 17)\n",
       "        value += \"ctrl+\";\n",
       "    if (event.altKey && event.which != 18)\n",
       "        value += \"alt+\";\n",
       "    if (event.shiftKey && event.which != 16)\n",
       "        value += \"shift+\";\n",
       "\n",
       "    value += 'k';\n",
       "    value += event.which.toString();\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, {key: value,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onclick = function(name) {\n",
       "    if (name == 'download') {\n",
       "        this.handle_save(this, null);\n",
       "    } else {\n",
       "        this.send_message(\"toolbar_button\", {name: name});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n",
       "    this.message.textContent = tooltip;\n",
       "};\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
       "\n",
       "mpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n",
       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
       "    // object with the appropriate methods. Currently this is a non binary\n",
       "    // socket, so there is still some room for performance tuning.\n",
       "    var ws = {};\n",
       "\n",
       "    ws.close = function() {\n",
       "        comm.close()\n",
       "    };\n",
       "    ws.send = function(m) {\n",
       "        //console.log('sending', m);\n",
       "        comm.send(m);\n",
       "    };\n",
       "    // Register the callback with on_msg.\n",
       "    comm.on_msg(function(msg) {\n",
       "        //console.log('receiving', msg['content']['data'], msg);\n",
       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(msg['content']['data'])\n",
       "    });\n",
       "    return ws;\n",
       "}\n",
       "\n",
       "mpl.mpl_figure_comm = function(comm, msg) {\n",
       "    // This is the function which gets called when the mpl process\n",
       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
       "\n",
       "    var id = msg.content.data.id;\n",
       "    // Get hold of the div created by the display call when the Comm\n",
       "    // socket was opened in Python.\n",
       "    var element = $(\"#\" + id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm)\n",
       "\n",
       "    function ondownload(figure, format) {\n",
       "        window.open(figure.imageObj.src);\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy,\n",
       "                           ondownload,\n",
       "                           element.get(0));\n",
       "\n",
       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
       "    // web socket which is closed, not our websocket->open comm proxy.\n",
       "    ws_proxy.onopen();\n",
       "\n",
       "    fig.parent_element = element.get(0);\n",
       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
       "    if (!fig.cell_info) {\n",
       "        console.error(\"Failed to find cell for figure\", id, fig);\n",
       "        return;\n",
       "    }\n",
       "\n",
       "    var output_index = fig.cell_info[2]\n",
       "    var cell = fig.cell_info[0];\n",
       "\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function(fig, msg) {\n",
       "    var width = fig.canvas.width/mpl.ratio\n",
       "    fig.root.unbind('remove')\n",
       "\n",
       "    // Update the output cell to use the data from the current canvas.\n",
       "    fig.push_to_output();\n",
       "    var dataURL = fig.canvas.toDataURL();\n",
       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
       "    // the notebook keyboard shortcuts fail.\n",
       "    IPython.keyboard_manager.enable()\n",
       "    $(fig.parent_element).html('<img src=\"' + dataURL + '\" width=\"' + width + '\">');\n",
       "    fig.close_ws(fig, msg);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.close_ws = function(fig, msg){\n",
       "    fig.send_message('closing', msg);\n",
       "    // fig.ws.close()\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.push_to_output = function(remove_interactive) {\n",
       "    // Turn the data on the canvas into data in the output cell.\n",
       "    var width = this.canvas.width/mpl.ratio\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Tell IPython that the notebook contents must change.\n",
       "    IPython.notebook.set_dirty(true);\n",
       "    this.send_message(\"ack\", {});\n",
       "    var fig = this;\n",
       "    // Wait a second, then push the new image to the DOM so\n",
       "    // that it is saved nicely (might be nice to debounce this).\n",
       "    setTimeout(function () { fig.push_to_output() }, 1000);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>');\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items){\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) { continue; };\n",
       "\n",
       "        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n",
       "    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n",
       "    button.click(function (evt) { fig.handle_close(fig, {}); } );\n",
       "    button.mouseover('Stop Interaction', toolbar_mouse_event);\n",
       "    buttongrp.append(button);\n",
       "    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n",
       "    titlebar.prepend(buttongrp);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(el){\n",
       "    var fig = this\n",
       "    el.on(\"remove\", function(){\n",
       "\tfig.close_ws(fig, {});\n",
       "    });\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(el){\n",
       "    // this is important to make the div 'focusable\n",
       "    el.attr('tabindex', 0)\n",
       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
       "    // off when our div gets focus\n",
       "\n",
       "    // location in version 3\n",
       "    if (IPython.notebook.keyboard_manager) {\n",
       "        IPython.notebook.keyboard_manager.register_events(el);\n",
       "    }\n",
       "    else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    var manager = IPython.notebook.keyboard_manager;\n",
       "    if (!manager)\n",
       "        manager = IPython.keyboard_manager;\n",
       "\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which == 13) {\n",
       "        this.canvas_div.blur();\n",
       "        // select the cell after this one\n",
       "        var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n",
       "        IPython.notebook.select(index + 1);\n",
       "    }\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.find_output_cell = function(html_output) {\n",
       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
       "    // IPython event is triggered only after the cells have been serialised, which for\n",
       "    // our purposes (turning an active figure into a static one), is too late.\n",
       "    var cells = IPython.notebook.get_cells();\n",
       "    var ncells = cells.length;\n",
       "    for (var i=0; i<ncells; i++) {\n",
       "        var cell = cells[i];\n",
       "        if (cell.cell_type === 'code'){\n",
       "            for (var j=0; j<cell.output_area.outputs.length; j++) {\n",
       "                var data = cell.output_area.outputs[j];\n",
       "                if (data.data) {\n",
       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
       "                    data = data.data;\n",
       "                }\n",
       "                if (data['text/html'] == html_output) {\n",
       "                    return [cell, data, j];\n",
       "                }\n",
       "            }\n",
       "        }\n",
       "    }\n",
       "}\n",
       "\n",
       "// Register the function which deals with the matplotlib target/channel.\n",
       "// The kernel may be null if the page has been refreshed.\n",
       "if (IPython.notebook.kernel != null) {\n",
       "    IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAyAAAAH0CAYAAADFQEl4AAAgAElEQVR4Xu2dB7gkRdW/q2du2rzLLrvsknMUhI9FMRHFCAbwD5gjJhDMqOiKCmIOmEAwCyjgp58YCAKCii4LggiiRMlh2Rxumun/U3PvVNXszuyt6pme7up653l4qJ17qurUe2pm6td9uiqK4zgWvCAAAQhAAAIQgAAEIAABCHSBQIQA6QJluoAABCAAAQhAAAIQgAAEagQQIEwECEAAAhCAAAQgAAEIQKBrBBAgXUNNRxCAAAQgAAEIQAACEIAAAoQ5AAEIQAACEIAABCAAAQh0jQACpGuo6QgCEIAABCAAAQhAAAIQQIAwByAAAQhAAAIQgAAEIACBrhFAgHQNNR1BAAIQgAAEIAABCEAAAggQ5gAEIAABCEAAAhCAAAQg0DUCCJCuoaYjCEAAAhCAAAQgAAEIQAABwhyAAAQgAAEIQAACEIAABLpGAAHSNdR0BAEIQAACEIAABCAAAQggQJgDEIAABCAAAQhAAAIQgEDXCCBAuoaajiAAAQhAAAIQgAAEIAABBAhzAAIQgAAEIAABCEAAAhDoGgEESNdQ0xEEIAABCEAAAhCAAAQggABhDkAAAhCAAAQgAAEIQAACXSOAAOkaajqCAAQgAAEIQAACEIAABBAgzAEIQAACEIAABCAAAQhAoGsEECBdQ01HEIAABCAAAQhAAAIQgAAChDkAAQhAAAIQgAAEIAABCHSNAAKka6jpCAIQgAAEIAABCEAAAhBAgDAHIAABCEAAAhCAAAQgAIGuEUCAdA01HUEAAhCAAAQgAAEIQAACCBDmAAQgAAEIQAACEIAABCDQNQIIkK6hpiMIQAACEIAABCAAAQhAAAHCHIAABCAAAQhAAAIQgAAEukYAAdI11HQEAQhAAAIQgAAEIAABCCBAmAMQgAAEIAABCEAAAhCAQNcIIEC6hpqOIAABCEAAAhCAAAQgAAEECHMAAhCAAAQgAAEIQAACEOgaAQRI11DTEQQgAAEIQAACEIAABCCAAGEOQAACEIAABCAAAQhAAAJdI4AA6RpqOoIABCAAAQhAAAIQgAAEECDMAQhAAAIQgAAEIAABCECgawQQIF1DTUcQgAAEIAABCEAAAhCAAAKEOQABCEAAAhCAAAQgAAEIdI0AAqRrqOkIAhCAAAQgAAEIQAACEECAMAcgAAEIQAACEIAABCAAga4RQIB0DTUdQQACEIAABCAAAQhAAAIIEOYABCAAAQhAAAIQgAAEINA1AgiQrqGmIwhAAAIQgAAEIAABCEAAAcIcsCbwt7/9TZx11lnipptuEo8//riYOXOm2GGHHcSznvUs8aUvfcm6nbrhJz/5SXH66aeLOI6d69pU+Mtf/iKuuOIKccopp9R87eTrjW98o7j22mvF/fff33aznWyrbWcsGnjsscfEZz7zGfHb3/5WPProo2Lu3Lni8MMPF4sWLRLbbLONRQuYQAACEOg+gb/+9a+136o//elP4qmnnhKbbbaZeM5zniPe//73iwMPPDCRQ2eeeabYY489xMtf/vJE9V0qPfLII+Lcc8+t9fX0pz/dpSq2EMgdAQRI7kKST4d+85vfiKOOOkocfPDB4m1ve5uYP39+bfG5ZMkScdFFF4mHHnrI2fG0BcgXv/hF8cEPflDcd999YrvttnP2b1MVOikaOtlWRwfZpLGhoSHxtKc9TSxfvrwmHuUP77///e+a+Ojp6RH/+te/xLRp09J2g/YhAAEIOBE4++yzaxejDjjgAPGud71LbLvttuKBBx4Q3/zmN8XixYvF1772NXHiiSc6tSmNp06dKo455hjxgx/8wLmuawX5e7tw4ULx/e9/X8jfDV4Q8JkAAsTn6HXR94MOOkg8/PDD4s4776wtNM1XtVoVpVLJ2RsEyBgynwTIVVddJZ7//OeL8847T7zlLW9RMb/wwgvFq1/9avGLX/xCvOIVr3CeC1SAAAQgkBaBP//5z+J5z3ueePGLXyz+93//t+E3bHR0tPadJe/oXnfddeLZz362kxsIECdcGENAEUCAMBmsCOy11161Kz3yFrbN62c/+5n46le/Kv7xj3+IKIpqX+oyfWvfffdV1VsJEJu6shGZEiZTgWSq1dq1a8WWW24pjjzyyFq/9bY39PWaa66p3cWRL9t+5JWtz372s7V0K5ly9pGPfERcffXV1ilYF1xwgZBX32677bZavzvttJM46aST1AK+mQCRV+Wkf1LwybHJfl/3uteJ9773vaK3t1cN6+9//7v4+Mc/XruCt3LlSjF79uwa43POOUdstdVWNbuLL75YyLtBsq2RkRGxxRZb1Bh873vfswllg80f//jHWl3ZprzqV3/97ne/q/24yx/xF73oRc7tUgECEIBAWgRe+tKXit///ve17/D696LZ14MPPii233772nfXr3/965YXhTb8zZK/bRu+5MU6mZ4rfzfe9KY31dKA5W/Ar371KyHvIB966KG1uy3yO73+knfo5ffqhndR6r9Vsj353yGHHLJRf/Lus/SLFwR8I4AA8S1iGfkr067kVW+5cH7Na14j9ttvv4aFsOmWzIk97bTTal++r3zlK8Xw8LD4whe+UBMjcqEs03bkq5kAsa17+eWX18TG7rvvLt73vvfVnj2QPy7yy15ejZcpYZ///OdrC395VV6mjMmX7Hv69OnCtp/6j8jLXvYy8da3vrW2yJd+yx8SeddnomdAPvGJT4hPf/rTNQ7HH3+8mDFjhvjnP/9ZS2H61Kc+VfOpmQCRY9pll11qP4p9fX3i1ltvFWeccUZtzHXhIIWJTCOQNh/60IfEvHnzhHw+Q4qs97znPTU2N9xwQ038HXvsseL1r3+9GBgYEP/9739rAupHP/qRCpv8oZPiYqLnceTVwmc+85licHCwlgYg+/jPf/5TE1OShxSopkDKaLrSLQQgAIEagUqlUvvOl6mjm7qA9oxnPKP23bxq1ara91mzZ/w2/M2S7UlBIYWBvBAkX7Iv+TtT/+3Yeuuta3eNX/WqVwkpdORvY39/f+33sP5soo0AkX7J3zL5uyrbeMlLXlLrTwqqZqKK8EMg7wQQIHmPUE78kw/syQff5MN78iUXmTIXVS6IZd6svDsiX/ILVl7Zeec73ym+/vWvK+/XrFkjdt5559ptcHllv5kAcakr7yLIl/zBkIvqZq9Wz4DY9iNTy+SPh7xjIHNv61e75AJejmXBggWbFCDy2RNpd9xxx4mf/OQnLSM5UQqW9EP+J4WV/PF58sknxaxZs2qbAey///7il7/8pZACqdlLPnD5gQ98QKxYsaImflq9DjvssJoAkQJjotfq1atrIlReKay/pIC59NJLaw918oIABCCQFwJywxT5HS6/h+V3aKuX/Lv8bZL28oKOjQCRbbVKwaoLEJneJYVD/SXv2MuLQvLu/cc+9rHa2zYCRNrxDEheZhV+dIIAAqQTFANqQ34B/uEPf6h9Ecov6KVLl9a+PG+88UYxZ86c2l0SebdE/nvDXTpe+9rX1q7Oyy94+drwapJtXXnFfdddd63dxZDpUK1erQSIbT/ygWp5JUu2I3dJMV9ywS3vfmzqDojcreTtb397LUVsUzusNBMgMrVK3lqXucvLli1r6FtedZNX6+TdGHn3Q+5CJR+ulOKufnepXkHmNMuUgCOOOKJ2VU/+8MlUtaQvmcIlxY4UfvLujoyDFFryx1Q+fC7vrGxK6CTtl3oQgAAEkhCwFSDyLvHPf/5z8cQTT9Q2L+mUALnkkkvE0Ucf3eC6/M2UF9HkM3UIkCRRpU4RCCBAihDFjMYgF6Mf/vCHxVe+8pXaF7ZMeZJpQvL2cKuXTNORt8SbCRDbunJRLrdOlClE8rkIVwFi24+82/Pc5z63dvdCXvE3X/JqmRQCmxIg9X7kTivyTkqr14YCRNpLISEX9/KZD/ljJe/yyPS1d7/73TURV88Nls+VyH5k6plM65KpZlIAyhjUU6Fk7rG8GyW5ydSxPffcs3blTaaEub6+853v1O5uSYEp777UX/fee6/Ycccda6JSCideEIAABPJAIM0ULDm+ie6AyN+RDR9sl2ms8rtYXmhCgORhluBDFgQQIFlQL1Cf8iq8zGOVD+/JB5Dlw8/veMc7hLzqI59PaPaqL1w3vANiW7fdOyC2/WR1B0SKhZNPPrkmbkyG9Ts3pgCp85XPbsicYnnbXz6ELx+aP/XUUxvwyx88KZrk3+QzNBPdmWkWOxlb6Yd8rmfDnc/kHTCZCy0fUOcFAQhAIC8EJnoIXT4zKC/01B9Cl99z8sKN3GrefMl0Y7lBiPms3EQCxOYOyG677VbbPGTDFDG5+Yv8XpV3Y+SLFKy8zCj86AQBBEgnKAbQhvwirj/IbQ5XLmhlepFM75ELU7lols89yKvyMo92U68NBYhLXXn7Wj6TIVOB5AN9zV7yAXT5MPYdd9xRe1i6/rLtRz53IR/uk896yCv+rs+AyH6kn3J7WvOB7w193fAOSN1vyVzmLsuX/MGTV83kXZBmAsRsUz4fIh96lOkEzV7ygXaZHid/SOV++C4v+eC8vMNRTwOr162LQpkKJu+I8YIABCCQFwL1bXjlg9tyG95yuaxck3dI5PON8gLa9ddfXztYV+7YKO8Sy4P/5OYe8iUvusi7x3fffXeDAJE7D8qDWOvPNtYbnugZELk5ST1b4IUvfGHt+cnbb79d+SW/U2V/8u5JXYDIO9577723+Na3vlW7E80LAj4TQID4HL0u+i6/9ORiXD50Lq/WyMX5LbfcUjtVVj6ULK+my11G5EteYZfPB0hRIr9Y5YJY5uHKxfOUKVNqB9jJV7NdsGzr1nfBkl/QMk1J7oIlU5fk+z/96U9r7de3LZTPYbzhDW+opSTJtCb5rIJtP+eff35t9yv53INMbZIPcyfZBUtuWVvfBUsKIvnsTJ3DhgJEbpe7zz771H4IpYiTO059+9vfrom7u+66SwmQyy67rPZDJH885YP/UqTIhx1lmpR8/kT6K+Mgr+7Jh8xl/KT/cgtImRYgb/9LfvJl+xC6/JGUc0HGUf54Sp4y/Uo+jyNjLB+Ml+/xggAEIJAnAuZBhPJORv03Q16IkVu6yzvHcpdH+ZLPtcnvMbn4l+nF8jtY3pmWvzHyb+YdEJkOK++Wywtw8iKd/H2Rdc1dsOQzePVdsKSwkb9F8o51fdMO+Zsln5GUokI+LyI3OpEpzXLnK7kbYl2ArFu3rnZHRO5CKZ+7k3df5AUy+R8vCPhGAAHiW8Qy8ldeTZe3pOWdAHllXqbzyC9b+YCzfBDcvMMgXZS2cqErF6TSVl7Jl7tmyVvbcrHbSoDY1pV28iq8vCIvxY/8gZALbHla+5e//GVF6aMf/aj44Q9/WNueVoom8+6BjY+yISlCPve5z9V+FORtepnaJHeMavaQYrPw/PjHP1bngMhDHOUdInlnpn6SbbOH0KW4kAt8ecq4vMIm76LI7R5likB9DPJvUsRIYSev1MnteuUPn7yrIQWXfMkT7OUPrPyxk7tnyXQ5mQInYyafo6m/bLfhlfbyCqDkLq8WyrkgrxDKu2BS7Gz4EHxG05VuIQABCGxEQP5myItm8rtL7uwoBYD8HpQ7BW64UYg820j+fsjvWflbJ7dGl9+h8jvXFCDyjrJ8Nk9e0JECodk5IPI3QO4aKH8LZZqq/G2s7+QonZTtyc1O5MUj+V0uU6+kAKlfpKoLEGl70UUX1d6/5557auc6cQ4IE91XAggQXyOH3xCAAAQgAAEI5JJA/Q7Ihht25NJZnIJABgQQIBlAp0sIQAACEIAABIpLAAFS3Ngyss4QQIB0hiOtQAACEIAABCAAgRoBBAgTAQKbJoAAYYZAAAIQgAAEIAABCEAAAl0jgADpGmo6ggAEIAABCEAAAhCAAAQQIMwBCEAAAhCAAAQgAAEIQKBrBBAgXUNNRxCAAAQgAAEIQAACEIBARwWIPGdB7mEtD+KpnxoNYghAAAIQ8JuAPKdAHjgqDzwrlUp+DwbvIQABCEAgcwIdFSDyxOWtt94680HhAAQgAAEIdJ7Agw8+WDvw0/XFBSlXYthDAAIQyDGBSPsWV+NEjnZUgKxcubJ20rL8kZo+fXoih6gEAQhAAAL5IrBq1araxaUVK1aIGTNmODuHAHFGRgUIQAAC+SWQNwEif6Tkj5MUIgiQ/M4bPIMABCDgQqDd73YEiAttbCEAAQjknAACJOcBwj0IQAACBSCAAClAEBkCBCAAgRQIyGcEk7w6moLV7o9UkgFQBwIQgAAE0iXQ7nc7d0DSjQ+tQwACEMiKAAIkK/L0CwEIQKDgBBAgBQ8ww4MABCCQkAACJCE4qkEAAhCAwKYJIEBCmSFGYrdRFAlTLEKhxjghEDIBBEjI0WfsEIAABFIkgABJEW6umkaA5CocOAMBDwggQDwIEi5CAAIQ8JEAAsTHqCXxGQGShBp1IBAyAQRIyNFn7BCAAARSJIAASREuTUMAAhDwmAACxOPg4ToEIACBPBNAgOQ5OvgGAQhAIDsCCJDs2NMzBCAAgUITQIAUOrwMDgIQgEBiAgiQxOioCAEIQAACmyKAACny/DC3u9LjbNgEa8LhJzuIbMJmMYAABHJPAAGS+xDhIAQgAAE/CSBA/IybndcIEDtOWEEAAs0IIECYFxCAAAQgkAoBBEgqWHPSKAIkJ4HADQh4SQAB4mXYcBoCEIBA/gkgQPIfIzyEAAQgkAUBBEgW1OkTAhCAQAAEECABBJkhQgACEEhAAAGSABpVIAABCEBgYgIIkIkZYQEBCEAgRAIIkBCjzpghAAEIdIEAAqQLkOkCAhCAgIcEECAeBg2XIQABCPhAAAHiQ5TwEQIQgED3CSBAus+cHiEAAQgEQQABEkSYGSQEIAABZwIIEGdkVIAABCAAARsCCBAbSthAAAIQCI8AAiS8mDNiCEAAAl0hgADpCmY6gQAEIOAdAQSIdyHDYQhAAAJ+EECA+BGnTXvZ/MDBxjpxEQbKGCAAgS4SQIB0ETZdQQACEAiJAAKkCNFGgBQhiowBAnkjgADJW0TwBwIQgEBBCCBAihBIBEgRosgYIJA3AgiQvEUEfyAAAQgUhAACZNOBNJf2JcO0bPzBTG7S5eaiIBbaInbMimrejxBmT45NFmQWMwwIQCANAgiQNKjSJgQgAAEICAQIAoSPAQQgAIFmBBAgzAsIQAACEEiFAAIEAZLKxKJRCEDAewIIEO9DyAAgAAEI5JMAAiSfccErCEAAAlkTQIBkHQH6hwAEIFBQAgiQggaWYUEAAhBokwACpE2AVIcABCAAgeYEECDMDAhAAAIQaEYAAcK8gAAEIACBVAggQFLBWvhGS1GLXb5abO3F7lyFnxIMsIAEECAFDCpDggAEIJAHAgiQPETBPx8QIP7FDI8h4EoAAeJKDHsIQAACELAigACxwoTRBgQQIEwJCBSfAAKk+DFmhBCAAAQyIYAAyQS7l51GJX0UY19Zl830qspoRY3NzNKKq82TsFodrkjKlpdTBKcLRgABUrCAMhwIQAACeSGAAMlLJPLvBwIk/zHCQwh0kgACpJM0aQsCEIAABBQBBAiTwZYAAsSWFHYQKAYBBEgx4sgoIAABCOSOAAIkdyHxwqGykY4ljA2xzFQrc58scyHTKr2KtCsvQo+TARFAgAQUbIYKAQhAoJsEECDdpF2cvhAgxYklI4FAKwIIEOYGBCAAAQikQgABkgrWwjeKACl8iBkgBAQCxGESRKc3PxzJoYmgTONF3PQOKuAMFgIbEAhBgFSrVeJuQaBkplVNYN/TU1YW5tmD1YreBctsgl8aiwBgAoGcEUCAOAQEAeIASwiBAHHjhTUEikYAAVK0iCYfDwIkOTtqQqCIBBAgDlFFgDjAQoC4wcIaAgUkgAApYFATDgkBkhAc1SBQUAIIEIfAmgKEq/utwdU5wchhcmEKgQISQIAUMKgJh+QiQMwteRu6a5HuRgpWwqBQDQIZEkCAOMBHgNjBQoDYccIKAkUngAApeoTtx4cAsWeFJQRCIIAAcYgyAsQOFgLEjhNWECg6AQRI0SNsPz4EiD0rLCEQAgEEiEOUESB2sBAgdpywgkDRCSBA8hhhM2Gpezs7OgmQckmDq2p/ky5Y8hgFfIJA6ASSfp6jOGnNJsTb/ZHqVhARIHakESB2nLCCQNEJtPvdHkXdWyAnjYV/2/AiQJLGmnoQgEDnCCSVEQgQzrhoOQsRIJ37gNISBHwmgADJY/QQIHmMCj5BIDQCCBCHiHMHxA4WAsSOE1YQKDoBBEg9wq32aerOHZ5Y6P5de4yNFKgoMlKjzBMCzYkcGWM17F1SsEo9ParF2Nj5yiwX/bPj8/jMOcYOZT5HMl3fESAOfBEgdrAQIHacsIJA0QkgQBAgdQIIkKJ/2vX4ECDhxLqdkSJAHOghQOxgIUDsOGEFgaITQIAgQBAgRf+Ubzw+BEh4MU8yYgSIAzUEiB0sBIgdJ6wgUHQCCJB0ItzdpzgseourxkCNcqRTqZzugPT3qfbikVFdbnEQYTqUaRUCEEiTAALEgS4CxA4WAsSOE1YQKDoBBEg6EbaQBB3s2KI3BEgHedMUBMIggABxiDMCxA4WAsSOE1YQKDoBBEg6EbaQBB3s2KI3BEgHedMUBMIggABxiDMCxA4WAsSOE1YQKDoBBEh2ER5cN6I6X/z3e1X59jseVeXlQ9q/uXMmqX9sPWuaKj/vWduq8sCUyaoclVrtp9VcsLikYPXM1P1XV69XfVYrOh0rO7L0DAEIdIIAAsSBIgLEDhYCxI4TVhAoOgEESHYRRoBkx56eIQCBiQkgQCZmpK/4nK6v+MQcRNiSHALEYVJhCoECE0CAZBdcBEh27OkZAhCYmAACZGJGCBAHRtIUAeIIDHMIFJQAAiRZYM2DA4UwLnyNH/73yKNLVcP3PrZClf/8B51edcujOhXq3ofXKptVK/UOU6UBbTNtQNsMCJ2bdcjhc1TdYw/YQZV321OnZpXKxgGFxqGHpu8uKVil6dNVP/Ea7VdcrSQDSi0IQCB3BBAgDiEhBcsOFgLEjhNWECg6AQRIsggjQBAgyWYOtSDgDwEEiEOsECB2sBAgdpywgkDRCSBAkkUYAYIASTZzqAUBfwggQBxihQCxg4UAseOEFQSKTgABkjDCxkZSS5euUo2c8/2ba+Wr71+m3ls+qFOqVj05Vb0f9fSrcv8UncY1UtHl+fN0R2uX6R2mesqG/chy1c6e2+oUqDNPOFC9v+W2m7cYqG7HJQVL9GrfxajezUs0bPebkG2Caq32+2rWlLkHWIKuklfh+PHk7KiZCQEEiAN2BIgdLASIHSesIFB0AgiQhBFGgGhwCBC7SYQAseOEVW4IIEAcQoEAsYOFALHjhBUEik4AAZIwwggQBIjr1EGAuBLDPmMCCBCHACBA7GAhQOw4YQWBohNAgGwc4YYf3fFdrcas9Apy6cN6l6uTvvlP1cg/7h+z6ZukDw3sjXW94apWLlNG16l6o/osPzHSq3esmtar6z62qlfZT51q7Lwl9Pu95aeUzXuPnaHKx71kYdOpHJV0X04pWFFZt9eQdpVNglNUx2HRvYVJ0T/2jA8CVgQQIFaYxn8eOAfEihYCxAoTRhAoPAEECAKkTgABUviPOwOEgBMBBIgDLu6A2MFCgNhxwgoCRSeAAEGAIECK/ilnfBBIRgAB4sANAWIHCwFixwkrCBSdAAJkLMLmD22lWlVhLwldvvcf96n3P3jh/ar8r8cHVDkeP5OvV+jdqOIRnfTT368PEHzurj2q3rQ+ndI0ME2nV03u1ylS66M1yv62B/TOU/c8rA8iHF6nd8o66f/pNt989LNU3SjSbUYqd0kIpzsgRjqaaDjYsDOfmIZdrVr8o2S8Xxr/R6WiWZtbJTdk0nXGRVqBQOEJIEAcQowAsYOFALHjhBUEik4AAYIAqc9xBEjRP+2MDwJuBBAgDrwQIHawECB2nLCCQNEJIEDGIjxa0XcszDsE6wf1g+JvOOlaNR1uXTtF373o0ed8zCqN3QI5cA9992H/ffQdit23m6nq7bj9fFXu7dcPkjfcXDAmoLkYWLvW8OsLNymrO28dVuVjDteVP/muQ9Q/opLxALnRvm8CpGzcAimPn4tSqeg7Vsbz/sK8A5J0UVX07wLGB4ENCST9rERx0ppNYtDuj1S3wooAsSONALHjhBUEik6g3e92M4Unr6yqRkpVKx8RIP6lYCFA8vqJw6+iEEgqIxAgi9hsr9WHAAFSlK8HxgGB9gggQMb4IUAQIO19kqgNgeIRQIA4xJQ7IHawECB2nLCCQNEJhC1AjAeWR3Tqzroh/SD3J75xpZoCl92o066qozplapd5+uHwj756i5r9Ac/cVdXr6e1XZfPcjVapVjZzbnCNTsF629lLVJWbrl+pyq8/Wp8Dcuqbn6ubNR48N/vKUwqWDQNs0iVgPvvfsA9Ak271pyddn/LWeisuNmdO+nCJHAHiMOMQIHawECB2nLCCQNEJIEDGIhwjQHK1C1bRP3c+jA8BMnGUECDNGZGCRQpWy08PAmTiLxYsIBACAQQIAqQ+z7kDEsIn3n6MCJCJWSFAECCKAHdAJv7ASAsEiB0nrCBQdAJhCxAd3ZH1+nyO085arP7wy3/q3bEmlfXuUZtNW6VsPn3Czqr8rIVj5cg8pKJTk8jI2bjqjzeqVk/72TJVHnlK74L1puNnq/dPfPmB2ouOpGB1alC0A4HiEbBJwfJh1KRgOUQJAWIHCwFixwkrCBSdAAJkLMIIENeH0Iv+yWB8EEhOAAGSVLo0Yd7uj1TyMLrVRIDY8UKA2HHCCgJFJ9Dud3tRtuFFgCBAiv5ZZ3zdI4AAQYB0b7Z51hMCxLOA4S4EUiKAABkDe9uNtyjCJ35Hp1c9ulSnXS2Yu17ZfOK4zVT50MP2Syk6jc3++lc67eozVzys/rhi+djOW/I1a0inY33towvU+z9LBE4AACAASURBVM/Y7+m6sRaJ627PgHRlyHQCAQhkSCCpjOAhdB5CbzltESAZfqLpGgI5IoAAQYDUpyMCJEcfTFyBQA4IIEAcgkAKlh0sBIgdJ6wgUHQCCBAECAKk6J9yxgeBZAQQIA7cECB2sBAgdpywgkDRCSBAxiJ8x99uVaE+6fv6ML+7H9OHEh6yl96G6pun6LSrKXNmdXSaLH9qtWrv7Et02tV1t+vDB+9/ZLKymSL0Tl3HvVDbfPgNL1E2pd6eCX3kDsiEiDAoPIGJjlyUAHw4QrAzgUKAOHBEgNjBQoDYccIKAkUngABBgHAHpOifcsZnTwABYrJCgNjPHHW+RU2j8gxIS3IIEIdJhSkECkwAAYIAQYAU+APO0BwJIEAQII5Tpm7OHRA7cAgQO05YQaDoBIIWIEYmxeAavcPVyWf9SYX9ssX9qrzPbrrCN07aTr2/wy7bJpom1Ypu73e/uU218e1r7lLlux+eqcqlaECVp/Y8pcrveqXu/g2vOFL9w/UwRFKwEoWRSoUiYCNAGpbohRr9hoPhDohDeBEgdrAQIHacsIJA0QkgQMYijADhHJCif9YZnw0BBEiDvIqTPe/CNrykYLX8tCFAbL6IsIFA8QkgQBAg9VnOHZDif94Z4UQEECAIkInmSIu/cwfEDhwCxI4TVhAoOoGgBUiL4F5+zZ/VXz5xzhpVfqo6TZUP3VfvVHXKS7ZW7++82461spFdJVau0LZX/FmnV/31dn1o4C3/1btUrR/UaV/D60uq7Xmz16ryW1+pfTn28APV+6UefXCi69xFgLgSw754BGwESLK7Aj6yIgXLIWoIEDtYCBA7TlhBoOgEECAbRxgBUvRZz/gg0IoAAsQkgwBx+KQgQOxgIUDsOGEFgaITQIAgQOoEuANS9E8745uYAAIEATLxLGlqgQCxA4cAseOEFQSKTgABsnGEh9bqHbG+89NrlMG3r9DpUJURffjgVlMfVDb77zZ2cOFDy/RC5p6R+ervK5ea6VUj6v2+kj7wcMFmq9T7e+2ofXn9UTvrfvbct/lvYGSzgGo+qxEgRf+0Mz4IuBHgDogDLwSIHSwEiB0nrCBQdAIIEARInQACpOifdsYHATcCCBAHXggQO1gIEDtOWEGg6AQQIAgQBEjRP+WMDwLJCCBAHLghQOxgIUDsOGEFgaITQIBsOsLV0aoyOO2zV6ryL/7aq8qD66eo8vC64Vq5NNlIr9J/FlvobCyx29yKqve8PWeo8l676B2u9tpbp10N9OqdstKYl9wBSYMqbULAXwIIEIfYIUDsYCFA7DhhBYGiE0CAIEC4A1L0Tznjg0AyAggQB24IEDtYCBA7TlhBoOgEECAIEARI0T/ljA8CyQggQBy4IUDsYCFA7DhhBYGiE0CA2Ee4Mqp3qvrJbxeriv98XB8QWB6aXHv/ObtOV3/fefu5qjx9ht4Fa95sbZOH9Kc8+GAfDSwhAIG0CSBAHAgjQOxgIUDsOGEFgaITQIDYRxgBYs8KSwhAwH8CCBCHGCJA7GAhQOw4YQWBohNAgNhHGAFizwpLCEDAfwIIEIcYIkDsYCFA7DhhBYGiE0CAFD3C9uMjBcueFZYQCIEAAsQhyggQO1gIEDtOWEGg6AQQIEWPsP34ECD2rLCEQAgEECAOUUaA2MFCgNhxwgoCRSeAACl6hO3HhwCxZ4UlBEIggABxiDICxA5WXgVIFJ1uNwCsNkkgjhdBCAJWBBAgVpiCMMqrADmiqncfq4hIxcI8ljEef7f+/w0DFgv9l9how7QvN7yv/6J7FMJsJzLsSw0dmn2ZfzBb2nhKmb6YZbNWq3KrcV9RKgcxdxlkOgQQIA5cESB2sBAgdpx8tUKA+Bq57vuNAOk+87z2iADRy/tGoaEjhgDJ6+zFrzQIIEAcqCJA7GAhQOw4+WqFAPE1ct33GwHSfeZ57REBggDJ69zEr2wIIEAcuCNA7GD5IEBYRNvFsm5lpq/Bzo1dyNYIkJCj3zj2vAqQF7RIwTKTi+opSI3pUnp8PUYKlpnGFRnvmylVrdOhWqVm6Z6rDVhNe1PgdHbeNfqu+yQFq7OcQ2sNAeIQcQSIHSwEiB0nn6wQID5FKz++IkDyE4usPUGANBcIjc9jIECynqf03z0CCBAH1ggQO1gIEDtOPlkhQHyKVn58RYDkJxZZe4IAQYBkPQfpP18EECAO8UCA2MFCgNhx8skKAeJTtPLjKwIkP7HI2pO8CpAXVisKjZneZKZPTZSCZdabbIAeMco9xq5Wa4z3zZQtc7erVjtSmbtzmT4mj2/zPbFa7ZRl9nM5u2Alx05NgQBxmAQIEDtYCBA7Tj5ZIUB8ilZ+fEWA5CcWWXuCANGSAgGS9Wyk/zwQQIA4RAEBYgcLAWLHyScrBIhP0cqPrwiQ/MQia08QIAiQrOcg/eeLAALEIR6mAHGoFqxpvKjV8UXZIGERnZw77JKzC7lmCAIkivpViON4KORwezl2cxcs8xDBxmP9Nk7CanWAoHng4FCsfwMrsU7UKt2zWLHa8i/3qfLKbeap8uDe26tyZeZWqjxqpD2NRnqvLjP1q/HgQpewNN/nqzEdS9tcXkrek4tX2BaTAALEIa4IEAdYQggEiBuvPFsjQPIcnfz6hgDJb2zwbIwAAsScCQgQPhfdI4AAcWCNAHGAhQBxg5VzawRIzgOUU/cQIDkNDG4pAggQBAgfh2wIIECy4U6vGRBgEZ0cOuySswu5ZggCJOT4FmHsL2pxEGGzFKyqsZOVuRuVmV41NKrT8Cb9/VqFaOvLb1blva69RpV3uec2Va706VafnDtVvb961gJVXrr3Abq8xXaqfOfxh6pyNGcXVe4d9zk2DkUUxjjcY0gKljszajQjgABhXgRDgEV08lDDLjm7kGsiQEKOvh9jR4C4xgkB4koM++YEECDMjGAIsIhOHmrYJWcXck0ESMjR92PsL6vq0zrWC/1QdbOnIep3E+TI1hvDm7r8KfWvvT/2aVXe7epLVHn2U9pmpE/3c98Ws5TNNsu0LzNGdA+9Rl/9Ze3ZqhWD6i8//+wiVb7nXaeqsnmnRjfT/FDE5k+ANN4vMR9Iv4JzQPyY5Dn1EgGS08DgVucJsIhOzhR2ydmFXBMBEnL0/Rg7AqSZLBEtE7YQIH7Max+8RID4ECV87AgBFtHJMcIuObuQayJAQo6+H2NHgCBA/JipxfMSAVK8mDKiFgRYRCefGrBLzi7kmgiQkKPvx9iPrg4rR1cKfa6GmbpUP8HDPPVi2LhHsP2PLlBtHH3yCap8/876XI+bDnuten/ZXjup8qM76fLs9cYD7INLlc28y25V5WMu+YYqTx/WKVhf+NTn1ftLT3yPKtdHFAl9Dol53kljlCZOzTLtLycFy49JnlMvESA5DQxudZ4Ai+jkTGGXnF3INREgIUffj7EjQMw4IUD8mLXF8BIBUow4MgoLAiyiLSC1MIFdcnYh10SAhBx9P8aOAEGA+DFTi+clAqR4MWVELKI7PgcQIB1HGkSDoQkQcxehVjsKmQ/xBjEJcj7IY6s67WmZ0IlXOhlLqOSlkpF2pWsJMevBJ9QoZy/W533c9kJ9ZoeYMlvZTDXO4TDTuszdrkaNvsrr16q6/++Fz1TlXf91nyqfde4PVXndy49R5Xrqlel7bOz2ZYbHnJut5q9pzy5YOZ/cOXcPAZLzAOFe5wiwiE7OEnbJ2YVcEwEyFn1zYYcAydcnAgGi44EAydfcLLo3CJCiR5jx6StB0emqHMd6z3QQTUwAATIxIyw2JoAAQYDk/XOBAEGA5H2OFtU/BEhRI8u4NiLAIjr5pIBdcnYh1wxNgIQca1/HfoyRgrVS6CSoZilI5h0CvaeUTNFqfo9rwEi1anWcn9mOiHU7UaQ92Pzxfyu8L37xUaq84MnHVfmsX/xClUv7H6bKo6JSKzemYJkJZskjxy5YydlRUwgECLMgGAIsopOHGnbJ2YVcEwEScvT9GDsCJHmcECDJ2VETAcIcCIgAi+jkwYZdcnYh10SAhBx9P8aOAEkeJwRIcnbURIAwBwIiwCI6ebBhl5xdyDURICFH34+xH2WkYA0au2CVjPSpemJU680E9F/MXa0iIzXL3Hmq8X3NyTwscMTof79ffFsZveSED6lyNEmnUn3vtM+p9x96+zu1zXipp8EXm9g03y7BPMSQXbBsOGLTigApWMyNYAiwiE4eatglZxdyTQRIyNH3Y+wIkJbLw6Z/QID4Ma998BIB4kOU8LEjBFhEJ8cIu+TsQq6JAAk5+n6MHQGCAPFjphbPSwRI8WLKiFoQYBGdfGrALjm7kGsiQEKOvh9jf2l1UDk6aqRgRQ07WNVNmp+UYe6CpY8yFGK0AYFOzmqVgmXaD4xov15z5CtUSzvccL0qPzljqio/sP9uqnzp+b9X5d4Z/bVybOSGNaaXuZ3+Ye4O9vuSOVo/4o2X+SGAAMlPLPAkZQIsopMDhl1ydiHXRICEHH0/xo4AQYD4MVOL5yUCpHgxZUTcAen4HECAdBxpEA0iQIIIs9eDRIAgQLyewB47jwDxOHi47kaARbQbL9MadsnZhVwTARJy9P0Y++HVEeVoJHSeknlUX/3BazN1yhydeZhgj9D/ajy40NwfS9c2bUaMurtfcoEyetW7PqDKt+67vSrfd8Sxqrz7/52nystnbKfKSz78pVp57cKd1Htxj/Yljpofftgqemb61u9KnTnQ0I+ZgpedJoAA6TRR2sstARbRyUMDu+TsQq6JAAk5+n6MHQGCAPFjphbPSwRI8WLKiFoQYBGdfGrALjm7kGsiQEKOvh9jR4AgQPyYqcXzEgFSvJgyIgRIx+cAAqTjSINoEAESRJi9HuTzjRSskpGC1Tiosb2fYuMwP3M3KDONqiQqRlXTykh7MiyGKuvUv3a7/UJV3vr0G7TV3EmqfNUZ79fvT9pclbc/4+OqPHmt3k/riecfV3t/+aH7qL9H/bo9YaZgGeNv9WSIMHYH4yR0r6d+5s4jQDIPAQ50iwCL6OSkYZecXcg1ESAhR9+PsSNAjDsgCBA/Jm1BvESAFCSQDGNiAiyiJ2bUygJ2ydmFXBMBEnL0/Rg7AgQB4sdMLZ6XCJDixZQRtSDAIjr51IBdcnYh10SAhBx9P8b+/KpOVzJ3eGruvZmYpC3qu2TJd0rGTlbmYYbD5t2F6rCqvMdF31PleUv+qMpXv+d9qjw8fw9V7umfrMql0TXa5l8r9Pvr1qty/95za+XSJKOe4YuZJCa9r79sUrDYBcuPOZ5XLxEgeY0MfnWcAIvo5Ehhl5xdyDURICFH34+xI0DMOCFA/Ji1xfASAVKMODIKCwIsoi0gtTCBXXJ2IddEgIQcfT/GjgBBgPgxU4vnJQKkeDFlRCyiOz4HECAdRxpEgwiQIMLs9SBtBEg9HakxXUknKZkpWKbNiLFj1ECsU722vuiritl+n/mpKl946gtUuee1Z6jySKQP/KsaKV5zjPLQkD4Asc84GbE00FtrpxJpz8wjEc09u3qMSJrjMG3WGDuBXcVBhF7P/aydR4BkHQH67xoBFtHJUcMuObuQayJAQo6+H2NHgOg4IUD8mLNF8RIBUpRIMo4JCbCInhBRSwPYJWcXck0ESMjR92PsCBAEiB8ztXheIkCKF1NG1IKAz4voQ4zb5wRYE7gmbr4rDYzyQQABko844EVrAkcYu2CZu1aZNZqnYGkLI+NJDBkpStuOrFJG8667QpUPfsd7VXlyuU+VT138W1WeNn13Ve41Uq1+PHiMen92/3TtxKjhRVWnbInSeMJVZNzfMHOwjDQxEQ0aw9ZtxNGAen8o1uVqPKLenxQNqXLjz5Xhi8Emjs3ELu27edhjvcE4NhLCYm1bNcbcP+lHTHPPCCBAPAsY7iYngABJzi6vNREgeY3MmF8IkHzHB++EQIAYi3sECB+JLhJAgHQRNl1lSwABki3/NHpHgKRBtXNtIkA6x5KW0iGAAEGApDOzaHUiAgiQiQjx98IQKIoAyWLRnSd2ZjpaFiwK84HowkAQIF2ATBdtETiiqlOBGne52rjZVjtc9VZ1+tFOf75cVdxh8aWqPOOPT6rynjder8pLFyxQ5c8uvkOVe8tGCpSRavS8J76kbE7f4m5VrlSaHx0YlzZOIGtINYua14uMdKmoIafKSPWKx3bYkq/ISBOLhT5oUTSUzcfcdd1qbJw/YqTV1t9tuVA1bHsHzm9rHlC5+wQQIN1nTo8ZEcjTItoVQdaL7jyxy5qFa+xCtkeAhBx9P8aOAEGA+DFTi+clAqR4MWVELQjkaRHtGqSsF915Ypc1C9fYhWyPAAk5+n6MHQGCAPFjphbPSwRI8WLKiBAgHZ8DCJCOIw2iQQRIEGH2epA2z4DUB2imLq01UpT2vf4GxeB5xx2lypuX9OGDUUnvHtVb0mlMg716J6vLjjpC1X3w1fup8lte9SlVvnmmxv3hGw9V/6iM6N20IiOtKm62g2JD2pVuz9xUMKrvniWEKBn2kZEuNbpurar8t99PUeVH/7S1Kj/5gG5/xu73qn8ct0inY+k0MWFQFSIaTx9r3OxQJ8KZ7/f3n+f1PAzReQRIiFEPdMx5WkS7hiDrq/55Ypc1C9fYhWyPAAk5+n6MHQGCAPFjphbPSwRI8WLKiFoQyNMi2jVIWS+688QuaxausQvZHgEScvT9GDsCBAHix0wtnpcIkOLFlBEhQDo+BxAgHUcaRIMIkCDC7PUg7Z4BGXtOwnxaYtA4wPCQb/5UMdjuonNVeWDNMlVeMXVzVV43a5Yq9/SuUeV7Zvar8pObb6vK7/j971T5X2t12tMJdx+p3q+MmilY5sF+9YMIWxza2vAIiN6NqqR2zxKiMqJ9vOoinQN21w1zVf+Lr52qytturtPK1s57VL1/8LH/UeUjX6FT0irGxlqR0Lt/RePutDpv1jygsL//u17PwxCdR4CEGPVAx5ynRbRrCLK+6p8ndlmzcI1dyPYIkJCj78fYESBGnOorfvncBwLEjwnssZcIEI+Dh+tuBPK0iHbzXIisF915Ypc1C9fYhWyPAAk5+n6MHQGCAPFjphbPSwRI8WLKiFoQyNMi2jVIWS+688QuaxausQvZHgEScvT9GPsLjFQqeZzepl5mtlLVOHhv0sgqVW1wUOcTlZfp1KXB8iRls27yZFWe/fA9qjy9R+8MNVrVvhz0rncrm6EnblPl99/xSlWu6A23GodQb6bVzldVnfJUMnbnqlTWqXauPn87Vf6/r26pyvO2GVTlLfd9RJWPOn65HtNO2rE+PWwhSjrdrFIxD4PUaWD1DbzMFCwzBqRg+fEZa+UlAsTv+OG9A4E8LaId3K6ZZr3ozhO7rFm4xi5kewRIyNH3Y+wIEASIHzO1eF4iQIoXU0bUgkCeFtGuQcp60Z0ndlmzcI1dyPYIkJCj78fYESAIED9mavG8RIAUL6aMCAHS8TmAAOk40iAaRIAEEWavB2kKkNhIwTKTsVrsH6XGbf5dL+eFkaQlRI9hZB4OaGwAJapGrtGoUe65VadpbX673nHrJ6/RaU+VEaMD4/DBeHwgDcllxt8jodOfRoZH1Jiu+vYMVb7mPH0o4twD9GGCrzn5v8pm/l46xSzqGdJzItZEqhXj8EFz1lQbtsFSf6mfhdiKvxmvvj52wfLtg4gA8S1i+JuYQJ4W0a6DyPqqf57YZc3CNXYh2yNAQo6+H2NHgCBA/JipxfMSAVK8mDIi7oB0fA4gQDqONIgGESBBhNnrQZq7YJkD2fTj6JsacsPBGsqw8Y6K/lfj1X39r8i4GzNqnEASC/1Q91XxCar9akXfaYiFfpC73mLJaMPYbbfhrsvl5+qzPP5yzj6q7Wm73arKp3zvKVXumzRPldeNLNV3LkSvKvfXb2PUzlHRZ5iUjIfyq1U97nhUs1F+GndsTPLcAfH6oycQIH7HD+8dCORpEe3gds0066v+eWKXNQvX2IVsjwAJOfp+jB0Bohf/CBA/5mxRvESAFCWSjGNCAnlaRE/o7AYGWS+688QuaxausQvZHgEScvT9GDsCBAHix0wtnpcIkOLFlBG1IJCnRbRrkLJedOeJXdYsXGMXsj0CJOTo+zH2VgKk0fuxtKDISGNyH13z1Ky44VF1s1WdRqVLQhiPd4s/mClYo/pB7liY6VhjbZopYKWyfhj8jhuWqU5//b4XqPLQjPtV+d0X6QfP52ymH06vVPU5J+Y5JNURnYJ1x03al6Fe/ZD7M56uzxmJ+mepviqx+QD7mJ+ReRCICcNA2tt3nntIqJEpAQRIpvjpvJsE8rSIdh131ovuPLHLmoVr7EK2R4CEHH0/xo4AQYD4MVOL5yUCpHgxZUQtCORpEe0apKwX3XlilzUL19iFbI8ACTn6fowdAYIA8WOmFs9LBEjxYsqIECAdnwMIkI4jDaJBBEgQYfZ6kC+o6m1o5T5N9VdjwlQ9gcn8u7FbU0O9Vvtn6RQp0yIyUrDM3asiYyer4Rbnk1xupGAJ4yyNODZ2wRp3uVTSvVaFToW69PQ5asyLf72tKv/PoutV+fgjJ6vymtH1qnzzbTrV6vHLt1HvP3LHFqr83zsGVLk6Sfu1/776DJPjP6vPEylP0zaV8dhEVSPvqgVeUrD8+xgiQPyLGR4nJJCnRbTrELK+6p8ndlmzcI1dyPYIkJCj78fYESAIED9mavG8RIAUL6aMiDsgHZ8DCJCOIw2iQQRIEGH2epAIEASI1xPYY+cRIB4HD9fdCJiLaLea2VsfLD6pnLjWKGfhWRwvyqJb1Sd3QDLF79Q5AsQJF8YQcCJQGX2Hso9jnVYVx0ae0ni5VNbvrR1cper96O0LVbk0+yH9PXvCw6q8+B59gOCjP9tJvf/UPbuo8uiI3r1q2wP1YYUvf6tuU5SmKfvvfnhLVd7tlQ+o8qvfoFO84nJf7f1qpXn6mrm1FylYTlMnF8YIkFyEASe6QQAB0hnKCJDOcAyhFQRICFFmjFkRQIBo8giQrGZh8n4RIMnZUdMzAgiQzgQMAdIZjiG0ggAJIcqMMSsCCBAESFZzrxP9IkA6QZE2IJAyAdKONGBYpDzZOtg8AqSDMGkKAhsQqAy/Vb0TmztxNaRgjZmUevROUo8v1Vvv/vBNz1VtTJ6kU7OWTdHl5XfqFKhd95ik7Oce9IQq73vQSlXear7uq6+/X71fqeg0se9/Rqdy3fxLvZvWyRffrOx33XHs+ZQR45TDUmSeRIgA8flDgQDxOXr4HgwBFt0IEB8nOwLEx6jhsy8EECAIEF/majM/ESA+Rw/fgyGAAEGA+DjZESA+Rg2ffSGAAEGA+DJXESA+RwrfgyaAAEGA+PgBQID4GDV89oXA6Oibtavm4YNV4xjF8WJkpGAtX6PTpb77Jr0LVjSq6x12ot69avv5K1Q/M7bXu12VB/Qhg1FJH1ZYibVNKSqrupWRNar8y6/tqsp/Om+2Kh//s+tU+Rn7bFYrjxqpW+YBjcLYBqu377u+hA0/xwlwB4SpAAEPCCBAECAeTNONXESA+Bg1fPaFAAJEby2MAPFl1hqaOTaEsoP7UZxUujTppN0fKQe/MYWAlwQQIAgQHyduu9/tUWScZ+AjAHyGQIoEECAIkBSnV+pNJ5URCJDUQ0MHEGDR3WwOIMb8+WQgQPyJFZ76R2DU2AXL9N7cEUuM74gVlfRivRrrgwW//NptVNXpk7TNG755r3p/YGCKKsdxjyoPD+tDA8tlvTtWSVSUTaln7DBB+brvlkFV/ukH9lDl0ZFhVX7dBbeo8k5bzayVR0eNgwgbLkogQPybtdpjBIjP0cP3YAiw6EaM+TjZESA+Rg2ffSGAAEGA+DJXm/mJAPE5evgeDAEECALEx8mOAPExavjsCwEECALEl7mKAPE5UvgeNAEECALExw8AAsTHqOGzLwRGh9+iXTUel2p4tLe+C5YwDvAr65Sn3144XbVx3flju07J13u//x9Vnr+VTruKY72rlRA67UpE+jBBEa9WdR9daex8derT1Ps3/23skEH5etre96vyO895TJV7B8bar1R0Slerx8J6+87zJWz4OU6AOyBMBQh4QAABggDxYJpu5CICxMeo4bMvBBAgOlIIEF9mrfYTAeJfzPA4QAIIEASIj9MeAeJj1PDZFwIIEASIL3O1mZ8IEJ+jh+/BEECAIEB8nOwIEB+jhs++EBgdNg4iNHKTzBSs+lEL5obWpbL+18igTpf6+ofGdp2Sr8HyDFV+zUGjqjx3L51SVZ49ot7/2839qrz2Xn0o4TUXaG+WP7ZA2Rz5iuWqfOCrH1TlzfecpcqVeHzXLONgxSjSO2I1HkRICpYv87buJwLEt4jhb5AEECAIEB8nPgLEx6jhsy8EECDmQ+gIEF/mLQLEt0jhb9AEECAIEB8/AAgQH6OGz74QQIAgQHyZq8385A6Iz9HD92AIIEAQID5OdgSIj1HDZ18IjBi7YJm7QzVLwRJGDlappHfEKkV6h6nVS1eood/612mqfPcfdGpWdYreEUtM1o3e90e9I9bkeXpHrNnb6sMKd91Pt7/wkCdU+wPTdLqXMFKsqtWxdiLRsK+XDo8xpt5etzsgZkpai9Z9mQbe+okA8TZ0OB4SAQQIAsTH+Y4A8TFq+OwLAQSIjhQCxJdZq/1EgPgXMzwOkAACBAHi47RHgPgYNXz2hQACBAHiy1xt5icCxOfo4XswBBAgCBAfJzsCxMeo4bMvBEZGTlCuRiW9I5WIdYpVfQOpuKITjSIjH6vcq1OqooZT/nR79Z20NsUlFtrebF8IneIlhJG+ZTRWNUziWKdvqdQrM1+qhROu54CQgpX9LEeAZB8DPIDAhAQQIAiQCSdJDg0QIDkMCi4VhgACxLgD4ngSOgIk+48BAiT7GOABBCYkgABBgEw4SXJogADJYVBwqTAEECAIEJ8nMwLE5+jhezAEECAIEB8nOwLEx6jhsy8EhkffTATCqwAAIABJREFUqFwtmelTDTlNYylN1RHjmr9Z1BlPolTSO1Y1pmOZh//p9C5hpF0JoXfNEmJI+RXHOjWrUtHtxKN9yiaK9A5aotSwh9e4jeGw8edY6Pb6B77nS9jwc5wAAoSpAAEPCCBAECAeTNONXESA+Bg1fPaFAAIEAeLLXG3mJwLE5+jhezAEECAIEB8nOwLEx6jhsy8EECAIEF/mKgLE50jhe9AEECAIEB8/AAgQH6OGz74QGKn8P+VqJPq121UjTWp8R6w4Luu/GxlN1ar+RxRtXK9WyXxfjKp2KqN6+6pS2dh5a1TblIxuhVHXPBixYXeshr7GuzKzsiL9DzMxbFL/+b6EDT/rYbXZXq0JrShOeu+kSWPt/kgRTQgUnQACBAHi4xxv97u9MQ/dRwL4DIH0CCBANFsESHrzLK2Wk8oIBEhaEaFdCDQhgABBgPj4wUCA+Bg1fPaFAAIEAeLLXG3mJwLE5+jhezAETAESzKAtBnpNwlu4Fk1j0gECCJAOQAymCU5mCCbUDBQCQggECNMAAh4QQIA0DxICJN+TFwGS7/jkyzsESL7igTcQSJcAAiRdvrQOgY4QQIAgQDoykbrcCAKky8C97g4B4nX4cB4CjgQQII7AMIcABCAAATsCCBA7TljJnZaM7ZJivbsSbCAAgWISQIAUM66MCgIQgEDmBBAgmYfAHwcQIP7ECk8h0AECCJAOQKQJCEAAAhDYmAAChFlhTQABYo0KQwgUgQACpAhRZAwQgAAEckgAAZLDoOASBCAAgRwQQIDkIAi4AAEIQKCIBBAgRYwqY4IABCDQPgEESPsM22ohOv30tur7WDletMhHt/EZAhBwJIAAcQQ2bl69+Gxd0Xg2W1TH3zY3jIrMf9QNhBCjRt89hk1s+mS+b/zBbFI0VDAqt9q1qkVfcd23EdVG6VUfSAaIWpkRaJgahhelkv6HeTxTq6OaWs2qzAZGx10ngADpOvLGDhEgGQeA7iEAgdQIIECSoUWAJONGrfQJIEDSZxxKDwiQjCONAMk4AHQPAQikRgABkgwtAiQZN2qlTwABkj7jUHpAgGQcaVOAFDk1KZRxZjyd6B4CuSKAAEkWjuol39QVIyOtqllzDZlQ+h/xqM6JiUpGG6UOJb/ErZaircY83q9RrXT0SckAUQsCEPCeAAIk4xCGsjAPZZwZTye6h0CuCCBAkoUDAZKMG7UgAAF/CCBAMo5VKAvzUMaZ8XSiewjkigACJFk4ECDJuFELAhDwhwACJONYhbIwD2WcGU8nuodArgggQJKFo3qpsQtWQ6rTeBqTkUWlNpcSQqxaNaw6vLXvKVV+9uT5qlwuGzlQDdldLXbBatZ/smGN19L9l445sa2WqAwBCPhLAAGScexCWZiHMs6MpxPdQyBXBBAgycKBAEnGjVoQgIA/BBAgGccqlIV5KOPMeDrRPQRyRQABkiwcCJBk3KgFAQj4QwABknGsQlmYhzLOjKcT3UMgVwQQIMnCUb30G7piw6ZVG+9gNWQcFHj2E5epev83dKcqnzvtzaq865x5uu2KzsGKGg4cNHKzjIMOY6HTp6KG1KyK4e9Eu2ORgpVsVlALAsUigADJOJ6hLMxDGWfG04nuIZArAgiQZOFAgCTjRi0IQMAfAgiQjGMVysI8lHFmPJ3oHgK5IoAASRYOBEgybtSCAAT8IYAAyThWoSzMQxlnxtOJ7iGQKwIIkGThqF5ipGA1pEZt3N4Tw6PqzdOevFCVHxCrVPnLk49X5Z3nLdCNjKxT5VKk07tKUZ96v7JeH2i4bGiZen9Kv7aZ3DdZvR+ZBydOkD5WOoaDCJPNEGpBwH8CCJCMYxjKwjyUcWY8negeArkigABJFg4ESDJu1IIABPwhgADJOFahLMxDGWfG04nuIZArAgiQZOFAgCTjRi0IQMAfAgiQjGMVysI8lHFmPJ3oHgK5IoAASRaO1gJkPKcp0mlRD61bqzr57LJfqPLBA1ur8stnH67KVaFTtoaFsQtWVedL9ZYGlP3Vj92ryuev+Zkq90yeospf7Dtalbecu4UqR1XzpMO675pJ6WhSsJLNEGpBwH8CCJCMYxjKwjyUcWY8negeArkigABJFg4ESDJu1IIABPwhgADJOFahLMxDGWfG04nuIZArAgiQZOFAgCTjRi0IQMAfAgiQjGMVysI8lHFmPJ3oHgK5IoAASRaO6iVnGxXNg/3GUpricln9/Y5lT6rye//7I1U+eUuddvWi+Qeo9yuDeuer0Ypuu2KkS1XLOsXr5Jt/oOreua3ua4/qdur9r016kSpP2Wy6KjemYI2/bQyHFKxk84NaECgCAQRIxlEMZWEeyjgznk50D4FcEUCAJAsHAiQZN2pBAAL+EECAZByrUBbmoYwz4+lE9xDIFQEESLJwIECScaMWBCDgDwEESMaxCmVhHso4M55OdA+BXBFAgCQLR/ViIwXLzMAaP5Qw7ulRDd++/AlV/sQ9F6jymTu8WpV3mT1fleOR9aoc9ejDBEW/3tVq3Tptc+Kt5yv7Oyc/rspfnPliVX72/IW6zaE1etDGbl2iYRxjJhxEmGx+UAsCRSCAAMk4iqEszEMZZ8bTie4hkCsCCJBk4UCAJONGLQhAwB8CCJCMYxXKwjyUcWY8negeArkigABJFg4ESDJu1IIABPwhgADJOFahLMxDGWfG04nuIZArAgiQZOGY6BmQakmnYF259N+qk988+mdVXrSVTsHabN5U7Uhc0eWhXlVeEw2q8mf+pQ8c/Mekper99av1IYbnbqkPH9x5zo7KJhoxUrAMP4XKwdIHHpKClWx+UAsCRSCAAMk4iqEszEMZZ8bTie4hkCsCCJBk4UCAJONGLQhAwB8CCJCMYxXKwjyUcWY8negeArkigABJFg4ESDJu1IIABPwhgADJOFahLMxDGWfG04nuIZArAgiQZOFoFCBGG6WxraRWrdGpUO95VO9SFa/Utqdt9kr1j1nTB1T5P8YuWA9XH1Pvf+Oxq1R59aQhVd565hxVHlyvDzHcp38r9f6nNz9Klft69XZXkdAHJjYjQQpWsvlBLQgUgQACJOMohrIwD2WcGU8nuodArgggQJKFAwGSjBu1IAABfwggQDKOVSgL81DGmfF0onsI5IoAAiRZOBAgybhRCwIQ8IcAAiTjWIWyMA9lnBlPJ7qHQK4IIECShaO1ABlLafrLQ4+qht+9+lxV3rx3uipvVtE7X62btEy9v3poWJVHR/RBhEPVfvX+F+a8TJX3mjFPlT/50K9V+ZHSKlX+7oLXqvKsSbrNqKJ3vBLRxicRkoKVbH5QCwJFIIAAyTiKoSzMQxlnxtOJ7iGQKwIIkGThQIAk40YtCEDAHwIIkIxjFcrCPJRxZjyd6B4CuSKAAEkWDgRIMm7UggAE/CGAAMk4VqEszEMZZ8bTqdb9IadvnOqQB7988eGaRUbaiC9O59RPBEiywDQIEHM6lku1Bu9d9oRq+B1P/FSVV/TpXaomx3rnq+mRTsfaa+02yn7hlO1V+TmzNlPlOUbalejXBxdecv8dyub8FVeq8gXbv1OVZ02epMpRxTj0UKVgcRBhslnRXq0bb/qraiCuGLuTmfOr4afD+Efk9p3Yqplq1RiD0WYsmqfqNVugRmMfgdrLdCsyU/yqrr+BxpxUB2YKsf/Che1Bp/YmCSBAMp4goSzMQxlnxtOp1j0CpL0oIEDa42fWRoAkY4kAScaNWq0JIEA2NTsQIFl8dhAgWVA3LwKcfrq+KrFoUcbepNc9AiQ9thu2jABpjzUCpD1+CJD2+SFA2mdIC40EECAIkLx9JhAgGUcklIV5KOPMeDptdAeExbRdREzRBjM7ZjZW3AGxobSxTctnQMbTTOJYH0T4yFK9G9UdQ4+rxraZtrkqb2ukRfX3TlbvNySqmFk2sf5H3KP9+/Uj96l/nPbkj1X5u/PfpMoL5+q0rmhUp2A12QRLsAtWsvmRpNbiGxeraqVmwRBCGGFv6MI1oamVf2aqVTXW+ViVqp5vFXPujZcjofOuymXtTdnIx2o1plbpYKaPrca9kBSsJFPNug4CxBpVOoahLMxDGWc6s8StVRbTbrykNczcmdnUQIDYUEKAJKNELRcCCBD50EhzYggQl5nUOVsESOdYJmoplIV5KONMNAk6XInFtDtQmLkzs6mBALGhhABJRolaLgQQIAgQl/nSDVsESDcob6KPUBbmoYwz4+lU657FtHsUYObOzKYGAsSGUhMBcvHZ+k3zqm09U8V8r2zuIGTkS5k7DhnpLg29tbr027ArkW7zD489oqq/57FvqfKJcw9R5RMWHKrKJSO1JmqykxIpWMnmR5JaS5bcaFTTE6hxehlzyW3jqw1c0q2aaVdmuVzW82pajz68cqBPl/unjKULjg4PqfbXDOqd3tYMj6j3KxU94SNjXtukj7UaKilYSWaafR0EiD2rVCxDWZiHMs5UJoljoyymHYEh2tyBWdZAgFiC2sCsigBJBo5aLQkgQFpPDgRINh8cBEg23FWvoSzMQxlnxtOp1j0CxD0KMHNnZlMDAWJDaWMbBEgybtRqTQABggDJ2+cDAZJxREJZmIcyzoynEwIkYQAQIAnBTVANAZKMa/WSrxsVJ0giaUijapavJXPfW7zf6qlcMwWrrA+tW71mrfLrxHt+oMoHTt1VlU/Y4XB9gc1Mi1EuaF9Kx5yYDBC1nAksWbLEuY6uYJGP1cokMlKjjPk2OqTTp/59y52qq3ufXKrKA+VptfLAFL2bWqlXe7XTdnrHtbmzZ6s/lPt1elfV6LNhwduQbab/EbV1EGGrz1kb6AtcFQGScXBDWZiHMs6MpxMCJGEAECAJwSFAUgGHAEkFa9CNIkDkNsPNn3Exn01BgHTvY4IA6R7rpj2FsjAPZZwZTycESMIAIEASgkOApAIOAZIK1qAbRYAgQPL2AUCAZByRUBbmoYwz4+mEAEkYAARIQnAIkFTAtT6IsMPdWWTWtErfGh4xDhk08mJ6jBSWqNSwFde486RgdTiKVs3ZCZAWE6LlPGmeclQ1U/iMOVAd1AcKXvq/lyu/z/7GmapcGR5U5f569p+RQjhY1X1On7mZsn3B4S9Q5ZcceYQq77KTTg8sl3T/jUMyDt405q/rLliR4WfSxbVVMAtilJRRFCet2QRcu3nCPscilIV5KOPMw1xkMe0eBZi5M7Op0e53u/mDbtNfUWwQIEWJZH7GgQARAgGSn/koPUkqIxAgHYpjKAvzUMbZoWnRVjMspt3xwcydmU0NBIgNpY1tqpcaD6HHEzyEbnPQQTI3Nq5l+lIyOjYfch8d1vVK+gF2fQy1vtpcOvqkTnlGOxMQcBYgFjdDzOclKrG+I9bbox8CX/nECuXZhRdeoso//vkPVbk8sl6VD9tK191hztgT56PGXFv6lL6rdsUDeq6tN84VOfFd71DtveZ1b1Hl4RF9nkjDvgwt2O2/cCHzKkUCCJAU4do0HcrCPJRx2sQ8bRsW0+6EYebOzKYGAsSGEgIkGSVquRBAgAiBAHGZMenbIkDSZ7zJHkJZmIcyzoynU617FtPuUYCZOzObGggQG0oIkGSUqOVCAAGCAHGZL92wRYB0g/Im+ghlYR7KODOeTgiQhAFAgCQEN0E1BEgyrg27YBkPzgrz4d560zYPkru6YXWcgdlxqzywJu8b1TgHxDUwye2X3GScA9JyzjTfpraxVx3TqpF2VTLS7davW62qnPaJz6jy3/5ynSrvvrlOz1s4t0+9P0U9eS5EaTyra2h0VP198lTtzfCQPhTknFtWqT8c9tLDVPn0j56hynHVGJ9F6iIpWMnnm01NBIgNpRRtQlmYhzLOFKeKddMspq1RKUOYuTOzqYEAsaG0sQ0CJBk3arUmgAARAgGSr08IAiTjeISyMA9lnBlPp1r3LKbdowAzd2Y2NRAgNpQQIMkoUcuFAAIEAeIyX7phiwDpBuVN9GEuzDN2pWvdx4sWda2vEDtiMe0edZi5M7OpgQCxoYQNBNIn0PoZEKvtrpSD5mHiPcZWUmvXrVU2Hz/jU6p83ZVXq/JBC/QOVwfsOEW9P6o3wRLDVb3L1ej4sR369A4hJpf133tj3d6P/6n73+Z/DlJtf/5z+oyRvrLOuzJPP9c7tMm9YXUsFh7ALlhpzkwESJp0LdpGgFhAwsSJAItpJ1w1Y5i5M7OpgQCxoYQNBNIngAARAgGS/jxz6QEB4kIrBVsESApQA2+SxbT7BICZOzObGggQG0rYQCB9AggQBEj6s8ytBwSIGy+sIZB7Aiym3UMEM3dmNjUQIDaUsIFA+gScnwExXWrY8Ez/Y91qnTt11rfOVjWu+dUvVfmgbfpVee+tB1R5cI0+uHCaccDlij7dfrV3LGWqd0inTvVG+vDBcqx30vrpP/TOW5O231X185NzzlPlvgHtS8XMJWuBfyEHEaY6MREgqeKlcQh0nwCLaXfmMHNnZlMDAWJDCRsIpE8AASIEAiT9eebSAwLEhRa2EPCAAItp9yDBzJ2ZTQ0EiA0lbCCQPgEECAIk/Vnm1gMCxI0X1hDIPQEW0+4hgpk7M5saCBAbSthAIH0CrQVI812wGrOumu8edcF3f6Yc//p5X1blI7fRBwQ+fZdp6v21q/WBgqPGbleTjcM2K8ZBhNXx1Kz+EZ2u9fCQ9uzuxwZV24sfG1LlZx9m7IJ15ucNuHo/LXPx2+pITQ4iTHdeIkDS5UvrEOg6ARbT7shh5s7MpgYCxIYSNhBInwACRDJGgKQ/0+x7QIDYs8ISAl4QYDHtHiaYuTOzqYEAsaGEDQTSJ4AAQYCkP8vcekCAuPHCGgK5J8Bi2j1EMHNnZlMDAWJDCRsIpE/A7hkQ7Ucc6wP/SkaK1JNPLlVGZ31Rp13ddN1V6v0PH7q5Kq9eq9OuyqM6lWqkpNO6Bso6CWpYm4i4MvaPK+5ep9r710ptO6lf72p12HNeqGzefMrbVXn+3LmqXK3oMTWcOGjgN9Ox2AUr3XmJAEmXL61DoOsEWEy7I4eZOzObGggQG0rYQCB9AggQIRAg6c8zlx4QIC60sIWABwRYTLsHCWbuzGxqIEBsKGEDgfQJIEAQIOnPMrceECBuvLCGQO4JsJh2DxHM3JnZ1ECA2FDCBgLpE2hHgAz06wME77jrP8rZ959wiipv1f+4Kh++4wxVXjc6osq9ZZ12FUX6EMEBIx3rKX3OoLjg78tqdVdV9cPjz3/+Iaq9Y17+SlXeebfdVXnKFL3zVrVqJFVZHD5oRoJdsNKdlwiQdPnSOgS6ToDFtDtymLkzs6mBALGhhA0E0ieAABFCIEDSn2gOPSBAHGBhCgEfCLCYdo8SzNyZ2dRAgNhQwgYC6RNAgCBA0p9lbj0gQNx4YQ2B3BNgMe0eIpi5M7OpgQCxoYQNBNIn4CpAqsZBgZP6+5SDd/33XlV+5xtPUOUtJ69Q5WN2maPKq40UrL5enQ5VNXa7Khvncyx5WB8uePn962vtnPXRT6n2XnLcsao8WtH5WsNDup7pux3Z5ocxkoJlRy+pFQIkKTnqQSCnBFhMuwcGZu7MbGogQGwoYQOB9AkgQDbFGAGS/gzcuAcESBbU6RMCKRJgMe0OF2buzGxqIEBsKGEDgfQJIEAQIOnPMrceECBuvLCGQO4JsJh2DxHM3JnZ1ECA2FDCBgLpE1iyZEmLTppf/Y+N3aN6e3tU3VVrdKrVKacsUu/f9c8/q/JpBy9Q5dVDQ6o80KP7Goz0zlarVumdsi6+Y62yXzswuVa++Jxz1Hsz526lymXjIMJWBwu2Q5YUrHboTVwXATIxIywg4BUBFtPu4YKZOzObGggQG0rYQCB9AggQd8YIEHdmLjUQIC60sIWABwRYTLsHCWbuzGxqIEBsKGEDgfQJIEDcGSNA3Jm51ECAuNDCFgIeEGAx7R4kmLkzs6mBALGhhA0E0ifQ+AxI87SrBi8MEyNbSvT06MMEz/72BarKD879siqf+cJ5qrxsmd6patIkncq1ZqRX2TywXO9gZfZ16R1j6V6bzdlM2X79G19U5Z133FuVq5Vqc4jaXc4BSX+aOfWAAHHChTEE8k+AxbR7jGDmzsymBgLEhhI2EEifAAKEc0DSn2VuPSBA3HhhDYHcE2Ax7R4imLkzs6mBALGhhA0E0ieAAEGApD/L3HpAgLjxwhoCuSfAYto9RDBzZ2ZTAwFiQwkbCKRPoB0BIow0pp4+/Y/L/u9q5fgnP/FhVT7j8LmqvH5Inzg4UtJ1b75Pp11tMV8fdLjf5pNU3etvX1orX3qv3knriIMPU3//4ld0Ota6kVEDos4fi0znhUXqmdEKz4CkOy8RIOnypXUIdJ0Ai2l35DBzZ2ZTAwFiQwkbCKRPAAEiGSNA0p9p9j0gQOxZYQkBLwiwmHYPE8zcmdnUQIDYUMIGAukTQIAgQNKfZW49IEDceGENgdwTYDHtHiKYuTOzqYEAsaGEDQTSJ7Dkpht1JxY3AszFYblcVnWHhtap8lc/d54q/+Y3P1Lljxykd61asVIfMnjX43pHrIdinY71zO2nqbrTBnVq1uYL+mvvf+uqZervfdvvo8o/Pu87qhyX9MGGsXGnw9wEy5UyKViuxNzsESBuvLCGQO4JsJh2DxHM3JnZ1ECA2FDCBgLpE0CAuDNGgLgzc6mBAHGhhS0EPCDAYto9SDBzZ2ZTAwFiQwkbCKRPAAHizhgB4s7MpQYCxIUWthDwgACLafcgwcydmU0NBIgNJWwgkD6BJUuMFCyL7lqlYA2vW6NqL/rUV1T5rr9epspvfe5sVV47qA8IHFyjc7/+8Kje2Wqecbjh0+bpVKqZk8cOK/z9LctVe8u22k+Vv/9tvQtW1KN3z4qMvCtSsCyCnZEJAiQj8HQLgbQIsJh2Jwszd2Y2NRAgNpSwgUD6BBAg7oy5A+LOzKUGAsSFFrYQ8IAAi2n3IMHMnZlNDQSIDSVsIJA+AQSIO2MEiDszlxoIEBda2ELAAwIspt2DBDN3ZjY1ECA2lLCBQPoEWgkQM0XJ3BwrNv5hbDAlymWdIvW5r5ytHL/yZz9W5Xc/T++CNTKod8FatkynY137lD44cK+ZParu9nN0+7N6x3bBuvYfT6m/PzTvf1T5u985U5V7+/VOWpGZg2WgdU3HQoCkOy8RIOnypXUIdJ0Ai2l35DBzZ2ZTAwFiQwkbCKRPAAEiGs5EtyGOALGhlNwGAZKcHTUhkEsCLKbdwwIzd2Y2NRAgNpSwgUD6BBAgCJD0Z5lbDwgQN15YQyD3BFhMu4cIZu7MbGogQGwoYQOB9Ak0ChCdjBQZh/aZaVetzirsG9CHEv71mhuU46d84GRVft2+k1V5dr9OqVqzXqddVQbGdriSr5nG8EdGtc30gYHaX678uz6I8PH5+yrrH33vq6pc6hlL19rw5Zp2ZdbnDki68xIBki5fWodA1wmwmHZHDjN3ZjY1ECA2lLCBQPoEECDujBEg7sxcaiBAXGhhCwEPCLCYdg8SzNyZ2dRAgNhQwgYC6RNAgLgzRoC4M3OpgQBxoYUtBDwgwGLaPUgwc2dmUwMBYkMJGwikT8AUIGaqVYsNo4SRmdXgXKmsk5rWDQ6qvx39qlep8vOmrlDlfRfMUuVlI8O6rbJO5Zpc1btjrTfKk6pj6Vs/uEXvgvW0F7xMtXHGxz6iyqORTvUyF7akYKU/t5L2gABJSo56EMgpARbT7oGBmTszmxoIEBtK2EAgfQIIEHfG3AFxZ+ZSAwHiQgtbCHhAgMW0e5Bg5s7MpgYCxIYSNhBIn8CNS5Y07cR8CL3VXQ+zYhzpx9N7+/rUnz546vtU+Z7r/6jKb1q4uSpXRiuqvG5YP2xe6tN3Q6qj+k7GpbeN3fl4dJ3u801ve5tu+6263OCjeYunBdpWd0bMh+8XLlyYfmAC7gEBEnDwGXoxCbCYdo8rzNyZ2dRAgNhQwgYC6RNAgDQyRoCkP+cm6gEBMhEh/g4BzwiwmHYPGMzcmdnUQIDYUMIGAukTQIAgQNKfZW49IEDceGENgdwTYDHtHiKYuTOzqYEAsaGEDQTSJ3DjkhuNTsxzQPTbNg9vV4V+YHzyJH32xnVXXqkaev+HTlXl43fUNvM3n6reX7FmnSr/e7lOx1ry0JB6f6R/7KyQE9/xfvXeYS89XJVnTZ2uyuYdjVZnmJiUbR5O5xmQdOclAiRdvrQOga4TYDHtjhxm7sxsaiBAbChhA4H0CSBAGhkjQNKfcxP1gACZiBB/h4BnBFhMuwcMZu7MbGogQGwoYQOB9AkgQBAg6c8ytx4QIG68sIZA7gmwmHYPEczcmdnUQIDYUMIGAukTWNJiFyybniNjV6mquQtWuUdVH1qrU6fef6ZOwfrHDXr3raq581VVJ0qVp+jUrGc9+1mqzde/8fW18o7b7KDdNM4JiYXeMathNy+bQVnYkIJlAakNEwRIG/CoCoE8EmAx7R4VmLkzs6mBALGhhA0E0ieAAHFnjABxZ+ZSAwHiQgtbCHhAgMW0e5Bg5s7MpgYCxIYSNhBInwACxJ0xAsSdmUsNBIgLLWwh4AEBFtPuQYKZOzObGggQG0rYQCB9Aotv1LtgmQ9gR8Y/zPP7Wp6TYfwhEsZuWmblXn2w4MpHH1eDGxlcq8rDI8OqvL40ttuVfC2Yt0CVJ43vstVqoWruduW6C1ZL4kZDC/fnIMI0ZyYCJE26tA2BDAiwmHaHDjN3ZjY1ECA2lLCBQPoEECCWjBEglqDaN0OAtM+QFiCQKwIspt3DATN3ZjY1ECA2lLCBQPoEECCWjBEglqDaN0OAtM+QFiCQKwIspt3DATN3ZjY1ECA2lLCBQPoE/n7j31Qno0beVUPqUmxzOob21dgQS5gZWOVI707V2693yor69KGEA/19qqE16wd1o4O6PFwZS7Iq6+Ya+mlIHxPNjx9sOFzRqNCw+DXTyozy/v9DClaaMxMBkiZd2oZABgSQuBM8AAALwUlEQVRYTLtDh5k7M5saCBAbSthAIH0CCBAhIgRI+hPNoQcEiAMsTCHgAwEW0+5Rgpk7M5saCBAbSthAIH0CCBAESPqzzK0HBIgbL6whkHsC5mI6987m0MFrFjW/lZ9DV3PvEgIk9yHCQQhAAAKZEECAZIKdTiGQHgEESHtsESDt8TNrI0A6x5KWIAABCBSJAAKkSNFkLBAQQiBA2psGCJD2+CFAOsePliAAAQgUlQACpKiRZVwQgAAEMibAHZCMA0D3EMiagLlVVcNOVRZHB5INm3X0Uu0fAZIqXhqHAAQgEC4BBEi4sWfkEKgRQIAwEVoQQIAwNSAAAQhAIBUCCJBUsNIoBPwhgADxJ1Zd9hQB0mXgdAcBCEAgFAIIkFAizTghAAEIuBFAgLjxwhoCEIAABCwJIEAsQWEGAQhAIDACCJDAAs5wIQABCHSLAAKkW6TpBwIQgIBfBBAgfsULbyEAAQh4QwAB4k2ocBQCXhAwHymJ2SXLi5i1chIB4nX4cB4CEIBAfgkgQPIbGzyDgI8EECA+Rq25zwiQ4sSSkUAAAhDIFQEESK7CgTMQ8J4AAsT7EKoBIECKE0tGAgEIQCBXBBAguQoHzkAAAhDIDQEESG5CgSMQgAAEikUAAVKseDIaCEAAAp0igADpFEnagQAEIACBBgIIECYEBCAAAQg0I4AAYV5AAAIQgEAqBBAgqWClUQ8JRC18ZiMnD4OJyx0hgADpCEYagQAEIACBDQkgQJgTEBgjgABhJkCgkQAChBkBAQhAAAKpEECApIKVRj0kgADxMGi4nCoBBEiqeGkcAhCAQLgEECDhxp6RQwACENgUAQQI8wMCEIAABFIhgABJBSuNQgACEPCeAALE+xAyAAhAAAL5JIAAyWdc8AoCEIBA1gQQIFlHgP4hAAEIFJQAAqSggWVYqRAwnxNhd6xUENNojgggQHIUDFyBAAQgUCQCCJAiRZOxpE0AAZI2YdrPEwEESJ6igS8QgAAECkQAAVKgYDKU1AkgQFJHTAc5IoAAyVEwcAUCEIBAkQggQIoUTcYCAQhAoHMEECCdY0lLEIAABCBgEECAMB0gAAEIQKAZAQQI8wICEIAABFIhgABJBSuNQgACEPCeAALE+xAyAAhAAAL5JIAAyWdc8AoCEIBA1gQQIFlHgP4hAAEIFJQAAqSggWVYEIAABNokgABpEyDVIQABCECgOQEECDMDAhCAAASaEUCAMC8gAAEIQCAVAgiQVLDSKAQgAAHvCSBAvA8hA4AABCCQTwIIkHzGBa8gAAEIZE0AAZJ1BOgfAhCAQEEJIEAKGliGBQEIQKBNAgiQNgFSHQIQgAAEmhNAgDAzIAABCECgGQEECPMCAhCAAARSIYAASQUrjUIAAhDwngACxPsQMgAIQAAC+SSAAMlnXPAKAhCAQNYEECBZR4D+IQABCBSUAAKkoIFlWBCAAATaJIAAaRMg1SEAAQhAoDkBBAgzAwIQgAAEmhFAgDAvIAABCEAgFQIIkFSw0igEIAAB7wkgQLwPIQOAAAQgkE8CCJB8xgWvIAABCGRNAAGSdQToHwIQgEBBCSBAChpYhgUBCECgTQIIkDYBUh0CEIAABJoTQIAwMyAAAQhAoBkBBAjzAgIQgAAEUiGAAEkFK41CAAIQ8J4AAsT7EDIACEAAAvkkgADJZ1zwCgIQgEDWBBAgWUeA/iEAAQgUlAACpKCBZVgQgAAE2iSAAGkTINUhAAEIQKA5AQQIMwMCEIAABJoRQIAwLyAAAQhAIBUCCJBUsNIoBCAAAe8JIEC8DyEDgAAEIJBPAgiQfMYFryAAAQhkTQABknUE6B8CEIBAQQkgQAoaWIYFAQhAoE0CCJA2AVIdAhCAAASaE0CAMDMgAAEIQKAZAQQI8wICEIAABFIhgABJBSuNQgACEPCeAALE+xAyAAhAAAL5JIAAyWdc8AoCEIBA1gQQIFlHgP4hAAEIFJQAAqSggWVYEIAABNokgABpEyDVIQABCECgOYF2BUipFKmG4xjKEIAABCDgNYHI+E6vVhMNJYqTSpcm3bX7I5VoBFSCAAQgAIFUCbT73Y4ASTU8NA4BCECguwTyJkBWrlwpZs6cKR588EExffr07sKgNwhAAAIQSIWAFCBbb721WLFihZgxY4ZzHwgQZ2RUgAAEIJBfAnkTIA899FDtR4oXBCAAAQgUj4C8uLTVVlsVb2CMCAIQgAAEukqgoylY1WpVPPLII2LatGkiMtRRV0dEZxCAAAQg0FECMlN39erVYsGCBaJUKnW0bRqDAAQgAIHwCHRUgISHjxFDAAIQgAAEIAABCEAAAi4EECAutLCFAAQgAAEIQAACEIAABNoigABpCx+VIQABCEAAAhCAAAQgAAEXAggQF1rYQgACEIAABCAAAQhAAAJtEUCAtIWPyhCAAAQgAAEIQAACEICACwEEiAstbCEAAQhAAAIQgAAEIACBtgggQNrCR2UIQAACEIAABCAAAQhAwIUAAsSFFrYQgAAEIAABCEAAAhCAQFsEECBt4aMyBCAAAQhAAAIQgAAEIOBCAAHiQgtbCEAAAhCAAAQgAAEIQKAtAgiQtvBRGQIQgAAEIAABCEAAAhBwIYAAcaGFLQQgAAEIQAACEIAABCDQFgEESFv4qAwBCEAAAhCAAAQgAAEIuBBAgLjQwhYCEIAABCAAAQhAAAIQaIsAAqQtfFSGAAQgAAEIQAACEIAABFwIIEBcaGELAQhAAAIQgAAEIAABCLRFAAHSFj4qQwACEIAABCAAAQhAAAIuBBAgLrSwhQAEIAABCEAAAhCAAATaIoAAaQsflSEAAQhAAAIQgAAEIAABFwIIEBda2EIAAhCAAAQgAAEIQAACbRFAgLSFj8oQgAAEIAABCEAAAhCAgAsBBIgLLWwhAAEIQAACEIAABCAAgbYIIEDawkdlCEAAAhCAAAQgAAEIQMCFAALEhRa2EIAABCAAAQhAAAIQgEBbBBAgbeGjMgQgAAEIQAACEIAABCDgQgAB4kILWwhAAAIQgAAEIAABCECgLQIIkLbwURkCEIAABCAAAQhAAAIQcCGAAHGhhS0EIAABCEAAAhCAAAQg0BYBBEhb+KgMAQhAAAIQgAAEIAABCLgQQIC40MIWAhCAAAQgAAEIQAACEGiLAAKkLXxUhgAEIAABCEAAAhCAAARcCCBAXGhhCwEIQAACEIAABCAAAQi0RQAB0hY+KkMAAhCAAAQgAAEIQAACLgQQIC60sIUABCAAAQhAAAIQgAAE2iKAAGkLH5UhAAEIQAACEIAABCAAARcCCBAXWthCAAIQgAAEIAABCEAAAm0RQIC0hY/KEIAABCAAAQhAAAIQgIALAQSICy1sIQABCEAAAhCAAAQgAIG2CCBA2sJHZQhAAAIQgAAEIAABCEDAhQACxIUWthCAAAQgAAEIQAACEIBAWwQQIG3hozIEIAABCEAAAhCAAAQg4EIAAeJCC1sIQAACEIAABCAAAQhAoC0CCJC28FEZAhCAAAQgAAEIQAACEHAhgABxoYUtBCAAAQhAAAIQgAAEINAWAQRIW/ioDAEIQAACEIAABCAAAQi4EECAuNDCFgIQgAAEIAABCEAAAhBoiwACpC18VIYABCAAAQhAAAIQgAAEXAggQFxoYQsBCEAAAhCAAAQgAAEItEUAAdIWPipDAAIQgAAEIAABCEAAAi4EECAutLCFAAQgAAEIQAACEIAABNoigABpCx+VIQABCEAAAhCAAAQgAAEXAggQF1rYQgACEIAABCAAAQhAAAJtEUCAtIWPyhCAAAQgAAEIQAACEICACwEEiAstbCEAAQhAAAIQgAAEIACBtgggQNrCR2UIQAACEIAABCAAAQhAwIUAAsSFFrYQgAAEIAABCEAAAhCAQFsEECBt4aMyBCAAAQhAAAIQgAAEIOBCAAHiQgtbCEAAAhCAAAQgAAEIQKAtAgiQtvBRGQIQgAAEIAABCEAAAhBwIYAAcaGFLQQgAAEIQAACEIAABCDQFoH/Dw4LA5kJQ4p1AAAAAElFTkSuQmCC\" width=\"800\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "853176327e934b179ede66948d7fa5cc",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "Button(description='Resample', style=ButtonStyle())"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "a2e4712184034d0686fc002b80548a50",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "HBox(children=(interactive(children=(Select(description='Class', options=('0', '1', '2', '3', '4', '5', '6', '…"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "89193b6ccc724124b2383113ff052cce",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "HBox(children=(interactive(children=(Select(description='Background Color', options=('white', 'cyan', 'salmon'…"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "31a760ed6b83407fb5e202975c5df81b",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "Button(description='Save as image', style=ButtonStyle())"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "e258054758da466ab4657d07f8803c23",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "Button(description='Reset', style=ButtonStyle())"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%matplotlib notebook\n",
    "\n",
    "# adapted from https://github.com/dariopavllo/style-semantics\n",
    "\n",
    "import matplotlib\n",
    "from matplotlib.widgets import TextBox, Button\n",
    "from ipywidgets import interact, interactive, fixed, interact_manual, HBox\n",
    "import ipywidgets as widgets\n",
    "import imageio\n",
    "\n",
    "from demo.dmnist import sample\n",
    "from utils.util_v2 import *\n",
    "\n",
    "label_nc = 10\n",
    "\n",
    "label_names = [i for i in range(label_nc)].copy()\n",
    "color_names = ['blue', 'red', 'green', 'violet', 'brown']\n",
    "bg_color_names = ['white', 'cyan', 'salmon', 'yellow', 'silver']\n",
    "style_names = ['thin', 'thick']\n",
    "\n",
    "cmap_arr = labelcolormap(label_nc + 1)\n",
    "cmap = matplotlib.colors.ListedColormap(cmap_arr/255)\n",
    "semantic_buffer = torch.FloatTensor(1, label_nc, 256, 256).zero_()\n",
    "if torch.cuda.is_available():\n",
    "    semantic_buffer = semantic_buffer.cuda()\n",
    "sem_buffer = torch.LongTensor(128, 128).zero_()\n",
    "\n",
    "im_buffer = np.zeros((128, 128, 3), dtype=np.float32)\n",
    "\n",
    "list_bbox = []\n",
    "list_label = []\n",
    "list_dcolor = []\n",
    "list_bgcolor = []\n",
    "list_dstyle = []\n",
    "\n",
    "im_size = 8\n",
    "fig = plt.figure(2, figsize=(im_size, im_size/2+1))\n",
    "plt.subplots_adjust(wspace=0, hspace=0.03)\n",
    "\n",
    "plt.subplot(1, 2, 1)\n",
    "plt.title('Mask')\n",
    "mask_im = plt.imshow(sem_buffer, cmap=cmap, vmin=0, vmax=label_nc + 1)\n",
    "plt.xticks([])\n",
    "plt.yticks([])\n",
    "ax = plt.gca()\n",
    "\n",
    "plt.subplot(1, 2, 2)\n",
    "plt.title('Output')\n",
    "output_im = plt.imshow(im_buffer)\n",
    "plt.axis('off')\n",
    "\n",
    "fig.tight_layout()\n",
    "plt.show()\n",
    "\n",
    "current_class = 1\n",
    "\n",
    "num_o = 8\n",
    "thres = 2.0\n",
    "z_obj = torch.from_numpy(truncted_random(num_o=num_o, thres=thres)).float()\n",
    "z_im = torch.from_numpy(truncted_random(num_o=1, thres=thres)).view(1, -1).float()\n",
    "\n",
    "def redraw():\n",
    "    with torch.no_grad():\n",
    "        global im_buffer\n",
    "        global list_bbox, list_label, list_dcolor, list_bgcolor, list_dstyle\n",
    "        global z_obj, z_im\n",
    "\n",
    "        im_buffer = sample(bbox=list_bbox,\n",
    "                           label=list_label,\n",
    "                           dcolor=list_dcolor,\n",
    "                           bgcolor=list_bgcolor,\n",
    "                           dstyle=list_dstyle,\n",
    "                           model_path='./pretrained/MNIST_Dialog_Generator.pth',\n",
    "                           z_im=z_im, z_obj=z_obj, resample=False)\n",
    "        im_buffer = (im_buffer * 0.5 + 0.5)\n",
    "        output_im.set_data(im_buffer)\n",
    "\n",
    "origin = None\n",
    "def on_click(event):\n",
    "    global origin\n",
    "    x = int(event.xdata)\n",
    "    y = int(event.ydata)\n",
    "    origin = (x, y) \n",
    "\n",
    "prev_coords = None\n",
    "def on_move(event):\n",
    "    global prev_coords\n",
    "    if event.inaxes and event.inaxes == ax and event.button in [1, 3]:\n",
    "        x = int(event.xdata)\n",
    "        y = int(event.ydata)\n",
    "        \n",
    "        temp_sem_buffer = sem_buffer.clone()\n",
    "        \n",
    "        temp_sem_buffer[origin[1]][origin[0]:x] = current_class\n",
    "        temp_sem_buffer[y][origin[0]:x] = current_class\n",
    "        for i in range(y - origin[1] + 1):\n",
    "            temp_sem_buffer[origin[1] + i][origin[0]] = current_class           \n",
    "            temp_sem_buffer[origin[1] + i][x] = current_class\n",
    "        \n",
    "        mask_im.set_data(temp_sem_buffer)\n",
    "        prev_coords = (x, y)\n",
    "        \n",
    "def on_release(event):\n",
    "    x, y = prev_coords\n",
    "    sem_buffer[origin[1]][origin[0]:x] = current_class\n",
    "    sem_buffer[y][origin[0]:x] = current_class\n",
    "    for i in range(y - origin[1] + 1):\n",
    "        sem_buffer[origin[1] + i][origin[0]] = current_class           \n",
    "        sem_buffer[origin[1] + i][x] = current_class\n",
    "    \n",
    "    mask_im.set_data(sem_buffer)\n",
    "    \n",
    "    global list_bbox, list_label, list_dcolor, list_bgcolor, list_dstyle\n",
    "    list_bbox.append([origin[0], origin[1], prev_coords[0]-origin[0], prev_coords[1]-origin[1]])\n",
    "    list_label.append(current_class - 1)\n",
    "    list_dcolor.append(current_digit_color)\n",
    "    list_bgcolor.append(current_bg_color)\n",
    "    list_dstyle.append(current_style)\n",
    "    \n",
    "    redraw()\n",
    "    \n",
    "def update_class(cl):\n",
    "    global current_class\n",
    "    current_class = int(cl)\n",
    "    ax.set_title(f'Selected class: {current_class}')\n",
    "    current_class += 1\n",
    "    \n",
    "def update_digit_color(dc):\n",
    "    global current_digit_color\n",
    "    current_digit_color = color_names.index(dc)\n",
    "    \n",
    "def update_bg_color(bc):\n",
    "    global current_bg_color\n",
    "    current_bg_color = bg_color_names.index(bc)\n",
    "    \n",
    "def update_style(s):\n",
    "    global current_style\n",
    "    current_style = style_names.index(s)\n",
    "    \n",
    "def save_to_file(_):\n",
    "    imageio.imwrite('saved_drawing_label.png', cmap_arr[sem_buffer])\n",
    "    imageio.imwrite('saved_drawing_image.png', (im_buffer * 255).astype(np.uint8))\n",
    "    ax.set_title(f'(!) Saved as saved_drawing_*.png')\n",
    "    \n",
    "def resample(_):\n",
    "    global z_obj, z_im\n",
    "    z_obj = torch.from_numpy(truncted_random(num_o=num_o, thres=thres)).float()\n",
    "    z_im = torch.from_numpy(truncted_random(num_o=1, thres=thres)).view(1, -1).float()\n",
    "    redraw()\n",
    "    \n",
    "def reset(_):\n",
    "    sem_buffer.zero_()\n",
    "    mask_im.set_data(sem_buffer)\n",
    "    \n",
    "    global list_bbox, list_label, list_dcolor, list_bgcolor, list_dstyle\n",
    "    list_bbox = []\n",
    "    list_label = []\n",
    "    list_dcolor = []\n",
    "    list_bgcolor = []\n",
    "    list_dstyle = []\n",
    "    \n",
    "    global im_buffer\n",
    "    im_buffer = np.zeros((128, 128, 3), dtype=np.float32)\n",
    "    output_im.set_data(im_buffer)\n",
    "\n",
    "fig.canvas.mpl_connect('button_press_event', on_click)\n",
    "fig.canvas.mpl_connect('motion_notify_event', on_move)\n",
    "fig.canvas.mpl_connect('button_release_event', on_release)\n",
    "\n",
    "class_box = widgets.Select(options=[str(label_names[x]) for x, y in enumerate(label_names)],\n",
    "                              description='Class')\n",
    "\n",
    "digit_color_box = widgets.Select(options=[color_names[x] for x, y in enumerate(color_names)],\n",
    "                              description='Digit Color')\n",
    "\n",
    "bg_color_box = widgets.Select(options=[bg_color_names[x] for x, y in enumerate(bg_color_names)],\n",
    "                              description='Background Color')\n",
    "\n",
    "style_box = widgets.Select(options=[style_names[x] for x, y in enumerate(style_names)],\n",
    "                              description='Style')\n",
    "\n",
    "\n",
    "resample_button = widgets.Button(description='Resample')\n",
    "resample_button.on_click(resample)\n",
    "save_button = widgets.Button(description='Save as image')\n",
    "save_button.on_click(save_to_file)\n",
    "reset_button = widgets.Button(description='Reset')\n",
    "reset_button.on_click(reset)\n",
    "\n",
    "i_class = interactive(update_class, cl=class_box)\n",
    "i_digit_color = interactive(update_digit_color, dc=digit_color_box)\n",
    "i_background_color = interactive(update_bg_color, bc=bg_color_box)\n",
    "i_style = interactive(update_style, s=style_box)\n",
    "\n",
    "display(resample_button)\n",
    "display(HBox([i_class, i_digit_color]))\n",
    "display(HBox([i_background_color, i_style]))\n",
    "\n",
    "display(save_button, reset_button)"
   ]
  },
  {
   "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.6"
  },
  "pycharm": {
   "stem_cell": {
    "cell_type": "raw",
    "metadata": {
     "collapsed": false
    },
    "source": []
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}