Spike_IO-2-Main.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# imports of both spynnaker and external device plugin.\n",
    "import spynnaker.pyNN as Frontend\n",
    "import spynnaker_external_devices_plugin.pyNN as ExternalDevices\n",
    "\n",
    "#######################\n",
    "# import to allow prefix type for the prefix eieio protocol\n",
    "######################\n",
    "from spynnaker_external_devices_plugin.pyNN.connections\\\n",
    "    .spynnaker_live_spikes_connection import SpynnakerLiveSpikesConnection\n",
    "\n",
    "# plotter in python\n",
    "%matplotlib notebook\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "import time\n",
    "import numpy"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "n_neurons = 3000\n",
    "\n",
    "run_time = 1000 # timeScaleFactor  = 10"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Allocating virtual identifiers\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Partitioning graph vertices\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Partitioning graph edges\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Generating partitioner report\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Filtering edges\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Placing graph vertices\n",
      "|0                           50%                         100%|\n",
      " ================================================================================================================================================================================================================================================================================================================================================================================================================================================================================================\n",
      "Generating placement report\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Generating placement by core report\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Generating SDRAM usage report\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Routing\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Allocating tags\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Reporting Tags\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Getting the number of keys required by each edge using theapplication graph\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Allocating routing keys\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Generating Routing info report\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Generating routing tables\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Generating Router table report\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Compressing routing Tables\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Generating compressed router table report\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Generating comparison of router table report\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Generating sPyNNaker data specifications\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Loading routing data onto the machine\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Clearing tags\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Loading Tags\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Executing data specifications and loading data\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Finding binaries\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Loading executables onto the machine\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Initialising buffers\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Updating run time\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Creating database\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Loading buffers (0 bytes)\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Getting spikes for pop_input\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Getting spikes for pop_output\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Getting provenance data from machine graph\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Getting provenance data\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Getting Router Provenance\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n"
     ]
    }
   ],
   "source": [
    "# initial call to set up the front end (pynn requirement)\n",
    "Frontend.setup(timestep=1.0) #, min_delay=1.0, max_delay=144.0)\n",
    "\n",
    "\n",
    "cell_params_lif = {'cm': 0.25,\n",
    "                   'i_offset': 0.0,\n",
    "                   'tau_m': 20.0,\n",
    "                   'tau_refrac': 2.0,\n",
    "                   'tau_syn_E': 5.0,\n",
    "                   'tau_syn_I': 5.0,\n",
    "                   'v_reset': -70.0,\n",
    "                   'v_rest': -65.0,\n",
    "                   'v_thresh': -50.0\n",
    "                   }\n",
    "\n",
    "\n",
    "weight_to_spike = 5.0\n",
    "\n",
    "# create ORDINARY populations\n",
    "pop_input = Frontend.Population(\n",
    "    n_neurons/10, Frontend.IF_curr_exp,\n",
    "    cell_params_lif, label='pop_input')\n",
    "\n",
    "pop_output = Frontend.Population(\n",
    "    n_neurons, Frontend.IF_curr_exp,\n",
    "    cell_params_lif, label='pop_output')\n",
    "\n",
    "# Create EXTERNAL INJECTION population\n",
    "spike_injector = Frontend.Population(\n",
    "    n_neurons, ExternalDevices.SpikeInjector,\n",
    "    {'port': 12345,\n",
    "     'database_notify_host':\"localhost\",\n",
    "     'database_notify_port_num':12346},\n",
    "    label='spike_injector')\n",
    "\n",
    "# port: internal port on SpiNNaker (you just can't use the same value to multiple SpikeInjectors)\n",
    "#       it's only really usefull if you are NOT using the spynnaker_live_spikes_connection.\n",
    "# database_notify_host: ip where the spikes are comming from (if using spynnaker_live_spikes_connection)\n",
    "# database_notify_port_num: port where the spikes are comming from (if using spynnaker_live_spikes_connection)\n",
    "\n",
    "# Create a connection from the injector into the input population\n",
    "Frontend.Projection(\n",
    "    spike_injector, pop_input,\n",
    "    Frontend.OneToOneConnector(weights=weight_to_spike),target='excitatory')\n",
    "\n",
    "\n",
    "# Propagates the spikes from pop_input to pop_output\n",
    "Frontend.Projection(\n",
    "    pop_input, pop_output,\n",
    "    Frontend.OneToOneConnector(weights=weight_to_spike))\n",
    "#     Frontend.AllToAllConnector(weights=weight_to_spike),target='excitatory')\n",
    "\n",
    "\n",
    "# record spikes so that we can read off valid\n",
    "# results in a safe way afterwards, and verify the behavior\n",
    "pop_input.record()\n",
    "pop_output.record()\n",
    "\n",
    "# Activate the SENDING of live spikes from the pop_output\n",
    "ExternalDevices.activate_live_output_for(\n",
    "    pop_output, \n",
    "    database_notify_host=\"localhost\",\n",
    "    database_notify_port_num=12346)\n",
    "\n",
    "\n",
    "# Run the simulation on spiNNaker\n",
    "Frontend.run(run_time)\n",
    "\n",
    "\n",
    "# Retrieve spikes from the synfire chain population\n",
    "spikes_input = pop_input.getSpikes()\n",
    "spikes_output = pop_output.getSpikes()\n",
    "\n",
    "# Clear data structures on spiNNaker to leave the machine in a clean\n",
    "# state for future executions\n",
    "Frontend.end()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "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,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4Xu3dDbhdVX3n8e+FGgZfgoqpUxAtWAEtKqhUGFEUldSqpdra1pdCq1StQVNobDGdBoioCJqpStACYkeKL9WpVutLKoPVqi0VtS1lWrRiFUVjRAgoI8HcO8//ZJ1hc3tvcnP+5667z97f/Tx5OOGetc9en/Xf5/6y9tsULgoooIACCiiggAK9EpjqVW/trAIKKKCAAgoooAAGQItAAQUUUEABBRTomYABsGcDbncVUEABBRRQQAEDoDWggAIKKKCAAgr0TMAA2LMBt7sKKKCAAgoooIAB0BpQQAEFFFBAAQV6JmAA7NmA210FFFBAAQUUUMAAaA0ooIACCiiggAI9EzAA9mzA7a4CCiiggAIKKGAAtAYUUEABBRRQQIGeCRgAezbgdlcBBRRQQAEFFDAAWgMKKKCAAgoooEDPBAyAPRtwu6uAAgoooIACChgArQEFFFBAAQUUUKBnAgbAng243VVAAQUUUEABBQyA1oACCiiggAIKKNAzAQNgzwbc7iqggAIKKKCAAgZAa0ABBRRQQAEFFOiZgAGwZwNudxVQQAEFFFBAAQOgNaCAAgoooIACCvRMwADYswG3uwoooIACCiiggAHQGlBAAQUUUEABBXomYADs2YDbXQUUUEABBRRQwABoDSiggAIKKKCAAj0TMAD2bMDtrgIKKKCAAgooYAC0BhRQQAEFFFBAgZ4JGAB7NuB2VwEFFFBAAQUUMABaAwoooIACCiigQM8EDIA9G3C7q4ACCiiggAIKGACtAQUUUGCHwDRwJrC+gMTrdcBPlJ/ppIACCnRGwADYmaG0IwookBT4OeCbwA1lPWeUAHg3A2BS1uYKKNA6AQNg64bEDVJAgZYIGABbMhBuhgIKjF/AADh+U9eogAJ1BR4CnAv8N2A58F3g74HnAo8HPgn8CvBM4ARgT+DDwCuA7zc2dfYh4LkC4M8D7wPeUdpH873LoePnAPsD3wIuBl4HzJT13wM4B/hF4P7AVuCfgFOAL9fl8tMUUEABMABaBQooMOkCXwFuLAEr/hsh7BeAFwKPKwHweuBy4L1ABMYIZ58HnrwbAfBE4KIS9qJ9LBEm/wY4tJw7+C/AUeXQ8fnAK8v7ot0zgFcB/w7sW7bt/cA/TPoAuP0KKDB5AgbAyRszt1gBBe4UiCC1pcys/dUcMMeWAPgx4OmNnz8PuBR4Svl5/GhnM4BrgFcDLy2zf8NV/Qbwp8ATgM821r+2hMAHAN8DrgY2AbEeFwUUUGDJBQyASz4EboACCiQFYkbtR8Afl9m4+PtwGQbAmA2MoDZc4sKO/1tC2mvL/5wvAG4ETgZ+FZgdMv+sHHqOWcXm8ijgykYwvaS8fiPw18CXvLAkOeo2V0CBlIABMMVnYwUUaIHAT5fDsnHY937A14DzgLcBwwAYs38xC9hcYubwz4FVOwmAcR5gzODFOp9YQmNzHRHmYhZxriXO/4vg+T+BOAcwZgXjPMEHAzcB7wT+cI51toDUTVBAga4LGAC7PsL2T4F+CTyiXFgRM3ZPKzODcRFIZgbwcOATwLVlnbc1SN8NHFmC3Vzfp/8x60KTaHpAuSjl9UDMCMZ5gS4KKKBAVQEDYFVuP0wBBSoI3KtcZRsXYFxVzvH7eLkwZPjxzy8zcAs9B/Bg4IpyAUcEyx+WFZ0EXAg8fISreb9Q7jkYVye7KKCAAlUFDIBVuf0wBRQYs0AErzeVq3vj3L+4Kve3gGeXq3HjtjAxAxhXAf9v4D3AIcDZ5Srg5uHbXd0GJs7ziyt+43Bw3A7mB+UpITE7GD+L2by4tcsy4Gcat52J8xM/B3yoXAwS7eJw8n8HTgXiamEXBRRQoKqAAbAqtx+mgAJjFlhR7gF4NBBX3EbYiituX1Nu+xLnAMbM3S+XizB+qYTECGOrZx2e3V7OJYyrfWOJ8//+qAS6CIexxPl7w0C5soTACHynA78OHFhmB79aLhiJ7Yi2cduY44GDSmi8rswcxgUmLgoooEB1AQNgdXI/UAEFKgoMA+BTSxCs+NF+lAIKKNBega4FwHhwe9ysNa4EjH/NX1Nuzho3gB0ucZL4W4BHAzeXG7ueNWuI4u9xEnkcPorzdOIqwViXiwIKTJbA8CrgONQbM4EuCiiggAJ070kgcR5OPAYqHrP0E+VRTXGPr7hZbJy0fc9yonbck2s9ECd2x60h3lDOI4qiiBPH4/FMcaJ3HMaJw0ARKuO9zav/LCAFFGi/gDOA7R8jt1ABBZZAoGszgE3CvYDfAX63nKB9BxBX7MWtF/Zr3IQ1ngf68vKeaB/n5mxonJgdJ5XfAJwGXLYEY+RHKqCAAgoooIACYxXoYgCMm8FGUNsH+Hp5/NP/KWoR7B5aZveGkHHy+GfK+/coh4Xj/8Vd/IdLPMIpTiz3MU5jLT9XpoACCiiggAJLIdDFADh0vHd5dmfcY+tnyyHgi8sd+Z/bwI6HuMf5fXFz1giA3yghMW76Olzi1hG3AC+eNUjhF7OJty7F4PmZCiiggAIKKDCyQNwzNI7wxVN7erd0OQDGYEb/4nzAOPT7gXJod5wzgPsD3+xd1dhhBRRQQAEFuiEQt4/6Vje6snu96HoAjAtBIgDGvb/iZq1xMce5s84BjHuBxUUfw4e5z3UO4LfLDVtnnwMYVwlvvf7661m+PF66LKXA2rVree1r45ofl6UWcCyWegTu+vmOR3vGw7Fox1jccsstHHBAHPgbnC4WR/h6t3QtAMYFHXG4Nq4EjhvExk1Yn1we0xRX8MZVwHFoN64Cjp9F6PtIuYN/PE0gljjPLwJhPDw+wmDcWuYF5ekBs68CHgTArVu3GgBbsOucdtppbNgQp3m6LLWAY7HUI3DXz3c82jMejkU7xiIC4D77RPYzALZjRPJb8WHgMSXo3QR8qgS4eHTTcDkMuKDcBzBmB99azhVsfvqZwEuAOD8gniU6330ADYD5MRvbGvxiHRtlekWORZpwrCtwPMbKmVqZY5HiG1tjA2D37gM4tuJY4IoMgAuEqvG2TZs2sXJlPJ3LZakFHIulHoG7fr7j0Z7xcCzaMRYGQANgthINgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJWcAzAraXgEFFFBAgcoCBkADYLbkDIBZQdsroIACCihQWcAAaADMlpwBMCtoewUUUEABBSoLGAANgNmSMwBmBW2vgAIKKKBAZQEDoAEwW3IGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkZALOCtldAAQUUUKCygAHQAJgtOQNgVtD2CiiggAIKVBYwABoAsyVnAMwK2l4BBRRQQIHKAgZAA2C25AyAWUHbK6CAAgooUFnAAGgAzJacATAraHsFFFBAAQUqCxgADYDZkjMAZgVtr4ACCiigQGUBA6ABMFtyBsCsoO0VUEABBRSoLGAANABmS84AmBW0vQIKKKCAApUFDIAGwGzJGQCzgrZXQAEFFFCgsoAB0ACYLTkDYFbQ9goooIACClQWMAAaALMlZwDMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJWcAzAraXgEFFFBAgcoCBkADYLbkDIBZQdsroIACCihQWcAAaADMlpwBMCtoewUUUEABBSoLGAANgNmSMwBmBW2vgAIKKKBAZQEDoAEwW3IGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkNAuCWLVtYvnw527ZtY9myZYN1dvl19K8P/Zyvj33sv33ux77drHnHvB9j3sdxjjqPALhixYro/j7ALdkwMIntpyZxo1u0zYMAeFyjevYFbiwb2NXX0b2u9m2+ftlnx7wPNW+dW+ddr/PoXyx3AFfseGkAbFGomqRNGQTAdwF7A1cBRwIzHX8dA9SHfjb7aJ8d8z7UvHVunXe9zmPW64QyaRPJzwA4SZGrXds6CIDPAy4FXlS27e0dfx0BN3airvczhnPYR/u8o7gd827v59a5dd71/fzk8j12a0l+BsB2hapJ2ppBADwVOLbM/F0JHNWD13s2ZgHtc3fHvI/j3Nx/+9h/+9zd/Xm+3019G/MfA8cbAAcTOS6jCwwC4GOY4jBmuBg4hj04lOnOvz6E6UHx2Oduj3kfx7m5//ax//a5H9/hfa7zx7MXH+V27uMh4NHTjy0ZBMB7s4Ij2ML92IPLWcHhbO7w6yku50AO57rB8Pejz3eOq33ucm0391/r3H3b7/Pu/j67P8fyXT7Ids8BNMiNLDAIgJs3b+bN69czPT3Nug0bOHvNmk6/Pv2cczh37VpmZmYGf/rQ52Eft2/fztTUlH22zju9n1vn3f4Ob35n9/X7fN3q1Zx34YUGwJHjjw29EbQ1oIACCiigwIQJeCNobwSdLVkDYFbQ9goooIACClQWMAAaALMlZwDMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAbB7AfB1wNOBBwE/AD4F/D7wzUZtTQM/AuJ50PE42xngaOCaxnvOAk6GwaPevgCsmvXz4VsNgJV3Wj9OAQUUUECBrIABsHsB8DXA+4GrgbsDbwUeBhwxKwA+GfjkPAX0SuAU4GnAV4EzgBOBg4HbZrUxAGb3QtsroIACCihQWcAA2L0AOLuEHgl8EbhvPLO3/DBmAJ8CXDFPvV0HbADOLz/fE7gBOA24zABYeS/14xRQQAEFFBizgAGw+wEwDv++FDho1gzgd4C7AV8H3gZcXH4eM3o3l0PCVzbabCqzimsMgGPeC12dAgoooIAClQUMgN0OgDHL9wHg2cAnGrX1JOBzwHbgqWVW71XAnwAPAL4BPBS4ttHmPcAtwIsNgJX3Uj9OAQUUUECBMQsYALsbAJ8BXAqcBHxoF3WzDjgeOKZc9LHbM4CrVq1i2bJlg49ZuXLl4I+LAgoooIACCrRHYNOmTcSfWLZt28bGjRvj5T5lgqc9G1ppS+Iq2K4tzy/n7z0HuHwBnYsAGIntceW9c50D+G3gVM8BXICmb1FAAQUUUKDlAs4Adm8GMK7eXQ88E/jsHPUXVwNH6I2rhONikLga+N1AhMDBPwWAOM8v1hO3k4kwGD97AXCIVwG3fI928xRQQAEFFFiAgAGwewEwQt0dwO1l/If3+YtbukQgjEPD55Zz/eI+gHERyAXARbPq5UzgJcC9gKu8D+AC9ibfooACCiigwIQIGAC7FwBrl573Aawt7ucpoIACCiiQFDAAGgCTJTR4UsjWrVu3snx5vHRRQAEFFFBAgbYLGAANgNkaNQBmBW2vgAIKKKBAZQEDoAEwW3IGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkZALOCtldAAQUUUKCygAHQAJgtOQNgVtD2CiiggAIKVBYwABoAsyVnAMwK2l4BBRRQQIHKAgZAA2C25AyAWUHbK6CAAgooUFnAAGgAzJacATAraHsFFFBAAQUqCxgADYDZkjMAZgVtr4ACCiigQGUBA6ABMFtyBsCsoO0VUEABBRSoLGAANABmS84AmBW0vQIKKKCAApUFDIAGwGzJGQCzgrZXQAEFFFCgsoAB0ACYLTkDYFbQ9goooIACClQWMAAaALMlZwDMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJWcAzAraXgEFFFBAgcoCBkADYLbkDIBZQdsroIACCihQWcAAaADMlpwBMCtoewUUUEABBSoLGAANgNmSMwBmBW2vgAIKKKBAZQEDoAEwW3IGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkZALOCtldAAQUUUKCygAHQAJgtOQNgVtD2CiiggAIKVBYwABoAsyVnAMwK2l4BBRRQQIHKAgZAA2C25AyAWUHbK6CAAgooUFnAAGgAzJacATAraHsFFFBAAQUqCxgADYDZkjMAZgVtr4ACCiigQGUBA6ABMFtyBsCsoO0VUEABBRSoLGAANABmS84AmBW0vQIKKKCAApUFDIAGwGzJGQCzgrZXQAEFFFCgsoAB0ACYLTkDYFbQ9goooIACClQWMAAaALMlZwDMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJWcAzAraXgEFFFBAgcoCBkADYLbkDIBZQdsroIACCihQWcAAaADMlpwBMCtoewUUUEABBSoLGAANgNmSGwTALVu2sHz5crZt28ayZcsG6+zy6+hfH/o5Xx/72H/73I99u1nzjnk/xryP4xx1HgFwxYoV0f19gFuyYWAS209N4ka3aJsHAfC4RvXsC9xYNrCrr6N7Xe3bfP2yz455H2reOrfOu17n0b9Y7gCu2PHSANiiUDVJmzIIgO8C9gauAo4EZjr+OgaoD/1s9tE+O+Z9qHnr3Drvep3HrNcJZdImkp8BcJIiV7u2dRAAnwdcCryobNvbO/46Am7sRF3vZwznsI/2eUdxO+bd3s+tc+u86/v5yeV77NaS/AyA7QpVk7Q1gwB4KnBsmfm7EjiqB6/3bMwC2ufujnkfx7m5//ax//a5u/vzfL+b+jbmPwaONwAOJnK6tLwOeDrwIOAHwKeA3we+2ejkI4C3AI8GbgYuAs6ahRB/j38oRMD7ArAKuGYOqEEAfAxTHMYMFwPHsAeHMt3514cwPSge+9ztMe/jODf33z723z734zu8z3X+ePbio9zOfTwE3KX8x2uA9wNXA3cH3go8DDii9PKewJeBS4D1wMHAx4A3AG8q73klcArwNOCrwBnAieW9t83SGgTAe7OCI9jC/diDy1nB4Wzu8OspLudADue6AUU/+nznuNrnLtd2c/+1zt23/T7v7u+z+3Ms3+WDbPccwE5FwLt25pHAF4H7RlADTgJeD+wHTJe3vgJ4OfCQ8vdINhuA88vfY3b8BuA04LK5AuDmzZt58/r1TE9Ps27DBs5es6bTr08/5xzOXbuWmZmZwZ8+9HnYx+3btzM1NWWfrfNO7+fWebe/w5vf2X39Pl+3ejXnXXihAbDDATAO/74UOKj0MYLdQ8vs3rDbRwOfKacD7FEOC8f/i1PbhsumMqu4Zq4AuHXr1sF9AF0UUEABBRRQoP0C3gi62zeCfgrwAeDZwCdKOcYpa/cAntsoz0PL+X0HABEAv1FC4rWN97ynXDX+YgNg+3dst1ABBRRQQIGdCRgAuxsAn1HuzBKHfD/UKIJFmQFctWrV/38yxsqVK4k/LgoooIACCijQHoFNmzYRf2KJp3Vt3LgxXnoj6PYMUXpLnl/O33sOcPmstcXFHOfOOgdwdbnoY2fnAH4biLu9zHkOoIeA02PmChRQQAEFFKgm4Axg92YA4+rduLr3mcBn56ikuAo4Du3GVcBxxXCEvo8Ab2xcBRzn+cV64nYycUHIOuAFwCHAnFcBGwCr7bN+kAIKKKCAAmkBA2D3AmBc2RuP+Lu9VEfcqi5ubh+3dBkGwsOAC8p9AOPK4LhVzKtnVdOZwEuAe5Wnuu30PoAGwPS+6AoUUEABBRSoJmAA7F4ArFY85YMG9wE0ANZm9/MUUEABBRQYXcAAaAAcvXp2tDQAZgVtr4ACCiigQGUBA6ABMFtyBsCsoO0VUEABBRSoLGAANABmS84AmBW0vQIKKKCAApUFDIAGwGzJGQCzgrZXQAEFFFCgsoAB0ACYLTkDYFbQ9goooIACClQWMAAaALMlZwDMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJWcAzAraXgEFFFBAgcoCBkADYLbkDIBZQdsroIACCihQWcAAaADMlpwBMCtoewUUUEABBSoLGAANgNmSMwBmBW2vgAIKKKBAZQEDoAEwW3IGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkZALOCtldAAQUUUKCygAHQAJgtuUEA3LJlC8uXL2fbtm0sW7ZssM4uv47+9aGf8/Wxj/23z/3Yt5s175j3Y8z7OM5R5xEAV6xYEd3fB7glGwYmsf3UJG50i7Z5EACPa1TPvsCNZQO7+jq619W+zdcv++yY96HmrXPrvOt1Hv2L5Q7gih0vDYAtClWTtCmDAPguYG/gKuBIYKbjr2OA+tDPZh/ts2Peh5q3zq3zrtd5zHqdUCZtIvkZACcpcrVrWwcB8HnApcCLyra9veOvI+DGTtT1fsZwDvton3cUt2Pe7f3cOrfOu76fn1y+x24tyc8A2K5QNUlbMwiApwLHlpm/K4GjevB6z8YsoH3u7pj3cZyb+28f+2+fu7s/z/e7qW9j/mPgeAPgYCLHZXSBQQB8DFMcxgwXA8ewB4cy3fnXhzA9KB773O0x7+M4N/ffPvbfPvfjO7zPdf549uKj3M59PAQ8evqxJYMAeG9WcARbuB97cDkrOJzNHX49xeUcyOFcNxj+fvT5znG1z12u7eb+a527b/t93t3fZ/fnWL7LB9nuOYAGuZEFBgFw8+bNvHn9eqanp1m3YQNnr1nT6denn3MO565dy8zMzOBPH/o87OP27duZmpqyz9Z5p/dz67zb3+HN7+y+fp+vW72a8y680AA4cvyxoTeCtgYUUEABBRSYMAFvBL30N4L+WrleYlelc9Cu3rBEPzcALhG8H6uAAgoooMCoAgbApQ+AwzunxBgeALwMeCcQwfBA4DeAC4CzRh3kRW5nAFxkYFevgAIKKKDAuAUMgEsfAJtjejnwR8DfNf5nXLV+NvCUcQ/+mNZnABwTpKtRQAEFFFCgloABsF0BMJ7FF1dlDy7LKUvcnugmGFxt28bFANjGUXGbFFBAAQUU2ImAAbBdAfCfgTeVm3QPh+2FQNxn+eEtrWQDYEsHxs1SQAEFFFBgPgEDYLsC4ErgL4GryzmAP12C37OAj7e0jA2ALR0YN0sBBRRQQAED4Pw10LYngcSFH88FHgB8E3h3CYNtrWIDYFtHxu1SQAEFFFBgHgFnANs1AziJhWoAnMRRc5sVUEABBXotYABc+gAYt3m5tFRhnO8333JJSyvVANjSgXGzFFBAAQUU8BBwew8B/wtwWNm8uPffXMsM4I2g3Y8VUEABBRRQYCwCzgAu/QzgWAZyCVfiDOAS4vvRCiiggAIKjCJgADQAjlI3zTYGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkZALOCtldAAQUUUKCygAHQAJgtOQNgVtD2CiiggAIKVBYwABoAsyVnAMwK2l4BBRRQQIHKAgbA9gXAY4GfA+41qxbWVa6NhX6cAXChUr5PAQUUUECBlggYANsVANcDpwP/CPywUSNxH8DjWlIzszfDANjSgXGzFFBAAQUUmE/AANiuAPht4FnA309QyRoAJ2iw3FQFFFBAAQVCwADYrgD4PWAFEDN+k7IYACdlpNxOBRRQQAEFioABsF0B8C3Ap4H3TVCFGgAnaLDcVAUUUEABBZwB3FEDUy0qhXcDvwR8Frhh1nad2KLtbG6KAbClA+NmKaCAAgooMJ+AM4DtCoDv2Emp/lZLy9gA2NKBcbMUUEABBRQwAM5fA22aAZzESjUATuKouc0KKKCAAr0WcAawXTOAw0PSRwIPBL4OXNXyi0IMgL3+CrHzCiiggAKTKGAAbFcA/CngQ8CjgJuA+wBfBE6Y45zA+ert14BVwCOBewJ3A6Ybb47XPwJ+XM5/jCuOjwauabznLOBkIMLdF8r6mj9vfrYBcBL3fLdZAQUUUKDXAgbAdgXA95RQ9jvA94H7AhuBPYAIdgtZnlra3R24eJ4A+GTgk/Os7JXAKcDTgK8CZwBxAcrBwG1ztDEALmRUfI8CCiiggAItEjAAtisAxo2gD4n7MzZqZB/gWuC/7mbdxCPlrpgnAD6l/GyuVV4HbADOLz/cs8w+ngZcZgDczVHw7QoooIACCrRQwADYrgD4HeBngB80aiWeCfzvwP13s352FgDjc+LQcJxj+LYyUxirj9m8m8sh4Ssbn7cJuBpYYwDczVHw7d76PbIAACAASURBVAoooIACCrRQwADYrgD4XuAO4GVlFjBm/2Imbi/gV3ezfuYLgE8CPgdsB+JwcczqvQr4E+ABwDeAh5ZZx+FHxqHpmJV8sQFwN0fBtyuggAIKKNBCAQNguwLgfsBHgIeXcwDjIpB/AZ4BfGs362e+ADh7NeuA44FjMjOAq1atYtmyZYN1r1y5cvDHRQEFFFBAAQXaI7Bp0ybiTyzbtm1j48a4zICYbGqeetaeDV7kLWnbfQDjgo+fAw4Argf+YdZVvAvl2J0AGGntcWXFc50DGOcmnuo5gAul930KKKCAAgq0W8AZwPbMAMb0Wdzy5THlNi2jVk4EyDi/LwLgx4A4hzAO924DDi9XGcf5fHE7mLgaOB4/F7OAg38GlPP84irgpwMRBuNnLygXp3gV8KijYjsFFFBAAQVaJGAAbE8AjLKIw7wHlrA2apmcBMQj5eL+frHEDGe8jnP/4iKPc8u5fnEfwLgI5ALgolkfdibwkhIe40bUcV9B7wM46ojYTgEFFFBAgZYJGADbFQDXA7cDr2lZnexsc7wP4AQNlpuqgAIKKKBACBgA2xUA/xZ4LLC5zM41n+DxhJaWrAGwpQPjZimggAIKKDCfgAGwXQEwnrox3xKPZ2vjYgBs46i4TQoooIACCuxEwADYrgA4icVqAJzEUXObFVBAAQV6LWAAbFcAjCt451uah4PbVLQGwDaNhtuigAIKKKDAAgQMgO0KgBHyhlfvzh6+eCZvGxcDYBtHxW1SQAEFFFBgJwIGwHYFwLh3X3PZH/i9cpuWeGZvGxcDYBtHxW1SQAEFFFDAALjTGmjbk0Bmb+xBwDvLo9raWMwGwDaOitukgAIKKKCAAXCiA2Ac+r2p3MS5jcVsAGzjqLhNCiiggAIKGAAnJgDGbF9zuQdwcnmKxyNaWskGwJYOjJulgAIKKKDAfAKeA9iucwBnXwQSh6fjebwnAp9raRkbAFs6MG6WAgoooIACBsD5a6BN5wA+aNZm3gp8v+XlawBs+QC5eQoooIACCswWcAawXTOAk1ihBsBJHDW3WQEFFFCg1wIGwHYFwLgR9B8ALwR+EtgHWAk8sNwKpo3FagBs46i4TQoooIACCuxEwADYrgC4HngmcB5wAXBv4MHA+4BHtbSSDYAtHRg3SwEFFFBAgfkEDIDtCoBfA54AXF/O/bsvELOC3wPidRsXA2AbR8VtUkABBRRQwBnAndZAmy4CiaC3ojwOLi7+iNB3N+Bb5ZBwG4vZANjGUXGbFFBAAQUUMABOTAD8ZDn0G4d8hwHwV4DfLucCtrGYDYBtHBW3SQEFFFBAAQPgxATAxwKXAx8CngW8C4gA+FTg8y2tZANgSwfGzVJAAQUUUGA+Ac8BbNc5gDFOhwAvAx4CfKfMCF7V4hI2ALZ4cNw0BRRQQAEF5hIwALYvAE5apRoAJ23E3F4FFFBAgd4LGADbFwDjecCPBu41qzovaWm1GgBbOjBulgIKKKCAAh4Cnr8G2nQV8EuB88sFID9sbPIMEMGwjYsBsI2j4jYpoIACCiiwEwFnANs1A/h14FTgLyaoag2AEzRYbqoCCiiggAIhYABsVwC8uTz9Y5Kq0wA4SaPltiqggAIKKGAAHNRAmw4BXwZcCHxqgqrTADhBg+WmKqCAAgoo4AzgjhpoUwD8Y+Ak4H8BN8wq0XUtLVkDYEsHxs1SQAEFFFBgPgEPAbcrAMaTQOZa4iKQ41paxgbAlg6Mm6WAAgoooIABcP4aaNMM4CRWqgFwEkfNbVZAAQUU6LWAM4DtmgGcxGI0AE7iqLnNCiiggAK9FjAAGgCzO4ABMCtoewUUUEABBSoLGAANgNmSMwBmBW2vgAIKKKBAZQEDoAEwW3IGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkZALOCtldAAQUUUKCygAHQAJgtOQNgVtD2CiiggAIKVBYwABoAsyVnAMwK2l4BBRRQQIHKAgZAA2C25AyAWUHbK6CAAgooUFnAAGgAzJacATAraHsFFFBAAQUqCxgADYDZkjMAZgVtr4ACCiigQGUBA6ABMFtyBsCsoO0VUEABBRSoLGAANABmS84AmBW0vQIKKKCAApUFDIAGwGzJGQCzgrZXQAEFFFCgsoAB0ACYLTkDYFbQ9goooIACClQWMAAaALMlZwDMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJWcAzAraXgEFFFBAgcoCBkADYLbkDIBZQdsroIACCihQWcAAaADMlpwBMCtoewUUUEABBSoLGAC7FwB/DVgFPBK4J3A3YLpRV48A3gI8GrgZuAg4a1bdxd9PBiLcfaGs75p5atMAWHmn9eMUUEABBRTIChgAuxcAnwrcF7g7cPGsABiB8MvAJcB64GDgY8AbgDeVYnolcArwNOCrwBnAieW9t81RcAbA7F5oewUUUEABBSoLGAC7FwCHJXQscMWsAHgS8Hpgv8as4CuAlwMPKQ2vAzYA55e/7wncAJwGXGYArLyH+nEKKKCAAgosgoABsF8BMILdQ8vs3rCcjgY+A+wD7FEOC8f/u7JRb5uAq4E1BsBF2AtdpQIKKKCAApUFDID9CoBxSPgewHMbdXYoEOf3HVAC4DdKSLy28Z73ALcALzYAVt5D/TgFFFBAAQUWQcAA2K8AuGgzgKtWrWLZsmWDEl25cuXgj4sCCiiggAIKtEdg06ZNxJ9Ytm3bxsaNG+NlHAGMSZ7eLVMd7fFc5wDGxRznzjoHcHW56GNn5wB+GzjVcwA7Wil2SwEFFFCgdwLOAHZvBjDO44tbv0QAjCt87wVsj7BfDv/God24Cvg15cKPjwBvbFwFHOf5xVXATwfigpB1wAuAQwCvAu7dV4QdVkABBRToooABsHsBMK70fQcwUwo2Zjjj9ZOATwOHAReU+wBuBd4KvHpWcZ8JvKSEx6u8D2AXd337pIACCijQZwEDYPcCYO16HtwHcMuWLSxfvnxwTsHwXMAuvw7kPvRzvj72sf/2uR/7drPmHfN+jHkfxznqPALgihUrovueA1g7OXXk8wYB8LjGGaT7AjeWznX1dXSvq32br1/22THvQ81b59Z51+s8+hfLHeVmwQbAjqSxJejGIAC+C9gbiOPFR5Zjzl1+Hc596Gezj/bZMe9DzVvn1nnX6zzOCzuhTNrE1J8BcAmSU0c+chAAnwdcCryodOrtHX8dJ1XGTtT1fsZwDvton3cUt2Pe7f3cOrfOu76fn1y+x24tyc8A2JE0tgTdGATAuEdMXHYcX57xCJGjevA6npE3/Jeife7umPdxnJv7bx/7b5+7uz/P97upb2P+Y+B4A+BgIsdldIFBAIxnx8V9YuJRI3G58YN78PqgMgton7s95n0c5+b+28f+2+d+fIf3uc6fDHwQuI+HgEdPP7ZkEADjLtIPLHeY/jywf09ex/jvB9jnbo95H8e5Wdd97L999ruty9/tcY7+EUCcv+8hYJPcqAKDALic+/IAvs/dga9wP/bnex1//UD2Jx6bTMf7+Z/H0j53vbabY26d9+P7bMeY93Hf3vE7q491/pMcxHf5kgFw1OxjO9gxA7h582bevH4909PTrNuwgbPXrOn069PPOYdz165lZmZm8KcPfR72cfv27UxNTdln67zT+7l13u3v8OZ3dl+/z9etXs15F17oDKBRbmSBQQDcunXr4EbQLgoooIACCijQfgGfBOKTQLJVagDMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJWcAzAraXgEFFFBAgcoCBkADYLbkDIBZQdsroIACCihQWcAAaADMlpwBMCtoewUUUEABBSoLGAANgNmSMwBmBW2vgAIKKKBAZQEDoAEwW3IGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkZALOCtldAAQUUUKCygAHQAJgtOQNgVtD2CiiggAIKVBYwABoAsyVnAMwK2l4BBRRQQIHKAgZAA2C25AyAWUHbK6CAAgooUFnAAGgAzJacATAraHsFFFBAAQUqCxgADYDZkjMAZgVtr4ACCiigQGUBA6ABMFtyBsCsoO0VUEABBRSoLGAANABmS84AmBW0vQIKKKCAApUFDIAGwGzJGQCzgrZXQAEFFFCgsoAB0ACYLTkDYFbQ9goooIACClQWMAAaALMlZwDMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJWcAzAraXgEFFFBAgcoCBkADYLbkDIBZQdsroIACCihQWcAAaADMlpwBMCtoewUUUEABBSoLGAANgNmSMwBmBW2vgAIKKKBAZQEDoAEwW3IGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkZALOCtldAAQUUUKCygAHQAJgtOQNgVtD2CiiggAIKVBYwABoAsyVnAMwK2l4BBRRQQIHKAgZAA2C25AyAWUHbK6CAAgooUFnAAGgAzJacATAraHsFFFBAAQUqCxgADYDZkjMAZgVtr4ACCiigQGUBA6ABMFtyBsCsoO0VUEABBRSoLGAA7GcAPAP4I+A2dvR/Bvgw8PxSf48A3gI8GrgZuAg4a57aNABW3mn9OAUUUEABBbICBsD+BsAnA0+Yo4DuCXwZuARYDxwMfAx4A/CmOd5vAMzuhbZXQAEFFFCgsoAB0AA4u+ROAl4P7AdMlx++Ang58BADYOU91I9TQAEFFFBgEQQMgP0NgGvKIeA4DPw54A+B/wA2AA8Fntaot6OBzwD7AD+YVYfOAC7CjukqFVBAAQUUWEwBA2A/A+DDgFuB64GfAs4DHgs8EngzcA/guY3COxS4BjgAuMEAuJi7pOtWQAEFFFBg8QUMgP0MgLMraxmwFXgm8AujzACuWrWKZctiNbBy5crBHxcFFFBAAQUUaI/Apk2biD+xbNu2jY0bN8bLOLp3S3u2st6WxFWwfV+GAfAXy4zgubPOAVwNnOI5gH0vE/uvgAIKKNAVAWcA+zkD+BzgCuBG4P7lEPAxwMPLbWGuLVcBv6aEvo8Ab/Qq4K7s9vZDAQUUUKDvAgbAfgbAvwSOKuf63QR8utwX8LqyQxwGXFDuAxiHht8KvHqencWLQPr+LWL/FVBAAQUmTsAA2M8AOM5CNQCOU9N1KaCAAgooUEHAAGgAzJaZATAraHsFFFBAAQUqCxgADYDZkjMAZgVtr4ACCiigQGUBA6ABMFtyBsCsoO0VUEABBRSoLGAANABmS84AmBW0vQIKKKCAApUFDIAGwGzJGQCzgrZXQAEFFFCgsoAB0ACYLTkDYFbQ9goooIACClQWMAAaALMlZwDMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJTcIgFu2bGH58uVs27aNZcuWDdbZ5dfRvz70c74+9rH/9rkf+3az5h3zfox5H8c56jwC4IoVK6L7+wC3ZMPAJLafmsSNbtE2DwLgcY3q2Re4sWxgV19H97rat/n6ZZ8d8z7UvHVunXe9zqN/sdwBXLHjpQGwRaFqkjZlEADfBewNXAUcCcx0/HUMUB/62eyjfXbM+1Dz1rl13vU6j1mvE8qkTSQ/A+AkRa52besgAD4PuBR4Udm2t3f8dQTc2Im63s8YzmEf7fOO4nbMu72fW+fWedf385PL99itJfkZANsVqiZpawYB8FTg2DLzdyVwVA9e79mYBbTP3R3zPo5zc//tY//tc3f35/l+N/VtzH8MHG8AHEzkuIwuMAiAj2GKw5jhYuAY9uBQpjv/+hCmB8Vjn7s95n0c5+b+28f+2+d+fIf3uc4fz158lNu5j4eAR08/tmQQAO/NCo5gC/djDy5nBYezucOvp7icAzmc6wbD348+3zmu9rnLtd3cf61z922/z7v7++z+HMt3+SDbPQfQIDeywCAAbt68mTevX8/09DTrNmzg7DVrOv369HPO4dy1a5mZmRn86UOfh33cvn07U1NT9tk67/R+bp13+zu8+Z3d1+/zdatXc96FFxoAR44/NvRG0NaAAgoooIACEybgjaC9EXS2ZA2AWUHbK6CAAgooUFnAAGgAzJacTwLpydNPfEJC/57+El8OPvGm+082ct/uZ537JBAD4FgCoE8C8eknXXz6S+wcXX8qwM762Mf+2+f+1XzfxtwngdwZe7wNTC4C+iSQnjwVJMqk63fIn/0kmz722ae/WOfu591+mpVPAjEA5mLfna19Ekix6PoTInxCwo6B7vo4N/vomPdjzPs4zn2uc58EYgAcawD0SSD9uHN+3+6WH0946WOffRLInbOAPuXH77YuPtnKJ4HsiEAeAs5FQZ8E4tNPOv30lz4+FaLPT0iIJxn1ccz72Oc+17lPAjEA5qLfjtY+CaQXTz/xSSB9eyrE/fBJIH0b8/hC71uf+1vnPgnEGcB8BPRJID4JpPNPBenjUyHiqT59fUJCPMmoj2Pexz73uc59EoiHgLMR0BtBZwVtr4ACCiigQGUBbwRtAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJWcAzAraXgEFFFBAgcoCBkADYLbkDIBZQdsroIACCihQWcAAaADMlpwBMCtoewUUUEABBSoLGAANgNmSMwBmBW2vgAIKKKBAZQEDoAEwW3IGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkZALOCtldAAQUUUKCygAHQAJgtOQNgVtD2CiiggAIKVBYwABoAsyVnAMwK2l4BBRRQQIHKAgZAA2C25AyAWUHbK6CAAgooUFnAAGgAzJacATAraHsFFFBAAQUqCxgADYDZkjMAZgVtr4ACCiigQGUBA6ABMFtyBsCsoO0VUEABBRSoLGAANABmS84AmBW0vQIKKKCAApUFDIAGwGzJGQCzgrZXQAEFFFCgsoAB0AC4s5I7CzgZiJD3BWAVcM2sBgbAyjutH6eAAgoooEBWwABoAJyvhl4JnAI8DfgqcAZwInAwcFujkQEwuxeOsf2mTZtYuXLlGNfoqkYVcCxGlVucdo7H4riOslbHYhS18bcxABoA56uq64ANwPnlDXsCNwCnAZcZAMe/M45jjaeddhobNsSwuSy1gGOx1CNw1893PNozHo5FO8bCAGgAnKsSY1bvZuBo4MrGGzYBVwNrDIDt2IFnb4VfrO0ZF8eiPWMRW+J4tGc8HIt2jIUB0AA4VyU+APgG8FDg2sYb3gPcArx4dgC8/vrrWb48cqPLUgqsXbuW1772tUu5CX52EXAs2lUKjkd7xsOxaMdYRAA84IADYmP2Kb/b27FhFbdiquJnTcpH7c4M4P7ANyelY26nAgoooIACCtxFICZ9vtVHEwPg3KM+1zmA3wZOnXUOYPjtB9zax+KxzwoooIACCkywwL3K+f0zE9yHkTfdADg3XZznF1cBPx2IMLgOeAFwyKyrgEeGt6ECCiiggAIKKLBUAgbA+eXPBF4CxL8QrprnPoBLNW5+rgIKKKCAAgooMLKAAXBkOhsqoIACCiiggAKTKWAAzI3bQp4WkvuEfrd+XTkM/yDgB8CngN+fdeHNI4C3AI8ut++5CIhxaS6O0/jr6APACcBTgCvK6h2L8Tvvao1xu6qzgSOB7eVpRcc4HrtiG/vP71PuHfvzwF7AvwKvAj7tWIzdurnCXytH5x4J3BO4GzDdeMM4vpMWso5F7eRirdwAOLrsQp8WMvon2PI1wPvL/RfvDrwVeBhwRKGJHf7LwCXA+vKklo8BbwDeVN7jOI2/juKpOM8Dnlr+RAB0LMbvvKs1Rvj7aDlfOfaTO8o/hD7veOyKbuw//zMg7inyLOCmcsFgnEb0QODHfk+N3Xu4wvgOui8Qvx8unhUAx/GdtJB1LFrnFnvFBsDRhRf6tJDRP8GWswXiX3lfLDv8VuAk4PXlSuzhv/peAbwceEhp7DiNt47ilgmfAWKWKe6XOZwBdCzG67yQtcXsUtysPv6RM3txPBYiOL73xEMCLixHI2Kt9yh3h3hsuafsuX5PjQ97jjUdW45ENGcAx7EPLGQdi9qxxVy5AXA03d25V+Bon2CruQTi8O9LgYPKD+O5b3HD7nhm83CJWZEIKHFzzz1246kuii9MIJ6I8+fA28uhlmEAdCwW5jeud+1dAsYbgScCDwa+BsRpE39RDke6b4xLe9fridNMjgOeA9xYHhv6QiAOH8aYOBa7Nsy8Y64AOI7vpF2tI05NmtjFADja0O3O00JG+wRbzRaIoBHnnT0b+ET5YUz5x7+0n9t486HlPKg4HBMBcKFPdVF81wIvK+f9rSxvjVnXYQB0LHbtN853xE3orwc2l/Nk/7GMTTyxKH4Znuy+MU7uXa7rv5Tvp9g34pDv98vh4L8rhyb9ntolYeoNcwXAcXwn7WodN6S2eokbGwBHGwBnAEdzG7XVM4BLyyHfDzVWsqt/nTkDOKr4f24Xs64xsxqHtCJ4xNIMgI7F+KwXsqbhd9A5wNpGg48DXyoXIjjrtBDJ8bznb8t5fqeVmdn4znon8ATgN50BHA/yTtbiDOAIxAbAEdBKk4U+LWT0T7BlCDwfOL8cWrl8FklcjDD73JrV5aT4nZ0DONdTXdTeuUCcC/Mn5ZmZw++NfYE4F/O9wOeA84CfalyF51gsblV9BXjfPAEwrkJ131hc/+HaYz/YUi5O+6fGR34BiBnZmKV1LBZ3LOYKgJnfD78LvAtYyDoWt2eLuHYD4Oi4Pi1kdLuFtoynscTVvc8EPjtHo7hC69pyFXBcMRyh7yNAnBc1vArYcVqo9s7fF4e44mq75hLPwY7bMMQh+Tjs5ViMx3qha4kLnk4v58D+c9lPInDErNO/OR4LZRzL+2JCIP6B+nvlllXxFKkI578AxFXZ7htjYf5PK4mjPHHhRwTAuANEPLghboe0rZwCkXVfyO+YxelZhbUaAHPIPi0k57er1nGIMW5tcXt5Y9RrPLMxLvoYBsLDgAvK7S9iNipuFfPqWSt2nHYlPdrP44s2bsMwvA+gYzGaY6bVH5T7oMVFTzEjGLX+V2WFjkdGdvfaxu2pYpYvTpGIQBKnSfxxuVgq1uRY7J7nQt8dRybeUX4vRJvh74gnlXswjsN9IetY6Pa26n0GwFYNhxujgAIKKKCAAgosvoABcPGN/QQFFFBAAQUUUKBVAgbAVg2HG6OAAgoooIACCiy+gAFw8Y39BAUUUEABBRRQoFUCBsBWDYcbo4ACCiiggAIKLL6AAXDxjf0EBRRQQAEFFFCgVQIGwFYNhxujgAIKKKCAAgosvoABcPGN/QQFFFBAAQUUUKBVAgbAVg2HG6OAAgoooIACCiy+gAFw8Y39BAUUqCfwqvJYqJ+v95F+kgIKKDB5AgbAyRszt1iBvgrc2njk017AnsBtjcc/NR8RWMPoZOCXgGfU+DA/QwEFFBingAFwnJquSwEFagnE854fBxxX6wPn+Jx4+Px7gT9dwm3woxVQQIGRBAyAI7HZSAEFllhgvgB4BvAU4PFl+z4JXA38JBAzhDcDrwS+DLwV+FngGuBE4CulzR7AqcCLgP3K//8D4IpGn/cBvg3sD9wEfK0EwaNLMP0W8NIyS7kB+Gng74HfAL5b1nMK8LvAijKTGYHyhUvs6scroEBPBAyAPRlou6lAxwR2FgCfDDyhEQAfCTwT+CzwCuC1wN8CLy4h7j3APUpAjGZnAr8I/Crw78AJwGXAw0vQi/e8APjNEjbj7xEAY4nPiUB5HvD8EvoiSN4O/DXwz8DvAD8D/BPwGOBfgbsDjwI+07FxsjsKKNBSAQNgSwfGzVJAgZ0K7E4AjBD322Vty8ss4K8Df17+37OBi4B9y99jlvA5wCcaWxDh7W9KeIz//YES6GIWcRgAYx0RLmN5BPAlIGYE/6H8v9NKKHx0mRGMoBgh8uNAnN/oooACClQTMABWo/aDFFBgjAK7EwBjtm9d+ey4cOQO4InAp8v/Wwl8GFhWDhV/B7gFmC4/j+/JnwDeCawqs3VxGDdm8eK9wwAY23RJ+fuDy2HmA4FvlP/3EuD3gIPL32O28GXAUeW9/wOI2UgXBRRQYNEFDICLTuwHKKDAIggsVgCMEBjn9EUonO9w7K+Uc/eOafQrDgHvbgAcNo9zDn+5hL9DymHnRSBzlQooVLFmHgAAAcVJREFUoMCdAgZAq0EBBSZRYLECYFi8sRy6jQsy/g3Yu5yrFxd9xOHkd5fDujFjN1zmCoBxUUlc/DHXDGDMAsYsYcxC/rAEzo8AD2mcTziJ4+I2K6DAhAgYACdkoNxMBRS4i8BCA2BcuRszec1DwNuAJ81zCDg+JL4X4wrduEjkAcCPgC8Ca4CvAlvKOX5fb2zRdcDZu3EI+DDgbeUq5Pi8CImvK+HSoVZAAQUWXcAAuOjEfoACCnRIIG76HFcJx9W7LgoooMDEChgAJ3bo3HAFFFgCgbjFTFxIElcFuyiggAITK2AAnNihc8MVUEABBRRQQIHRBAyAo7nZSgEFFFBAAQUUmFgBA+DEDp0broACCiiggAIKjCZgABzNzVYKKKCAAgoooMDEChgAJ3bo3HAFFFBAAQUUUGA0AQPgaG62UkABBRRQQAEFJlbAADixQ+eGK6CAAgoooIACowkYAEdzs5UCCiiggAIKKDCxAgbAiR06N1wBBRRQQAEFFBhNwAA4mputFFBAAQUUUECBiRUwAE7s0LnhCiiggAIKKKDAaAIGwNHcbKWAAgoooIACCkysgAFwYofODVdAAQUUUEABBUYT+H/ajTonW1zQLwAAAABJRU5ErkJggg==\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "\n",
    "# If there are spikes, plot using matplotlib\n",
    "if len(spikes_input) != 0 or len(spikes_output) != 0:\n",
    "    plt.figure()\n",
    "    if len(spikes_input) != 0:\n",
    "        plt.plot([i[1] for i in spikes_input],\n",
    "                   [i[0] for i in spikes_input], \"bo\", markersize=5, label=\"input\")\n",
    "    if len(spikes_output) != 0:\n",
    "        plt.plot([i[1] for i in spikes_output],\n",
    "                   [i[0] for i in spikes_output], \"rv\", label=\"output\")\n",
    "    plt.ylabel('neuron id')\n",
    "    plt.xlabel('Time/ms')\n",
    "    plt.title('spikes')\n",
    "#     plt.xlim(450,650)\n",
    "    plt.show()\n",
    "else:\n",
    "    print \"No spikes received\"\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.11"
  },
  "widgets": {
   "state": {},
   "version": "1.1.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}