Autoencoders.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "%matplotlib notebook\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import matplotlib\n",
    "\n",
    "matplotlib.rc(\"image\", cmap=\"gray\", interpolation=\"nearest\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "First load the patterns, and display them. There are a few numbers, plus \"HI\", plus a smiley face:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\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",
       "            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",
       "        this.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 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);\n",
       "        canvas.attr('height', height);\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'];\n",
       "    var y0 = fig.canvas.height - msg['y0'];\n",
       "    var x1 = msg['x1'];\n",
       "    var y1 = fig.canvas.height - msg['y1'];\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, fig.canvas.height);\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;\n",
       "    var y = canvas_pos.y;\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 overriden (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",
       "    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 + '\">');\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 dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\">';\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",
       "        event.shiftKey = false;\n",
       "        // Send a \"J\" for go to next cell\n",
       "        event.which = 74;\n",
       "        event.keyCode = 74;\n",
       "        manager.command_mode();\n",
       "        manager.handle_keydown(event);\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,iVBORw0KGgoAAAANSUhEUgAAAyAAAACgCAYAAAD0KAgtAAAfuklEQVR4Xu2dUXIctw5F2wuI/J/E+1+Yk/yPsgC/akWqZ1ka4YINECB19JWqoAnyXBDknR5ZXw5+IAABCEAAAhCAAAQgAAEITCLwZVIe0kAAAhCAAAQgAAEIQAACEDgwIBQBBCAAAQhAAAIQgAAEIDCNAAZkGmoSQQACEIAABCAAAQhAAAIYEGoAAhCAAAQgAAEIQAACEJhGAAMyDTWJIAABCEAAAhCAAAQgAAEMCDUAAQhAAAIQgAAEIAABCEwjgAGZhppEEIAABCAAAQhAAAIQgAAGhBqAAAQgAAEIQAACEIAABKYRwIBMQ00iCEAAAhCAAAQgAAEIQAADQg1AAAIQgAAEIAABCEAAAtMIYECmoSYRBCAAAQhAAAIQgAAEIIABoQYgAAEIQAACEIAABCAAgWkEMCDTUJMIAhCAAAQgAAEIQAACEMCAUAMQgAAEIAABCEAAAhCAwDQCGJBpqEkEAQhAAAIQgAAEIAABCGBAqAEIQAACEIAABCAAAQhAYBoBDMg01CSCAAQgAAEIQAACEIAABDAg1AAEIAABCEAAAhCAAAQgMI0ABmQaahJBAAIQgAAEIAABCEAAAhgQagACEIAABCAAAQhAAAIQmEYAAzINNYkgAAEIQAACEIAABCAAAQwINQABCEAAAhCAAAQgAAEITCOAAZmGmkQQgAAEIAABCEAAAhCAAAaEGoAABCAAAQhAAAIQgAAEphHAgExDTSIIQAACEIAABCAAAQhAAANCDUAAAhCAAAQgAAEIQAAC0whgQKahJhEEIAABCEAAAhCAAAQggAGhBiAAAQhAAAIQgAAEIACBaQQwINNQkwgCEIAABCAAAQhAAAIQGDUg53O/H8fxLwjLCfx2HMc/x3H8cM4EDZ3AksLRLwnsxGFHNGT/TRTISIV+fbQYmcmIfmce9uAI7fhn0C+e6ewRhzQcNSB/HMfx1+wVku8ugT+P4/jbyQcNncASw9EvEe6kob0asv8mCSOmQT8RVNMwr37nMtiDfcREvz5ajM7EreGoAXk4juP2/fv34+Hh/M/7P1+/fh1dzKd/7na7fcjg8fHx+Pbt2xlzQn50AnvS0PnMh+HWfCNzVY6VUNOp+qm6RK9LyVuR82nDxPclr4bh+69yT2yQG/3WFtGr37la9mAfzdGvjxajM3FreMmAnBcMy4B8+TKaYpTBPs/9+PHxt6pOA/J8kXILn9F8rfnuokxCTafqp+oSvS4lb0XOsw6j8w58CMDlp1dD8O5B9FtbPwwI+vUisP5svD306TuQIz9PzRcDMoJOf8a6wGFAdJaRkQ0ur67D06qjFzbR61LyVuTEgETuhm3G8h6eGJBe0nv1c/XQXkvdcjbot76sbg0xII1Fty5wGJAa8aIvzQOfnrsOT6uOMCAhdeRtvlxgQ7CHDYJ+YShLBvLq5+qhJSv6XEnRb3293RpiQBqLbl0cMSA14mFANO5W/Wa8iVByZuQdMJEYEK2MZkV5D0/0m6WMlserHwZE4zorCv1mkc7L49YQA5InxuWRrcsUBuQy4qEBMCAaNqt+M4yAkjMjLwZEq4nGUd7DEwPSS0yvfhgQ9OtFYP3ZuPcgBqSx6NZlCgNSIx4GRONu1W+GEVByZuTFgGg10TjKe3hiQHqJ6dUPA4J+vQisPxv3HsSANBbdukxhQGrEw4Bo3K36zTACSs6MvBgQrSYaR3kPTwxILzG9+mFA0K8XgfVn496DlwyIwqvwQqBMryQmgYlbeM8/w5sw3xLuUUkTeGypn8Ip2swpOT+7AVEZRe2XJN6R0zvH8u5BDEi0AtfG8+pXZkAq9t81tG+fju7bA/svRT9Vm4T1h0lUuAb3HsSAhMmuD5RQIG7hMSC6Xr9Gop/GTuEU3ciVnEkXYu8eLLvAqow0lbWoaJ21rK6oZfRzrerzBHv1S7nAKrgr9p8yL09Mwn5uoZ+qTcL6Pfg/jC1cg1tDDEiY7PpACQXiFh4DouuFARljpdR5dCNXcmJAPv4Dp2Nqf/xUtM4Jc/T20DIDmbD2HYb06ocBuaB6wn5uoV/h+XFBjdePFq7BrSEGJEx2faCEAnELjwHR9cKAjLFS6jz6IFNyYkAwIO9UtLeHYkDG2kLWU179MCAXlIju23wF64IYvzy60Bl47S+hK8gKYSjTK4lJYJLafBPmW8I9KmkCjy31UzhFH2RKTgwIBgQDEtUN24yT2kMjV6n2qMic0WNF920MSJxCan110JA3IHG6yyMlFEhq802Yr8yqY2ACjy31UzhFN0ElJwYEA4IB6dhZL80ptYdemtngJ9SROaPHiu7bGJA4hRY6A3kDEie7PlJCgaQ234T56rAaRibw2FI/hVP0QabkxIBgQDAgDRvrtSml9tBrU3v9tNqjInNGjxXdtzEgcQqp9dVBQ96AxOkuj5RQIKnNN2G+MquOgQk8ttRP4RTdBJWcGBAMCAakY2e9NKfUHnppZrwBUfC10K/w/FAYSTGFa3BriAGRJI0NSigQt/D8Evq4puinsVM4YUA0lpFRii6R+ZIMX/QUvT2UX0KPVuDaeF79zmwlGlbsv2to3z4d3bd5AxKnkFpfHTTEgMTpLo+UUCCpzTdhvjKrjoEJPLbUT+EU3QSVnEkXYq+GJZefc+0qo8i9F61z5Nyex1pGv4S17zCkVz8MyAXVE/ZzC/3U3piw/gtqvH60cA1uDdMNiEq1EJo6RTOucA1u4V8+/bndbsfDw3kXuv/TebOZohQGWPXw+Ph4fP16Suf+K8yuw9OaxwuiaJ2VvBU5MSDaV7AUbRSNk3hH72xvDw03kCrL6IVHjqfUTGS+n8by6ufqoUlzbjmsUocJOrfQT1l7935WuAa3hhiQwBawkvAYkEDh7wxl1QMG5DiiDzKLeZbZGjCR4RdYtaIjGUWOpc4/Kc57eIbrp7JMWn/IsNH72TEpr34YkMFzK+kC3kI/dQ8W1rm5JQrX4NYQA2LKqQesJDwGRNd1NNKqBwwIBmS0tq48Z9Wlx6RFjnVlTQHPeg9PDMg70AsvZl79MCAYkDcEduhnhWtw70EMSMDJ9TLESsJjQAKFH2zkGBAMSH4Vvs0Q2acix6pg8VNO7+GJAcGAFJdsTnplTycYTe/+SzGQytqT3gCFiVm4BreGGJAw2fVf7uy0efkdkMAC+GUoqxFgQDAgedV3f2SrLnkDIqmCAcGASIWyWpDSHzrdYSL5KmvHgNwljgGJLEbvWIXF6xaeNyBedf3xVj1gQDAg/qq6/oRVlxgQiTEGBAMiFcpqQUp/wICMfnafXw2Kfkkmyn0PHaVY1nwTCj+sIlYSHgMSJvvdgax6wIBgQPKr8G0Gqy4xIJIqZWegNLuioMLz2X35qfo7IEXSyGmV/pCgcwv9lLUnXd5lfazAwjW4NcSAWGo6/v9KwmNAHMIOhlr1gAHBgAyW1qXHrLrEgEh4MSC8AZEKZbUgpT9gQEavzvnVoOiXZKIwIPny3s+wkvAYkPxKseoBA4IBya9C3oCIjL2HJwYEAyKW1lph1rnV5fKa8QZLWXvS+sOKpHAN3h56jNq4suab4Lw/pfAYkDDZ7w5kNQIMCAYkvwoxICJj7+FZdgaK6ykJKzyfvfqdfMI1LIEenNQ6t5Iu4C30U9aetP4wFQvX4NawjQFR6VfArcip8niOcwufYUAUTuoBpYzlaQTR4yn6WDkxIBgQpY46x1g1/jJ3dd8XrtXbQz/V5XUBnb36YUDubDZF64T93EI/Ze2ee0dFPytcg1tDDIhQIYWCCrN7CnELjwF5H21kY7XqBgOCAVE3eNc4q8YxIF2V881rAZ2Hz0Afif2jFa0jz8mrH6JGKqKsHQNyl7h7D2JAhOpdoCjdwmNAMCC/Eog+VJR9U5Ez6QDx7sEtPkFXNE7iLXRuV8in1E8ltIDOXv14A8IbkDcEFqhzc8sWrsG9BzEgppylf2BQmB1vQBRIFZvSyskbEN6AKLXbOcaqcd6AdFZPn9sCOrsvP/wOyPv6K1pHf3B05VscehXbkcrau3+gUrgG9x7EgNg1eRQKKswOA6JAqtDQyokBwYAotds5xqpxDEhn9fS5LaCz+/KDAcGA/EpggTo3N23hGtx7EANiyskbEPXTDqXwI8fyfBKhzM0znlA2pnHFgGBAlDrqHFOxr5J4eA/PLb5Cp7JcQGevfufSP5WGkVqr57iakzcgDlJGaOFede9BDIige6Ggwux4A6JAqtDQyokBwYAotds5xqpx3oB0Vk+f2wI6uy8/GBDegPAGRO8BQqR7D2JABKo7N9/b7XY8PJwfBN3/UT/tUDhFjuV5Y6HMzTOeUDa8Afnxw8Sk1oM50HNAhc7Pqb3Nd4tPXwt5qyWhxn1K/VQ4C+js1e9c+hZ7UNVQjVO0ju7bvAFR1bHjFP2i7zqDZ2CfP0RoY/0vogJuRU6Vx6jw/CtY7xOObKxW3fAGhDcgzn3eLtyq8ZcJR+6rJAjeC+ynurwuoLNXPwzInY2kaJ2wn1vop6w96fIe1tYK1+DWkDcgguyFggqzewpxC48BwYD8SiD6UFH2TUXOpAPEuwe3uMAqGifxVnujGvcp9VPhLKCzVz8MCAbkDYEF6tzcsoVrcO/B5QyISf85QBVBGS/6kqTkdMa4hc8wIM45fxiu6tdZG2sN3d6AROqnjhWtn8U88RN57x7EgKhFMifuU+qnoi3cV+oUvfq1NyAqcxVQZFx0377yIWrkulTmCesPW0bhGtx7EAMiyN652J6n7xYeAyIIfzHEagQYEL6CdbHEyh+3ajzR8EWv3dtDtzCQKsQFdPbqhwFRxX8nLuFO1EK/BercVK1wDW4NMSCmnPGXJCGlN8QtPAbEi9gfbzUCDEj83rKYJ16IvXtwiwtsIW//hvz4iU+pnwpxAZ29+mFAVPExIG8IJBiwC2q8frRwr7r3IAZEkL1zsfEGZLSEBeEvhliNAAOCAblYYuWPWzWeaPii1+49PLcwkCrEBXT26ocBUcXHgGBAtFpx78HR21v75qs2TIUrBmS0TBS678eo+nXWxloDBgQDMr5Dejxp1TgGpIdOV2exgM7uy0/3f4ZXZX5V25HnE87dFvqpzBPWPyLDu88UrsGt4ejNEgMSVi4hA7mF5ytYIdw/HMRqBBgQDEh+FeZmsGocA5LLf9boC+g8fAbOYujNozL3jhsRn3ABb6Gfyjxh/RGyPI1RuAa3hhgQQfbOxfY8fbfwGBBB+IshViPAgGBALpZY+eNWjWNAyiUKmcACOg+fgSGAEgZRmSekNodMuBO10E9lnrB+k7kaULgGt4YYEEHVzsWGARktYUH4iyFWI8CAYEAullj541aNY0DKJQqZwAI6uy8/fAVrvDQS7kQt9Fugzk3RCtfg1nD09sZXsMwymBrgFp43IPn6WI0AA4IBya/C3AxWjWNAcvnPGn0BnYfPwFkMvXlU5t5xI+IxIKNX5wj6H4+h1k0HDUcpYkDy68iTYbj53m634+HhlPP+T0Khmmsr3ETm3NQAaw0YEAyIWktd46wax4B0Vc43rwV0Hj4DfSTmRavM583o/5kS7gQt9FOZJ6w/TMbCNbg1xIAIsncutufpu4XnDYgg/MUQqxFgQDAgF0us/HGrxjEgukQqS33E+ZGFZ+XwGTifUu+MSh0m6NxCP2XtvdXTZ9dBQwyIoFeCUEJWV8jw5uUNiIuzK9hqZhgQDIiroBoGWzWOAdFFU1nqI86PLDwrh8/A+ZR6Z1TqMEHnFvopa++tnj67DhpiQAS9EoQSsrpChjcvBsTF2RVsNTMMCAbEVVANg60ax4Dooqks9RHnRxaelcNn4HxKvTMqdZigcwv9lLX3Vk+fXQcNMSCCXglCCVldIcObFwPi4uwKtpoZBgQD4iqohsFWjWNAdNFUlvqI8yMLz8rhM3A+pd4ZlTpM0LmFfsrae6unz66DhhgQQa8EoYSsrpDhzYsBcXF2BVvNDAOCAXEVVMNgq8YxILpoKkt9xPmRhWfl8Bk4n1LvjEodJujcQj9l7b3V02fXQUMMiKBXglBCVlfI8ObFgLg4u4KtZoYBwYC4CqphsFXjGBBdNJWlPuL8yMKzcvgMnE+pd0alDhN0bqGfsvbe6umz66AhBkTQK0EoIasrZHjzYkBcnF3BVjPDgGBAXAXVMNiqcQyILprKUh9xfmThWTl8Bs6n1DujUocJOrfQT1l7b/X02XXQEAMi6JUglJDVFTK8eTEgLs6uYKuZYUAwIK6Cahhs1TgGRBdNZamPOD+y8KwcPgPnU+qdUanDBJ1b6Kesvbd6+uw6aIgBEfRKEErI6goZ3rwYEBdnV7DVzDAgGBBXQTUMtmocA6KLprLUR5wfWXhWDp+B8yn1zqjUYYLOLfRT1t5bPX12HTRczoCoBRIJtyKnXkZPkcObt6sBca7fDK/Q0MrZzYBE7plTEGv9Z0xFzoy8A3vw4TiOm1m4zQMUjZN4R5Px9tBw/TZiGa2NMp5Xv3PMcA2ViXaPUeowum8P9E/0u1NIin5JPdm9BzEgQjcoFFSYHQZEgVShoZUTA4IBUWq3c4xV4y9zT7iwRGPxHp7hl9eNWEZro4zn1Y8L7IULbMJ+Rj+lyoWYwj7i1hAD0ltQYXYYEAVSxaa0cmJAMCBK7XaOsWocA6KrtxFLfdFxke7LD29A3oev1CEGJK5wo0dS9OMNyCD1CrgVOZ14hpsvX8F6TTqysVp1gwHBgDj3ebtwq8YxILpkG7HUFx0XOXwGxk1hj5GUOow8J5+poV9Q+Sj6YUAGYVfArcjpxDO8eTEgGJCsS6Kyb6IPMiVnk+Yb/hUeZ88ICS/kHTL/nwbx9tBw/TZiGa2NMp5Xv3PMcA2ViXaPUeowum/zOyBxVaHo1+QMPPgKlqB7oaDC7J5ChpsvBgQDggFRt9mHcd49uMXlZ4HeqIpbrt9GLFXmkXFe/TAgd+grdYgBiSzd2LEU/TAgg8wr4FbkdOIZbr4YEAwIBsS5294P9+5BDEgI9rBByvVb4JwJg50wkFc/DAgGJKEM64cs7CPuPcgbEKFeCgUVZscbEAVShYZWTn4HhN8BUWq3c4xV41nmNoGJ9/AMN5AbsUyQxxzSqx8GBANiFtWKAYV9xL0HMSBChRUKKswOA6JAqtDQyokBwYAotds5xqpxDIiu3kYs9UXHRbovP/wOyPvwlTrkK1hxhRs9kqLfmbODhhgQQf1CQYXZYUAUSBUaWjkxIPFN0GKeeCH2XoDCP0FX9kF0TCHv6KWU67cRy2htlPG8+vEGhDcgSl0tF1PYR9x7sI0BKYQWVmCFa3AL//Lpz+zfAVEZhYniHCjyUwFrrRgQDIizPNuFWzWeaPiiWXh7aLiB3IhltDbKeF79ygyIqrOy6KqYyHPyeQ3L6FfFXM2r1lcHDTEgqqpC3ErCY0DeFzRyU1r1gAHBgAhtpXWIVeMYEF2+jVjqi46LXOYCq+ochyZ+pMhzEgMSq49aXx00xIAEar+S8BgQDMivBKIbkrIfKnKe647OO/BPYYd/gh7YyuShFI2TeMtzFAO9F9hw/TZiKSIPDfPqxxuQC/gb9M8y/S5gm/JoYR9x70EMSGBJrCQ8BgQDggEJ3Pz+v8UTfoENXY04WGHPE2coh3kPz3D9NmIpQw8M9OpXdoFVdQ5kEz4UBiQcadiAan110BADEib7cawkPAYEA4IBCdz8GJAPYSYcdqHidXiDVXh+RLOsGA8DMpF6wn5eRr+JmIdSFfYRt4YYkCGJ339oJeExIBgQDEjg5seAYEAullPh+XFx5i0ed19+qv4ZXlXnFlTvTAID0lcdtb46aIgBCayjlYTHgGBAMCCBmx8DggG5WE6F58fFmbd4HAMyUYYOl9cqAzkR81Cqwj7i3oMYkCGJeQMyik3dHKPjX30usrFaa+VfwYr/ZXCL+Ut9ROr8PKa3+Yb/DsHV2h95vpD3yHQ/eqZcv41YRmujjOfV7xyzZA+qOiuLropp0D/L9KtiruZV66uDhhgQVVUhbiXheQPCGxDegAibWg/xXoBKLj/6crTIwp6nTVCPKtdvI5Y69bhIr35lF1hV5zg08SN1uLxWGch4mrEjqvXVQUMMSKD2KwmPAcGAYEACNz9fwfoQZsJhFyoev4QejXP6eBiQicgT9vMy+k3EPJRqpXtougEphDEk3oyHEpikbt6E+c7AnJYjgceW+imcog8yJedZGNF5O1xg1YJXGanjVcTtqB+6XKqk1B56aWY8rBBAP4WSEKP2kQY99MCACIJGhyQUSOrmTZhvNNKp4yXw2FI/hVN0E1RyYkB+TN0vGcmi66aDgVRrN4Nn1JgJuqhTS+2h6iSIGyaAfsPoXj+o9pGEverWEAMSJLpnmIQCcQvv+f5kwnw9uNrFJvDYUj+FU3QTVHJiQDAg7zQV7x4M/x0etXbbNcSfJhS9nx1r9ep3Dh2uoWO+hL4mgH5BFaH2kYS96tYQAxIkumeYhAJxC+9pvgnz9eBqF5vAY0v9FE7RTVDJiQHBgGBActpq9H52zDK1hzrmQegYAfQb4/bmqYXOQL6CFaS5a5iEAkndvAnzdfHqFpzAY0v9FE7RFxYlJwYEA4IByemq0fvZMcvUHuqYB6FjBNBvjBsG5CNuhReCIDnjh0lgkrp5E+YbD3XiiAk8ttRP4RR9YVFyYkAwIBiQnIYZvZ8ds0ztoY55EDpGAP3GuGFAMCC+ykm4JKVu3oT5+oA1i07gsaV+CqfoC4uSEwOCAcGA5DTV6P3smGVqD3XMg9AxAug3xg0DggHxVU7CJSl18ybM1wesWXQCjy31UzhFX1iUnBgQDAgGJKepRu9nxyxTe6hjHoSOEUC/MW4YEAyIr3ISLkmpmzdhvj5gzaITeGypn8Ip+sKi5MSAYEAwIDlNNXo/O2aZ2kMd8yB0jAD6jXHDgGBAfJWTcElK3bwJ8/UBaxadwGNL/RRO0RcWJScGBAOCAclpqtH72THL1B7qmAehYwTQb4zb5zUgt9vteHg4/ynt+z+FDSlIzrphrMvU4+Pj8fXruW+PFpvXmm8dydjMCTWdqp+qS/S6lLwVOT+7AYndDduM5t2D/A2JXtJ79Ttnj4Z9NES/IC2Uc7fJGXjtn+HFgARVzJ1hrELCgOTyvzd69KU520BadfSyzuh1KXkrcjZpvlx+arbvvazeCxD6ra0fBgT9ehEImo1y7jY5AzEgQZqnDGMVEgYkBbs5aPSlGQNiIpcDrD2TZbYGNOQCK6s6JRADMgVzWhKvfhiQNCmGBka/IWxvH1roDMSABGmeMoxVSBiQFOzmoBgQE9FTgFW/GZ/CKDkz8mJAtJpoHOW9AGEge4np1Q8Dgn69CATNZqEzEAMSpHnKMFYhYUBSsJuDYkBMRBgQGxEXWJvRzAjvBRb9Zqpj5/LqhwGxmc6MQL8g2ta98SVNh3vMl8E1PzVffgdkkJ74mFVIGBARZHBYh43r+QVKq46yGpKSN5qlkpM3IMEbYo/hvBcgDEgv3b36YUDQrxeBoNksdAbyBiRI85RhrELCgKRgNweNvjQPfH3HdXhadYQBMSVXArwXIC6wCtV5Meg3j3VGJq9+rh6aMWHGfEUA/YIKouq8H7nH8AYkSPSMYaxCwoBkULfHxIDYjM4Iq34z3kQoOTPyDjRfDIhWRrOivBcg9JuljJbHqx8GROM6Kwr9gkgvdAbyBiRI85RhrELCgKRgNwfFgJiIngKs+s0wAkrOjLwYEK0mGkd5L0AYkF5ievXDgKBfLwJBs1noDLxmQL5//27+IcLnP5QXhPZzDXP+js1HP6cB+fbt2xnSovla891FvYSaTtVP1SV6XUreipxPG+a/P+AZ+ePVkAtsJP3rY6HfdYaVI3j1w4BUqvU2N/oF6aGcu03OwGED8sdxHH8F8WKY6wT+PI7jb+cwaOgElhiOfolwJw3t1ZD9N0kYMQ36iaCahnn1O5fBHuwjJvr10WJ0Jm4NR38H5Hzu9+M4/h2dKc+FEfjtOI5/zm+8OEdEQyewpHD0SwI7cdgRDdl/EwUyUqFfHy1GZjKi35mHPThCO/4Z9ItnOnvEIQ1HDcjsxZEPAhCAAAQgAAEIQAACENiAAAZkAxFZAgQgAAEIQAACEIAABFYhgAFZRSnmCQEIQAACEIAABCAAgQ0IYEA2EJElQAACEIAABCAAAQhAYBUCGJBVlGKeEIAABCAAAQhAAAIQ2IAABmQDEVkCBCAAAQhAAAIQgAAEViGAAVlFKeYJAQhAAAIQgAAEIACBDQhgQDYQkSVAAAIQgAAEIAABCEBgFQIYkFWUYp4QgAAEIAABCEAAAhDYgAAGZAMRWQIEIAABCEAAAhCAAARWIYABWUUp5gkBCEAAAhCAAAQgAIENCGBANhCRJUAAAhCAAAQgAAEIQGAVAhiQVZRinhCAAAQgAAEIQAACENiAAAZkAxFZAgQgAAEIQAACEIAABFYhgAFZRSnmCQEIQAACEIAABCAAgQ0IYEA2EJElQAACEIAABCAAAQhAYBUCGJBVlGKeEIAABCAAAQhAAAIQ2IAABmQDEVkCBCAAAQhAAAIQgAAEViGAAVlFKeYJAQhAAAIQgAAEIACBDQhgQDYQkSVAAAIQgAAEIAABCEBgFQIYkFWUYp4QgAAEIAABCEAAAhDYgAAGZAMRWQIEIAABCEAAAhCAAARWIYABWUUp5gkBCEAAAhCAAAQgAIENCGBANhCRJUAAAhCAAAQgAAEIQGAVAhiQVZRinhCAAAQgAAEIQAACENiAAAZkAxFZAgQgAAEIQAACEIAABFYhgAFZRSnmCQEIQAACEIAABCAAgQ0IYEA2EJElQAACEIAABCAAAQhAYBUCGJBVlGKeEIAABCAAAQhAAAIQ2IDA/wBvLMSvqG46WAAAAABJRU5ErkJggg==\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "f = np.load(\"patterns.npz\")\n",
    "\n",
    "# extract the labels and data\n",
    "labels = sorted(f.keys())\n",
    "data = np.array([f[key].reshape((10, 10)).T.ravel() for key in labels])\n",
    "\n",
    "# reverse numbers\n",
    "data[labels.index('num1')] = 1 - data[labels.index('num1')]\n",
    "data[labels.index('num2')] = 1 - data[labels.index('num2')]\n",
    "data[labels.index('num3')] = 1 - data[labels.index('num3')]\n",
    "data[labels.index('num4')] = 1 - data[labels.index('num4')]\n",
    "\n",
    "# display them\n",
    "fig, axes = plt.subplots(1, 7)\n",
    "for j in range(7):\n",
    "    axes[j].imshow(data[j].reshape((10, 10)))\n",
    "    axes[j].set_xticks([])\n",
    "    axes[j].set_yticks([])\n",
    "fig.set_size_inches(10, 2, forward=True)\n",
    "plt.tight_layout()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Define some helper functions:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "def sigmoid(x):\n",
    "    \"\"\"Sigmoid function\"\"\"\n",
    "    return 1.0 / (1.0 + np.exp(-x))\n",
    "\n",
    "def dsigmoid(x):\n",
    "    \"\"\"Derivative of the sigmoid function\"\"\"\n",
    "    return sigmoid(x) * (1 - sigmoid(x))\n",
    "\n",
    "def add_noise(X0, p=0.5):\n",
    "    \"\"\"Adds random zero-masking noise\"\"\"\n",
    "    X1 = X0.copy()\n",
    "    idx = np.random.rand(*X1.shape) < 0.5\n",
    "    X1[idx] = 0\n",
    "    return X1\n",
    "\n",
    "def init_plotting(data, M, N):\n",
    "    \"\"\"Initialize figure and subplots\"\"\"\n",
    "    fig, axes = plt.subplots(4, 7)\n",
    "    axes = axes.T\n",
    "    ims = np.empty(axes.shape, dtype=object)\n",
    "    for j in range(k):\n",
    "        ims[j, 0] = axes[j, 0].imshow(data[j].reshape((10, 10)))\n",
    "        ims[j, 1] = axes[j, 1].imshow(M[j].reshape((10, 10)))\n",
    "        ims[j, 2] = axes[j, 2].imshow(N[:, j].reshape((10, 10)))\n",
    "        ims[j, 3] = axes[j, 3].imshow(data[j].reshape((10, 10)))\n",
    "\n",
    "    for ax in axes.flat:\n",
    "        ax.set_xticks([])\n",
    "        ax.set_yticks([])\n",
    "        \n",
    "    axes[0, 0].set_ylabel(\"Corrupted input\")\n",
    "    axes[0, 1].set_ylabel(\"Encoder weights\")\n",
    "    axes[0, 2].set_ylabel(\"Decoder weights\")\n",
    "    axes[0, 3].set_ylabel(\"Decoded image\")\n",
    "\n",
    "    fig.set_size_inches(12, 7, forward=True)\n",
    "    plt.tight_layout()\n",
    "    return fig, axes, ims\n",
    "\n",
    "def draw(fig, data, M, bM, N, bN):\n",
    "    \"\"\"Draw the current weights, as well as an example\n",
    "    corrupted input and reconstruction from that input.\n",
    "    \n",
    "    \"\"\"\n",
    "    k = data.shape[0]\n",
    "    for j in range(k):        \n",
    "        X1 = add_noise(data[[j]])\n",
    "\n",
    "        # compute hidden representation\n",
    "        Y = sigmoid(np.dot(M, X1.T) + bM)\n",
    "        # reconstruct the input\n",
    "        X2 = sigmoid(np.dot(N, Y) + bN).T\n",
    "\n",
    "        ims[j, 0].set_data(X1.reshape((10, 10)))\n",
    "        ims[j, 1].set_data(M[j].reshape((10, 10)))\n",
    "        ims[j, 2].set_data(N[:, j].reshape((10, 10)))\n",
    "        ims[j, 3].set_data(X2.reshape((10, 10)))\n",
    "\n",
    "    fig.canvas.draw()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false,
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\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",
       "            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",
       "        this.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 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);\n",
       "        canvas.attr('height', height);\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'];\n",
       "    var y0 = fig.canvas.height - msg['y0'];\n",
       "    var x1 = msg['x1'];\n",
       "    var y1 = fig.canvas.height - msg['y1'];\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, fig.canvas.height);\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;\n",
       "    var y = canvas_pos.y;\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 overriden (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",
       "    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 + '\">');\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 dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\">';\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",
       "        event.shiftKey = false;\n",
       "        // Send a \"J\" for go to next cell\n",
       "        event.which = 74;\n",
       "        event.keyCode = 74;\n",
       "        manager.command_mode();\n",
       "        manager.handle_keydown(event);\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,iVBORw0KGgoAAAANSUhEUgAAA8AAAAIwCAYAAABJMwI2AAAgAElEQVR4Xuy9CZg0ZXX2f78oKtv74hJDBDSKURJjXDDGuGsw4gJGDCjGoBJFRRFFVOAfwR31L1HADQU3YvzA5XNBDAbjEjTGJR/G5QNcoiyKGpEZQBQRvqume2b6nXdm6tzVdbqfqv71danInHrq1H3u5zzPr6u6e4N4oQAKoAAKoAAKoAAKoAAKoAAKoMAMKLBhBq6RS0QBFEABFEABFEABFEABFEABFEABAcCYAAVQAAVQAAVQAAVQAAVQAAVQYCYUAIBnosxcJAqgAAqgAAqgAAqgAAqgAAqgAACMB1AABVAABVAABVAABVAABVAABWZCAQB4JsrMRaIACqAACqAACqAACqAACqAACgDAeAAFUAAFUAAFUAAFUAAFUAAFUGAmFACAZ6LMXCQKoAAKoAAKoAAKoAAKoAAKoAAAjAdQAAVQAAVQAAVQAAVQAAVQAAVmQgEAeCbKzEWiAAqgAAqgAAqgAAqgAAqgAAoAwHgABVAABVAABVAABVAABVAABVBgJhQAgGeizFwkCqAACqAACqAACqAACqAACqAAAIwHUAAFUAAFUAAFUAAFUAAFUAAFZkIBAHgmysxFogAKoAAKoAAKoAAKoAAKoAAKAMB4AAVQAAVQAAVQAAVQAAVQAAVQYCYUAIBnosxcJAqgAAqgAAqgAAqgAAqgAAqgAACMB1AABVAABVAABVAABVAABVAABWZCAQB4JsrMRaIACqAACqAACqAACqAACqAACpQKwFVet5F0JSXqjAI7SPqRpBsSMsYPCaJOYMgsT+CHCRQv4RRZfqhSxRMJBZvAkFmewA8TKF7CKbL8QI9IKNYEhsQPExC5Y6dozROlAvDOki7pWFFIV9pF0qUJQuCHBFEnNGSGJ/DDhIqXcJoMP1Rp4omEYk1oyAxP4IcJFS/hNBl+oEckFGpCQ+KHCQndodO04olSAXijpLkOFaO3qc7N1Zdhfn5eu+66a6XBJknzCWKk+SFyfU2uZ9OmSorYy8nBGTd29vSoDE9Yfui5vuECFqJDhh8qDSxPhEUjcBIKZHgCP0yicjnnyPADPSKnVpMYFT9MQuVunaMVT2QB8JGSXrOKni+S9LqAzixeAZEmEXLDDfVPNFcAPASzVky5ynWl+SFyfU103rAhPrWcHJxxm+SdcEyGJyw/9FzfcMkK0SHDD2xuwy4oMjDDE1aPKFKV2U0qww/0iO76CT90t3ZZmbfiifgu3buM6i5gtQCtfF0u6RaBoVi8AiJNIiSyaQaAt6yEA6oRjRfP4Iw7CX8EztFKo1pxHqs/9FzfQAkGIYXokOEHNrdhFxQZmOEJq0cUqcrsJpXhB3pEd/2EH7pbu6zMW/FE2wC81fDLSK4YAvDo+LtL+oyk3w0owuIVEGkSIZFNMwAMAK/jxVYaFQA8/myPzOUJvNGS4Qc2t+PbY5ojZHiCPcQ0KzreuTP8QI8YrybTPBo/TFP9Ms/diifaBuDra74F+CRJzwvoyeIVEGkSIZFNMwAMAAPAk5iN450jMpcB4PE05uhGCrSymRnnTbJGWXNQlgIZfgCAs6qVPy5+yNe4a2doxRNtA/CDhneAz5L0iBFFKzC+TNJ3gioDwEGhssMim2YAGAAGgLNn4vjjR+YyADy+zoxgK9DKZgYAtnUv9YAMPwDApVa7Pi/8UK/RrEW04om2AXixCNVXAl88RkUA4DHEa/PQyKYZAAaAAeA2Z13OWJG5DADnaM+o6yrQymYGAO6NyzL8AAB31x74obu1y8q8FU9kAfBD17nqfw0oAgAHRJpESGTTDAADwADwJGbjeOeIzGUAeDyNObqRAq1sZgDgRtqXeFCGHwDgEisdywk/xHSapahWPJEFwNUjzytfi7+nc6NAlQDggEiTCIlsmgFgABgAnsRsHO8ckbkMAI+nMUc3UqCVzQwA3Ej7Eg/K8AMAXGKlYznhh5hOsxTViieyAHhlIW4j6ThJH5b00UCVLAAuZGMXuKwyQhL1asWUq6hk+aEMlbuVRcc80Tk/JOobNloJOUjqdY/I0jhr3LB5cgMzPNG5HpErcadGz/BDMQDcpblcSK699kOnZmY5ybbiiUkB8GLz+YqkOwc0tBavQiZp4LLKCEnUqxVTAsCT90nHPGH1h8mrueUZE/UNX14JOQDAy+VyftO7kNqFvWYGZqwbnesRpmZ9Ds/wAwDcwDGF9J1e+6FBWTikpTfSJwnAlYmrL8aqFqa6l7V4FTJJ666pmL8n6kWjKqbKXiId84TVHzwlcqIT9Q0nXEIOADAAvIphM9aNzvWI8ETuf2CGHwDgBr5hzWggGodMQoFWekQWAB+0QoHtJB0gaV7SXgF1rMWrkEkauKwyQhL1asWUq6hk+aEMlbuVRcc80Tk/JOobNloJOQDAADAAHJ6ysxrY631EIX045K1Ccu21H0KFIGilAq14IguA/3tFtldK+qqkv5f0o0AtrQ1uIZM0cFllhCTq1YopAeDJ+6RjnrD6w+TV3PKMifqGL6+EHABgABgADk/ZWQ3s9T6ikD4c8lYhufbaD6FCENQpAB63XNYGt5BJOu41T+z4RL1oVBOrYrsn6pgnrP7QrlLNRkvUN5xQCTkAwAAwAByesrMa2Ot9RCF9OOStQnLttR9ChSCokwC8s6Rdh5/9vdSoobXBLWSSGpc33dBEvWhU0y1t47N3zBNWf2gsSosHJuobzrKEHABgABgADk/ZWQ3s9T6ikD4c8lYhufbaD6FCENQpAP4dSf8o6WHDrKvfAP60pCdJ+mmgltYGt5BJGrisMkIS9aJRlVFiO4uOecLqD7YYCQck6hvOtoQcAGAAGAAOT9lZDez1PqKQPhzyViG59toPoUIQ1CkA/pCkm0o6QlL1eeDbS3qdpOsk7RuopbXBLWSSBi6rjJBEvWhUZZTYzqJjnrD6gy1GwgGJ+oazLSEHABgABoDDU3ZWA3u9jyikD4e8VUiuvfZDqBAEdQqALx9C79xI1jsOYfjmgVpaG9xCJmngssoISdSLRlVGie0sOuYJqz/YYiQckKhvONsScgCAAWAAODxlZzWw1/uIQvpwyFuF5NprP4QKQVCnAPh7ku4paRSAK/Ctvgl6t0AtrQ1uIZM0cFllhCTqRaMqo8R2Fh3zhNUfbDESDkjUN5xtCTkAwAAwAByesrMa2Ot9RCF9OOStQnLttR9ChSCoUwD8RElPkHSkpB9I+n1Jr5Z0+vA/ixdz/Rp1TdvgFjKh+2pnGlVfK9v8ujI8YfWHEuY8OSwZKMMP1eCWJ5rbmSMTFMjwRJofnLmcoNXCkBs2xH/B0snXGTfr2vr+Jlmibn0dOqM/sGZ02y2teCLeRT2xFsG2+vKrxVd1rtH/X/37GwHAnrCFR7diylWuMW0zU7iefUgvwxOWH0rYAJIDANyHyZx0DVPvEc51OXPZGdeJdUDVydcZ18nXjM3wA8BjFqGgcPxQUDEKSaUVT2QB8IOCIn0OAA4q1Y2wVkwJAHej2MEsMzwBAAfFHw0rZCOc4Qc2tw38UNAhGZ6weoSjhTOPnHGdWAdUnXydcZ18zdgMP9AjzCIUFI4fCipGIam04oksAB5XoyIWr0IWg3G1nOTxrZgSAJ5kydLPleEJqz+UsAEkhyWfZfiBzW36NE49QYYnrB7hXJ0zl51xnVhnb+Lk64zr5GvGZviBHmEWoaBw/FBQMQpJpRVPZALwHSTtIWmHFYK9MyBgEYtXIYtBQK5iQloxJQBcTD3bSCTDE1Z/KGEDSA4AcBuTqadjTL1HOLo6c9kZ14l19iZOvs64Tr5mbIYfAGCzCAWF44eCilFIKq14IguAnynpTZKqn0O6ekSw6jPAFRjXvawNbt1go3/v4GLgXN60Y1sxJQA87TK2ev4MT1j9oYQ5Tw4AcKuzql+DTb1HOHI6c9kZ14l1QNXJ1xnXydeMzfADAGwWoaBw/FBQMQpJpRVPZAHwDyU9X9KHG4plbXCdc3RwMXAub9qxrZgSAJ52GVs9f4YnrP5QwpwnBwC41VnVr8Gm3iMcOZ257IzrxDqg6uTrjOvka8Zm+AEANotQUDh+KKgYhaTSiieyAPgKSTuOIZS1wXXO08HFwLm8ace2YkoAeNplbPX8GZ6w+kMJc54cAOBWZ1W/Bpt6j3DkdOayM64T64Cqk68zrpOvGZvhBwDYLEJB4fihoGIUkkornsgC4PdJeruktb7luU5Da4NbN9jo3zu4GDiXN+3YVkwJAE+7jK2eP8MTVn8oYc6TAwDc6qzq12BT7xGOnM5cdsZ1Yh1QdfJ1xnXyNWMz/AAAm0UoKBw/FFSMQlJpxRNZAPxGSU+W9CFJP1oh2DEBAa0NbmC8pZAOLgbO5U07thVTAsDTLmOr58/whNUfSpjz5AAAtzqr+jXY1HuEI6czl51xnVgHVJ18nXGdfM3YDD8AwGYRCgrHDwUVo5BUWvFEFgB/Zg2Rqi/BemhAQGuDGxgPAHZEah7biikB4OYFKPDIDE9Y/aGEDSA5AMAFzs1SUpp6j3CEcOayM64T64Cqk68zrpOvGZvhBwDYLEJB4fihoGIUkkornsgC4HE1sja4456s1ONZuJYqs+CHubk5bdxY/eP6r0IW8bo0Z+XvrTSqFWJ1rj90cC5n+TPDD73f3PbcPxme6FyPyJpwHRw3ww+97xEdrHM0ZfwQVWp24lrxBABcsGE6uOlpxZSrlAQALtinNalleKJzm9sOzuUsx2X4ofeb2577J8MTnesRWROug+Nm+KH3PaKDdY6mjB+iSs1OXCueaBOA/1nSXkP9/01S9bjzaq8HBmrE4lUJeMNaEm6pYCF3PVsxJQAcmCHdCcnwROf6QwfncpbDMvzQ+81tz/2T4YnO9YisCdfBcTP80Pse0cE6R1PGD1GlZieuFU+0CcBHSTpuqP+x69ThZYEasXgBwKM24Q5wYNIUGtJKo1pxbZ3rDz0HGMd6GX7o/ea25/7J8ETneoQziXoem+GH3veIHnsCP/S4uA0vrRVPtAnADa9j1cNYvABgALjNGTW9sVppVADw9ArY8pkz/ND7zS0AbLuQPYQtWTEH0COKKUURieCHIspQVBKteAIALqqmmyfTwU1PK6ZcpSTcAS7YpzWpZXiic5vbDs7lLMdl+AEAHqlWIR+HcfyT4YnO9QhHsJ7HZvih9z2ix57ADz0ubsNLa8UTAHBD9SdxWAc3za2YEgCehLsmdo4MT3Ruc9vBuZxlkAw/9H5z23P/ZHiicz0ia8J1cNwMP/S+R3SwztGU8UNUqdmJa8UTAHDBhungpqcVUwLABZvSTy3DE53b3HZwLvuVjh2R4Yfeb2577p8MT3SuR8Smz0xEZfih9z2ix87ADz0ubsNLa8UTAHBD9SdxWAc3Pa2YEgCehLsmdo4MT3Ruc9vBuZxlkAw/9H5z23P/ZHiicz0ia8J1cNwMP/S+R3SwztGU8UNUqdmJa8UTAHDBhungpqcVUwLABZvSTy3DE53b3HZwLvuVjh2R4Yfeb2577p8MT3SuR8Smz0xEZfih9z2ix87ADz0ubsNLa8UTbQLwZ9b57d/Ra3xo4IJZvPgW6FGb8CVYgUlTaEgrjWrFtXWuP/QcYBzrZfih95vbnvsnwxOd6xHOJOp5bIYfet8jeuwJ/NDj4ja8tFY80SYAv2LkQnaU9HeSzpL035J+X9IjJZ0q6bmBC+7c4uVsUALXvxDCt3kuKdU5P0RrPANxrTSqcQDYmZsdnHNds1CGH9jcds0Fm+eb4QnWjKHGHex/GX7odY9wauy0ikLWw875walHIRo7tighthVPtAnAo6J8UNK7JZ058i8fJempkv46oF7nFi/H8IHrB4A3F6lzfojWeAbiWmlUAHBvnJLhh15vbntT+bUvJMMTrBkA8ErH9dYTGfvPgm7CZPSH1DXDqQcA3GiFa8UTWQA8J+nmkq4fubStJP1CUpV43atzjcoxfN3FL/69gxOjFVOuok/n/BCt8QzEZXjC8oMzNzs457pmoQw/pG5muiZwB/PN8ITVIzqoWTjlDva/DD/0ukc4NQ4bp5ynEDvnB6ce7DkcRy7FtuKJLAC+UNLzho9AL2b8CEknSvqDwOV2bvFyDB+4/oWQDk6MVkwJAEcd0om4DE9Y/cGZmx2cc50wwUiSGX7o9ea2awVukG+GJ6we0SDnzhzSwf6X4Yde9winxo5xC1kPO+cHpx6FaOzYooTYVjyRBcB/O/y870dHPgO8j6SDJb03oF7nFi/H8IHrB4A3F6lzfojWeAbiWmlUK3Sy/ODMTRajdEdm+KHXm9v0ikz/BBmesHrE9CXIy6CD/S/DD73uEU6NHacVsh52zg9OPQrR2LFFCbGteCILgCuB7ivpKZJ2kXSJpPdI+kJQuc4tXo7hgxpwB3hZqM75IVrjGYhrpVEBwL1xSoYfer257U3l176QDE+wZgz1dvYmhWzGM/zQ6x7h1NjpJ/jBUWs51qlHIRo3u9DpHdVKj8gE4HGk6dzi5Rg+KkwHJ0YrplxFn875IVrjGYjL8ITlB2dudnDOdc1CGX7o9ea2awVukG+GJ6we0SDnzhzSwf6X4Yde9winxo5xC1kPO+cHpx6FaOzYooTYVjyRCcCLd4B/T9Leku4paVtJ5wbU69zi5Rg+cP0LIR2cGK2YEgCOOqQTcRmesPqDMzc7OOc6YYKRJDP80OvNbdcK3CDfDE9YPaJBzp05pIP9L8MPve4RTo0d4xayHnbOD049CtHYsUUJsa14IguAHy/p7ZJOl/QESdVidC9Jr5P00IB6nVu8HMMHrh8A3lykzvkhWuMZiGulUa3QyfKDMzdZjNIdmeGHXm9u0ysy/RNkeMLqEdOXIC+DDva/DD/0ukc4NXacVsh62Dk/OPUoRGPHFiXEtuKJLAD+hqRnSPri8KePqp9EuqmkiyXdOqBe5xYvx/CB6weAAeCoTUqPa6VRAcCllzmcX4Yfer25DSvb3cAMT3RuD5FVPmdvUshmPMMPve4RTo0dn+EHR63lWKcehWjc7EKnd1QrPSILgK+QtONQm8sl3aJ6olfSz4f/XCdb5xYvx/B1F7/49w5OjFZMuYo+nfNDtMYzEJfhCcsPztzs4JzrmoUy/NDrzW3XCtwg3wxPWD2iQc6dOaSD/S/DD73uEU6NHeMWsh52zg9OPQrR2LFFCbGteCILgM+T9FRJ/0fSIgBXj0C/bfgodJ2AvV28ej4xWjFlHwC4z3WOXNv8/Lw2barsoAxPWP0hkm+TN52ccesa3ujfe74gZvih15vb6uIcr3XQPxmesHqEMz+dWKduzrgdrLFzeRl+6H2PcASOxjr+TfQkfhgWLKsezrhR71RxpXsiC4APlPRySa+SdLykQyS9RNJLJb0/IGARi1cgTzvEMVqieey8gwfQqJIbVbAOqWERDwPAzUvQwXnvXCw9wlFrBvpJCW+SNShJ6JBIrwwNtCKIHtFEtYXvoplrdOSMHuT4N9GTrBnJ64BTZ2cqlO6JLACuNKp+A/j5kv5A0o8lnSDpxKB4vW1UjtESzRMsgx1Go0puVHZFEg6IeBgAbi58B+e9c7H0CEetGegnALBvCHqEr9nwy1gBYEO6yFq/OFyiJ1kzktcBp86GfWb2DrCj0WqxAHDu4wPj1met42lUyY0qq3DOuJFmCQA7im4em7iRaJ5Ue0fSIxpoGZlzE9iINsg8dEiGJ4rYQzh1Cyk1DKJHOGotxRbhiUaZT+kgx7+JnszoD5WinfNDVj2ccR0rlu6JrDvA1bdA33UVoarPBt89IGDnjBm4poUQx2iJ5omm68bRqADgBQUAYHfqLMd3cN47F0uPcNSagX7CHWDfEPQIX7MuAk+jq2zxoEL2q6wZyeuAU2fHXol9qhVPZAHwlZJ2WEWoX0iqfhKp7gUAcwd41COd84PTUBKbRN08a/T3yLUBwI2kXTioa34wr7SVhWuVc3auRzi6Rebc4ngd9E+GJ4rwg1M3xw8drLFzeRl+qM5fhCccIaYd6/g30ZP4AQBeORVa8UTbAHzQMMs3SXr28KePFhO/s6THSNo9MKl726gKaSiBEjQKacWUfdjc9rnOkWsDgBvNHwC4uWy9XTMqSSJzDgDezDxF+MGpm2P9RNhw0siKZR+Rpaw5ruPfRE/iBwC4EwD838MsbyvpopGMr5d02fCboc8OzMEiFq9AnnZIIQ3Fzjt4AI0quVEF65AaFvEwANy8BIkbieZJtXckPaKBlpE5BwADwA2sVeIh9IhCqlJI38EPyftKp86ONRP3Mq14ou07wIvanCXpkY5QK2IB4G4+CtmKKVfxTef84DSUxCYxxhRc+9DItQHAzaXvmh/MK6VHmIJV4ZE5BwADwA2sVeIh9IhCqlJI38EPAPDKGdGKJ7IAeGktlvS7kn5SreHGnO4c8ESvrZCGEk3XjWvFlACwK/tk4yMeBoCb1wQAbqRdb9cMALi7foj0yiZXR49oohqfAXZVc/yb6En2lQBwpwB4W0lvlHSgpK0l/UbSe4e/C3x1YBL2djNTSEMJlKBRCI0quVE1qkrLB0U8DAA3Fz1xI9E8qfaOpEc00DIy55bedd6Q/b52gwtY/5AMTxSxh3Dq5qhKj3DUWootwhONMp/SQY5/Ez2Z0R8qRTvnh6x6OOM6VizdE1kr5anDn0E6RlL1ueDbS3qppG9J+ruAgEUY0zFFYqEDchUTQqNKLkUJnozkMAsA7JQ6olk2wJSQQ9JP3nRyM+P4p+exGetGEXuIntct6/Iy/ECPaFAt1owGovX8kD55IguAfybpTyT9eMQLt5H0X5JuFfBHEYtXIYUOyFVMCAtXcilK8GQkBwB4cyNENAOAx5o8RawZY13B7B6csW7gh+76KcMPAHADP5SwbvGmaYPCJR7SJ09kAfDFkv5A0q9G6lA9Fn2hpF0CtSli8Sqk0AG5iglh4UouRQmejOQAAAPAa0wFekRyj+jg8BmeKGIP0cFalJByhh8A4AaVjaz12W/cAsANCpd4SJ88kQXAz5V0N0mHSbpK0g6Sjpf0DUknBWpTxOJVSKEDchUTwsKVXIoSPBnJAQAGgAHg5GbQn+Ez1o0i9hD9KdFEryTDDwBwgxJG1noAuIGwHT6kT57IAuDqDvBOwxr/QtLNh/88+kh09a+q3wte7VXE4lVIobs0VVi4kqtVgicjOQDAADAAnNwM+jN8xrpRxB6iPyWa6JVk+AEAblDCyFoPADcQtsOH9MkTWQD85GB93wMAB5XqRhgLV3KdSmg+kRwAYAAYAE5uBv0ZPmPdAIC7648MPwDADfwQWesB4AbCdviQPnkiC4DHLW8Ri1chhR5Xy0kez8KVrHYJnozkAAADwABwcjPoz/AZ60YRe4j+lGiiV5LhBwC4QQkjaz0A3EDYDh/SJ09kAfAd1qnv9wO1L2LxKqTQAbmKCWHhSi5FCZ6M5AAAA8AAcHIz6M/wGetGEXuI/pRooleS4QcAuEEJI2s9ANxA2A4f0idPZAHw9ZJukLQ4fvXPi68bBWpfxOJVSKEDchUTwsKVXIoSPBnJAQAGgAHg5GbQn+Ez1o0i9hD9KdFEryTDDwBwgxJG1noAuIGwHT6kT57IAuDbrajvzpKOkfReSf8UqH0Ri1chhQ7IVUwIC1dyKUrwZCQHABgABoCTm0F/hs9YN4rYQ/SnRBO9kgw/AMANShhZ6wHgBsJ2+JA+eSILgFcr7+9KOlvS3QO1L2LxKqTQAbmKCWHhSi5FCZ6M5AAAA8AAcHIz6M/wGetGEXuI/pRooleS4QcAuEEJI2s9ANxA2A4f0idPTBKAbyrpZ5KqhanulbZ4FVK8uuvv6t97vXBleccZ1zHGhg050zuSb0kA7GjW59hI3djMlOmAQmqXJU7GusEeIqta+eNm+KEYAHbmcr7U658haw9hXlev/WBqQfhAgVY8kbNDlh66okrbSXrK8LeB7xeoIItXQKQCQ1ox5SrXleYHR0Nn4XIWDmdcJ18nB2fcSL4AsKPoZGIjdQOAJ1ML9yyF1M5NOxqfsW6krRk9r0W0ZplxGX4AgBtULGsPYabSaz+YWhDeAQCuvgRr9HWVpK9Kerak/xuoIItXQKQCQ3rdqLI2Ps64Ts2zFq9IvgCwU6nJxEbqBgBPphbuWQqpnZt2ND5j3WAPEVW/vLgMPwDADeqctYcwU+m1H0wtCO8AAI9bJBavcRWczvG9blRZm1BnXKesWYtXJF8A2KnUZGIjdQOAJ1ML9yyF1M5NOxqfsW6wh4iqX15chh8A4AZ1ztpDmKn02g+mFoQXDsA3kfSfku4l6VcNq8Xi1VC4KR/W60aVtQl1xnXqm7V4RfIFgJ1KTSY2UjcAeDK1cM9SSO3ctKPxGesGe4io+uXFZfgBAG5Q56w9hJlKr/1gakF44QBcpXeppNtLurZhtVi8Ggo35cN63aiyNqHOuE59sxavSL4AsFOpycRG6gYAT6YW7lkKqZ2bdjQ+Y91gDxFVv7y4DD8AwA3qnLWHMFPptR9MLQjvAAC/XNKvJb2qYbVYvBoKN+XDet2osjahzrhOfbMWr0i+ALBTqcnERuoGAE+mFu5ZCqmdm3Y0PmPdYA8RVb+8uAw/AMAN6py1hzBT6bUfTC0I7wAA/5uk+0i6TNIPJY1+KdYDAxVk8QqIVGBIrxtV1ibUGdepedbiFckXAHYqNZnYSN0A4MnUwj1LIbVz047GZ6wb7CGi6pcXl+EHALhBnbP2EGYqvfaDqQXhHQDgY9ep0ssCFWTxCohUYEivG1XWJtQZ16l51uIVyRcAdio1mdhI3QDgydTCPUshtXPTjsZnrBvsIaLqlxeX4QcAuEGds/YQZiq99oOpBeGFA/CNJO0h6Tw+AzxzXu11o8rahDrjOo7KWrwi+QLATqUmExupGwA8mVq4Zymkdm7a0fiMdQMAjqpfXlyGHwDgBnXO2kOYqfTaD6YWhBcOwBskVb/7u72kGxpWi8WroXBTPjkSHKEAACAASURBVKzXjSprE+qM69Q3a/GK5AsAO5WaTGykbgDwZGrhnqWQ2rlpR+Mz1g32EFH1y4vL8AMA3KDOWXsIM5Ve+8HUgvDCAbhKr7r7+0hJP2pYLWvx6vnmoKGEUzmsiEaFH6ZS+7VOmuGJzvUHPLlkjww/FLO5LWrmdSeZDE9YPcKRKmsud21cRzMzNsMP9AizCAWF44eCilFIKq14orpbm/E6SNJTJVWf9/3Bii/B+n7ghNbilbVwBPIkZHMFWjHlKqLih+46LcMTnfMDPQoA7u4UTs986j3CucKsudy1cR3NzNgMPwDAZhEKCscPBRWjkFRa8UQWAI9+6/PiY9DVuap/rj4jXPfq3Aa37oJm5O+tmBIA7pVbMjzRuf6QtbntoFMy/MDmtoNGGEk5wxNWj3Dky5rLXRvX0cyMzfADPcIsQkHh+KGgYhSSSiueyALg260jUvWzSHUva/HKWjjqkuTvWyjQiikB4F45K8MTnesP9KglT2f4gc1tt1tGhiesHuHIlzWXuzauo5kZm+EHeoRZhILC8UNBxSgklVY8kQXA42pkLV5ZC8e4FzGDx7diSgC4V87J8ETn+gM9CgDu1axu92Km3iOcy8may10b19HMjM3wAwBsFqGgcPxQUDEKSaUVT2QC8D6SninptpIuknSypI8GxevcBjd4XX0Pa8WUAHCvbJLhic71h6zNbQedkuEHNrcdNMJIyhmesHqEI1/WXO7auI5mZmyGH+gRZhEKCscPBRWjkFRa8UQWAB8o6c2STpX0XUl3HH4p1nMlvScgoLV4ZS0cgTwJ2VyBVkwJAPfKVhme6Fx/oEcteTrDD2xuu90yMjxh9QhHvqy53LVxHc3M2Aw/0CPMIhQUjh8KKkYhqbTiiSwA/i9Jz5f06RGxHirpBEl3DQhoLV5ZC0cgT0IAYDywvgKtNKoVp+hcf6BHAcA0ijUVmHqPcGqTNZe7Nq6jmRmb4QcA2CxCQeH4oaBiFJJKK57IAuArJN18+K3Pi3ptJekXkqrE616d2+DWXdCM/L0VU66iFX7oroEyPNE5P2Rtbjtoiww/sLntoBFGUs7whNUjHPmy5nLXxnU0M2Mz/ECPMItQUDh+KKgYhaTSiieyAPgbkg6V9NkRsR4k6S2S7hIQ0Fq8shaOQJ6EbK5AK6YEgHtlqwxPdK4/0KOWPJ3hBza33W4ZGZ6weoQjX9Zc7tq4jmZmbIYf6BFmEQoKxw8FFaOQVFrxRBYAP2X4uPMpkr4naTdJB0l6gaR3BgS0Fq+shSOQJyEAMB5YX4FWGtWKU3SuP9CjAGAaxZoKTL1HOLXJmstdG9fRzIzN8AMAbBahoHD8UFAxCkmlFU9kAXCl0eMkPV3SrpIullTB8AeD4i1scOfm5rRxY/WP6782bMi8jLqz8/cRBVox5SqKWsBDRYpSIMMTlh+yNpaOyuQAADt+mbHYqfeIGdO79MvN8AMAXHrV184PPxRUuxL2MsOP0s6PK0up5AgAj1vZ6RxPo5qO7iWfNcMTAHCDivdp4eJNsgYGKPeQqfeIcqWZycwy/AAAd9dK+KGg2vVpH9E2AN9jeOf371ep1yuGd4C/HqglABwQqcAQGlWBRZlyShmeAIAbFLVPCxcA3MAA5R4y9R5RrjQzmVmGHwDg7loJPxRUuz7tI9oG4PdJ+tQav/X7t5IeLulJgVoCwAGRCgyhURVYlCmnlOEJALhBUfu0cAHADQxQ7iFT7xHlSjOTmWX4AQDurpXwQ0G169M+om0A/r6ku0ta7dnsHSRVd3/vEKglABwQqcAQGlWBRZlyShmeAIAbFLVPCxcA3MAA5R4y9R5RrjQzmVmGHwDg7loJPxRUuz7tI9oG4Ap81/vWqrq/L5YZAC7I8EYqNCpDrBkJzfAEANzAPH1auADgBgYo95Cp94hypZnJzDL8AAB310r4oaDa9Wkf0TYA/0jSfSRdtEq9bivpy5J2CtQSAA6IVGAIjarAokw5pQxPAMANitqnhQsAbmCAcg+Zeo8oV5qZzCzDDwBwd62EHwqqXZ/2EW0D8Psl/VjS4avU6/WSdpZ0QKCWAHBApAJDaFQFFmXKKWV4AgBuUNQ+LVwAcAMDlHvI1HtEudLMZGYZfgCAu2sl/FBQ7fq0j2gbgO8q6T8knSHpNEmXSNpFUvUFWPtJ+jNJ3wzUEgAOiFRgCI2qwKJMOaUMTwDADYrap4ULAG5ggHIPmXqPKFeamcwsww8AcHethB8Kql2f9hFtA3BVpgdLepukO0m6QVJ1jgslPUPS54J1BICDQhUWRqMqrCAFpJPhCQC4QWH7tHABwA0MUO4hU+8R5Uozk5ll+AEA7q6V8ENBtevTPiIDgBdLdUdJt5b0U0nfNeu3sMG9+OKLtXHjet+pNRh106ZqfvAqQAEaVQFFKCyFDE9YADw3NxeWJKuXkMNSCTL8wOY27PAiAzM8YfWIIlWZ3aQy/ECP6K6f8ENBtSthL1Nh3xq/NmQplQnAViIrgqvPClePT/PqlgLV4+6XJqSMHxJEndCQGZ7ADxMqXsJpMvxQpYknEoo1oSEzPIEfJlS8hNNk+IEekVCoCQ2JHyYkdIdO04onSgXgKq/bSLqyQwWZ9VSr33muvgW8euy97Rd+aFvRyYyX5Qn8MJn6tX2WLD9UeeKJtqs1mfGyPIEfJlO/ts+S5Qd6RNuVmsx4+GEyOnfpLK15olQA7lIxyBUFUAAFUAAFUAAFUAAFUAAFUKADCgDAHSgSKaIACqAACqAACqAACqAACqAACoyvAAA8voaMgAIogAIogAIogAIogAIogAIo0AEFAOAOFIkUUQAFUAAFUAAFUAAFUAAFUAAFxlcAAB5fQ0ZAARRAARRAARRAARRAARRAARTogAIAcAeKRIoogAIogAIogAIogAIogAIogALjKwAAj68hI6AACqAACqAACqAACqAACqAACnRAAQC4A0UiRRRAARRAARRAARRAARRAARRAgfEVAIDH15ARUAAFUAAFUAAFUAAFUAAFUAAFOqAAANyBIpEiCqAACqAACqAACqAACqAACqDA+AoAwONryAgogAIogAIogAIogAIogAIogAIdUAAA7kCRSBEFUAAFUAAFUAAFUAAFUAAFUGB8BQDg8TVkBBRAARRAARRAARRAARRAARRAgQ4oAAB3oEikiAIogAIogAIogAIogAIogAIoML4CAPD4GjICCqAACqAACqAACqAACqAACqBABxQAgDtQJFJEARRAARRAARRAARRAARRAARQYXwEAeHwNGQEFUAAFUAAFUAAFUAAFUAAFUKADCgDAHSgSKaIACqAACqAACqAACqAACqAACoyvAAA8voaMgAIogAIogAIogAIogAIogAIo0AEFAOAOFIkUUQAFUAAFUAAFUAAFUAAFUAAFxlcAAB5fQ0ZAARRAARRAARRAARRAARRAARTogAIAcAeKRIoogAIogAIogAIogAIogAIogALjKwAAj68hI6AACqAACqAACqAACqAACqAACnRAAQC4A0UiRRRAARRAARRAARRAARRAARRAgfEVAIDH15ARUAAFUAAFUAAFUAAFUAAFUAAFOqAAANyBIpEiCqAACqAACqAACqAACqAACqDA+AqUCsBVXreRdOX4l8gIE1JgB0k/knRDwvnwQ4KoExgyyxP4YQLFSzhFlh+qVPFEQsEmMGSWJ/DDBIqXcIosP9AjEoo1gSHxwwRE7tgpWvNEFgDfT9Ilkn4o6daSjpd0naQjJP08IPbOw+MDoYQUpMAuki5NyAc/JIg6oSEzPIEfJlS8hNNk+KFKE08kFGtCQ2Z4Aj9MqHgJp8nwAz0ioVATGhI/TEjoDp2mFU9kAfB/SXqcpO9IOlXS7SX9WtIvJD0xIPJGSXMXX3yxNm6s/pFXmwocdthh4eFOOOGE2tj5+XntuuuuVdwmSfO1B/gBaX7YtKlKmddLXvKSsAhHHFG9j7X+K9kTlh+ocV21tvz7kUceGT7oqKOOqo1N9kN1/gVP1CbSIGBuLmXYBplM95CDDz44nMDpp58ejk1aNyw/ODWmnwxK+7rXvS5c4xe96EXh2CQ/LPWI6L6SOjslk5w5v9dee9UOzppRK1HxAYccckg4x/e9733h2LZ6RBYAXyFpx+HVXCbpHpKuHgLx7waucmHxqhYlADiglhnyjGc8I3zEySefXBtbNarhYpEKwBl+2LAhawrUylZUwKte9apwPkcffXRtbLInrP5AjWvLtUXAS1/60vBBxx57bG1ssh+WNre1iTQIuOGGjE91NEhkyocceOCB4QxOO+20cGxbm5kVJ7QA2Kkx/WSg9EknnRSu8aGHHhqOTfLDUo+I7iOos1My6eMf/3j4gEc/+tG1sawZtRIVH3DQQQeFc3zXu94Vjm2rR2Tt/qvHnHeS9AeSPixp9+FntKq7g9Xz23Uva4NbNxh/31wBAHhZDxa5gRYAMF1iVAEAeFkNB4767CIAeFBd1gwAuM/zvOm1AcCsGSu9M6sA/InhZ3grCP6upBdIuoOkf5X0+4EJBgAHRGoaAgADwCu9AwA3nU39PA4AZjOz0tkAMAA86gnuAPez9ze9KgCYNQMAHihQfSD0OEnXDuG3+uzv4yXdTVL985PDz3NFH1VpOmFn9TgAGAAGgGd19seuGwBmMwMArz5XuAPMHeBYF52tKACYNQMAHihQfUNX9S3QK19r/fuVcdwBTuydADAADAAnTrAeDA0As5kBgAHg9VoZd4B70OhbvAQAmDUDAB4oUH3Wd7Wvb75c0i0Ccw4ADojUNAQABoAB4KazZzaOA4DZzADAADAAPBv9vo2rBIBZMwDggQJXrvJlV1tJ+h8AuI1WM94YADAADACPN4f6fjQAzGYGAAaAAeC+d/r2rg8AZs2YdQB+71CA6vO+K38IcLfh3+4XmHLcAQ6I1DQEAAaAAeCms2c2jgOA2cwAwAAwADwb/b6NqwSAWTNmHYAXf8jpbySN/qrx9ZKq3wN+u6QfBibbAgBXX6F9k5vcpDb8rW99a23MYgBfYCE5P+vxkY98pFbbX/7yl/qbv6lKrtTfAf7whz+s7bbbrjafhz/84bUxiwGOFl3zTta1RX6wvPLE05/+9CxPLPSHY445Rje72c1qax353eImfqg98UjA/e9//3D4ueeeG47NCnS8vscee9Sm8dvf/lbnnXdelh+qcYvwhKNbrWgTCMjqES972ctqs//Vr36l446rviszZd1Y8MOZZ54ZWjMe8pCH1ObbJMDRt8n40zzG8foHP/jB2lSrNWP4TeOp+4i3vOUt2mabbWrzeepTn1ob02TdcHQLJ9CxwGc961m1GV977bU69dRTs/rD0prx7ne/W9tuu21tPvvtt19tzGJA12pcQp86+OCDa/WtPPGe97ynNU9k/Q7wUcNvga69oDUCAOCmygWOc8wOAC8L2uem5lwbAByYZCtCAGAA2HdN/hHOWuD0CAB4UDtH3/xqt3sGxw8AcHf3Ee26ZjAaAJyhavMxS+hTfQLg0UpUn/0dfVV3g+teAHCdQmP83TE7ANzdhcups7OZAYD9yQcAA8C+a/KPyOoRADAAPOpeALi7+4iMLgQAZ6jafExnHWh+lvWP7BMA7y6pei75TyWtfNbkRgEBAeCASE1DHLMDwN1duJw6A8C5G1YAGABu2q8zj8vqEQBwbj/J9ER0bGfNAIC7u4+I+sGJA4AdtfJjnXUgK5s+AfAXJV0qqfpM8NUrBPtcQEAAOCBS0xDH7ABwdxcup87OZoY7wP7MA4ABYN81+Udk9QgAGADmDvDq89dZa/M7wHTOAABPR/e1zuqsA1mZ9wmAq98BvpWkaxuKBQA3FC5ymGN2ABgAXukpADgyyzaPAYABYN81+Uc4a4GzcQeAAWAAGABeq4MBwPm93TmDsw444zqxfQLg6us+/1LSTx0BRmIB4IbCRQ5zzA4AA8AAcGRWrR8DAAPA47uo/RGctQAA9vV39PVHn+4Rjh94BLq7+4gMlwHAGao2H7OEPtV1AL7DiPzVb348SVL1Wwg/XlGW7wfKBAAHRGoa4pgdAO7uwuXU2dnMcAfYn3kAMADsuyb/iKwewR3gQe0cffOr3e4ZnDUDAO7uPqJd1wxGA4AzVG0+Zgl9qusAXH278w3DEoz+vNLov6v+mS/Bau7TVo50zA4Ad3fhcursbGYAYH8aAsAAsO+a/COyegQADACPuhcA7u4+IqMLAcAZqjYf01kHmp9l/SO7DsC3Cwrzw0Acd4ADIjUNccwOAHd34XLqDADnblgBYAC4ab/OPC6rRwDAuf0k0xPRsZ01AwDu7j4i6gcnDgB21MqPddaBrGy6DsBt6gIAt6nmirEcswPA3V24nDo7mxnuAPuTEwAGgH3X5B+R1SMAYACYO8Crz19nrc3vANM5AwA8Hd3XOquzDmRl3icAPnANkX4lqboD/BVJ1SPTa70WAHhubk4bN1b/yKtNBfbZZ5/wcB/72MdqY+fn57Vp06Yqrvqv6hvA234t+CE6aNZkdhYuJwdn3KgGVZyTwze/+c3w0He9613DsUmesPzgJOto5ozr1NjJIWvchz3sYeHLO+ecc8KxSX6ozm+tGY5uzsU5tXPGzYp1dPjkJz8ZTmOvvfaqjU1eN6we4dTN0cwZt1awwgLe8pa3hDN69rOfHY7N7hFOItHYrDpHNT7kkEOiqabFPe1pTwuPfeqpp4Zjs/0Q5Qxn3jsX17VYx+vHH398+PJe8IIX1Ma2vWaMfla39uRGwMWSbj38vO8vJN1c0m8lVf/8O5IukPQYSd9ZY0xrM2PkRagkALiZDZwG6DQJZ1wncycHAHigrKOZUwunxk4OWeMCwE51l2Od2jU7Q7tHOf4BgAfaO5p1zQ+Ou6JwVo0JADvKLsdGNQaAG+lrcYYz7xtl05GDnJ42qwD8XEnVbaLDJV05fHf+dZK+Jek0SW8d/rtHAcCTdz0A3ExzpwE6TcIZ18ncyQEABoBHvQUAOzMNAF6pFneAc/tJM3e2e1QUzgDg5rpHNQaAG2kMADeQzdlXzioAXyTpTpKqR54XX9sO7/zuKumWw3++FQDcwIFjHgIANxPQAVWnSTjjOpk7OQDAuRtWp8ZO3bLGBYCdmQYAA8Cr+8WZy80cN72jonAGADevUVRjALiRxgBwA9mcnjarAPw/ku4i6Scj+u40vANcwe9Ww8ehFz44usrLMmaDGs70IQBws/JnwYYzrpO506gAYAB41FsAsDPTAGAAGABeb8bwCHSzfgIAL+iW+t0yfAbY86azr5xVAH63pN0lHSupuhtc/UTSMZIulPQUSQ+UdJKkuwHAnvnaiAaAm6nogKrTJJxxncydHABgABgAdmZXP4DH6T18BnhQc0czpweP777JjhCFsyorALhZbaIacwe4kb7WjTZn3jfKpiMHOT1tVgF4e0knSnqipJtIulbSP0mqPht8laTqMejq338PAJ686wHgZpo7DdBpEs64TuZODgAwAAwAO7MLAF5PLT4DnNtPxnfq+CNE4QwAbq51VGMAuJHGAHAD2Zx95awC8KKs1aPO1ed8q0ei1/vZo5VlsIzZoIYzfQgA3Kz8Dqg6TcIZ18ncyQEAzt2wOjV26pY1Lo9AOzNtOdapXbMztHuU4x/uAA+0dzTrmh8cd0XhDAB2VN08NqoxANxIY4sznHnfKJuOHOT0tFkH4KYltYzZ9CSzehwA3KzyTgN0moQzrpO5kwMADACPegsAdmYaALxSLe4A5/aTZu5s96gonAHAzXWPagwAN9LY4oysfVqjzKd4kLOvnCUA/r+S/nBYl+p3gG9Yo0a3DdTOMmZgPEJGFACAm9nBaYBOk3DGdTJ3cgCAczesTo2dumWNCwA7Mw0ABoBX94szl5s5bnpHReEMAG5eo6jGAHAjjS3OcNbaRtl05CCnp80SAFef960+51u9nrxOLd8TqLNlzMB4hADAY3vAaYBOk3DGdS7CyQEABoBHvQUAOzMNAAaAAeD1ZgxfgtWsnwDAC7rxLdDN7JNylLOvnCUAblPsNAB2YMMpdPTis87vjPv5z38+mq4e8IAH1MbOz89r06aFX7RKbVS1iQwDnLo5umWNG70uN87J98wzzwwPv/fee4djkzyx0B+cJDJiHX1L8JmT7x577BGW7D//8z/DsUl+qM5vrRldq4eTr1MMxxP7779/eOgzzjijNjZ53bB6hKNDVi1qBRsJyMrXGfdJT3pSOOX3ve994djsHuEkkhHraBz12gUXXBBO9U53ulM4Nnr+asD99tsvPO4HPvCBcGy2HzJ+BimjxpVgJYzrFO7LX/5yOPze9753bWzba8aG2jM2D7i9pAMk3UbScyTdUdLWkqpHpete1mambrDRvzsT2jFbNIes8zvjAsDL1XJ0c/zgjBv1jhvn5AsAu+qWsRg5PnP8AAAv+8HRLasezriOk51rA4AHymbVIqtuTr6OHwBgp2K5/QQAblQLizOy5lHXxnWUnlUAfqikj0n6oqQ/l7SDpOpW4pGSHhUQ0DJmYLylkCyzRXPIOr8zLgAMAK/0KwAcncG5G5lqdGcT6sx7Z1wAuLt1dpzseAIABoBHvQUAOzMtt58AwI1qYXFG1lrbtXEdpWcVgL8i6ZWSPirpF5JuLmkbSd+X9HsBAS1jBsYDgEdEAoABYADY6Rqrxzrw0LVFDgDO3bBmvtHhONvxMAAMAAPAzuya3LoBADeqi8UZXVvDs/J1lJ5VAL5C0o5DoS6XdIvhPy/CcJ2GljHrBhv9+7RNkXV+Z1wAGAAGgJ2uMbmNTCYYObADAAPAK10PAAPAAHCZ6wYA3KguFmc4e2xnre3auI7SswrA35a0j6TvSloE4N0lVZ98v2tAQMuYgfGWQrLMFs0h6/zOuAAwAAwAR2fs2nF9XuQAYAAYAF597jtr7fhdptw333gEull1M9YNALhRLSzOcOZ9Ro1LeXPcUXpWAfgFkv5G0ouH0PtoSa+T9H5JJwUEtIwZGA8AHhEJAAaAAWCna5S7Cc1alAFgABgABoDX65IAcLM1JAOOAOBGtbA4I2ut7dq4jtKzCsBbSTpG0vOGP09xjaS3STqi+o6XgICWMQPjAcAA8Ko2KaH5OP51Yp2Fli/BcpQdxDr6luAzJ18AGAAGgAFgANhfF+qOcPpwdN0AgOtUX/XvFmdEa9H3vYGj9KwC8KhGt5ZUfSb4WkM4y5jGuNZPGDiNKppDCZOIO8DL1SqhHlHvuHGOfwFgV10AeFExfge4u/3E6RF8BnhQZ2fN8LtK7Ainbk6+zrjcAY7VamWUo3G0dgBwo1pYnBGtBQC8XItZBeDq8edPSfpZI1tKljGdc2SZOJpD1vmdcQHg7m5Yoz5zmzAA7Cg7iM3YyJQyLneAl/3QtTo7TnauDQAGgEe9BQA7My23nwDAjWphcYazx3b6atfGdZSeVQCuvgTrTpK+MQThf5H0eeMusGVMpyBZZovmkHV+Z1wAGABe6VcAODqDczcyALBfh+ER1prh9MuubWYcBZ1rA4ABYADYmV2rxzpzLtqnAOBGdWHNGMrmeNJRelYBuNJoF0l/KWnP4X+2G0LwIwICWsYMjLcUEm0o7kY0moNz/uiYbtzxxx8fPuTwww+vjZ2fn9emTZuquOq/5msP8AMW/BA9LGsyR89fxZVQZ0eHQw89NHx5b3rTm8KxSZ6w/OAkW0KsU7esfN/xjneEhz744IPDsUl+qM5vrRnO/HTq4YzriObEOvk6437gA9WPOMRe++23X21g8rph9YgszWpFGAlwvFNCvuedd1748u5xj3uEY7N7hJNINNaph1PnjPNHx3T3MQ984APDQzs3YbL9MDc3p40bq3ax/iujbnXnXPn3afusysfJ4e1vf3v4EiP7iLbXjA3h7JoFbi3pAZIq6H1axQUjvw+83ojWZsZJzTGxU+hoDs75o2O6cQCwq5gfX0KdHf8CwH6Ns45w6paVAwC8rKxTj67Ne8c/ALCjlh/reMfxpJ9J7AgAuJwekeUHx5MAcGzejBPl1NmpnZOTk8OsAnB1O6m6+1u9JXT+8DHo6jPB/y7puoDYAHBApKYhAHBT5eLHZTWfeAbeO3UAsKNsbqyzwGRlAgCXs7l1a5zlHwDYrYQX76wZWTV2MgaAy+kRWX5wPAkAO7OnWaxTZ6d2TjZODrMKwNdLulDS0ZI+FoTe0RoAwI4jzVgA2BSsQXhW83FScRoVAOwomxvr1C0rEwC4nM2tW+Ms/wDAbiW8eGfNyKqxkzEAXE6PyPKD40kA2Jk9zWKdOju1c7JxcphVAH7Q8A5wdRf4jsPP/lZ3gKv/fCcgNgAcEKlpCADcVLn4cVnNJ54Bd4AdrUqKdRaYrLwB4HI2t26Ns/wDALuV8OKdNSOrxk7GAHA5PSLLD44nAWBn9jSLders1M7JxslhVgF4VM9bSHq6pCOHX1Ryo4DYAHBApKYhAHBT5eLHZTWfeAYAsKNVSbHOApOVNwBczubWrXGWfwBgtxJevLNmZNXYyRgALqdHZPnB8SQA7MyeZrFOnZ3aOdk4OcwqAG8r6SEjd4FHfxLpRQGxAeCASE1DAOCmysWPy2o+8QwAYEerkmKdBSYrbwC4nM2tW+Ms/wDAbiW8eGfNyKqxkzEAXE6PyPKD40kA2Jk9zWKdOju1c7JxcphVAP61pP+RdM7wsefqf39iiAwAG2K5oQCwq5gfn9V8nEycRsVngB1lc2OdumVlAgCXs7l1a5zlHwDYrYQX76wZWTV2MgaAy+kRWX5wPAkAO7OnWaxTZ6d2TjZODrMKwH8s6ZuOqCtiAeAxxKs7FACuU2j8v2c1Hyczp1EBwI6yubFO3bIyAYDL2dy6Nc7yDwDsVsKLd9aMrBo7GQPA5fSILD84ngSAndnTLNaps1M7Jxsnh1kFYEfP1WIB4HEVXOd4ADhR3OHQWc3HydxpVACwo2xurFO3rEwABWZlMgAAIABJREFU4HI2t26Ns/wDALuV8OKdNSOrxk7GAHA5PSLLD44nAWBn9jSLders1M7JxskBAHaUXY4FgJvpFjoKAA7JNFZQVvNxknIaFQDsKJsb69QtKxMAuJzNrVvjLP8AwG4lvHhnzciqsZMxAFxOj8jyg+NJANiZPc1inTo7tXOycXIAgB1lJwDAzdKZzlGvfOUrwyd+yUteEo59wxveEI593vOeVxs7Pz+vTZs2VXHVf83XHuAHLLwhEj3MmaDRMTPjnEaVdW2HHXZY+BJPPPHEcGySJyw/OMmWEJtVY+fa7na3u4XD/+u//iscm+SH6vzWm6YlzDlHtBLyPfXUU8Mp/93f/V1tbPK6YfWIEuZcVo2zxr3wwgtra7wYcOc73zkcm90jnESisSX4J5qrE+d459WvfnV46KOPPjocm+2Hubk5bdxYtYv1X44WdWM1/bvjs6x8nRzOPPPM8KU++tGPro1te83YUHvG6QRYm5nppJh/VgB4SePObWYcdziNymk+Tg4AsKNWbmxWjZ2sAeBltUqoRwk9AgB2ZpAfm1XjrHEB4LJ6hO+4+iMc7wDA9XqOG+GsRU7tnLycHGYRgG8i6T8l3UvSrxxhR2IBYEkAMAC8cv44zceZewCwo1ZubFaNnawB4LI2t85mJss/ALAzg/zYrBpnjQsAl9UjfMfVH+F4BwCu13PcCKe3O7Vz8nJymEUArrS8VNLtJV3rCAsAb64WAAwAA8CrdhDriYCGPWhqhzkLTFaSAHBZm1tnM5PlHwA4a7YNxs2qcda4AHBZPSLDnY53AOCMCmw+ptPbndo5mTs5zCoAv1xS9VvAr3KEBYAB4DX8YgGPM0Eb+rPVw5xGlXVt3AFutaRjDZZVYycpALiszW0JPQIAdmaQH5tV46xxAeCyeoTvuPojHO8AwPV6jhvh7A2c2jl5OTnMKgD/m6Q/k/QTST+UdP2IwA8MiM0j0DwCPWoTAHiohtN8AvNsKQQAdtTKjc2qsZM1AFzW5tbZzGT5BwB2ZpAfm1XjrHEB4LJ6hO+4+iMc7wDA9XqOG+H0dqd2Tl5ODrMKwMeuI+jLAmIDwAAwALzKRHGaT2CeAcCOSBOKzaqxkz4AXNbm1tnMZPkHAHZmkB+bVeOscQHgsnqE77j6IxzvAMD1eo4b4fR2p3ZOXk4OswrAjp6rxQLAADAADACv1UesJwLGbUaTPt5ZYLJyA4DL2tw6m5ks/wDAWbNtMG5WjbPGBYDL6hEZ7nS8AwBnVGDzMZ3e7tTOydzJYZYBuNqkVj/stIuk10n63arHS7osIDYADAADwAAwABxolhkhAHBZm1tnM+NsUBzvAMCOWn5sVo2zxgWAy+oRvuPqj3C8AwDX6zluhNPbndo5eTk5zCoA313S2ZLmJN1G0vaS9pL0NEl/HRAbAAaAAWAAGAAONMuMEAC4rM2ts5lxNiiOdwBgRy0/NqvGWeMCwGX1CN9x9Uc43gGA6/UcN8Lp7U7tnLycHGYVgD8r6QxJb5H0C0k3l7SDpP87vCNcpzcADAADwAAwAFzXKZP+DgCXtbl1NjPOBsWxDwDsqOXHZtU4a1wAuKwe4Tuu/gjHOwBwvZ7jRji93amdk5eTw6wC8M8l/c7w258vl3SLocDVHeFNAbF7C8COKR2jZY0bqJXm5+e1adNCWav/mo8cY8ZYn/nsim6LGji1c3RzdDjqqKPCQ7/mNa8JxyZ5wvKDk6wT6+jr1LiEcffaq3pgJ/Y6++zqYZ/wK7VHzM3NaePGyh7rv0qoR12Ok/i747X3v//94ZQOOOCA2tjkdcPqEY4OtRc2EpDlsz7nkLRmVJJZnnA0nrZ/Lr744nC6u+66azjWCXzOc54TDn/Tm94Ujs32Q8aa4VxcCbFZ/m37TZG214zqM7kZr+9KeoCkH0taBODbSvoXSXcOnBAAluSYcpoLbdumXMUf1sLVFd0Wr9OpXWDuLIU4OgDAjrKDWEdfp8YljAsAL/shqx6+49o/wrk2ANjXP2veO5l0LYds4HG0i8Y68yg6ZhUXrR0A7Ki6FGtxRrQWjTKZ8kFZ/p1VAH61pHtJerak/5D0J5LeLOlrkl4eqLVlzMB4xYQ4k8gxZda4EeEA4IhKa8c4tXPO5PgHAHaUBYBH1eIO8LIaWXPZd2fsCKdHAMAxTUejHD84tXAy6VoOALDfTwBgZ0YAwCvVyuo9swrAN5V0sqQDh0LfIOkjkp4o6dcBqwLAhdxhCtSKR6AjIq0T42xQnFM5TQ0AdpQFgAHg1f2SNZd9d8aOcHoEABzTFADeUidzXqR+TMKvYv0RzjyqHw0AXqFRqh94BNp7ms3x76wC8KJGt5S02/Cnjy4yhAOAAeBRu/AItDF5FkOdRRkA9gV29HU2gCWMyyPQy37IqofvuPaPcK4NAPb1z5r3TiZdy4E7wAAwAOzM8PFjnXXAOdusA7Cj1RbAE31npulJpnFc1mKUNW5EIx6Bjqi0doxTO+dMTlMDgB1lB7GOvk6NSxgXAAaAV84IANjvEVnz3smkazkAwAAwAOzM8PFjnT2Hc7ZZAuDPVHvCgDgPDcRwB7iQDXagVjwCHRFpnRhng+KcymlqALCjLAA8qhafAfY3rL7bco5wegQA7NfA6e1OLZxMupYDAOz3Ez4D7MyIpViLM5x51CibKR6U1XtmCYBfMVK/HSX9naSzJP23pN+X9EhJp0p6bqDOljED4xUT4kwix5RZ40aE4w5wRKW1Y5zaOWdy/AMAO8oCwADw6n7Jmsu+O2NHOD0CAI5pOhrl+MGphZNJ13IAgAHgFf7mM8DOhG8Qm9V7ZgmAR2X/oKR3Szpz5F8+StJTJf11oD4AMHeAR23CZ4ADk2ZliNPUAGBfYEffrE1o1rg8Ar3sh6w6+45r/wjn2gBgX/+s+elk0rUcAGAAGAB2Zvj4sc464JxtVgF4TtLNJV0/ItZWkn4xbG51GgLAADAAXDdLav7uNDUA2Bfb0TdrE5o1LgAMAK+cEQCw3yOy5qeTSddyAIABYADYmeHjxzp7GedsswrAF0p63vAR6EW9HiHpREl/EBAQAAaAAeDARFkvxGlqALAvtqNv1iY0a1wAGAAGgP2esPKIrPnpZNa1HABgABgAdmb4+LHOXsY526wC8N8OP+/70ZHPAO8j6WBJ7w0IWAQAOwtH4JpSQxwDf+ADHwjnst9++9XGzspngGuFKCzA8YSTujkvMj6/Yz0S71ybo5mpQziNrByccQ8//PBwvm94wxvCsdmb2+gvB5RQO0e05z2vej859nrjG98YC5Tk6HDooYeGxz3xxOq97vVfyeuG1SOcueFoVqdBaX8/66zqa1tir2c+85mxQEkXXeT8CqYy1owqV8sT4YtLvFnh5JAR68yLo48+OpzCcccdF45lzXCkKiv2kY+svvop9or0nrbXjA2x1BpF3VfSUyTtIukSSe+R9IXgSABwUKjFMKdRAcDL4vZ5M+N4wrGbqVnGZqbXGxmnbk4tnHEBYGdGLMc6GjtnAIAdtRZirR7h1M2Zc3bWUz4gsgldTBEAbjbvu+QfZ14AwM0mr6NxszNM76hZBuBxVAeATfWcSQQAA8CmvTYLNxdwANgU25nLTi2ccQFgs2jDcEdj5wwAsKMWAGyrNTwAAG6mnDPvnZ7dLJv2jnKuCwBuprujcbMzTO+oWQbguw8feb5t9QSMpHdI+j/BUgDAQaEWw5xJBAADwKa9AOBVBMvayDhz2cnBGRcAbjZDHI2dMwDAjloAsK0WANxUsoXjnHnv9OyxkmrhYOe6AOBmgjsaNzvD9I6aVQDeS1L1+d9PSvqupN0kVV+C9djhv6urCABcp9CKvzuTCAAGgE17AcAA8JqW4TPAy9I4fdiZgwCwoxYAbKsFADeVDAAeKgcAN7NQ1prRLJt2j5pVAP6ypOMlnT4i5/6SXijpTwMSA8ABkUZDnEkEAAPApr0AYAAYAA5MGqcPB4ZbCgGAHbUAYFstALipZAAwADwx74x1oikcPKsAXP3e7y1X/A7wjST9XNKOgToAwAGRAODVRXI2oV16HMm0hPVYljO2qRmfAXbETXyczpkXPAJtFm0Y7mjsnAEAdtQCgG21AOCmkgHAAPDEvDPWiaZw8KwC8AWSqju+Xx/R/G6SPtil3wE2N/pTsNfyKZ2NF3eAl3XrUo1dgzmecMY2NQOAHXEBYFOtpXDrTVPTw+GcsuYcABwuwWIg3wJtSybxJVgNREvs2c2yae8op5/xCHQz3R2Nm51hekfNKgBXjzpXPxr4OknfH34G+AhJbx7+u7qKWJuZusGa/j1rk9Q0n/WOcyYRAAwAj+NBc14AwKbYzlx2auGMyx1gs2jDcEdj5wwAsKPWQiwAbEsGADeQbOEQZ947PbtpPm0d51wXANxMdUfjZmeY3lGzCsBbSXq+pKdL2lXSxZJOkfQPKx6LXqsyALDpWWcSAcAAsGmvzcLNBRwANsV25rJTC2dcANgsGgC8mWAnnnhirYDz8/PatKlqD5p6j3DmhjPnakUoLIA7wM0K0lf/ONcFAOd7p9kZpnfUrALwuIoDwKaCTqMCgAFg014A8CqCZW2Enbns5OCMCwA3myGOxs4ZuAPsqLUQyx1gWzLuADeQbOEQZ947PbtpPm0d51wXANxMdUfjZmeY3lGzCsAPG971PX9E+t0l7SLpnEA5Fhavf//3f9f2229fG/7Hf/zHtTF9D3CaqjPhTjrppFrprrnmGr34xS+u4jLeyV/azOy5557aeuuta/P55CerX9/i5dTZ8c/97ne/WnGvu+46/cd//EeWJxb6QwUFN73pTWtzee1rX1sb0yQgS19n3PPOOy+c+j3ucY9wrBP4sIdV7X79V+WHz3zmM1l+WOoRc3Nz2rixssf6L8fvTj3qzjv69xJyyMr34x//eO3Qv/zlL/X4xz8+yxMLPWL33XfXjW5Uff/m+q9vfetbdSFLf8/yQziBQgId/15wQfW1MOu/rrrqKu2xxx5ZfljqEc961rNC68Yb3/jGupTTPfHNb34zlEMJe2DHD7vuWj0Yuv7r+uuv16WXXpruh7e//e3aZptt6tLRk570pNqYxQBHi/CgHQx0euV3vvOd2iusesQ973nP1jyxofaMzQK+LemRkn4wcvjvS/qEpLsEhgSAAyKNhjgTzjElAGwWoqBwp86OfwDgQZGz9HXGBYCXJpz11JDjd6cezvQvIYesfAFgR9luxjr+BYCb1RgABoCbOaeMo5y1s08APD98BGllFeaGdwnrqgMA1ym04u/OYuSYEgA2C1FQuFNnxz8AMAA8anPuADeb9M6cc+Zys2zqj3LyBYDr9ex6hOMHALhZtQFgALiZc8o4ylm3+gTA1Tc/P0DSgnuHr50l/buk2wZKAwAHRBoNcRYjx5QAsFmIgsKdOjv+AYABYAB4/InuzDlnLo+f2eojOPkCwFlVKGdcxw8AcLO6AcAAcDPnlHGUs271CYDfMgTdAyVdLukWkt4l6ceSnhkoDQAcEAkAXl0kPgOcC2gAcK6+zqLBI9BLPYBHoM01ww13gAcAdtXtXrzjBwC4WX0BYAC4mXPKOMrZy/QJgG8u6aOS7ivpF5Kq//9FSY8Z/v+66gDAdQqt+LuzGDmm5A6wWYiCwp06O/4BgAFg7gCPP9GdOefM5fEz4w7wehqWUIusGjvjOv4FgB1ll2MBYAC4mXPKOMrplX0C4EX17yWp+vKr6suwvmqUBAA2xKpCncXIMSUAbBaioHCnzo5/AGAAGAAef6I7c86Zy+NnBgADwPUucvwLANfruVoEAAwAN3NOGUc561YfAbhpFQBgUzlnMXJMCQCbhSgo3Kmz4x8AGAAGgMef6M6cc+by+JkBwABwvYsc/wLA9XoCwJsrwM8gNfNMSUc561afALj68d4XSLq3pB1WFOSBgQIBwAGRRkOcxcgxJQBsFqKgcKfOjn8AYAAYAB5/ojtzzpnL42cGAAPA9S5y/AsA1+sJAAPAzVxS7lHOutUnAD5D0t0lfUTS1SvK87JAuQDggEgA8Ooi8SVYuYAGAOfq6ywafAnWUg/gS7DMNcMNd4CHL8Fy1e1evOMHALhZfXkEmkegmzmnjKOcvUyfALj64qvdJf2kYRkAYFM4ZzFyTMkdYLMQBYU7dXb8AwADwNwBHn+iO3POmcvjZ8YdYO4A17vI8S8AXK8nd4C5A9zMJeUe5axbfQLg/5Z0Z0nXNiyN9W5+w3P06rCrrroqfD1HHnlkOPZNb3pTbez8/Lw2bdpUxVX/NV97gB+Q5gdnEffTnu4RTvM5+uijw8m++tWvro1N9oTlhz7XuLYQDQNe9rLIgzqDwY855pjasyT7oTr/gidqE2kQ4MyjBsN35pC99947nOuZZ54Zjk1aNyw/ODWmnwxKe/bZZ4dr/PCHPzwcm+SHpR4xNzenjRsre6z/os51Cm3+94MPPjh8wMknn1wby5pRK1HxAaeffno4xyc84Qnh2LZ6xAbnjEZsNRMqAH6hpOuN4xZDrQ1ug/F7dwgA3KykfV7knE0dANzMP309CgBerqwzj/rqh+q6AOBBdfu8Zjj+BYAdtfofCwCzZqx0+awC8MWSdpL0G0n/s0KU2wZaAQAcEGk0BAA2BRuG93kz42zcAeBm/unrUQAwm5mV3gaAAeBRTwDAfe3+za4LAGbNAIAHCjx5nSn0nsD0AoADIgHApkirhAPAA1EA4PG91KcRAGA2MwDw6jO6z2uG08MAYEet/scCwKwZAHA78xwANnXkDrAp2DC8z5sZ7gBzx6bZrJAAYDYzADAAvF7/AICbdtd+HgcAs2bMOgDvJ+kDIyLsJul7I///CEmvD0x/ADgg0mgIAGwKBgBvJhh3gJv5p69HAcBsZgBgABgA7muHb/+6AGDWjFkH4OobgEe/Xu9ySbcYEWXl39eahQCw2Z8AYFMwABgAbmaZmTgKAGYzAwADwADwTLT7Vi4SAGbNmHUAvlLSDiMiVL8HfPOR/7/y7wBwK61HAoCbCckj0APduAPczD99PQoAZjMDAAPAAHBfO3z71wUAs2bMOgBzB7j9vhIaEQAOybRFEAAMADdzTr+PAoDZzADAADAA3O8+3+bVAcCsGQAwj0C32VPCYwHAYak2CwSAAeBmzun3UQAwmxkAGAAGgPvd59u8OgCYNWPWAfjXkl47IsILJf3/I///RZJuFph0C58Bvvjii7Vx4+hHigNHzmiIA8AvfelLwyq9/vX131k2Pz+vXXfdtRpzk6TqKYC2X2l+2LSpSrmfr7m5ufCFOcBz7LHH1o6b7AnLD32ucW0hGgYcddRR4SOPPPLI2thkP1TnX/BEbSINApx51GD4zhyy//77h3N1viE4ad2w/ODUmH4ysMGHP/zhsB/23XffcGySH5Z6RHRfSZ2dkklPecpTwgeccMIJtbGsGbUSFR/woQ99KJzjQQcdFI5tq0dscM4YiP2spBtq4h4SGGdnSZcE4ggpS4FdJF2akBJ+SBB1QkNmeAI/TKh4CafJ8EOVJp5IKNaEhszwBH6YUPESTpPhB3pEQqEmNCR+mJDQHTpNK55oG4Db0q/K6zaSqi/N4tUNBaovP/tR4A2QJleDH5qoNv1jsjyBH6Zf2yYZZPmhygVPNKnI9I/J8gR+mH5tm2SQ5Qd6RJNqTP8Y/DD9GpSWQWueKBWASxOcfFAABVAABVAABVAABVAABVAABTquAADc8QKSPgqgAAqgAAqgAAqgAAqgAAqgQEwBADimE1EogAIogAIogAIogAIogAIogAIdVwAA7ngBSR8FUAAFUAAFUAAFUAAFUAAFUCCmAAAc04koFEABFEABFEABFEABFEABFECBjisAAHe8gKSPAiiAAiiAAiiAAiiAAiiAAigQUwAAjulEFAqgAAqgAAqgAAqgAAqgAAqgQMcVAIA7XkDSRwEUQAEUQAEUQAEUQAEUQAEUiCkAAMd0IgoFUAAFUAAFUAAFUAAFUAAFUKDjCgDAHS8g6aMACqAACqAACqAACqAACqAACsQUAIBjOhGFAiiAAiiAAiiAAiiAAiiAAijQcQUA4I4XkPRRAAVQAAVQAAVQAAVQAAVQAAViCgDAMZ2IQgEUQAEUQAEUQAEUQAEUQAEU6LgCAHDHC0j6KIACKIACKIACKIACKIACKIACMQUA4JhORKEACqAACqAACqAACqAACqAACnRcAQC44wUkfRRAARRAARRAARRAARRAARRAgZgCAHBMJ6JQAAVQAAVQAAVQAAVQAAVQAAU6rgAA3PECkj4KoAAKoAAKoAAKoAAKoAAKoEBMAQA4phNRKIACKIACKIACKIACKIACKIACHVcAAO54AUkfBVAABVAABVAABVAABVAABVAgpgAAHNOJKBRAARRAARRAARRAARRAARRAgY4rAAB3vICkjwIogAIogAIogAIogAIogAIoEFMAAI7pRBQKoAAKoAAKoAAKoAAKoAAKoEDHFQCAO15A0kcBFEABFEABFEABFEABFEABFIgpAADHdCIKBVAABVAABVAABVAABVAABVCg4woAwB0vIOmjAAqgAAqgAAqgAAqgAAqgAArEFJgUAP+FpOskfS6Wlqq8biPpymA8YdNXYAdJP5J0Q0Iq+CFB1AkMmeUJ/DCB4iWcIssPVap4IqFgExgyyxP4YQLFSzhFlh/oEQnFmsCQ+GECInfsFK15IguAPyXp1ZI+K+m5kl4r6XpJR0s6ISD2zpIuCcQRUpYCu0i6NCEl/JAg6oSGzPAEfphQ8RJOk+GHKk08kVCsCQ2Z4Qn8MKHiJZwmww/0iIRCTWhI/DAhoTt0mlY8kQXAPxtuSK6V9PUhBM9JOkPSnQIib5Q0d/HFF2vjxuof13995StfqQtZ+vv3vve9cOw222wTjn3MYx4Tip2fnw/FVUGXXXZZOPZOd4rIOhju4Q9/eHjcd77znbWxV111le5973tXcZskxS+wduSlgAU/3O52t9NWW21Ve9R5551XG7MYsGlTlXL7r3e9613hQZ/61KeGY29xi1uEY0877bRw7Bve8IZw7IUXXlgbe/311+uSSxbew8rwxIIfzj77bG233Xa1udz//vevjSkp4MY3vnE4HadHXXll/IGaQw45JJxDRN9f/vKXetrTnpblh2rcBU9Ekz7ggAOioXrb294WjnUCzznnnHD4nnvuGY51Aj/xiU+Ew7/73e+GY4855phwbGaPiCbxta99LRqqPfbYIxzrBD7jGc8Ih5988snhWCfwAQ94QDj8Zz+rtnmx1/nnnx8LHERlrBlLPeJDH/pQaN3Ya6+9nJzDsbe+9a3DsT/96U/DsdMOvN/97hdO4eUvf3lt7NVXX6199tkn3Q+HH364bnrTm9bmc+aZZ9bGLAaccsop4Vgn8GY3u1k4/Ac/+EE4dqeddgrHDvd1ofhDDz00FFcF7b333rWx1157rYZ761Z6RBYAXyFpR0m3lHSBpFsNr6yCo3qiHW5m5ubmQgD8pS99qVa4xYDvfOc74dhtt902HPu4xz0uFOsA8I9+VD1RHHvtvvvusUBJziL3/ve/v3bcamP9R3/0R+mN6va3v30IgJ2N2oYNOVPg9NNPr9VtMeDxj398OPZWt1qcSvWHfPCDH6wPGka85jWvCcdGNjMVAF900UVZnliAnXPPPVfbb799bd53v/vda2NKCth6663D6TgA7PSe5z//+eEcHvSgB9XGVgD8xCc+McsP1bgWAD/5yU+uzXkx4N3vfnc41gms3sCJvpw3LaNjVnEf/ehHw+GRN74WB3vRi14UHjcJeCw/XHBBtU2Jve585zvHAs2o5z63elgu9jrxxBNjgWbUQx7ykPARDpx961vfCo+b5IelHvHP//zPIQB29knOxTmw4dwEcXLIiI2sA4vnfe1rq4dC139VAPwXf1F9ejL3DZGjjjpKEbD83//7f9elvPT3973vfeFYJ9BZ77///e+Hh/693/u9cGx1YzL6Gr7pHQp/7GMfWxtXAfDwzb+iAfibkp4tqaKih0nad2ji6vZrZAe/sHgBwADwcEYs+AEAlgDgBUcAwMOJ4SyIAPDy+goAD7QAgAc6AMADHQDg5R4BANfyyGYBAPCyHADwQItZBeD9JVXPX1aPQFfPkXxB0l9LeqakyLNcAHD1jVLcAV7sKADwUAkAGAAe3XUAwEtqWHf8AGAAeHQeAcAA8ErcA4AB4FEFuAM8UIM7wLF5sfgB2muG4b87/N+fBA4HgAHgUZsAwADwFn7gEWgJAAaAA+vpmiHcAR5IAwADwADwOJ1E4g7wsn7cAR5oMat3gKtv1njUKtPpY5IWPtVe8wKAAWAAeJVJwh3gBVF4BHroDQAYAK5bTNf7OwAMAI/6g0egl9XgDrDXWQBgAHilY2YVgNf6sqvLJUW+xhYABoABYAB4rRUYAAaAV3qDR6C9/epCNAAMAAPAq08cANhrKAAwADzrAPzQoQAfl/RoSaNfsVt9deKRkm4XmFYAMAAMAAPAAHBNs+QOMHeAA+vpmiEAMAAMAAPA4/SQxWMBYAB41gH4+qEAN6yA3+r//1jSUcMvx6qbbwAwAAwAA8AAMABct1Ys/p07wFGlRuIAYAAYAAaAG7SOLQ4BgAHgWQfgxeuvfgbpj8eYVAubmfvc5z668Y1vXDvMv/3bv9XGLAY4d0x+9atfhcc99thjQ7H/8A//EIqrgqrf18143XBD9X5E7LXvvtUvWK3/+s1vfqPhj4S38ttcq5xtwQ/RHyyvfiss+nr9618fDZXzm8F77VV9+XnsVf0uYfT1pCc9KRqq006rvog99nKu7VnPelbtoFUNTj311CouwxMLfqh+uy+S92677Vab72LAAQccEI49+uijw7HHHXdcOLb6XcLo69Of/nQ0VN/+9rfDsb/+9a/DsV/72tdqY6se8aEPfSjLD9W4C564973vHVoznJ+Eqr5sLfr6nd/5nWioKk2ir5tS5lgoAAAgAElEQVTe9KbRUO26667h2GOOOSYc+/nPfz4cG8m36hHveMc7sjyx4Ic//MM/1I1udKPavL/5zWrL0v7LWWvbP/tgRKef7LNP5CtaBuPe9773zUo5Y81Y6hH777+/bnKTm9Tm/tWvfrU2ZjHA6cN3uctdwuNeccUVodhNmyrJYq/q50Wjrz/6o+qXTGMv57fKzz///NpBq/54zjnnZPWHJT/UJjIMuP76xft79UdstdVW9UHDiD322CMc+/KXvzwce8opp4RjzzjjjHDsJZdcEo695prF70CuP8TxWlv7ytFHlOsznFwEAAwAj7oNAB6qAQAvCAEAD/3gbLwA4OWWAgAPtACABzoAwAMdAODlHgEASw6UAMDL3gGAB1rMKgBXb3/8raR7S9phBTcfGOBoABgABoBXmSgAMAA8agsAeEkN7gAPpeAO8HKP4A4wd4BH+uVCj+AOMHeAh56wPjbDHeCBatwBrifYt0naT1L1fN7VK8KfWn/44A4Pj0DzCPRoo+IRaAkABoAB4FVXEAAYAN7iTVMAGAAGgDfvlzwCvfmbpgEeWQgBgAHgqFf+R9KfS/pO9IAVcQAwd4C5A8wd4LXaB49AD5XhDjB3gFdOEu4Acwd41BN8BnjzHsEdYO4Acwd481WDzwA3JNU1DvuRpNtKuq7hsAAwAAwAA8AAcE0DBYABYAB41UnCl2ANZQGAAeDRGcIdYO4Ar+yYAHBDUl3jsOqrJavnd9/QcFgAGAAGgAFgABgAji4hPAI9VIo7wNwB5g7w2m+KcAeYO8DcAeYOcKVAm98CXf0W0eLv61Tj/pmki6qfs13Rih4Y2NEAwAAwAAwAA8AAcGC5WAYefgZJ/AzSwDHcAeYO8MrewZdgVV+uw88gLfqCL8EaKsEd4Og2Y+242A/hSi8LnAoABoABYAAYAAaAA8sFADwqEneAAWDuAHMHeK3GCQAvKQMAA8DR/cVE4wBgABgABoABYAA4uvDwCPRQKQAYAAaAAWAAuHbpAIAB4FqTTCMAAAaAAWAAGAAGgKPrDwAMAG+xZvAzSPwM0ogpeASaR6C36BHRBYafQRooxe8A1zvm+pHPA49G/1rSDyX9o6TXSfrNGkMtNKr3vOc92nbbbWvPdu2119bGLAZUX4AQfR1//PHRUL34xS8OxVa/ZRt9Pec5z4mG6utf/3o4dt999w3HnnbaabWx11xzjQ4++OAqbpOk+doD/ADrnbobblj8KHr9ia666qr6oGHEDjvsEI79+7//+3DsK1/5ynCsE+jo8PSnPz089CmnnBKOTfKE5Qcn2S984Qvh8P/1v/5XOPakk04Kxz760Y8Ox37lK18Jx1522WXh2Ic85CHh2M9+9rPh2CQ/VOdf8MQVV1yhjRurf1z/9YQnPKEuZOnvr3tdtVTFXp/73OdigZJ+/vOfh2O33nrrcOwhhxwSjt1pp53CsY4v3/nOd9aOOz8/r02bqiUjZd2wesQ3vvGN2nwXA+5617uGY50eHB7UDHR6z5vf/Obw6BdccEE41gxM30ds2FD/9TdO7ZxYpw8fe2zsk4XOXnX77bcPl8PpfR//+MfD45577rnh2KT+sLRmnHPOOdpuu+1q83nve99bG7MY8JnPfCYce/7554djSwh0vB7ht8VrivyixdVXX62HPexhra0Z9V2gmeLPl/RMSSdK+oGk35d0qKRTh9D7AkmVm/6/NYYHgCUBwEvusDYzzgQFgAcaA8ADHQDggQ4A8PLK5GwCAeCBbgDwQAdnLWq21ao/CgDech8BANf7popweh8AvKwpADzQYlYB+MuSniTpwpFpdqfhnd97S7qXpDMk3QEAXrsRAcAAcGyZWj3K2XwBwADwqIsAYAB4ZVfhDvBAEe4AD3TgDvDyDHHWWu4AS9wBXvYOd4AHWvTpDnD1Q2O3lHTdyCJaPcNVPfO1+Hxa9ezpWs9icAeYO8Cj+y/uADcgYWdRBoABYAB49Unm3AXhDvBAQ+4AD3RwenCDFh86hDvAW76Rzh3gkHW4AzwiE49A+z1tlu8Af1TSq0b8c5Sk6sOnfyppF0n/IWnnNaYhAAwAA8CxNWrNKGfzBQADwAAwALxey+EO8EAd7gAPdOAO8PJscdZa7gBzB3i0z3IHeKBGn+4AV485nyWp+naqiyXtKukmkh4pqXo8+qGSdpf0FgB47S0Hj0Bv+c5thAmdxYjPAA8UBYABYAAYAAaA61cYABgAXukSZ88BAAPAAPCWfbZPAFxdXfVNfntLuo2kSyWdWX1LZ/3yshDBHWDuAI9ahUeggxNnNMxZlAFgABgABoAB4PpGCwADwADwlvOEL8Fa1oQvwRpoMauPQNevIutHAMAAMAA85iwCgH0B+RbogWZ8Cdayd/gM8EALHoEe6AAAA8AAMAC83u4CAJ49AK5+0mjxM78vX8ccxwS2pQAwAAwABybKeiEAsC8gAAwAr3QNAAwAj3oCAAaAAWAAGACu31/N0h3g6jO/1Wd8q9davwJ9w/Dzv3XKAcAAMABcN0tq/g4A+wICwAAwALz6vOEO8EAXABgABoABYAC4fn81SwBcr0Y8AgAGgAHg+HxZNRIA9gUEgAFgABgAXq9zAMAAMAAMAAPA9furWQfg6guwbivpS/VSbRaxAMD/+q//qu23X+ungpfjf/KTn4SHP/vss8Oxzu/nRQfdcccdo6G64oorwrGR37VbHGzXXasv5Y693va2t9UG/vKXv9R+++1XxVVffDZfe4AfkPYlWNdcc004m1vesvpp69jrFa94RSxQ0hFHHBGOdQIdAH79618fHvqFL3xhODbJE5YfnGQ//OEPh8PvfOc7h2Pvcpe7hGM///nPh2OPO+64cOwjH7n4gE79Iaeffnp90DDi3HPPDccm+aE6/4In5ubmtHHj4k/Nr53W/vvvH875jDPOCMc6ffiLX/xieNxHPOIR4dhPfvKT4dgLL7wwHHvRRReFY1/ykpfUxs7Pz2vTpmrJSFk3rB5xwgkn1Oa7GHDYYYeFY51Ap24PfvCDw0Nvs8024Vjni4wOPfTQ8Lg/+MEPwrFJfljqEU4i0di3v/3t0VA5+6/ovK/2ytHXn/5p9Uuksdduu+0WC5R0wAEHhGNPPvnk2thq//LrX/86qz8s+aFabyOcseeee9bmvBhw5JFHhmMvueSScKzTp5y1yNmDHnTQQeF8DzzwwHDs+9///trYK6+8Uve85z1b88SG2jM2C7iVpH+SVDmmIoztJFU7jgdIinRNAFgCgJe9Z21mHPADgAciA8ADHQDggQ4A8HLzAYAHWgDAAx0A4IEOAPByjwCAAeBRVAKAB2rMKgC/T9J1ko6W9E1JN5d0a0nVrYI7BZgaAAaAR20CAAcmzcoQ540AABgAHvUPAAwAr+wnADAAPOoJABgAHvUDd4CX1QCAZxuAL5N0B0m/lHS5pFsMrVE90xt5BhgABoAB4AbQO3oIAOwLyB3ggWYAMAAMAK/eP7gDPNAFAAaAAeDVewQAPNsAXD3UfkdJvxoB4ApqvyHpdoFtKQAMAAPAgYmyXggA7AsIAAPAK13DI9ADRbgDPNABAAaAV/YIHoHmEehRTwDAsw3A1SPQFQS/eASAq98Grr4Q6ymBbSkADAADwIGJAgCPKdKKwwFgABgAXn1OAcAA8KgzuAPMHWDuAHMHeL0d2Kx+Brj69udPS6q+wnknSdVXAG4t6c8l/TiwZQWAAWAAODBRAOAxRQKAVxWQR6CXZeEO8EALABgABoBXX2+4A8wdYO4Abzk3ZhWAKyVuIunRw0ehq88EV78vclVwuwoAA8AAcHCyrBXGI9C+gNwB5g4wd4C5A7xe5+AR6IE63AFedgkADAADwADwogLV3d7f+NvPpSMAYAAYAB5jAlWHAsC+gAAwAAwAA8AAcH3vBIAB4FGX8C3Qy2rwGeCBFrN6B7i60/s5SZ+SdLak8+vb6WYRADAADACbk2ZlOADsCwgAA8AAMAAMANf3TgAYAAaAV58nAPBsA/CfSdpT0l9Kuo+kn0r6lyEMn17fWgUAA8AAcGCirBcCAPsCAsAAMAAMAAPA9b0TAAaAAWAAeL1OMat3gEc1qb4I6xBJR0kLYHuj+tY6AOCvfvWr2n776vD1X1//+tfrQpb+/tOfViwee93rXveKBVaUf5+K8+tf73jHO+qDhhHf/va3w7EnnXRSOPbP/7z6LrLY67jjjqsNvPrqq7XXXntVcZskzdce4Acs+CF6mAN+0THduAc/+MHhQz73uephidjrxje+cSyw+gzCb+KfQthvv/3C437wgx8MxyZ5wvKDk+yOO0Z+pnww4hVXVD9rHnudddZZsUBJj3jEI8KxL3zhC8Oxd73rXcOxT37yk8OxZmBqj5ibm9PGjZU91n/tu+++dSFLfz/88MPDsc5nAf/lX6r3hGOvyy6rvkYj9vrsZz8bC5R06aWXhmPPPffccOxb3vKW2tj5+Xlt2lTZIWXdsHrEN75R/UJj7OXc3XEg0fHDX/3VX8WSlXT/+98/HPulL30pHPua17wmHPvRj340HJvkh+r8liechF//+teHw4844ohwbDTwrW99azRUTn9w/PDYxz42nMMb3/jGcGy2Hz71qU9pu+22q83ntNNOq41ZDHC+RHLnnXcOj+v0qfCgkpw17nnPe1546NvdLvKrt4PhTjnllNpxr7nmmsXvHmhlH7Gh9ozNAqrkqjvADx/+Z6uRx6HPCAwJAEsCgJecYi1cAPBANwA40GlWhADAA0EA4GVjOJsDAHigGwA80AEAHugAAC/3EwBYAoCX/QAA9w+Ar5P0XUnV21LV29zxW5kDLQBgAHgUTQDgoRrcAV7uDz7e1h8BAAPAK10CAA8U4Q7wQAfuAA904A7wcqfgDrDEHeBlP3AHeKDFrN4Bfpekh0m6duTOb/W7wNHHYwFgABgAXoXXAGAAeNQWPAK9pMbCmsEj0LIeceQR6IF/nEcLAWAAeOXSDAADwKOeAIBnG4AXvXCX4RdhVV+Gdd9qnak+klJ/H4Y7wJVGPAK9+eY24JuFEB6BHijFI9BRxyzHcQd4oAWPQC97gjvAAy24AzzQAQAGgAHgLddW7gAvawIAA8CVAvcYAnD1WeAKgH8h6fcC21LuAAPAozbhEeihGtwBXhDC8kOg3yyFAMAA8Eq/AMAA8KgnAGAAGAAGgNfbVwDAsw3A/yTpLyTtIOkLw88BV78JfF5wMwoAA8AA8CqTBQAGgEdtwSPQS2rwCPRQCudbXnkEeiAaj0APdHC+9ZfPAC93Yh6B5hHo0XUZAJ5tAK6+E7768qvq911+FYTeLYCHn0GKf3cYP4O0bB8egR5owSPQfufhDjB3gLkDvPq84RHogS7cAR7oAAADwKOdgkegl9UAgGcbgP2d5+ZHcAeYO8BbvCESNRUADABHvbIyDgAGgAFgAHi9/gEAA8Ar/cEdYO4Aj3oCAAaAm+5Bq+MAYAAYAF5lBvEI9IIofAZ46A0egV6aJDwCPZSCR6D9HsEj0APz8Ah0s20rAAwAA8Bbzp1Z/RmkZl1k+SgAGAAGgAHgtfoIAAwAr/QGAAwAN14zAGAAeJxNKwAMAAPAAPA4PWSLxev888/XDjtU36O1/uu886LfrSX9yZ/8Sd1wjf6+yy67hI6rfqcy+jr++OOjoXrFK14Rjt1qq63CsWeddVZt7NVXX63HPe5xVdwm47eea8cdCbCAp4RHoJ2LO+GEE8Lhhx12WDjWCTzkkEPC4W9961vDsUmeWPBD5eMNGzbU5vLb3/62NmYxwLnD7swjp+985StfCed76KGHhmN32223cOwPf/jDcGzk+weqOXn99den94jo7wDvv//+4es744wzwrFZgc7nb+9//8gvDQ4yffe7352S8lOe8pTacefn57VpU7VkpKwb1prxsY99rDbfxQBH38svvzw87h3veMdw7N577x2O3WmnncKxj3rUo8Kxf/VXfxWONQOL2Ec4OV933XXh8J///Ofh2Fvf+tbh2IzAnXfeOTysk6uzZ0/qD9V1LfSIf/zHf9S2225be52nnnpqbcxiwCc+8YlwrBP46U9/Ohz+F39RfQ9x+y+ndo985CPDCZx22mm1sRVr7LPPPq2tGfW7x9qUtgjYWtLzJZ3Y8AuwlowJAAPAQ3dZmxkA2J+0APBAMwB4oAMAvDyHAGC/nwDAA80AYN872cDTKKOagwBgCQBeNgkAPNBiFgG4uu7qNufCW7sNXwvAAwADwADw5jOIO8ALenAHeGgL7gAvzQ/rEWjuAA904w7wQAfuAA904A5wsx0rAAwAjzoHAJ5tAP6kpGMlfblZOxlscAFgABgABoBX6SEAMAC80hYA8FAR5xFdABgAHp1IAHCzHSsADAADwFvOnVm9A/xSSQdLeqekH0ha+PDX8FX9u7oXACyJzwAv2YRHoIdScAd4QQgAGAAGgNdYRQHg5R5Rt9FY/Dt3gAdKAMBRx2weBwADwAAwALyowH+v0UZukHSHQIsBgAHgUZsAwADwFn7gS7AkHoHe/E0yvgRLAoAB4NFmyZdgNXsjPbBPXQoBgAFgABgAdnrGerEAMAAMAK8yQ7gDvLy5BYAB4JEpwiPQQzEAYAAYAF51e2m9ke5sZgFgABgABoBXKlB9y3T1Hfw/dprJ4iOOfAaYzwAPfWMtXHwLtDnbJPEt0APN+BbogQ58C/TyHOJboP1+wrdADzTjW6B97/At0MuaOd+s3EjpmoP4GaRlgfgZpIEW/AxS/UyrflTrjZIOlFT9COd2kh4j6Y8lvar+cL4Eq9KIzwAvOQUAHkrBHeAFIfgM8NAPPAK9eY/gEWgegW7ypimfAR6oxmeAA7vTVUK4A8wd4FFb8C3QAzVm9Uuw3izpDyS9TFL1i9A7StpVUvXt0BUE1714BBoAHvUIAAwAb+EHHoHmEegRU/AI9FAMHoFefpOsbqOx+HcAGACOemW1OAAYAAaAt5wZswrAF0u6W/UE0PA/txhK8wtJNw80GgAYAAaAV5ko3AFe3twCwAAwALxlkwCAAeBRV/AlWEtqWG+kB/apSyEAMAAMAAPAiwpcJmkXSdeNAPA2kr4n6TaBxrLQqM466yxtt1319PT6r2ojHH05m4Of/exn0WF1wQUXhGL32WefUFwV5Hx+aMOG6uPWsdeee+4ZC5T0whe+sDb26quv1r777lvFbZI0X3uAH2AtXM5ngB3d7nnPe4Yz/9rXvhaOdXIIDyrJ0eGxj31seOiPfOQj4dgkT1h++H/tnQvUXVV5rl9yTyAX7gkhiWCgCCIIkVY7UEEZYAu1ldN6O3jwWFEc1mJPa+pQEHTYUWy9gLRQldNSbPGCVRGOqCloRArUgFZ0oJZbYgAJkCsQIAlnzL3W3nv9+7+s793Z88+/2c8aoyXD/1tzzf3Ob35zPuubcy6nskuXLg2bf+ITnwjbOv0+GkvSw6+++upwHfbee++wrVOHT37yk+FyM/lDen7DJ1LMnDMn/XPs64wzzqgzaf39iiuuCNtOmzYtbPvmN785bHvTTTeFbY8//viw7cknnxy2XblyZdj2H/7hH2ptN23apHnz0uKwLOOGFSNuuOGG2vo2DXbsqH7VcezbnLH2Jz/5SbgOS5YsCdtG+kOzsC9/+cvhcj/96U+HbR3fyeQPrRgRrrRh6Iy1d9xxR7jk6Jzjq1/9arjMNK+OXp///Oejpnr5y18eto34Q9J069atueJDyx/+5V/+RbNmpV2bY1+f+9zn6kxaf3eWQH/gAx8Il/u9730vbOuMGTfffHO4XGe+evrpp4fLvfbaa2ttt2zZ0vSznrBGnJpqqzbE4BuSviPp4goAny3pJEkNUqq5AGAAuOoi1mTGGYyczhwdjFLFAeBW8/UkUHXEC8sf6oJN9e8AcKEGANz2CgC40CIyaW2qBgAXSgDAvu8AwO3YE51zAMDOKN+ybcwjAGAJAO7Kf0a96bA0Xkr6paRlkr4v6cWSXirpF4FHAcAAMAAc6ChjmTgvAsgAF0oCwABwZ58CgH2IAYAB4Go/cl6eAMAAcFKADPDw2R0Z4P7IAKeWS+vt0inQ6TCstCT6HyWlvcGRCwAGgAHgSE8ZwwYA9gUEgAFgAHjkfuNADAAMAAPAI/cjlkDHVpMAwABwpwL9sgTan3kOvQMABoAB4J3sRQCwLyAADAADwADwWJGDPcCFOuwBbnuJM9YCwABwNb6wB7hQo9/3AKdsb+T654ARAAwAA8CBjjKWiTMoswS6UBIABoABYAAYAK4ffABgALjqJRyC1VaDQ7AKLQbpEKzO5c37S5pcHoKVPn2Ujk5MS6EX14fW4kRPToFOX5GKXc5hTs7BHJwC3dY/eiBFuoNDsFq6cQiWJE6BzuoPqXBOgS4l5hTohhDWQXmcAl04D6dAx+ZbnVbOy2YywGSAq/5DBrhQo98zwNU2fY+koyWdU34WJ02C/1bSf6WVM4EQAwCTAa66iTWZcQYj58UBAFw0CZ9BKnTgM0iFDnwGqR2q+AxSoQV7gAsdnJfNfAap1Y9yvDS1X4oE5qktE2fOAQADwADw8N71XALglA0+VNKTlZ+ZPuh7l6RFgcACAAPAAHCgo4xl4gzKLIEulGQJdKEDn0Fq9yxOgS604BCsQgf2ABc6sAS6HSOcsRYAjsUSDsEaPrvjFOj+OAX6EUlHSnqw0oQHSEpfek+nQ9ddADAADADX9ZKavzuDMgAMAFfdCQAGgDvDCwAMAFd9AgAGgKv+wB7gthrsAS60GKQ9wNW+cLmkoySdJ+k+Sc+TdH4JwG8LzOsBYAAYAA50lLFMAGBfQDLAhWYAMAAMAI8cP8gAF7oAwAAwADxyjACABxuA03LnT0o6Q9J0SU9J+ryk90raEpiWNgD4tttu0x577FFr/sgjKeEcu5wDQg45JH3COHa95z1p23P99alPfareqLS45557wraf/exnw7Zvf/vbw7ZXXnllre2TTz6ps846K9lNiL07Dvg5e4BrhagYvPGNbwybX3XVVWFbx9DR4V3vele46EsvvTRsm8knrD3hTmXf8IY3hM3f/OY3Z7H9z//8z3C5X/3qV8O2y5cvD9vuv386wzB2PfzwwzHDwiprjFi7dq3mzEnuMfb11re+tc6k9fdnnnkmbPvhD384bHv99deHbV/2speFbZ96Kg23sctZ+XHmmWfGCpV08cUX19pu2rRJc+cmd8jiE1aMcLLbJ598cu1vaxq8+93vDtu+971pehS7FixYEDOUdOKJJ4ZtncMbTznllHC5X/rSl8K2mfwhPd/yCafCzliba87h1Ddq+6EPfShqqgsuuCBsu99++9XaphdN5dw+65ixYsUK7b57Qpaxr3PPPbfOpPX3W2+9NWyb4uBz9Zo9e3b4p9188821tuk7wOVY2BOf2K32iTtnkMrfV9I6Sc8aRQHAkgDglsdYA9dEGIwA4Fbb9SRQdcQOyx+MuCMAuFALAG57DQBcaAEAFzoAwIUOAHA7RkyEOYczzkVtAeC2UgBw1GvadoMMwAl+X1J+9uh+ST80IBgABoCrvc0CnokwGAHAAHDVgZ1sMRlgf6BtZnfIAEtkgBv+Y40ZZICLPkcGuKvYo4kw5+iu5mPfBQADwDvjV4MKwGl9zjWSjpG0XlL6DvDtkl4r6YGAoAAwAAwABzrKWCbOoMwS6EJJMsCFDmSA2z2LDHChBRngQgcywIUOZIDbMcIZa1kCLbEEuu07LIEutHguLYH+gqSUAT47neUkaS9JfydpkqTXB+b1ADAADAAHOgoAvJMiddwOAAPAnR4FAAPAVZ8AgAHgzhgBAIs9wBWnYA9wIcagZoDT549+Q1J1d3faC/hzSfMDU1YAGAAGgAMdBQDeSZEA4BEFJAPclgUABoAB4OFhggwwGeCqV3AIVlsNAHiwAfghSUs7TnxOx4H9d1pdF5iyAsAAMAAc6CgA8E6KBAADwDUuBAADwAAwANyrsZYl0CyBrvoSS6ALNZ5LS6C/KCl9OyJ9WyVlgVP295Lyk0h/FJiyAsAAMAAc6Ci9GpTZA1woyRLoQgcywO2eBQADwAAwANyrsRYABoAB4OG96bkEwAdIuk7SkeUe4HQI1p2STpW0NjCvB4ABYAA40FF6NSgDwABw1ZcAYAC4M7ZwCFahCHuACx1YAt3uIewBZg9wNV6yBLpQY1D3AKffng68Ok7SIklrJN0maUdwTg8AA8AAcLCzjGbmDMoAMAAMAI/ck8gAF7oAwABwtYcAwABw1R/YA9xWAwAebACeXsJuWgbdvKaWUPxUYF4PAAPAAHCgo4xlAgD7ArIEutCMDHDbdwBgALgaScgAF2oAwAAwADzyHAMAHmwA/ndJF0haWXGPl0s6V9JJgWmp9RF7Z6K/efPmwOMLk0suSduWY9f73//+kOHHP/7xkF0yWrduXdjW0eDCCy8Ml/vKV76y1nbbtm36wQ9+kOzSXu/qyd+19wYNGv7wvve9T9Onp3crY19PP/10nUnr744W4UIlHXBA2gUQux54IPJp7FhZVSvHJ5x9SW984xtrK5Og4eqrr87lEw1/+N3f/V1NnZreq4193X57+gR576/0/Oh16aWXRk21fPnysK1zgMbSpelcwtj1xBNPxAwlfeEL6at3Y1/bt2/XXXfdlcsfUrkNn/jYxz6mmTNn1lVHF110Ua1N0+CnP/1p2Paqq64K2379618P2/7hH/5h2NYZ42655ZZwuZG+1izsxz/+cW25adxYtWpVLp+w5hC1la0YOO32e7/3e07Ru9z20UcfDddhn332Cdu+6EUvqrVNMaLsa1nnEekl57Rp02rrs3Jldfo6tvm9995bW17TYMOGDWHba6+9NmT76le/OmSXjObPj3yIpSjuu9/9brjcsi+H7P/jP/6j1i7NIa655ppc8SGVa8WIRx55pLbOTQOnb4QL7UPDNWvS4t/Yddhhh9UapjltOTfpSYxI3+rNcSVPSeStD1sAACAASURBVL1sW6XwKZLS6dCRqGk5pjPRdyYHALAEAHfXPQBgALjqOQAwANwZSRyQAoDtOGzNIZzSnXYDgAtlAeC2hwHAEgDsRJz+tR1UAE6py7T3d2ul6WaUB2DtHWhOa/ACgCVHAyfrCQAHvHUEEwAYAAaA2wqQAR4eJByQAoDtOGzNIZzSnXYDgAHgTt8CgAFgJ970s+2gAnBaAv1lSZdVGu8sSW+SVL+m1lya4MAfGWAJAG57paOFE4gAYAAYAAaAx4oZDkgBwE70bdgCwLZkEkug26KxBJol0NUuxBJoP6AMKgC/VFKC4BWSfi7p0HLvb9qocHNARmvwAoDJAFd9ij3AhRpOv2APcCAqdZiwB7gQhD3AbcdgD3ChBXuACx3IABc6sAS6HSPIAJMB9mcb/XnHoAJwaq0jJL1T0kGS7pOUToGJniYCAHMIVrXHcwhWF/EPAC5E4xCsQgcOwWp3Ig7BKrTgECw/sDqZewAYAO70MAAYAPajTn/eMcgAvDMtBgADwADwzvQgMsAt9QBgALizKwHAAHC34RUALpRzTrolA0wGuNrfOASr2+jTX/cNMgC/TNKZkhZIOk3SMZJmSbop0IQAMAAMAAc6ylgmZIALdQBgABgAHjlSkAH2gywADAB3eg2fQVLzk2ahDgUAh2Tqe6NBBeDXS/qMpC9KekN5IMUySR+TdGKgVQFgABgADnQUAJjvAEfdhCXQbaXIAJMBjvabTjsAGAAGgIf3Hr4D3NbEWR3RbRzqh/sGFYB/Iukd5YFX6yXtKWm6pPRV5P0CDQcAA8AAcKCjAMAAcNRNAGAAuNNXyABHe0/bDgAGgAFgAHisyAEAF+oMKgBvkDSvdJDHJO0laTdJj5b/rht1AGAAGACu6yU1f2cJdCEQS6ALHQBgABgA3smgKgkABoABYAAYAK6PpYMKwD+S9FZJd0hqAnBaAp2+C5z+W3c1AHjjxo2aMyf9c+zrs5/9bJ1J6++/9Vu/FbY98sgjw7ZRw5//PH0VKnbdeOONMUNJxx13XNj22GOPDdtGIGrTpk2aO3duKjP9v03hwuOGDX9YuHChJk2aVHvX6tWra22aBs7nf8KFSrr++uvD5qecckrY1jGMtF2zvKOPPjpc9JYtW2ptd+zYoXJfVA6faPjDV77yFe2+++61dcmlb+2DJ5jBgQceGK7R/Pnzw7bvfve7a22ffPJJnX322dljRG1FSgPndN4vfelL0WI1efLksK0TI0499dRwuZH+2Sxsjz32CJd75pnpSI/YdcUVV8QMC6tsMSJaiRRLotfpp58eNbXsnHKd+lqVmBjGOfwh/bLGuHHbbbcp4veHH354WI3Zs2eHbWfMmBG2XbduXcjWef7mzZtDZSajk08+OWw7fXpa5Bm7IvEhzSuXLFmSKz60/OFDH/qQIm3ivPj667/+65gQaXnsnmmBbOyK1LNZUtIvejmg6tj+6Z/+abQKuvzyy2tt0zyinG/0JEakrGyO6y2SPizpo5I+Luldks6VdL6kqwIPBIAlAcAtTwGAA52m0wQALhQBgAsdAOB2DwGACy0iINBUDQAulHBA1QnbTrkAsKPs0HkEAAwAlx7RmFcCwN5SZQA4FnvS6+L3SjpE0kOSPiXp4titxZs6MsBkgKuBigxwsPeUZgAwAFz1GAAYAO6MIABwoYgDlA6oOhHbKdepr1OHCWLbk+zOCL+FDLAkMsBDX4gAwADwBIl7Qx0TAAaAAeDuuyYADAADwCP3HzLAhS4AMADc/QiT7U4AuJSWJdAsga72MpZA98cS6IMlvUrS3pIekXSDpHuMcEkGmCXQVXdhCbTReZqmADAADAADwGOFDgAYAO5iaMl9CwAMAIs9wMO7GQA88QH4I5LeXzZdgt99Je0ovwH8gWDkBIABYAA42FlGMwOAAWAAGAAGgOsDqbOk2FmqXP/ktoVTrlNfpw4TxBYABoAB4BE6IwA8sQH4dyR9QdI5kv5V0lZJ6ci7N0n6ZPnf6wJBFgAGgAHgQEcZywQABoABYAAYAK4PpA5QOqBa/2QAeASNAGAAGAAGgIcpMNFPgf43Sd+RdOkIbfeOdCCrpD8IDAoAMAAMAAc6CgDMZ5CibsIhWG2l2ANcaMES6EIHADgaRcbFDgAGgAFgALjvAPh+Sekjs2npc+e1j6TbJS0OhFAAGAAGgAMdBQAGgKNuAgADwJ2+AgADwNH4MY52ADAADAADwH0HwOvTd53HCJQbJM0LBFIAGAAGgAMdBQAGgKNuAgADwADwyL2FDHA0ioyLHQAMAAPAAHDfAfBGSSl4jXZtkhrf+K27GgCcPrg8Z069+T/90z/Vldf6+0te8pKw7RFHHBG2jRr+8pe/jJpq5cqVYdtjj02J99j1ile8ImZYfIy51jad1rdo0aJkl3XgWrBggSZNmlRbn5/97Ge1Nk2DuXPHctdwMcMMJ8KEKtJ2zYr/9m//dvjHPv7447W2O3bs0P33pwUhWXyiER+uvPJKzZo1q7Yuufbs1T54ghkccMAB4Rrtv//+Yduzzjqr1jbt3fmzP/uzXP6Qym34RG1FSoPXvOY1UVNdccUVYdvJkyeHbVesWBG2PeWUtHsodm3ZsiVmaC6Bfuc73xku96qrrgrb5owR0UqkWBK9zjjjjKipZecsy7/mmmussvvMOOs84oYbbggt/T/uuOPCsjkrKaZPnx4u99FHHw3ZOs934sOrXpU+5hK7pk2bFjOUdNlll9Xapu8Vv/CFL8w+ZvzlX/6lIm1y3XWR44uKn3XBBRfU/r6mwbx5kZxgYZ3rEKy1a9eG6+vYLl++PFzuJZdcUmub5hF/8Rd/0TOf2K32iZ7BU5IuHOOW95WHYtWVulDSr+qM+PuEU+BASfGeFK8+/hDXaqJZ5vAJ/GGitXK8Pjn8IT0dn4i3wUSzzOET+MNEa+V4fXL4AzEirv9Es8QfJlqL7Pr69MQneg3A35X0bI02JwS0S/VK6YrNAVtMJoYCsyU9EGj/bmqLP3Sj2q6/J5dP4A+7vm27qUEuf0h1wSe6aZFdf08un8Afdn3bdlODXP5AjOimNXb9PfjDrm+DiVaDnvlErwF4oglFfVAABVAABVAABVAABVAABVAABVCgoQAAjCOgAAqgAAqgAAqgAAqgAAqgAAoMhAIA8EA0Mz8SBVAABVAABVAABVAABVAABVAAAMYHUAAFUAAFUAAFUAAFUAAFUAAFBkIBAHggmpkfiQIogAIogAIogAIogAIogAIoAADjAyiAAiiAAiiAAiiAAiiAAiiAAgOhAAA8EM3Mj0QBFEABFEABFEABFEABFEABFACA8QEUQAEUQAEUQAEUQAEUQAEUQIGBUAAAHohm5keiAAqgAAqgAAqgAAqgAAqgAAoAwPgACqAACqAACqAACqAACqAACqDAQCgAAA9EM/MjUQAFUAAFUAAFUAAFUAAFUAAFAGB8AAVQAAVQAAVQAAVQAAVQAAVQYCAUAIAHopn5kSiAAiiAAiiAAiiAAiiAAiiAAgAwPoACKIACKIACKIACKIACKIACKDAQCgDAA9HM/EgUQAEUQAEUQAEUQAEUQAEUQAEAGB9AARRAARRAARRAARRAARRAARQYCAUA4IFoZn4kCqAACqAACqAACqAACqAACqAAAIwPoAAKoAAKoAAKoAAKoAAKoAAKDIQCAPBANDM/EgVQAAVQAAVQAAVQAAVQAAVQAADGB1AABVAABVAABVAABVAABVAABQZCAQB4IJqZH4kCKIACKIACKIACKIACKIACKAAA4wMogAIogAIogAIogAIogAIogAIDoQAAPBDNzI9EARRAARRAARRAARRAARRAARQAgPEBFEABFEABFEABFEABFEABFECBgVAAAB6IZuZHogAKoAAKoAAKoAAKoAAKoAAKAMD4AAqgAAqgAAqgAAqgAAqgAAqgwEAoMFEBONXrAEmbB6IVnhs/crakByQ9m+Hn4A8ZRB2HInP5BP4wDo2X4RG5/CFVFZ/I0GDjUGQun8AfxqHxMjwilz8QIzI01jgUiT+Mg8h99oie+cREBeCFkn7VZ41CdaUDJa3NIAT+kEHUcSoyh0/gD+PUeBkek8MfUjXxiQyNNU5F5vAJ/GGcGi/DY3L4AzEiQ0ONU5H4wzgJ3UeP6YlPTFQAniNp41133aXZsxPsj30tXJjGOq4cCqxZs6a22M2bN+vwww9PdnMlbaq9wTdo+MPkyZO12271Lvvss/Ek9H333efXJnDHjh07AlaFSeQ3NQtzflu4ApKWLFkSNp8yZUqtbarn9u3bc/lEwx9qK1EaRHy4WdaiRYuixVp2v/pV/H3egQem2N77a/Xq1eFCFy9eHLY1DbPGCLMumE8MBXL4hBUjJoIMaXyLXmVsjZqH7SZNmhS2dca4cKGFYQ5/SOX2nU+Yuj1XzSeEPxx77LFhfY844oiwrTMun3feeeFy3/SmN4Vtnd/2N3/zN+FySy4I2xuGPfGJepowatRD00agWrt2rebMSf8c+4pAcl0ZI/09B5g4ZTp1zgVGGzfWc8amTZtUgkNPnHKE393whwReEf0cLR599FFH5rCtMzmI/Kbmg53fFq6spD333DNsPnXq1FrbVM9t27blmsxYE5mIDzd/0Ny5yYV7f6U+Er0iMS9aVtVuw4YN4dvmzZsXtjUNs8YIsy6YGwo4cSpSbCWW5fAJK0Y4v82JwU65uQDYqW+uOkT8oWKTwx9S8S2fiLSLo5v5+8LmkXqmwiZCXcM/KvjSP3N8GOIPkbofd9xxEbOGzYte9KKwrZOE+au/+qtwub//+78ftnV+26c//elwubkSCr16SZYTgF8m6UxJCySdJukYSbMk3RRQDwAOiFQ1yRUAI/AAAA9vLAAYAK56BQDcUiP75NYMnZgHFYhOxoPFVSftOXwCAC4bwpkbAMBt73V0i/q8axftcxOhrs5vi/wuAHi4ogBwb+cRuQD49ZI+I+mLkt5QvnlbJuljkk4MdBQAOCASADyySM5gQAa40JAMcKEDGeBCBzLAZgAeAPPIpNWRIfMEFwAGgDvdkQyw00Ez2kZiSeb4kH6dFSOcLCkZ4MJ5BjUD/BNJ75B0s6T1aX4tabqktKF0v0C/AoADIgHAALDpJqOaA8AAcNU5AOBe9aznTjmRSavzazNPcK3JrfPbnBesTrm5sq9OfXPVwfGLXi1vHOGZALDZELnMI/0ic3wAgCuN68A9S6Dre0XabNbcRPaYpL3Kz1SkDZfp33UXAFynUMffnUHOKZol0I5abVuWQLMEuuo5LIFuqZFjuas9memuVw/2XZFJq6NQ5gkuAFw2hjM3AIDbHuzo5vi9YxvtcxOhrr3+XZnjgz1mOJBIBrjwhkHNAP9I0lsl3SGpCcBpCfRlktJ/6y4AuE4hAHhUhZzBgCXQhYxkgAsdWAJd6EAG2AzAA2AenYxHpcg8wQWAAeBOVyQDHO2cme0isSRzfACAK23swD0Z4PrO8RZJH5b0UUkfl/QuSedKOl/SVfW3F2vzOQU6oFQXg1y8VIkMsKNW25YMMBngqueQAW6pQQa4u5Cyy++KTFqdSmae4ALAXcwNyAC3Pdh5ke74vWMb7XMToa69/l2Z4wMADAAr1yFYSdp0AvR7JR0i6UFJF0m6ONhJAOCgUE2zXAEQADYbojQHgAFgAHjEvgMAdxdSdvld0cl4tKKZJ7gAMADc6YpkgKOdM7NdJJZkjg8AMACcFYB3pgtZg1cE0pqVcZb2OR+Fj/5YB1QdiHLe3G7dujVaXUW++VopbEJMbp966qnw75s5c2bY1mmPcKEZDffYY49w6Q888EDY1vxObQ6fsOLD008/Hf5t06ens/piV2QQb5bk+KRTB6ffO3VI39yOXma/yOEP9mQm+tuSndPO27dvDxft+OW0adPC5R5//PFhW6e+t9xyS7hc0zCHT1gxYsaMGeEqO5o988wz4XJnzUpfiYxdTzzxRMwwo5XTL5x5z3gcghWRxZkbPPnkk5EiGzaO/zgaRyuQVlZGr/Xr0xm2seuoo46KGfrfLc4RH+wxY+nSpeHft2rVqrBtrm1X4Qqkg5n2ihzNVJTobBt0GGpXxIicGWBH/05ba/ACgCVnIgwAt93NGeTMif7O+H9P7gWACxkd0HDg05mcOPDp1MHp904dAOB2F3Ta2ZncOn4JANsh0ZpDAMC2vtaLoV0xuR3hF1k+4cwNAGAJAG57HABcaDGoALwjveQZIQCl1Nz9kj5ffhN4tNejVqACgAHgqq85E31nkAOAC5XJABc6OGDk+CQA7E/G3W86Ok9w2hkAdpRt2ebI8FhzCADYbzenXwDAbX2dGOFoHG1BMsAtpawYQQa40I0McH1PS3t/31nu+b1P0vMk/YmkyyUl6P0/kv5Z0gdGKcpyTAAYAAaAh/ckMsCFJk6mzYFPZ3ICAGeFncZ7mXRwYv3Q5Fs47exMbh2/JANst5vlDwCwra/1AhAABoCrHtaP/gAAA8DRKHmbpP8p6ReVGw4tM7/HlZ9C+pKkgwHg0SV1Mo7OUkiWQLc1JwNcaMEe4EIHALjQgSXQ7RgBABdasAe40MF5ycEe4EKzfgQeZ27AEmiWQFdn8iyBLtQY1CXQ6U383pK2VZxiasqel2/q0/+8RdJop/RYb2/JAJMBrgYfJ9vmDHLOC4nom6KcdmSAC3WdTBsADAB39kkAGACu+gQAXKjh9AsAuO1Bjv84GkfnEiyBbillcQYZ4EI3lkDX97SUAf56+R3gpvX7Jb1O0kskHSjpVkkLRynKckwAGAAGgIf3JAAYAO62X5ABbivnTEKdya3zYoYl0PWTjg4Law7BEmhbXwC4IhkZYDLA1R5EBrhQY1AzwGmZ8/9LyRdJayQtkpS+4/A7khIcnyjpMEl/DwCPPvA4GUeWQLd1JANcaAEAA8AAsD+x77wDAC4UYQl0oYPzkoMl0IVmZIDbUcXxHyf2RCMdGeCWUtZLMjLAhW5kgGM9LZ3seJqkAySlD49daxxSYjkmGWAywN1O9FkCXSjHHuBCB5ZAFzqQAW5HFGcS6kxuyQC3NOYUaEl8BzirP6TCrXmlMzcgA0wGuDoHJQNcqDGoGeAYIo9uZQUqABgABoCHdyYywIUmDmgAwABwZ08CgAtFyAAXOjgvOcgAF5qRAW5HFcd/nNgTnXSTAW4pZXEGGeBCNzLAsZ72gnKp837pvITKLecFbrccEwAGgAFgAHi0uAIAF8o4WwPIALe9yZmEOpNbxy/ZAxyYNQw1seYQ7AG29WUPcEUyMsBkgKs9iAxwocagZoD/h6R/lfQzSYeX/z1C0k2STgiE2tbgFZl8OG8XnbcXe++dDrLujyvXZMrUP8dSttQADX+4//77NWdO+ufY17777ltn0vr7YYelreix684774wZSnLqsG7dunC5U6emw9Rj17x582KGkpw6rF+/vrbcTZs2acmSJckuh09Yk1vnBdn8+fNrf1vTwJn0PPjgg+FyFyxYELZ1DJ04Gen3zrMrtjn8oRUjuqxTz2774Q9/GC7rpJNOCtuuWLEibLts2bKw7QEHpB1KscvZJuH4Ws4Ykfw44stOfX/wgx/EBJN07LHHhm0j9QwXVjF0zhJxXgQ4dYn8tkobPKdjhKNx9KWls3LJ+Qzm6tWrw828ePHisK1pmNUfjj/++NCWn7lzUzVi17Zt1Q/gjH3Ptdem3aGxK9KPmiU5MS329MLKWdXibCN4+9vfXluNxDmXX355suuJT1Qzs7UPNwx+JOkiSf8oKc2W95T0J4kLJFkZ4EiDOw0NAEtONsHUvydOOYKfAcClKABwQwgA2AjG3QyIkX7fRRV6NnCNFiO6rFPPbgOACymdMblXk5mORmzECABYAoBbnmGNGz0LCh0FAcC2slnnlQCw1x4AcL1em0ro3S5pg6SUikqnQN9TfgKprgQywHUKdfydDHBbECf7Sga40I0McKEDGeBCBwDYDMClOQAMAFc9hwxwPJ6QAR4ec8gANzQBgEvXcMZl8yVkeMADgOulekjSQZKeLKH3+DIT/GtJs+tvb2d4Ig3uNDQZYDLAVf8DgAHgqj8AwPEJayCOj2SSdTLTZZ16dhsADAADwMO7kzmPe07HCDLAdrjN6g9kgL32AIDr9fqGpLRQ+2uSPiPpNyQ9kU76l/SK+tsB4IBGQ0zIALflIANcaMEe4EIH9gD7UBKZsLoxqrTPOpnpsk49uw0A9n0tU4aHJdClV7MEutW9WQItiT3AQ/0BAPaGPwC4Xq90kszk8vu/afnzheW+vQ9Kurv+dgA4oBEAPIpIADAAXHUNANiHEgDYjcCFPQDs+xoAXGiWq88BwABwNZoBwABwd6NbcRcAvDPqxe5lD3BMp5YVGeC2YAAwAAwADw8gzlaRXJPxTLCTfuyEyO4AwABwteexBzgO9+wBHh6z2QPc0CTrqiEywB5sAMAxvRZKOmaEPb/p80h1FwBcp1DH3wFgALjTZVgCXShCBtiHEgDYDMClOQDs+1qmCS5LoEufJAPc6ssT4iUZe4Dt2AoAl5I547LzwttpEQC4Xq2zJF1SHoKV9v42r2clRT48CADXazzEAgAGgAHgkTsNAOxDiTPQmqEq62TGrEvPzQFg39cA4EKzXH0OAAaAq4GOJdBD/YEMsDcMAsD1eq2V9B5JX6k3HdECADaFA4ABYAAYAB4rbDhvhHNNxjPBTvrZEyK7AwADwNU+yBLoONyzBHp49GYJdEOTrC9NAWAPNgDger0ekbSvpJTx7eYCgE3VAGAAGAAGgAFgM3D22BwABoAB4OGdKvJCDQAGgEcJxwBwKUykHzU1dF54O8MgAFyv1t9L+pakr9ebjmiRDYD33HPPcJVWr14dtp09O/J543BxDcPo279kO3PmzHDhzpKoSIcbr4Hr/vvv15w5yTXGvpxDsFatWlVXXOvvy5YtC9uuW7cubOvs1XWCj1Pu448/Hq7vhg0bam03bdqkxYsXJ7scg5eV7YvUt/mDnPjgDDBPPFHdCTK2fLNmpa/Fxa5I/2yWtH379lihkiZNmhS2NQ1z+EOqguUTZp3D5k57OP7zi1/8IlyHQw89NGx7wgknhG1vvPHGsK1pmMMnss0hdt999/DPc+Lq1KlTw+U6Y/j69evD5UbG13BhFcNIvxiveUQ39e/lPREt3Oc5Prlly5Zw8T/+8Y/DtkcddVTY1jTMER9aY8bRRx+tyZPTR2vGvk477bQ6k9bfzzvvvLCtM9YuX748XO6FF6aP7/T+csatiK7NGr7mNa+prWya/377299Odj3xid1qn9idwR6SbpG0RtKDHUX870CR2QYvZ4ILAMf2JY3XwAUAe0fQA8BFpAGACx0A4MDIs5MmzuTWmUgAwHbDZJtDOLABABftFukX4zWPsD2pxzdEtHAf6fgkANxStxEjAGDP25xxa1AB+P9Keq2k70nqTC2dEZA72+AFAEvO2+NIsB6vgQsABoDL2GFl+wBgADgw5vTEJBIvmw9yJhIAsN082eYQDmwAwABwp+c6MSLq9Y5PAsAAcNSvRrJzxq1BBeDNko6UdF+XQmcbvABgALjqkyyBLtRwJmoRoGQJ9NDIxxLolh49Wbo0wrhivRTpclyqvc2Z3DoTCQC4VvpOg2xzCAc2nLjKEujWkTHECNPdHZ8EgAFg072GmDvj1qACcALftBHp6S6FzjZ4AcAAMAA8vFc6EzUAuNDPGQgAYAC4s9c5/gMA2zOJbHMIBzacuAoAA8C2l5c3OD4JAAPA3fqZO+8ZVAA+W9JBkt6ftp51IXa2wQsABoABYAB4tJjkxAcHYABgABgAHnUmkCPjl20O4cAGAFy0eWRlxHhtpepiPtrTWyJauA90fBIABoBd/6raO/OeQQXgdPjVfEnPSEqfRKpejaNha65sg5czweUQrIk1cLEHmD3AZdywlrtGMtbNeOTEB2cgAIABYAAYAB5r3kMGmAxw3cR4tL8DwF0pxyFYXcjmzHsGFYD/1xi6XhHQHADmM0hVN2n4AwAMAAPAQ6Onk03gFOjAyLOTJk57OBMJlkDbDZNtDuHABhngot0i/YIMsO3jrRscnyQD3JINAO7C5Zxxa1ABuAtZh9ySbfByMjxkgCfWwAUAA8AAMAA8yuBirQrY2QFqtPsjE/3mvc5EAgC2WyzbHMKBDQAYAO70XCdGRL3e8UkAGACO+tVIds64NUgA/LzKqc8HjyHwPQHxsw1eADB7gKv+xynQhRrORC2ypJhToIdGOZZAt/TIsd8zFQ4AlxIfemg6fzJ2nXDCCTFDSTfeeGPY1jTM4RPZ5hAObDhxlSXQLIE2+03L3PFJABgA7tbP0n0A8MjqbSonIemvO5JOaeVLadr8d/rv5ID42QYvABgABoCH90BnogYAF/o5AwEADAB39jrHf8gAB2YNQ02yzSEc2HDiKgAMANteXt7g+CQADAB362fuvGeQMsCLJKXDr9K1ZAyB7w+In23wCjy7K5Po/jpn+cv06dPDdYk+v/F2Ykd6P5HlyvEmP1W04Q9pSfqcOemfY1/77LNPnUnr71OmTAnbOhpv3pw+hR27ZsyYETOUtHTp0rDtgw8+GLZ1Jmr9BsApGx29Hnmk88y+0e/ca6+9osWG/LZZ2MyZM8PlLlqUwm7sciBq0qRJsUJ9q6wxwq/Orrtj7733Dj/8zjvvDNt+61vfCttedtllYdtbbrklbGsa5vAJaw7h1Ndpt0cffdQpOmzrvDy5/PLLw+W+7W1vC9s6hpF5z6DsAXZ0i744dcaMlStXhqvgjHFHHnlkuFzTMEd8SFVoxIhly5YpMg90kmdf+9rXwj/x9NNPD9vefvvtYdsHHnggbOuwxtatW8PlRvp9s7BTTjmlttxt27ZpxYoVya4nPtHM0NY+eJwNrMHLGQxy/Y4oHDkO4Thl9Pnp9wPAbS+IBL6mtaMxALxJixc3DnzvSaDq6LfWclcAuFAPtfX8pwAAG+9JREFUAM4V/bsr1wEpANjW2JpDOKU77QYAF8pG5j0A8HAvBIAbmuSYQwDAFXdzWAMAdkaL7mytwQsAlhw4A4AB4M5uSQa4UIQMcKEDGeDuBi7nLgekAGBH2YatNYdwSnfaDQAGgB3f6rQFgAHgqk+QASYDPCyeAMAAcNUpWAJdqMES6EIHMsCFDmSAd2Yq2vt7HZACgG39AeBSMpZAt3zHWjlke1yGGwBgABgAbivAEugRggwADAADwMM7BgAMAFe9AgDOMEPdiSIB4JZ4OZY4AsAAcGfvBIAlsQd46AsR9gBLLIHeiYE8w63W4AUAA8AAMAA8WhwiA0wGOMMYtdNFAsAAcLdO5Mx5yACTAa76GQAMAHfGHQC420jcvq9x4k3gWh2wAYDlvZVhD3Dbq1gCXWhBBrjQAQAGgANjzribAMAAcLdOBwB3pRwZYDLAVcfhFOhSDQC4q3gy5Kbmt3/rStql3wGuq1y3f48CaOQ0xGYdHKeMPj+VzSFY7VbmFOhCCw7BKnTgEKxCBw7B6nYkiN8HAAPAcW8ZagkAd6UcAAwAA8AjdB2HNTgFeuTY8/zK/3yCpLdL+oikeyUdJOkDkj4nKfJROjLAZICHBSq+A8x3gEunsCYyZIAL1dgD3NWkOdtNADAA3K1zAcBdKWeNG109occ3cQhWQ9AcZwSkcskAl/4KAPe24/5U0kmSql9iXijp25KOCDwKAAaAAeAROsrSpUsD3acwYQl0oQMADACHO804GgLAAHC37gYAd6UcAEwGeNi8kkOwvO2WZIDrY89GSftL2loxnZXm5OXbnLoSAGAAGAAGgEeLE9ZEBgAGgOsGnF3xdwAYAO7W7wDgrpSzxo2untDjm8gANwQlA1z6Fd8B7o/vAK+QdJ+kcyRtkTRb0ickHSzpVYEY0XeBKn2fKnI5e+umTp0aKbJh4wyIzh7gyJ7lyrOzBqocS6Bnz06u2fvrscceCxc6Z05y99j1zDPPxAwlOfubt2xJ3TR2bdiwodYwQefixY1z8XL4hBUfIvVt/iCnbzh9bu7cJEPsmjFjRszQ3M8fjVHp4U6cCle2MMzhD6lcyyecOjtaROJl89nOuQ3f/nZaPBW7Tj755JhhGoxfFRmOi+JWrEjDepYrh080/CG1R6RNnL7sxGBHLcfPnPquX78+XA3nAMlwoZLbBjn8IWuMcLRwbCO+m8rbb7/9wsX++te/DtvecccdYdsXv/jFYVvTMKs/LFiwIDTenXjiieFq//mf/3nYdvfddw/bPv/51Z2mY9/mzP8++MEPhutw/vnnh22nTZsWto2MRf3yHeAEut+QdKikFH33lPRLSadJujugSLbJTODZXZlEJ5fOIAcAt5qi4Q8AsORMvpwACAAXvgYAFzo4ccoMmFknM2ZdQuaOFtEJa3owANySP4dPAMClvADw0HlEqNNPEKNoPAGAu2qwRowAgCUAuCv/GfOmSZJeKint/f2VpFvTmB98DAAsCQAGgDv7CwDcUMSKD2SACy+KvqQDgIf2OgC40IMMcKGDE4OD8x37pRMZYEdZALiqFhngof4AAAPAXUWTwE27SZpf7v0NmPdvoIpOLp3JFAAMAAPAI4YNALiUxclYR2MUAAwAj9TrAGAA2JnEVW0jmczx2krV7W/YFfdFdEv1IgPcVeuQAS5lIwPclf+MelM68OpTkt5SZn3TIvfXSnqhpI8GHmVNcAPlZTeJTi4B4K6agiXQpWxO9oEl0IVoZIALHaIxCgAGgAHg0ccpJwY7o50zNyAD7Cjbsu27eSUA3Gi7HFskUrkAMADcVSCpu+nvJB0i6QJJ10maJ2mRpG+WEFx3f98Fqujk0hnkyAAPHbjYA+wtvwOAAeBqoI3GKAAYAAaAAeC6SZrz9wjIkQEermhEt3QXGWDHG4fOK1kCzRLorrxnjJvWSDpKUjoKN/3fXqVt80CsuucBwOwBrvoIGeBSDSf7AAADwABw3VBT/3fnpWV0wpqeyiFYLe1zZHg4BKs56eIU6Kaj9d28MhpPAOD6OD6CBRngUhSWQHflP6Pe9JCkA9OquwoAzyxPgD4g8Ki+C1TR7IozmSIDPPRNHRlgMsClR1jxgSXQhWrRGJVsnTgViOdVkxywk8q3fMKps6NFdMIKAA9pgRw+AQADwJ3dPFuMcOKJYxuNJwCwo+rQeSUZYDLAXXnPGDelTyB9R9LFFQA+W9JJkl4XeFjfBaro5NKZTAHAAHBnXyED3FDEig8AMAAcGHNGNXFidnTCCgADwJ0O5/gZe4C76tHWuNHVE3p8UzSeAMBdCU8GuJSNDHBX/jPqTYdJWll++3eZpO9LSl/JTp9F+kXgUX0XqADgRqvmeJPfAh4ywGSAy9hhxQcAGAAOjDkAcI1InAJdCOS8hHT8DgB+tilX1nmE0ya72hYAzj+vJANMBjhHP9+7PAU6HYaVlkT/o6S0NzhyWRPcSIG5bQDg/IEKAPYmX+wBLno9AAwA70z8d8AkOmFN9WEPcKtVcgAPS6BLedezB7jpaH03r4zGEzLAXUV4MsClbGSAu/KfbDe1AlUkADjLgZwab9q0KWw+Z06qcm+vp59+OlzgjBkzwrbOt0NN/XNMZNLvynYIluM7ES2ajeDYOhkFxycXLlwY9onNmzeHbSMTqlTPJUuWpDJz+IQ1kXEAeN68dGB976+NGzeGC507N0kWuxw4c/xs8uTJsQr4Vjn8oRUj/Or0xx333ntvuKIHHXRQ2PbUU08N2153XfqgQ+xy4mruGOHE4sive/LJJyNmDRtnXHbq6dg6beHYhkWQFKnvoJwCHdGiqW20PfbYY49wczz++ONh21WrVoVtjznmmLCtaZh1zDjqqKMUGe/233//cLW/+c30wZveX7/5m78ZLvTWW28N2zo+2Wt+aFbyhBNOqK1vSjR+//tpQXFv5pW71T4xbpC++Ru5/jlgBABLAoBbngIAl1IAwA0hAODSHwDgoTEiMLb0pQkAbDebNYdwSgeAHbUK28gEGwAerisA3NAEAC5dAwCeuADcubw5vS5JaYT0GaQ9Je0ol0IvDoRPa/CKBonAc4eYOLBBBjh/oMqxBNrxncgg3nQgx9bJzDk+SQa4aA0ywIUOjp9F3oi78bS0zzqZ6bJOE/42ANhuImsO4ZQOADtqAcCdajlzg+j8hAyw75PNF+lkgGMvqJoKkwGu97X3SDpa0jmS0jriNOn5W0n/JenT9be3MzyRYBENEoHnAsAdIpn6Z53cAsASANxwUDLAZT8lA9wKWJZPuGPBrrYHgO0WAIBLyZz5kWPrtMhEmkc49c5hG9Gi+dxoewDAXbVUI0YAwABwV94zxk0pG3yopOpmmd0l3SVpUeBh1uAVDRKB5wLAAPCobuIMXI6tk5kDgAHgqoMCwABwZ8BiD/DQGOHE4sgcgQxwRKWhNpE2YAn0cF2jc1sA2PdJMsBtzSL9s2lNBrje1x6RdKSkByumB0j6iaR0OnTdBQCzB7jqI+wBLtUAgAFgAHjE4YMMcCkLAAwAV3tIFKLSPY5t3SSu+vfIBBsABoBH8amsKwvJAJMBdmJZxPZySUdJOk/SfZKeJ+n8EoDfFigAAAaAAeAROgoADAADwADwWGMoAAwAA8ATN0ZEXgY0ax99IUEGOEAVw01YAl1q4vgkGeB6X0vLnT8p6QxJ0yU9Jenzkt4raUv97ewBThpxCnTLU8gAkwEe9kIkEEcaJhyCVSjlLLXnEKyod42PHXuAbZ2tl+hO6SyBdtQqbCMTbDLAw3UFgBuakAEuXYNToCfuKdAjRcX0maV9Ja1Lq2uMsGkNXtEgYTy/Yepk2zgFOn+g4hAszyc5Bbro8QAwAOzG/olmDwDbLWLNIZzSAWBHLQC4U63Iy4DmPdG5LRlg3yfZA9zWzPFJMsAxX0vw+xJJ6bNH90v6oQHB1uAVDRKxaretAOCJ9eYWAAaAy95p7fcEgAFgN/ZPNHsA2G4Raw7hlA4AO2oBwACw7y8dd5ABLgUhA9wfGeAFkq6RdIyk9eV3gG+X9FpJDwS6gzV4AcDSjBkzArIWJr1+gzNeS5cAYAAYAB7azTkFuqWH9VIkHCwniCEAbDeENYdwSgeAHbUAYADY9xcAeGTNAOD+AOAvpG0fks6W9JikvST9naRJkl4f6A7ZJjNONmj//fcPVLUwie7XPfjgg8NlrlmTviYVuzZv3hwzTJuyp6dt2bErsjSinwF48eK0QCF2feMb34gZpnX/+6aV/7HrpptuihlKOvzww8O2L3/5y8O2Dz/8cNh2/fr0TmvsK62eWLJkSTLK8fbWig+5VnLsvns66iB2OXFn7twkWeyK9M9mSVu2RI5fKKydcmM1bVnl8IdUuOUTZp3D5s5WmFmzZoXLdcYC5xAsZyy45557wvU1X0rn8IlsAPyzn/0srMMLXvCCsK3T5/bZZ59wuevWpR1ou/aK/LbxmkfsWiWkrVu3hqsQGWtTYfPnzw+X6bzA+fWvfx0u14k74UILwxzxoTVmpHpHXiQ7v2/FihXhn3j55enM4NjlzOmc2OPEtOc///mxykqaMmVK2PaP//iPa20TZ11xxRU984kEqTmu9Pmj30jbaCuFJyf+eeqrgQdmm8w4E1EAODYRHq+BK0cGGAAueiMAXOjgAAwAHIjkw02yTma6qlEPb3L8BwBuCZ/DJwDgUl4AuOVn2eaVTggBgB21GrY54kNjuJe0EQCWAGDbJ8e84SFJSztOfJ4t6b8lRdKq2QIVAEwGuOq5ADAAXPUHB2AA4K4GjayTma5q1MObHP8BgAHgqutFsqRNezLAXXXabPNKpzYAsKMWAFxViwxwf2SAv5i+uiHpXWUWOE16Lik/ifRHAffPFqgAYAAYAB7eA8kAF5o4AAMAByL5cBMAuNQEAAaAAeC2AuO1kqyrqNXDmwBgW8ysYwYZYDLAtkfW3HCApOskHVnuAd5T0p2STpW0NvAwAFiSs++LPcBtr3LejpMBJgNMBnh4RHayUYF4XjXJOpkx69Jzc+cFCgAMAAPAAPBYQYg9wGSAyQC3FeiXPcCpxunAq+MkLUosJ+m2dABxcMYBAAPAQ7gk7dVgDzCHYJVOYcUHDsEqVOMQrODosxNmAHAhHodgFTo4B9E4L52cl7zsAW51aGvc2IkwMOatZIBtZbO+NCUDTAbY9siaG9Ixwwl20zLo5jW1hOKnAg/LFqhYAs0S6Kr/kQEu1GAJdKGDAzAsgQ5E8uEmWSczXdWohzc5/kMGuCV8Dp/gEKxSXgAYAK6GOE6BHuoPADAA3MMpQKOof5d0gaSVlYLTN1nOlXRS4GEAMBngqps0/IEMMBng0ims+EAGuFCNDHBg5NlJEwC4EJAMcKEDGeBCh0h2mz3Aw4MPS6AbmuR4QZbK5RTo0uU4BXonB/6O2x8pP3e0rfK/pw9CpdOhIx+wsya4TtXJAJMBrvoLGeBCDTLAhQ4OwJABdiJvyzbrZKarGvXwJsd/yABn9QkywKW8ZIBbfpZtXumEEJZAO2oBwFW1OAW6P06BTl9eT3t/q1/8nlEegLV3wP2zBSoAGAAGgIf3QAAYAK56RSRjE4jjI5kAwKUqADAA3G2fYw9wV9En27zSqQ0A7KgFAAPAbQX65RCstAT6y5IuqzTeWZLeJOmVAffPFqgAYAAYAAaAR4tBTgaPDHAgkg83AYAB4E6vyOETZIDJAHf6WbZ5pRMJAWBHLQAYAO4/AH5puQ94haSfSzq03Pv7akk3B9zfClSPPfZYoMjCZK+99grb9pPh5MmTw9V96qnIOWRFcVOmpJXr4SvHRCY93PKH9JYoejkQ88wz1TPdxn7CpEnpEPTYtWNH9HD0WHndWOXaJ5tp/042f5g2bVo38tXe89BDafdH7Jo/f37M0LS64447wnccc8wxYdsJsN9zSIyIZK+dOkfKa4p16623hnVbuzbyRcCiuAMPPDBcrtOXr7nmmnC5F110UdjWNMwxblgxYsaMtEAtdkX3ZabSnHFg6tR0Tmjscsa4WImF1Yknnhg2v+WWW8K25hiXwx+GxIhIxWfOnBkxa9g4B0s55Ub3kDt7OB2fnD49nWcbuzZu3BgzTCfkenOerP6Q9IjE+Je+NGFN7Lrppptihn1otWzZsnCtV61aFbZ1xuRezSt3C9fONzxC0jslHSTpPkmXSvppsBhr8AKAJQC47VnO5AAALnRzJs1OlrRXgaojbljxwfEHALhQGgBue1xkctS0BoCDI/xQsxwTXCtGAMBFgwDAbcd0QBUAlgDgtu8AwIUWgwzAXY2E5U3W4AUAA8BVZ3OABwAGgKu+AwADwJ0DFwBcKEIGuNCBDHChAxngdqQAgAHg6rgBAAPAL5N0pqQFkk5LSQVJsyRF1gYAwObrAzLAbcEAYNN5yAC3BAOAAWAAeOT4AQADwFXPAIAB4Ko/kAFuqwEADzYAv17SZyR9UdIbyj2caeH4x9Iqm8D0HAAOiFQ1AYABYNNlhpizBLqQAwAGgAFgAHisWEoGmAxwp3+QASYDTAZ4eNQc1CXQP5H0jvLAq/WS9pSUdtOvkbRfYKIOAAdEAoBHFokMsOk8ZIDJAHe4DHuA24KwBLrQggxwoQMADAADwMPnGGSA25qQAS60GFQA3iBpXukO6YjmdPRyOnDr0fLfdTN0ALhOoY6/kwFuCwIAm84DAAPAAPConQYABoCrzgEAA8AAMAA81iwLAB5sAP6RpLdKSt/daAJwWgKdvgscOUMbADYZBgAGgE2XGWLOEuhCDpZAFzqQAW53DwAYAAaAh48u7AFua8ISaJZAV3sIADzYAPwWSR+W9FFJH5f0LknnSjpf0lWBiToAHBCpagIAA8CmywDAIwgGAAPAnW4BAAPAADAAPNb4CgADwADw8B4yqEugkxLpBOj3SjpE0kOSPiXp4uAkHQAOCtU0A4ABYNNlAGAAeFSXIQNMBrjTOdgDXCjCEuhCBzLA7R4CAAPAADAAvDNz8Oq9ALCpJAAMAJsuAwADwABwoNOQAS5EAoAB4Gp3AYAB4Ko/cAhWWw2WQBdaDGoG+GBJr5K0t6RHJN0g6Z7AXKNp0gDgtWvXas6c9M+xr9mzZ9eZtP7uTGbChUp69tlnQ+bO86NlpgdPmjQp9PxktGPHjrDtE088UWub9o/Onz8/2c1N5ynV3uAbNPxhypQpcvSLPObhhx+OmNk2udrOKXfbtm3heu+7775h28jLllTP0s9y+IT1gmzDhnQmX+yaN695dl/MPmrlTA7mzk2S9f5at25duFDHHyJ9suK3Ofwh/S7LJ8JCYNjV2OnEqUzjhuUPER9uCuHYOjo4Y7hTruPCBx10UNj87rvvDtuahtljRKQNc2lsavGcNDf1z+4PEZGXL18eMWvYXHvttWFbZ36ycuXKcLknn3xy2Paoo44K21599dVh24wxrSc+kU5m7vX1EUnvLwtN8Jtm1om40jeAPxB8GABsQHXS1HE0ALjthQBwoYUDPABwMIpVzADg1gvCngxcI7SABTx+C3JHZNLaVMmEhxw+YfmD89scW0cHZwx3ynU8FwBuq5VLY6c9nqu2kT400V6aAsCFNwLAo/fK35H0BUnnSPpXSVslzZD0JkmfLP97XaBTA8AAcNVNyACXajiDMhngQjTnDSsZYP+FSD9OZgJjECYdCkTaGQAeKpoTrwHglnY5XoikwlsvRSK+7LQdwcJTwNQ/uz9Eag8AA8B1fvJvkr4j6dIRDN8h6RRJf1BXSDNQsQQ6tqw66ekMnmSA2x5IBtgHHjLAgQjWYUIGmAyw7zUT647IpBUABoB74LXZgSfiywBwD1pylCJM/bP7Q+SXAsAAcJ2f3C/p2HLfb6ftPpJul7S4rhAAuFDICcAAcMCrRjABgAHgqluQAfb9oR8nM91Fi8G+K9LOADAA3INekh14Ir7szL968JsHqghT/+z+EBEfAAaA6/xkvaQ9xzBKp9FETplhCTQAXHUjlkCXajiDMkugC9FYAl3owCFYdcMXfx9LgcikFQAGgHvQi7IDT8SXnbG2B795oIow9c/uDxHxAWAAuM5PNpYnOo5ml04Irj/WudyrwRJolkCXjgQAA8DDXojUBaPm3wFgADjqK9iNrkBk0goAA8A96EPZgSfiywBwD1pylCJM/bP7Q+SXAsAAcJ2fPCXpwjGM3lceilVXDhlgMsDDgIfPIHnL4skAFy4EAAPAdQMOf69XIDJpBYAB4HpPqrXIDjwRXwaAa9upawNT/+z+EPkhADAAXOcn303cVmN0Ql0h7AEuFHICMHuAA141ggl7gAtR+AxSoQN7gH1/6MfJTHfRYrDvirQzAAwA96CXZAeeiC87868e/OaBKsLUP7s/RMQHgAHgiJ/0woYMMABc9SOWQJdqOIMyGeBCNDLAhQ7sAe7F0DS4ZUQmrQAwANyDHpIdeCK+7Iy1PfjNA1WEqX92f4iIDwADwBE/6YVNA3juuusuzZ49u7a8hQsX1to81w0iAaWbyclDDz1UK93mzZt1yCGHJLusgSp9fsf5nbUVl3T33XdHzGwbZ/B0svdOudu3bw/X++CDDw7bRuqb6lnWNYdPtL7nGKn06tWrI2YNm8WLI4fUh4trGa5ZsyZ806JFi8K2juE999wTNnf8IVxoYZjDH1K5lk+YdcY8rwI5fCKbPzhjkBOvI3G1mzHcabrnPe95YfN77703bGsa5vAHYoTZCBPIfEL4wznnnBOW5Prrrw/bOp9I/OY3vxku93Wve13Y9oUvfGHY9sorrwzbzp2bmi7L1ROf2C1L1Xa+0ES0v9r5YihhnBU4UNLaDM/EHzKIOk5F5vAJ/GGcGi/DY3L4Q6omPpGhscapyBw+gT+MU+NleEwOfyBGZGiocSoSfxgnofvoMT3xiYkKwKleB0ja3EcNMuhVTan6BwJ7wLvRCX/oRrVdf08un8Afdn3bdlODXP6Q6oJPdNMiu/6eXD6BP+z6tu2mBrn8gRjRTWvs+nvwh13fBhOtBj3ziYkKwBNNcOqDAiiAAiiAAiiAAiiAAiiAAijQ5woAwH3egFQfBVAABVAABVAABVAABVAABVAgpgAAHNMJKxRAARRAARRAARRAARRAARRAgT5XAADu8wak+iiAAiiAAiiAAiiAAiiAAiiAAjEFAOCYTlihAAqgAAqgAAqgAAqgAAqgAAr0uQIAcJ83INVHARRAARRAARRAARRAARRAARSIKQAAx3TCCgVQAAVQAAVQAAVQAAVQAAVQoM8VAID7vAGpPgqgAAqgAAqgAAqgAAqgAAqgQEwBADimE1YogAIogAIogAIogAIogAIogAJ9rgAA3OcNSPVRAAVQAAVQAAVQAAVQAAVQAAViCgDAMZ2wQgEUQAEUQAEUQAEUQAEUQAEU6HMFAOA+b0CqjwIogAIogAIogAIogAIogAIoEFMAAI7phBUKoAAKoAAKoAAKoAAKoAAKoECfKwAA93kDUn0UQAEUQAEUQAEUQAEUQAEUQIGYAgBwTCesUAAFUAAFUAAFUAAFUAAFUAAF+lwBALjPG5DqowAKoAAKoAAKoAAKoAAKoAAKxBT4/2w4ZzdA2wjEAAAAAElFTkSuQmCC\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/jhamrick/miniconda3/lib/python3.5/site-packages/ipykernel/__main__.py:3: RuntimeWarning: overflow encountered in exp\n",
      "  app.launch_new_instance()\n"
     ]
    }
   ],
   "source": [
    "alpha = 0.01\n",
    "M = np.random.randn(7, 100)\n",
    "bM = np.random.randn(7, 1)\n",
    "N = np.random.randn(100, 7)\n",
    "bN = np.random.randn(100, 1)\n",
    "k = data.shape[0]\n",
    "\n",
    "fig, axes, ims = init_plotting(data, M, N)\n",
    "\n",
    "for i in range(10000):\n",
    "    for j in range(k):\n",
    "        X0 = data[[j]]\n",
    "        X1 = add_noise(data[[j]])\n",
    "\n",
    "        # compute hidden representation\n",
    "        Y = sigmoid(np.dot(M, X1.T) + bM)\n",
    "        # reconstruct the input\n",
    "        X2 = sigmoid(np.dot(N, Y) + bN).T\n",
    "\n",
    "        # do backpropagation to compute weight\n",
    "        # updates\n",
    "        d1 = 2 * (X0 - X2) * dsigmoid(X2)\n",
    "        dN = alpha * np.dot(d1.T, Y.T)\n",
    "        dbN = alpha * d1.sum(axis=1)[:, None]\n",
    "        d0 = np.dot(d1, N).T * dsigmoid(Y)\n",
    "        dM = alpha * np.dot(d0, X0)\n",
    "        dbM = alpha * d0.sum(axis=1)[:, None]\n",
    "\n",
    "        # update weights\n",
    "        N += dN\n",
    "        bN += dbN\n",
    "        M += dM\n",
    "        bM += dbM\n",
    "\n",
    "    # update plots every 100 steps\n",
    "    if (i % 100) == 0:\n",
    "        draw(fig, data, M, bM, N, bN)"
   ]
  }
 ],
 "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.5.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}