Spike_IO-BEE-Main.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Liquid State Machine with SpiNNaker\n",
    "\n",
    "- I've been tweaking the .spynnaker.cfg file to slow down SpiNNaker by increasing the timeScaleFactor to 10. This way it will have more resources to deal with the input/output spikes... I hope!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# plotter in python\n",
    "%matplotlib notebook\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "import time\n",
    "import numpy\n",
    "\n",
    "import save_load_file as slf"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "import sys\n",
    "import BEE # see https://github.com/ricardodeazambuja/LiquidStateMachine-Python"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "net_shape = numpy.array([20,5,6])\n",
    "\n",
    "membrane_rand=[13.5E-3,14.9E-3]\n",
    "current_rand=[13.5E-09, 14.5E-09]\n",
    "\n",
    "noisy_current_rand=1E-18 #1E-9 # NO NOISY CURRENT TO MAKE IT COMPATIBLE WITH SPINNAKER\n",
    "\n",
    "factor=10E6 #2*10.0 # VIRTUALLY NO NOISY RESET VOLTAGES TO MAKE IT COMPATIBLE WITH SPINNAKER\n",
    "\n",
    "SpkLiq_vresets_rand=[13.8E-3,13.8E-3*(1+1/factor)] # the membrane reset values are drawn from a uniform distribution\n",
    "\n",
    "lbd_value = 1.2\n",
    "\n",
    "parameter_list=\\\n",
    "[[[0.1  ,  0.32 ,  0.144,  0.06 ,  -47.  ,  0.8],\\\n",
    "  [0.4  ,  0.25 ,  0.7  ,  0.02 ,  -47  ,  0.8]],\\\n",
    " [[0.2  ,  0.05 ,  0.125,  1.2  ,  150  ,  0.8],\\\n",
    "  [0.3  ,  0.5  ,  1.1  ,  0.05 ,  70  ,  1.5]]]\n",
    "\n",
    "\n",
    "# Neuron parameters\n",
    "taum=0.03 # membrane time constant\n",
    "cm=3e-08 # membrane capacitance\n",
    "taue=0.003 # excitatory synapse time constant\n",
    "taui=0.006 # inhibitory synapse time constant\n",
    "\n",
    "\n",
    "time_step = 2/1000.\n",
    "\n",
    "number_of_threads = 5"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Liquid 0: [8504 6373 6528 9820 4728 5408]\n",
      "Liquid 1: [6878 1704 9127 4863 8948 3511]\n",
      "Liquid 2: [ 214 9852 3943 8654 8283 5606]\n",
      "Liquid 3: [4730  341 4093 7541  811  467]\n",
      "Liquid 4: [3202  913 4191 4309 9239 8001]\n"
     ]
    }
   ],
   "source": [
    "# \n",
    "# Defines how many liquids\n",
    "# \n",
    "number_of_liquids = 5\n",
    "\n",
    "# \n",
    "# Generates the seeds used to generate the liquids\n",
    "# \n",
    "# rng_seed=123\n",
    "# rng = numpy.random.RandomState(rng_seed)\n",
    "# liqseeds = []\n",
    "# for i in range(number_of_liquids):\n",
    "#     liqseeds.append(numpy.array(rng.randint(0,10000,6),dtype=numpy.uint32))\n",
    "#     print \"Liquid %d:\" % i, liqseeds[i]\n",
    "\n",
    "liqseeds = slf.load_from_file_gz(\"liqseeds.gzpickle\")\n",
    "number_of_liquids = len(liqseeds)\n",
    "for i in range(number_of_liquids):\n",
    "    print \"Liquid %d:\" % i, liqseeds[i]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Processing liquid number: 0\n",
      "initializing BEE...\n",
      "BEE_setup done!\n",
      "Liquid->Liquid connections and parameters\n",
      "Processing liquid number: 1\n",
      "initializing BEE...\n",
      "BEE_setup done!\n",
      "Liquid->Liquid connections and parameters\n",
      "Processing liquid number: 2\n",
      "initializing BEE...\n",
      "BEE_setup done!\n",
      "Liquid->Liquid connections and parameters\n",
      "Processing liquid number: 3\n",
      "initializing BEE...\n",
      "BEE_setup done!\n",
      "Liquid->Liquid connections and parameters\n",
      "Processing liquid number: 4\n",
      "initializing BEE...\n",
      "BEE_setup done!\n",
      "Liquid->Liquid connections and parameters\n",
      "\n",
      "Total number of neurons: 3000\n",
      "Total number of neurons LSM_EXC population: 2400\n",
      "Total number of neurons LSM_INH population: 600\n",
      "Total number of EXC2ANY connections: 4537\n",
      "Total number of INH2ANY connections: 1284\n"
     ]
    }
   ],
   "source": [
    "pop_lsm_exc = []\n",
    "pop_lsm_inh = []\n",
    "\n",
    "input_synapses_exc = []\n",
    "input_synapses_inh = []\n",
    "\n",
    "exc_neuron_idx = []\n",
    "\n",
    "inh_neuron_idx = []\n",
    "\n",
    "total_size_exc = 0\n",
    "total_size_inh = 0\n",
    "\n",
    "i_offset_exc = []\n",
    "i_offset_inh = []\n",
    "\n",
    "init_membrane_v_exc = []\n",
    "init_membrane_v_inh = []\n",
    "\n",
    "connections_pre_e = []\n",
    "connections_pos_e = []\n",
    "weights_pre_e = []\n",
    "\n",
    "connections_pre_i = []\n",
    "connections_pos_i = []\n",
    "weights_pre_i = []\n",
    "\n",
    "total_exc2_conn = 0\n",
    "total_inh2_conn = 0\n",
    "\n",
    "\n",
    "for lsm_i in range(number_of_liquids):\n",
    "    \n",
    "    print \"Processing liquid number:\", lsm_i\n",
    "\n",
    "    if BEE.BEE_initialized():\n",
    "        # Resets the simulator\n",
    "        BEE.BEE_free()\n",
    "\n",
    "    if not BEE.BEE_initialized():\n",
    "        print \"initializing BEE...\"\n",
    "        BEE.initialize_sim(my_net_shape = net_shape, \\\n",
    "                       my_lbd_value = lbd_value, \\\n",
    "                       my_seeds=liqseeds[lsm_i], \\\n",
    "                       SpkLiq_step=time_step,\\\n",
    "                       SpkLiq_threads_N = number_of_threads,\\\n",
    "                       SpkLiq_vresets = SpkLiq_vresets_rand,\\\n",
    "                       SpkLiq_membrane_rand=membrane_rand,\\\n",
    "                       SpkLiq_current_rand=current_rand,\\\n",
    "                       SpkLiq_taum=taum, SpkLiq_cm=cm, SpkLiq_taue=taue, SpkLiq_taui=taui, \\\n",
    "                       SpkLiq_noisy_current_rand=noisy_current_rand)\n",
    "\n",
    "\n",
    "    if parameter_list!=[]:\n",
    "        BEE.change_parameters(parameter_list)\n",
    "\n",
    "    if BEE.BEE_initialized() and (not BEE.BEE_connected()):\n",
    "        BEE.generate_connections()\n",
    "\n",
    "    if BEE.BEE_initialized() and (not BEE.BEE_connected()):\n",
    "        BEE.process_connections()\n",
    "\n",
    "    NofN,NoINHN,NoEXCN,NoINHC,NoEXCC=BEE.output_stats(stats=0)\n",
    "    # NofN: total number of neurons\n",
    "    # NoINHN: total number of inhibitory neurons\n",
    "    # NoEXCN: total number of excitatory neurons\n",
    "    # NoINHC: total number of inhibitory connections\n",
    "    # NoEXCC: total number of excitatory connections\n",
    "    \n",
    "    print \"Liquid->Liquid connections and parameters\"\n",
    "    \n",
    "    #\n",
    "    # The populations (exc and inh) are ordered as the arrays (exc_neuron_idx and inh_neuron_idx) below:\n",
    "    #\n",
    "    \n",
    "    exc_neuron_idx.append(BEE.output_exc_indices(NoEXCN))\n",
    "    # exc_neuron_idx[spinnaker_index]=>BEE_index\n",
    "    total_size_exc += len(exc_neuron_idx[-1]) # [-1] because I'm appending each time the loop runs...\n",
    "    # SpiNNaker will have two populations: one with excitatory and another with inhibitory neurons.\n",
    "    # Since I have more than one liquid, I need to add to obtain the total number of neurons.\n",
    "\n",
    "    inh_neuron_idx.append(BEE.output_inh_indices(NoINHN))\n",
    "    # inh_neuron_idx[spinnaker_index]=>BEE_index\n",
    "    total_size_inh += len(inh_neuron_idx[-1])\n",
    "\n",
    "    i_offset = BEE.output_noisy_offsets(NofN)*1E9 # SpiNNaker uses nA\n",
    "    i_offset_exc.append(i_offset[exc_neuron_idx[-1]])\n",
    "    i_offset_inh.append(i_offset[inh_neuron_idx[-1]])\n",
    "\n",
    "\n",
    "    init_membrane_v = BEE.output_initial_voltages(NofN)*1E3 # SpiNNaker uses mV\n",
    "    init_membrane_v_exc.append(init_membrane_v[exc_neuron_idx[-1]])\n",
    "    init_membrane_v_inh.append(init_membrane_v[inh_neuron_idx[-1]])\n",
    "    \n",
    "    connections_pre_e.append(BEE.output_pre_e_connections(NoEXCC))\n",
    "    connections_pos_e.append(BEE.output_pos_e_connections(NoEXCC))\n",
    "    weights_pre_e.append(BEE.output_pre_e_weights(NoEXCC))    \n",
    "    total_exc2_conn+=len(connections_pre_e[-1])\n",
    "    \n",
    "    connections_pre_i.append(BEE.output_pre_i_connections(NoINHC))\n",
    "    connections_pos_i.append(BEE.output_pos_i_connections(NoINHC))\n",
    "    weights_pre_i.append(BEE.output_pre_i_weights(NoINHC))\n",
    "    total_inh2_conn+=len(connections_pre_i[-1])\n",
    "    \n",
    "    # Resets the simulator\n",
    "    BEE.BEE_free()\n",
    "\n",
    "total_size=total_size_exc+total_size_inh\n",
    "\n",
    "print\n",
    "print \"Total number of neurons:\",  total_size\n",
    "print \"Total number of neurons LSM_EXC population:\",  total_size_exc\n",
    "print \"Total number of neurons LSM_INH population:\",  total_size_inh\n",
    "\n",
    "print \"Total number of EXC2ANY connections:\",  total_exc2_conn\n",
    "print \"Total number of INH2ANY connections:\",  total_inh2_conn"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[0, 14]"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#\n",
    "# Converts from the SpiNNaker internal populations indices to BEE individual liquids/indices.\n",
    "#\n",
    "# SpiNNaker will always have two only two populations for all liquids: \n",
    "# - one excitatory and one inhibitory.\n",
    "from_spinnaker_exc2bee = [[lsm_i, spnkr_exc_idx] for lsm_i,exc_lsm in enumerate(exc_neuron_idx) \n",
    "                                                      for spnkr_exc_idx in exc_lsm]\n",
    "from_spinnaker_inh2bee = [[lsm_i, spnkr_inh_idx] for lsm_i,inh_lsm in enumerate(inh_neuron_idx) \n",
    "                                                      for spnkr_inh_idx in inh_lsm]\n",
    "\n",
    "# Example:\n",
    "# SpiNNaker excitatory neuron with index 10 => [BEE_lsm_i,BEE_index]\n",
    "from_spinnaker_exc2bee[10]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['inh', 81]"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#\n",
    "# Converts from lsm_i/BEE to SpiNNaker internal populations\n",
    "#\n",
    "\n",
    "# from_bee2spinnaker[lsm_i][bee_idx]=>['exc/inh',idx]\n",
    "from_bee2spinnaker = [range(NofN) for i in range(number_of_liquids)]\n",
    "\n",
    "for i,(lsm_i,idx) in enumerate(from_spinnaker_exc2bee):\n",
    "    from_bee2spinnaker[lsm_i][idx]=['exc',i]\n",
    "\n",
    "for i,(lsm_i,idx) in enumerate(from_spinnaker_inh2bee):\n",
    "    from_bee2spinnaker[lsm_i][idx]=['inh',i]\n",
    "    \n",
    "# Example:\n",
    "# BEE liquid number 0 and neuron index 440 => ['exc' or 'inh', neuron_index]\n",
    "from_bee2spinnaker[0][440]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "# Saves the arrays to be used by the receiver/sender\n",
    "numpy.save('from_spinnaker_exc2bee.npy',from_spinnaker_exc2bee)\n",
    "numpy.save('from_spinnaker_inh2bee.npy',from_spinnaker_inh2bee)\n",
    "\n",
    "# Saves the arrays to be used by the receiver/sender\n",
    "numpy.save('from_bee2spinnaker.npy',from_bee2spinnaker)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "#\n",
    "# These are the cell (neuron) parameters\n",
    "#\n",
    "cell_params_lsm = {'cm'          : 30,    # Capacitance of the membrane\n",
    "                   'i_offset'    : 0.0,   # Offset current\n",
    "                   'tau_m'       : 30.0,  # Membrane time constant\n",
    "                   'tau_refrac_E': 3.0,   # Duration of refractory period - EXCITATORY\n",
    "                   'tau_refrac_I': 2.0,   # Duration of refractory period - INHIBITORY\n",
    "                   'tau_syn_E'   : 3.0,   # Decay time of excitatory synaptic current\n",
    "                   'tau_syn_I'   : 6.0,   # Decay time of inhibitory synaptic current\n",
    "                   'v_reset'     : 13.5,  # Reset potential after a spike\n",
    "                   'v_rest'      : 0.0,   # Resting membrane potential\n",
    "                   'v_thresh'    : 15.0,  # Spike threshold\n",
    "                   'i_noise'     : 0.2    # mean 0 and SD=0.2nA\n",
    "                }\n",
    "\n",
    "\n",
    "# Parameter of the neuron model LIF with exponential currents \n",
    "# (SpiNNaker does NOT accept to individualy vary the refractory period)\n",
    "cell_params_lif_exc = {\n",
    "                         'tau_m': cell_params_lsm['tau_m'],\n",
    "                         'cm': cell_params_lsm['cm'],\n",
    "                         'v_rest': cell_params_lsm['v_rest'],\n",
    "                         'v_reset': cell_params_lsm['v_reset'],\n",
    "                         'v_thresh': cell_params_lsm['v_thresh'],\n",
    "                         'tau_syn_E': cell_params_lsm['tau_syn_E'],\n",
    "                         'tau_syn_I': cell_params_lsm['tau_syn_I'],\n",
    "                         'tau_refrac': cell_params_lsm['tau_refrac_E'],\n",
    "                         'i_offset': cell_params_lsm['i_offset'],\n",
    "                      }\n",
    "# The only difference is the refractory period\n",
    "cell_params_lif_inh = {\n",
    "                         'tau_m': cell_params_lsm['tau_m'],\n",
    "                         'cm': cell_params_lsm['cm'],\n",
    "                         'v_rest': cell_params_lsm['v_rest'],\n",
    "                         'v_reset': cell_params_lsm['v_reset'],\n",
    "                         'v_thresh': cell_params_lsm['v_thresh'],\n",
    "                         'tau_syn_E': cell_params_lsm['tau_syn_E'],\n",
    "                         'tau_syn_I': cell_params_lsm['tau_syn_I'],\n",
    "                         'tau_refrac': cell_params_lsm['tau_refrac_I'],\n",
    "                         'i_offset': cell_params_lsm['i_offset'],\n",
    "                      }"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# imports of both spynnaker and external device plugin.\n",
    "import spynnaker.pyNN as spnkr\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"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# initial call to set up the front end (pynn requirement)\n",
    "spnkr_timestep = time_step*1E3 # SpiNNaker uses ms\n",
    "\n",
    "spnkr.setup(timestep=spnkr_timestep) #, min_delay=1.0, max_delay=144.0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "# Creates two separated populations: one for excitatory neurons and other for inhibitory\n",
    "# It's necessary to divide excitatory and inhibitory neurons because it's not possible to set\n",
    "# different values for the refractory period inside the same population \n",
    "# (Is that still true for 3.0.0???)\n",
    "pop_lsm_exc = spnkr.Population(total_size_exc, spnkr.IF_curr_exp, cell_params_lif_exc, label='LSM_EXC')\n",
    "\n",
    "pop_lsm_inh = spnkr.Population(total_size_inh, spnkr.IF_curr_exp, cell_params_lif_inh, label='LSM_INH')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "# Initialises some parameters individually for each neuron.\n",
    "pop_lsm_exc.tset('i_offset',[i for j in i_offset_exc for i in j])\n",
    "pop_lsm_inh.tset('i_offset',[i for j in i_offset_inh for i in j])\n",
    "\n",
    "pop_lsm_exc.initialize('v',[i for j in init_membrane_v_exc for i in j])\n",
    "pop_lsm_inh.initialize('v',[i for j in init_membrane_v_inh for i in j])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Activate the SENDING of live spikes from the pop_lsm_exc\n",
    "ExternalDevices.activate_live_output_for(\n",
    "    pop_lsm_exc, \n",
    "    database_notify_host=\"localhost\",\n",
    "    database_notify_port_num=12346)\n",
    "\n",
    "# Activate the SENDING of live spikes from the pop_lsm_inh\n",
    "ExternalDevices.activate_live_output_for(\n",
    "    pop_lsm_inh, \n",
    "    database_notify_host=\"localhost\",\n",
    "    database_notify_port_num=12346)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "LSM number:  0\n",
      "Liquid->Liquid connections... EXC to ???\n",
      "Original number of EXC=>??? connections: 913\n",
      "Filterer number of EXC=>??? connections: 827\n",
      "Liquid->Liquid connections... INH to ???\n",
      "Original number of INH=>??? connections: 255\n",
      "Filterer number of INH=>??? connections: 199\n",
      "LSM number:  1\n",
      "Liquid->Liquid connections... EXC to ???\n",
      "Original number of EXC=>??? connections: 936\n",
      "Filterer number of EXC=>??? connections: 826\n",
      "Liquid->Liquid connections... INH to ???\n",
      "Original number of INH=>??? connections: 259\n",
      "Filterer number of INH=>??? connections: 199\n",
      "LSM number:  2\n",
      "Liquid->Liquid connections... EXC to ???\n",
      "Original number of EXC=>??? connections: 888\n",
      "Filterer number of EXC=>??? connections: 790\n",
      "Liquid->Liquid connections... INH to ???\n",
      "Original number of INH=>??? connections: 267\n",
      "Filterer number of INH=>??? connections: 203\n",
      "LSM number:  3\n",
      "Liquid->Liquid connections... EXC to ???\n",
      "Original number of EXC=>??? connections: 895\n",
      "Filterer number of EXC=>??? connections: 805\n",
      "Liquid->Liquid connections... INH to ???\n",
      "Original number of INH=>??? connections: 268\n",
      "Filterer number of INH=>??? connections: 207\n",
      "LSM number:  4\n",
      "Liquid->Liquid connections... EXC to ???\n",
      "Original number of EXC=>??? connections: 905\n",
      "Filterer number of EXC=>??? connections: 812\n",
      "Liquid->Liquid connections... INH to ???\n",
      "Original number of INH=>??? connections: 235\n",
      "Filterer number of INH=>??? connections: 189\n"
     ]
    }
   ],
   "source": [
    "connected_internally = True\n",
    "\n",
    "internal_delay = spnkr_timestep\n",
    "\n",
    "liquid_proj = []\n",
    "connections_exc2exc = []\n",
    "connections_exc2inh = []\n",
    "connections_inh2inh = []\n",
    "connections_inh2exc = []\n",
    "for lsm_i in range(number_of_liquids):\n",
    "\n",
    "    print \"LSM number: \",lsm_i\n",
    "\n",
    "    lsm_proj = []\n",
    "\n",
    "    #\n",
    "    # Creating the EXC to ??? connections\n",
    "    #\n",
    "    print \"Liquid->Liquid connections... EXC to ???\"\n",
    "\n",
    "    indices_pre_exc_temp = connections_pre_e[lsm_i]\n",
    "    indices_pos_exc_temp = connections_pos_e[lsm_i]\n",
    "    weights_pre_exc_temp = weights_pre_e[lsm_i]\n",
    "\n",
    "    weights_pre_exc_temp = weights_pre_exc_temp.astype(numpy.float16) # This is what actually SpiNNaker uses.\n",
    "\n",
    "    # Now I need to verify which connections went to zero after changing to float16 and cut them out.\n",
    "    weights_pre_exc_temp_filter = numpy.arange(len(weights_pre_exc_temp))[weights_pre_exc_temp>0]\n",
    "\n",
    "    indices_pre_exc=indices_pre_exc_temp[weights_pre_exc_temp_filter]\n",
    "    indices_pos_exc=indices_pos_exc_temp[weights_pre_exc_temp_filter]\n",
    "    weights_pre_exc=weights_pre_exc_temp[weights_pre_exc_temp_filter]*1E9 # SpiNNaker uses nA!\n",
    "\n",
    "    print \"Original number of EXC=>??? connections:\",len(indices_pre_exc_temp)\n",
    "    print \"Filterer number of EXC=>??? connections:\",len(indices_pre_exc)\n",
    "\n",
    "\n",
    "    # EXC2EXC\n",
    "    indices_pre_exc_e2e=[from_bee2spinnaker[lsm_i][i][1] for i,j in zip(indices_pre_exc,indices_pos_exc) \n",
    "                                                             if from_bee2spinnaker[lsm_i][j][0]=='exc']\n",
    "    indices_pos_exc_e2e=[from_bee2spinnaker[lsm_i][i][1] for i in indices_pos_exc \n",
    "                                                             if from_bee2spinnaker[lsm_i][i][0]=='exc']\n",
    "    weights_pre_exc_e2e=[i for i in indices_pos_exc if from_bee2spinnaker[lsm_i][i][0]=='exc']\n",
    "\n",
    "    delays_exc2exc=[internal_delay]*len(weights_pre_exc_e2e)\n",
    "    if connected_internally:\n",
    "        # Generates the internal connections exc=>exc\n",
    "        connections_exc2exc.append(zip(indices_pre_exc_e2e,indices_pos_exc_e2e,weights_pre_exc_e2e,delays_exc2exc))\n",
    "        if len(connections_exc2exc[-1]):\n",
    "            lsm_proj.append(spnkr.Projection(pop_lsm_exc,pop_lsm_exc,\n",
    "                               spnkr.FromListConnector(conn_list=connections_exc2exc[-1]),\n",
    "                               label=\"EXC2EXC_conn_\"+str(lsm_i), target='excitatory'))\n",
    "\n",
    "\n",
    "    # EXC2INH\n",
    "    indices_pre_exc_e2i=[from_bee2spinnaker[lsm_i][i][1] for i,j in zip(indices_pre_exc,indices_pos_exc) \n",
    "                                                             if from_bee2spinnaker[lsm_i][j][0]=='inh']\n",
    "    indices_pos_exc_e2i=[from_bee2spinnaker[lsm_i][i][1] for i in indices_pos_exc \n",
    "                                                             if from_bee2spinnaker[lsm_i][i][0]=='inh']\n",
    "    weights_pre_exc_e2i=[i for i in indices_pos_exc if from_bee2spinnaker[lsm_i][i][0]=='inh']\n",
    "    \n",
    "    delays_exc2inh=[internal_delay]*len(weights_pre_exc_e2i)\n",
    "    if connected_internally:\n",
    "        # Generates the internal connections exc=>inh\n",
    "        connections_exc2inh.append(zip(indices_pre_exc_e2i,indices_pos_exc_e2i,weights_pre_exc_e2i,delays_exc2inh))\n",
    "        if len(connections_exc2inh[-1]):\n",
    "            lsm_proj.append(spnkr.Projection(pop_lsm_exc,pop_lsm_inh,\n",
    "                               spnkr.FromListConnector(conn_list=connections_exc2inh[-1]),\n",
    "                               label=\"EXC2INH_conn_\"+str(lsm_i), target='excitatory'))\n",
    "\n",
    "\n",
    "    #\n",
    "    # Creating the INH to ??? connections\n",
    "    #\n",
    "    print \"Liquid->Liquid connections... INH to ???\"\n",
    "    indices_pre_inh_temp = connections_pre_i[lsm_i]\n",
    "    indices_pos_inh_temp = connections_pos_i[lsm_i]\n",
    "    weights_pre_inh_temp = weights_pre_i[lsm_i]\n",
    "\n",
    "    weights_pre_inh_temp = weights_pre_inh_temp.astype(numpy.float16) # This is what actually SpiNNaker uses.\n",
    "\n",
    "    # Now I need to verify which connections went to zero after changing to float16 and cut them out.\n",
    "    weights_pre_inh_temp_filter = numpy.arange(len(weights_pre_inh_temp))[weights_pre_inh_temp<0]\n",
    "\n",
    "    indices_pre_inh=indices_pre_inh_temp[weights_pre_inh_temp_filter]\n",
    "    indices_pos_inh=indices_pos_inh_temp[weights_pre_inh_temp_filter]\n",
    "    weights_pre_inh=weights_pre_inh_temp[weights_pre_inh_temp_filter]*1E9 # SpiNNaker uses nA!\n",
    "\n",
    "    print \"Original number of INH=>??? connections:\",len(indices_pre_inh_temp)\n",
    "    print \"Filterer number of INH=>??? connections:\",len(indices_pre_inh)\n",
    "\n",
    "\n",
    "    # INH2EXC\n",
    "    indices_pre_inh_i2e=[from_bee2spinnaker[lsm_i][i][1] for i,j in zip(indices_pre_inh,indices_pos_inh) \n",
    "                                                             if from_bee2spinnaker[lsm_i][j][0]=='exc']\n",
    "    indices_pos_inh_i2e=[from_bee2spinnaker[lsm_i][i][1] for i in indices_pre_inh \n",
    "                                                             if from_bee2spinnaker[lsm_i][i][0]=='exc']\n",
    "    weights_pre_exc_i2e=[i for i in indices_pre_inh if from_bee2spinnaker[lsm_i][i][0]=='exc']\n",
    "\n",
    "    delays_inh2exc=[internal_delay]*len(weights_pre_exc_i2e)\n",
    "    if connected_internally:\n",
    "        # Generates the internal connections inh=>exc\n",
    "        connections_inh2exc.append(zip(indices_pre_inh_i2e,indices_pos_inh_i2e,weights_pre_exc_i2e,delays_inh2exc))\n",
    "        if len(connections_inh2exc[-1]):\n",
    "            lsm_proj.append(spnkr.Projection(pop_lsm_inh,pop_lsm_exc,\n",
    "                               spnkr.FromListConnector(conn_list=connections_inh2exc[-1]),\n",
    "                               label=\"INH2EXC_conn_\"+str(lsm_i), target='inhibitory'))\n",
    "\n",
    "    # INH2INH\n",
    "    indices_pre_inh_i2i=[from_bee2spinnaker[lsm_i][i][1] for i,j in zip(indices_pre_inh,indices_pos_inh) \n",
    "                                                             if from_bee2spinnaker[lsm_i][j][0]=='inh']\n",
    "    indices_pos_inh_i2i=[from_bee2spinnaker[lsm_i][i][1] for i in indices_pre_inh \n",
    "                                                             if from_bee2spinnaker[lsm_i][i][0]=='inh']\n",
    "    weights_pre_exc_i2i=[i for i in indices_pre_inh if from_bee2spinnaker[lsm_i][i][0]=='inh']\n",
    "    \n",
    "    delays_inh2inh=[internal_delay]*len(weights_pre_exc_i2i)\n",
    "    if connected_internally:\n",
    "        # Generates the internal connections exc=>inh\n",
    "        connections_inh2inh.append(zip(indices_pre_inh_i2i,indices_pos_inh_i2i,weights_pre_exc_i2i,delays_inh2inh))\n",
    "        if len(connections_inh2inh[-1]):\n",
    "            lsm_proj.append(spnkr.Projection(pop_lsm_inh,pop_lsm_inh,\n",
    "                               spnkr.FromListConnector(conn_list=connections_inh2inh[-1]),\n",
    "                               label=\"INH2INH_conn_\"+str(lsm_i), target='inhibitory'))\n",
    "\n",
    "    liquid_proj.append(lsm_proj)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Creates ONE BIG population that will receive spikes from the ethernet\n",
    "# and inject them into the liquid.\n",
    "# This population will be equivalent to all liquids concatenated.\n",
    "spike_injector = spnkr.Population(\n",
    "                            total_size, ExternalDevices.SpikeInjector,\n",
    "                            {'port': 12345,\n",
    "                             'database_notify_host':\"localhost\",\n",
    "                             'database_notify_port_num':12346},\n",
    "                            label='spike_injector')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "# # Generates the connections from the input (as if it was all LSMs concatenated as inside BEE) to the \n",
    "# # excitatory and inhibitory populations.\n",
    "\n",
    "# input2liquid_delays = spnkr_timestep\n",
    "# input2liquid_weights = 5.0\n",
    "\n",
    "# input_synapses_exc = []\n",
    "# exc_offset = 0\n",
    "# for lsm_i in range(number_of_liquids):\n",
    "#     input_offset = NofN*lsm_i\n",
    "#     exc_offset += (lsm_i>0)*len(exc_neuron_idx[lsm_i])\n",
    "#     # exc_neuron_idx has the BEE indices and the ones inside SpiNNaker \n",
    "#     # were simply concatenated for form a big excitatory population with all liquids.\n",
    "#     for i,e in enumerate(exc_neuron_idx[lsm_i]):\n",
    "#         input_synapses_exc.append((e+input_offset, \n",
    "#                                    i+exc_offset,\n",
    "#                                    input2liquid_weights,\n",
    "#                                    input2liquid_delays))\n",
    "        \n",
    "# input_synapses_inh = []\n",
    "# inh_offset = 0\n",
    "# for lsm_i in range(number_of_liquids):\n",
    "#     input_offset = NofN*lsm_i\n",
    "#     inh_offset += (lsm_i>0)*len(inh_neuron_idx[lsm_i])\n",
    "#     for i,e in enumerate(inh_neuron_idx[lsm_i]):\n",
    "#         input_synapses_inh.append((e+input_offset, \n",
    "#                                    i+inh_offset,\n",
    "#                                    input2liquid_weights,\n",
    "#                                    input2liquid_delays))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Generates the connections from the input (as if it was all LSMs concatenated as inside BEE) to the \n",
    "# excitatory and inhibitory populations.\n",
    "\n",
    "input2liquid_delays = spnkr_timestep\n",
    "input2liquid_weights = 5.0\n",
    "input_synapses_inh = []\n",
    "input_synapses_exc = []\n",
    "for lsm_i in range(number_of_liquids):\n",
    "    for idx in range(NofN):\n",
    "        spnkr_idx = from_bee2spinnaker[lsm_i][idx]\n",
    "        input_offset = NofN*lsm_i\n",
    "        if spnkr_idx[0]=='exc':\n",
    "            input_synapses_exc.append((idx+input_offset, \n",
    "                                       spnkr_idx[1],\n",
    "                                       input2liquid_weights,\n",
    "                                       input2liquid_delays))\n",
    "        else:\n",
    "            input_synapses_inh.append((idx+input_offset, \n",
    "                                       spnkr_idx[1],\n",
    "                                       input2liquid_weights,\n",
    "                                       input2liquid_delays))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Total number of synapses INPUT=>LSM_EXC population: 2400\n",
      "Total number of synapses INPUT=>LSM_INH population: 600\n"
     ]
    }
   ],
   "source": [
    "input_proj = []\n",
    "\n",
    "input_proj.append(spnkr.Projection(spike_injector, pop_lsm_exc,\n",
    "                            spnkr.FromListConnector(conn_list=input_synapses_exc), target=\"excitatory\"))\n",
    "\n",
    "input_proj.append(spnkr.Projection(spike_injector, pop_lsm_inh,\n",
    "                            spnkr.FromListConnector(conn_list=input_synapses_inh), target=\"excitatory\"))\n",
    "\n",
    "print \"Total number of synapses INPUT=>LSM_EXC population:\",  len(input_synapses_exc)\n",
    "print \"Total number of synapses INPUT=>LSM_INH population:\",  len(input_synapses_inh)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "poisson_rate = 10\n",
    "poisson_weight = 1.0\n",
    "pop_poisson_exc_exc = spnkr.Population(total_size_exc, \n",
    "                                   spnkr.SpikeSourcePoisson, {\"rate\": poisson_rate}, label='POISSON_NOISE_EXC')\n",
    "pop_poisson_exc_inh = spnkr.Population(total_size_exc, \n",
    "                                   spnkr.SpikeSourcePoisson, {\"rate\": poisson_rate}, label='POISSON_NOISE_EXC')\n",
    "\n",
    "pop_poisson_inh_exc = spnkr.Population(total_size_inh, \n",
    "                                   spnkr.SpikeSourcePoisson, {\"rate\": poisson_rate}, label='POISSON_NOISE_EXC')\n",
    "pop_poisson_inh_inh = spnkr.Population(total_size_inh, \n",
    "                                   spnkr.SpikeSourcePoisson, {\"rate\": poisson_rate}, label='POISSON_NOISE_EXC')\n",
    "\n",
    "\n",
    "poisson_proj = []\n",
    "poisson_proj.append(spnkr.Projection(pop_poisson_exc_exc, pop_lsm_exc, \n",
    "                                     spnkr.OneToOneConnector(weights=poisson_weight, delays=spnkr_timestep), \n",
    "                                     target=\"excitatory\"))\n",
    "poisson_proj.append(spnkr.Projection(pop_poisson_exc_inh, pop_lsm_exc, \n",
    "                                     spnkr.OneToOneConnector(weights=poisson_weight, delays=spnkr_timestep), \n",
    "                                     target=\"inhibitory\"))\n",
    "\n",
    "poisson_proj.append(spnkr.Projection(pop_poisson_inh_exc, pop_lsm_inh, \n",
    "                                     spnkr.OneToOneConnector(weights=poisson_weight, delays=spnkr_timestep), \n",
    "                                     target=\"excitatory\"))\n",
    "poisson_proj.append(spnkr.Projection(pop_poisson_inh_inh, pop_lsm_inh, \n",
    "                                     spnkr.OneToOneConnector(weights=poisson_weight, delays=spnkr_timestep), \n",
    "                                     target=\"inhibitory\"))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "# record spikes so that we can read off valid\n",
    "# results in a safe way afterwards, and verify the behavior\n",
    "pop_lsm_exc.record()\n",
    "pop_lsm_inh.record()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Launch the Sender/Receiver before starting the simulation!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {
    "collapsed": false,
    "scrolled": 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 LSM_EXC\n",
      "|0                           50%                         100%|\n",
      " ============================================================\n",
      "Getting spikes for LSM_INH\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"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "CPU times: user 6.13 s, sys: 607 ms, total: 6.74 s\n",
      "Wall time: 22.1 s\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "\n",
    "run_time = 1000\n",
    "\n",
    "# Run the simulation on spiNNaker\n",
    "spnkr.run(run_time)\n",
    "\n",
    "\n",
    "# Retrieve spikes from the synfire chain population\n",
    "spikes_exc = pop_lsm_exc.getSpikes()\n",
    "spikes_inh = pop_lsm_inh.getSpikes()\n",
    "\n",
    "# Clear data structures on spiNNaker to leave the machine in a clean\n",
    "# state for future executions\n",
    "spnkr.end()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "times = []\n",
    "spikes = []\n",
    "\n",
    "for lsm_i in range(number_of_liquids):\n",
    "    times.append([])\n",
    "    spikes.append([])\n",
    "    \n",
    "# If there are spikes, plot using matplotlib\n",
    "if len(spikes_exc) != 0 or len(spikes_inh) != 0:\n",
    "    if len(spikes_exc) != 0:\n",
    "        for i,j in spikes_exc:\n",
    "            lsm_i, idx = from_spinnaker_exc2bee[int(i)]\n",
    "            (times[lsm_i]).append(j)\n",
    "            (spikes[lsm_i]).append(idx)\n",
    "    if len(spikes_inh) != 0:\n",
    "        for i,j in spikes_inh:\n",
    "            lsm_i, idx = from_spinnaker_inh2bee[int(i)]\n",
    "            (times[lsm_i]).append(j)\n",
    "            (spikes[lsm_i]).append(idx)\n",
    "else:\n",
    "    print \"No spikes received\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "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,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4Xu3dCZhkZX3v8e+ADgoyE/AqVxaNCzgI4lUiBKOiRjMSNQZQ0ahADHFhEGVcQJLLpqCgTPSBgUQUjMaIC5AYCU70mohC4r4gV8AtiqK4wQzKDYPAff7Tb4VDUd1d1f+qt09Vfc/zoNPd9Z7l8/6r+tfvOe85S3BRQAEFFFBAAQUUmCqBJVN1tB6sAgoooIACCiigAAZAi0ABBRRQQAEFFJgyAQPglHW4h6uAAgoooIACChgArQEFFFBAAQUUUGDKBAyAU9bhHq4CCiiggAIKKGAAtAYUUEABBRRQQIEpEzAATlmHe7gKKKCAAgoooIAB0BpQQAEFFFBAAQWmTMAAOGUd7uEqoIACCiiggAIGQGtAAQUUUEABBRSYMgED4JR1uIergAIKKKCAAgoYAK0BBRRQQAEFFFBgygQMgFPW4R6uAgoooIACCihgALQGFFBAAQUUUECBKRMwAE5Zh3u4CiiggAIKKKCAAdAaUEABBRRQQAEFpkzAADhlHe7hKqCAAgoooIACBkBrQAEFFFBAAQUUmDIBA+CUdbiHq4ACCiiggAIKGACtAQUUUEABBRRQYMoEDIBT1uEergIKKKCAAgooYAC0BhRQQAEFFFBAgSkTMABOWYd7uAoooIACCiiggAHQGlBAAQUUUEABBaZMwAA4ZR3u4SqggAIKKKCAAgZAa0ABBRRQQAEFFJgyAQPglHW4h6uAAgoooIACChgArQEFFFBAAQUUUGDKBAyAU9bhHq4CCiiggAIKKGAAtAYUUEABBRRQQIEpEzAATlmHe7gKtETgduAE4KSyP/Hv44B7APGzti/d+7/Y+3t88du8jx15BPB2YB/gFuCjwGuAG/po60sUUGBCBAyAE9KRHoYCYyawF/BD4Lqy350Ac08D4IJ6cntgR+Dz87R+APA14P8CJwPbAG8DfgA8fkFbtpECCoylgAFwLLvNnVZg4gQMgHW69K+AQ4AHATeVTT4B+DRwAPAPdXbDrSigwGILGAAXuwfcvgLjI7AzcBrwOGAZ8FPgP4AXlFG7fYF/BZ4DPAt4NhCnJP8JOBL4ZeNQu0+h9gqATwc+DJxX2kfze5dTx88FdgB+BLwLeDNwR1n/VsBbgD8CtgPWl1GvI4BrhsTdvf/92uwP/GExis/fvy2nX/cGTgceBXwPWA38ywD72jmFvtk8bb4N/Dvw4q7X/Wfpuz8dYJu+VAEFxljAADjGneeuK1BZ4FvAL0q4iv+PABZh5iXAb4BOALwW+CTwQSCCUYSzLwC/P0AAPBg4p4S9aB9LhMl/A1aUawe/AfxuufbtTOB15XXR7pnAG4AIPPcFfg/4SB+nSPsl7Q6A/dpE0LoQuAR4IvC/gb8Bnlyc4pR4fO/RZZSuGZrn2rde1wC+BwjHTii8F/Br4Fjg1K6VXQz8DyCCqIsCCkyBgAFwCjrZQ1RgCAIRon5WRtU+Nsv6OgEwws0zGq/5E+B9wFPLKFP8aK4RwNcCbwReXkb/OquKUasINRGcLmusPwJNTCCJa+B+DlwBrANiPaNamvs/iM27gT9v7NSXgP9Vrr+LkblYHllGLONUbbj1s/QKgDEyGmZblBXE9X8xYhqu7+xaaWwnwnQEdhcFFJgCAQPgFHSyh6jAkARiNO2/ygzSGImLr5tLJwDGiGAEtc4SEzv+Xwlpp5RvzhYA1wKHAc8DuoPm35XTz90h5THA5xrh9Nzy7zilGqdRv9LHxJL4LGyePo3965xS7sXXvf/92sSp6wsaK3x/OV0ep9SbXjE79y+Bjtd8XdjPLGAD4HyK/lyBKRIwAE5RZ3uoCiQFfrucko3TvnG6MK5Veyvw12W9nQAYo38xCthcYvTwQ8CqOQJghJgYwYv1PqmExuY6IszFKGKvJcJaBM+4pi6uAYxRwQhbDy23N3kv8Bc91tlZV1xnGCNuscS6Yj2xvtmW7gDYr03s/6caK43txqnxB3ZtKNb/phKa++m2fgJgXD/5K08B98PpaxSYfAED4OT3sUeowCgE9gBiUkWM1u1XTrkOYwQwTod+Ari6rPfmxs5/AHhsCXa9Prvi+rrua+Z2KhMu4pq3GBGM6wJ7LRHAItR2lgiicWuUfgNg83Vz2SxmAIx9/E45fR7XBjaXCN0xquskkFG8W1ynAi0UMAC2sFPcJQXGRGDrMsM2Jl9EuOoEwI+XySGdw3ghECNw/V4DuEsZJYvTqhEuY+JCLDFCF9euxTVyg87mjWvtYoJFzE4exjLfjaB72cTI39MWcQQwjjtuAxPhL0YsO7eBifv/XeptYIZRFq5DgfERMACOT1+5pwospkCErneUmb0RzGJGbowWxb3jYvJAXGfXnAX8f4DzgYeXU5kxC7h5+na+28DEdX4xIhUjU3E7mDh1GU8JidHB+FkEzrih8VLgYY3bzsQ1ipeXp1vEZJBoF6eT43q6o4CYLTyMpbn/g9jUHAGMCScR9uIazM4SN4z+arkRdMyu/q0yIzgmh8RMaRcFFJgSAQPglHS0h6lAUuB+5R6A8fiwmG0bQSsCVjxNIm75EksEwBjlOrBMwvjjEhTjUWOv6jo9e1u5njBm+8YS17DF7U8i0HUeBRfX78V9BeO2MitLmIufHwM8H3hwGR2M05oxYST2JdpGsPkD4CElNH63jBzGBJNhLc39H8Sm1wjgU8otX5r7FuuPawDDpZ+l4xchubPE9YUxC7j5vfjZbsCaMqFmY7n5c8yY9lFw/Uj7GgUmRGDSAuBB5SLzuJnqfcpfvr2eK7pnuRlq3MQ2binRXE4s1zXFrLw4bRQXrV85If3tYSgwSoFOAOwOOaPcputWQAEFFFiAwKQFwPjFsy2wZXk6QK/nisY9seJ01E+AuDFqMwDGtUxxYXtcdxSjCvFXdZxCiWuSmhejL4DaJgpMvEDnFHD3ac6JP3APUAEFFBg3gUkLgB3/zkhErwAYDz6P+33F46Hi9gvNABiniuLUSOc6objOKS4cj8cyxf26XBRQYHYBRwBHUx3xOTTXEqeLXRRQQIGBBKYtAEbYi8cuxWOWju4KgHHK90YgrnGKm8p2lniiQFzrNMqnCgzUab5YAQWmRqAzqjrbAcc9C+NayLluWTM1WB6oAgr0LzBNATBuDhuz3w4t98GK07vNEcC4sD0+RHct9yDrKMZMxg3AS3uwhl/MquvcTqF/eV+pgAIKzC8Ql7PM93i2eCayo4DzW/oKBboF4nZNcZZvrqf+TKzaNAXAeFpBzFx8denN7gC4kBHAHYAfTmx1eGAKKKCAAgpMtkAM/sRtkKZumaYAGPcTW974Szn+so5rBOO0b9zHLK7/63UN4I/L/cN6XQMYoXH9tddey7JlzUd5Tl0dteKAjz32WE45pd9Hp7Zilyd2J+yLdnWt/dGe/rAv2tEXGzZsYKed4kFBm3JBnOWbumXSAmBM7ohQF9fNxLNIY3g3To3Eva7iXl3N+2G9ptz4dH/g+nL/sLjOL2YBx7NMIwweB7yo3My21yzgTQFw/fr1BsAWvHVWr17NmjUxh8dlsQXsi8Xugbtu3/5oT3/YF+3oiwiAy5dH9jMAtqNH8nsRj4qKm592zudHwI1/P7k86qi5he5TwJ2fnQC8rITHL85zH0ADYL7PhrYGP1iHRplekX2RJhzqCuyPoXKmVmZfpPiG1tgACJM2Aji04uhzRQbAPqFqvGzdunWsXBkPjHBZbAH7YrF74K7btz/a0x/2RTv6wgBoAMxWogEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJWcAzAraXgEFFFBAgcoCBkADYLbkDIBZQdsroIACCihQWcAAaADMlpwBMCtoewUUUEABBSoLGAANgNmSMwBmBW2vgAIKKKBAZQEDoAEwW3IGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkZALOCtldAAQUUUKCygAHQAJgtOQNgVtD2CiiggAIKVBYwABoAsyVnAMwK2l4BBRRQQIHKAgZAA2C25AyAWUHbK6CAAgooUFnAAGgAzJacATAraHsFFFBAAQUqCxgADYDZkjMAZgVtr4ACCiigQGUBA6ABMFtyBsCsoO0VUEABBRSoLGAANABmS84AmBW0vQIKKKCAApUFDIAGwGzJGQCzgrZXQAEFFFCgsoAB0ACYLTkDYFbQ9goooIACClQWMAAaALMlZwDMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJWcAzAraXgEFFFBAgcoCBkADYLbkDIBZQdsroIACCihQWcAAaADMlpwBMCtoewUUUEABBSoLGAANgNmSMwBmBW2vgAIKKKBAZQEDoAEwW3IGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkZALOCtldAAQUUUKCygAHQAJgtOQNgVtD2CiiggAIKVBYwABoAsyVnAMwK2l4BBRRQQIHKAgZAA2C25AyAWUHbK6CAAgooUFnAAGgAzJZczwD4sY/BMcfMrPqAA+DCCyfr3xs2wJIlsPXWk3l8k9pvmeOyzyfrPTxbLdjP09HPzf6f1j6/7bYNXHXV8qCI/9mQDQPj2H7JOO50i/a5ZwDcfXe48sqZvdxiC7jllsn7d6cPJvX4PK6716x9Prnv52a928/T0c/2eWQ+A2CL8tTY7YoBcEIDrgHQADiJf7j1U9cGQAPgJA9e3PkeMAA6ApjLnPOeAj7wQLjggpmNTMq/m6cMJuWYJql/RnEs9vlkvYdnqxH7eTr6udn/09rnngKevGsADwJWAY8C7gPcE7i9ZLy9gL8EHgvcG/g+8FfAe7oy4InAYUCEuy+V9ZUTundLi04CyQVoWyuggAIKKFBdwEkgkxcAnwZsC2wJvKsrAO4H3A/4Z+DnwJOAfwReDHy0VN/rgCOAeO13gOOBg4FdgJt7VKgBsPrb1g0qoIACCiiQEzAATl4A7FTEvsCnugJgr2q5CPhP4Kjyw+8Ca4Azy9ebA9cBq4H3GwBzbzhbK6CAAgoo0AYBA+B0B8AYvbsKOBp4XznleyOwD/C5RoGuA64AXmsAbMPb1n1QQAEFFFAgJ2AAnN4AGNcGxmnf+P8/KNcJ7gj8ANgVuLpRWueXewS91ACYe8PZWgEFFFBAgTYIGACnMwDGBJA49XsP4I8a1/bFiKAjgG14Z7oPCiiggAIKjFDAADh9AfC3gIvLJJDnALd21VevawB/XK4RnPUawFWrVrF06dJNq1q5cuWm/1wUUEABBRRQoD0C69atI/6LZePGjaxduzb+6ZNA2tNFqT3ZrJzWjUkglwBbA7dFXwP3Bz4BxC1dXlS+372xuM4vZgE/A4gweFx57cOdBZzqFxsroIACCijQGgFHACdvBPAQ4DzgjlJlcaPr+PeTy21f4rYundu5dF7zmRL4OoV5AvCyEh6/6H0AW/N+dUcUUEABBRQYioABcPIC4FAKY4CVeB/AAbB8qQIKKKCAAm0QMAAaALN1aADMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJWcAzAraXgEFFFBAgcoCBkADYLbkDIBZQdsroIACCihQWcAAaADMlpwBMCtoewUUUEABBSoLGAANgNmSMwBmBW2vgAIKKKBAZQEDoAEwW3IGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkZALOCtldAAQUUUKCygAHQAJgtOQNgVtD2CiiggAIKVBYwABoAsyVnAMwK2l4BBRRQQIHKAgZAA2C25AyAWUHbK6CAAgooUFnAAGgAzJacATAraHsFFFBAAQUqCxgADYDZkjMAZgVtr4ACCiigQGUBA6ABMFtyBsCsoO0VUEABBRSoLGAANABmS84AmBW0vQIKKKCAApUFDIAGwGzJGQCzgrZXQAEFFFCgsoAB0ACYLTkDYFbQ9goooIACClQWMAAaALMlZwDMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJWcAzAraXgEFFFBAgcoCBkADYLbkDIBZQdsroIACCihQWcAAaADMlpwBMCtoewUUUEABBSoLGAANgNmSMwBmBW2vgAIKKKBAZQEDoAEwW3IGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkZALOCtldAAQUUUKCygAHQAJgtOQNgVtD2CiiggAIKVBYwABoAsyV3twD4sY/BMcfMrPaAA+DCCyfr3xs2wJIlsPXWk3l8k9pvmeOyzyfrPdxPLdjn09Hn09zPxx23gYMOWh4dHf+zIRsGxrH9knHc6Rbt890C4O67w5VXzuzhFlvALbdM3r87/pN6fB7X3WvWPp/c9/Ns9W6fT0efT2s/r1ixgauuMgC2KE+N3a4YACcw4BoADYCT+IfboHU9rcFgUKdxf/209rMB0FPA2cR5lwAYp38PPxx++lPYbjs49FB4z3vg+utnTpluueXM5m6+GW66aeY1j3scXHTRzPf33x8uv7zdr499v/FG2Hxz2Gsv+Pznx2ffa9v30//9vGYxayRO999ww0wfb799++q1jX5XXDHz/h7X97l93v7P4exnWed9c911M+/tbbaBO+4Y399Lg3rE7+Pjj9/A857nCGA2BE1z+7sEwObp3912g298A7q/F1idU8TxmnH+uvmX77gfyyj2v5/+7+c1i1kjbe/jNvq1cZ8GqSH7fLw/l/v5LOuu0bb3+SD12+/xOwlk8kYADwJWAY8C7gPcE7i9kVD3AM4A9gRuBM4BTuxKsPH1YUCEuy+V9ZWr+u6WdQ2A5RRwP2+6Yb+J276+foJAP69ZzONs+y+GNvq1cZ8GqSH73AA4SL2M42d/vEcNgJMXAJ8GbAvEydZ3dQXACITXAOcCJwG7AJcAbwPeUaLd64AjgP2A7wDHAweX197cY6jzvwPgpZcuu8vp37POmhlS7z4lfN55dz1FPK5fL1t211MGcUphXI+lc7p+mPvfT//385rFdu0+/XzZZe2p3zb6tXGfBqmhtr+vx913FJ81g/RvbL/bsO19Pujx9fP6Cy6A225zEsikzgLeF/hUVwA8BDg1LmVqjAoeCbwS2LmEu+8Ca4Azy9ebA3GVxGrg/XMFwMc9btldTu2O+yjAfH/VtX2UYLH/gu2n//t5zWIfR3P7bevzNvq1cZ8GqaG29XH3vo+773yfqzV+Pu2nf++sqbjzi9cATuI1fL0CYAS7XcvoXueY9wE+W6pgs3JaOL73uQbKOuAK4LUGwN63uKnxoTXIL7E27E8/v6j6eU2bjrtt4aCNfm3cp0FqqG19bAAc/uloA2Dn95gBcJpGAOOU8FbACxpBbkXMyQB2AiIA/qCExKsbrzm/3CTypYOcAj7kEPjbv4XOLKuYQbnPPvAP/xBDzzOzru5977vOFByHn88Mnc/Mal669M5ZY+Ow7zXt73//mRtlxyzwzqzwbqN+XrPYNRKXMXTP+m5L/bbRb9zf93Gngpi9HEsbP7Pa2Ofj9tnXbRihv9Pn4/p7adDP9vhc22orTwFPUwAc6QjgsmXL7jLjt/mXdCc4zve9cfv5Qm8g28tj3I49jmG+4+8eiex1jP28Zr7tjPrnba7fcfAbx9q2z3vfxH9SPru63zfdgxuj/kxZ7PXH8V9++QaWL/cU8LScAo7JHKd1XQP4qjLpY65rAH8MHDXXNYCrVq1i6dKlvPe98ItfrARW3iUctPnDNPNGzLTtJ0CN+/r7CSf9vGaxHdpcv+PgZwC866+YbD2PW5+3sf+nNwDGFV3ruO994fnP38jatWujOH0U3ISkwDiNG7d+iWsAY4bv1sBtwMZy+jdO7cYs4JPLxI+LgdMbs4DjOr+YBfwMICaEHAe8CHh43L95rlPAMQLYfA7wgQdCnC5tPmtxru/Fusfh590zZeMYx2XfZ+uPUe3/qafOzJTuPBu6V//285pOu1Ht51zrj/3vPPu5jfXZdr+4VKKNbvb5zOfzYrynmjWxWNtvvm86v6N6vdcXa/9G/ZkXx/+EJzgCOGmngGOm73nAHSWsxfHFv58MXArsDpxV7gO4HjgbeGNXsDsBeFkJj18c5D6AExKiPQwFFFBAAQUmWsD7AE7efQBrF+zdngVcewfcngIKKKCAAgoMJmAANAAOVjF3f7UBMCtoewUUUEABBSoLGAANgNmSMwBmBW2vgAIKKKBAZQEDoAEwW3IGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkZALOCtldAAQUUUKCygAHQAJgtOQNgVtD2CiiggAIKVBYwABoAsyVnAMwK2l4BBRRQQIHKAgZAA2C25AyAWUHbK6CAAgooUFnAAGgAzJacATAraHsFFFBAAQUqCxgADYDZkjMAZgVtr4ACCiigQGUBA6ABMFtyBsCsoO0VUEABBRSoLGAANABmS84AmBW0vQIKKKCAApUFDIAGwGzJGQCzgrZXQAEFFFCgsoAB0ACYLTkDYFbQ9goooIACClQWMAAaALMlZwDMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJWcAzAraXgEFFFBAgcoCBkADYLbkDIBZQdsroIACCihQWcAAaADMlpwBMCtoewUUUEABBSoLGAANgNmSMwBmBW2vgAIKKKBAZQEDoAEwW3IGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkZALOCtldAAQUUUKCygAHQAJgtOQNgVtD2CiiggAIKVBYwABoAsyVnAMwK2l4BBRRQQIHKAgZAA2C25AyAWUHbK6CAAgooUFnAAGgAzJacATAraHsFFFBAAQUqCxgADYDZkjMAZgVtr4ACCiigQGUBA6ABMFtyBsCsoO0VUEABBRSoLGAANABmS84AmBW0vQIKKKCAApUFDIAGwGzJGQCzgrZXQAEFFFCgsoAB0ACYLTkDYFbQ9goooIACClQWMAAaALMlZwDMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAXA6A+A2wBrg6cAWwDeBNwCXlvrbAzgD2BO4ETgHOHGW2jQAVn7TujkFFFBAAQWyAgbA6QyAfwfsBOwP3AAcBZwAPBD4DXANcC5wErALcAnwNuAdPQrOAJh9F9peAQUUUECBygIGwOkMgFcA7yyjfFFyWwE3AXsDuwKnAdsDt5d6PBJ4JbCzAbDyO9TNKaCAAgooMAIBA+DiB8DvAXf00bcP6eM1/b4kTuc+BXgu8AtgNfASIE79vrmEwP0aK9sH+CywHPhV10YcAexX3dcpoIACCijQEgED4OIHwD9r1EKclj0ceC8QwfDBwIuBs+a4Bm8hpXQv4CJgZTnl+8tyOvjfgXeVEcEXNFa8AriynDa+zgC4EHLbKKCAAgoo0B4BA+DiB8BmNXwS+N9ABLHO8rvAm4CnDrFsPlOu84uRvzj1+8wSOp8IHOoI4BClXZUCCiiggAItFDAAtisAbgBihu5tjVrZvEzUiFOtw1juC/wMeDTwtcYKvwScD1zf4xrAVwFHzHUN4KpVq1i6dOmm1a1cuXLTfy4KKKCAAgoo0B6BdevWEf/FsnHjRtauXRv/jMu7In9M3bKkRUf89TLT9t2NfYpr82KW7iOHuJ/fBWK08TXlmr5nAB8G/hD4AnB1mQV8cgl9FwOnOwt4iD3gqhRQQAEFFFhEAUcA2zUCGMNm/wjELN24BvC3S/CL27V8fIh18ogyyhezfu8JXAu8HegEz93LdYdxH8D1wNnAG2fZvpNAhtgxrkoBBRRQQIEaAgbAdgXA6POY+BETMHYEfgh8oITBGvWwkG0YABeiZhsFFFBAAQUWUcAA2L4AuIjlsKBNGwAXxGYjBRRQQAEFFk/AALj4ATBu8/K+UgJxvd9sSzyZo42LAbCNveI+KaCAAgooMIeAAXDxA+A3gLjmLpa47q/XEjeKHuaNoIf5pjAADlPTdSmggAIKKFBBwAC4+AGwQjePdBMGwJHyunIFFFBAAQWGL2AANABmq8oAmBW0vYY/emsAACAASURBVAIKKKCAApUFDIAGwGzJGQCzgrZXQAEFFFCgsoAB0ACYLTkDYFbQ9goooIACClQWMAAaALMlZwDMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwADYvgC4L7AXsHVXLRxXuTb63ZwBsF8pX6eAAgoooEBLBAyA7QqAJwHHAF8Fft2okbgP4FNaUjPdu2EAbGnHuFsKKKCAAgrMJmAAbFcA/DGwP/AfY1SyBsAx6ix3VQEFFFBAgRAwALYrAP4cuB8QI37jshgAx6Wn3E8FFFBAAQWKgAGwXQHwDOBS4MNjVKEGwDHqLHdVAQUUUEABRwBnamBJi0rhA8AfA5cB13Xt18Et2s/mrhgAW9ox7pYCCiiggAKzCTgC2K4AeN4cpfqnLS1jA2BLO8bdUkABBRRQwAA4ew20aQRwHCvVADiOveY+K6CAAgpMtYAjgO0aAeyckn4s8EDg+8AXWz4pxAA41R8hHrwCCiigwDgKGADbFQAfAHwUeAxwA7AN8GXg2T2uCWxLvRkA29IT7ocCCiiggAJ9ChgA2xUAzy+TUl4B/BLYFlgLbAYc1Gef1n6ZAbC2uNtTQAEFFFAgKWAAbFcAjBtBPzzuz9jo1+XA1cD/TPb1qJobAEcl63oVUEABBRQYkYABsF0B8CfAw4BfNfo7ngn8bWC7EdVAdrUGwKyg7RVQQAEFFKgsYABsVwD8IHArcHgZBYzRvzOBLYDnVa6NfjdnAOxXytcpoIACCijQEgEDYLsC4PbAxcAjyzWAMQnkG8AzgR+1pGa6d8MA2NKOcbcUUEABBRSYTcAA2K4AGP0UEz72AnYCrgU+D9ze4hI2ALa4c9w1BRRQQAEFegkYANsTAJeWW778DvBfY1SuBsAx6ix3VQEFFFBAgRAwALYnAEZ/xGneBwMbx6g8DYBj1FnuqgIKKKCAAgbAmRpo06PgTgJuAU4eo/I0AI5RZ7mrCiiggAIKGADbFwA/A+wNXF8eA9e89u+JLS1ZA2BLO8bdUkABBRRQYDYBTwG3awTw+DlK9cSWlrEBsKUd424poIACCihgAJy9Btp0CngcK9UAOI695j4roIACCky1gCOA7RoBjFvAzLa09VYwBsCp/gjx4BVQQAEFxlHAANiuABgh745ZCmnzlhaYAbClHeNuKaCAAgoo4Cng8TgFvG/Xbu4AvAY4B/jrlpaxAbClHeNuKaCAAgooYAAcjwDYay8fArwXeHxLy9gA2NKOcbcUUEABBRQwAI5vAIxTvzcAEbTauBgA29gr7pMCCiiggAJzCHgNYLuuAYzRvuayFXAY8GRgj5ZWsgGwpR3jbimggAIKKOAI4HiMAHZPAolb1HwXOBi4vKVlbABsace4WwoooIACChgAxyMAPqhrN28Cftny8jUAtryD3D0FFFBAAQW6BTwF3K5TwONYoQbAcew191kBBRRQYKoFDIDtCoBxI+ijgZcA9weWAyuBB5ZbwbSxWA2AbewV90kBBRRQQIE5BAyA7QqAJwHPAt4KnAX8FvBQ4MPAY1payQbAlnaMu6WAAgoooMBsAgbAdgXA7wFPBK4t1/5tC8So4M+B+HcbFwNgG3vFfVJAAQUUUMARwDlrIGbatmWJoHe/8ji4mPwRoe+ewI/KKeG27GdzPwyAbewV90kBBRRQQAED4NgEwH8tp37jlG8nAD4H+PNyLWAbi9kA2MZecZ8UUEABBRQwAI5NANwb+CTwUWB/4O+BCIBPA74w5EreB3gT8FjgNuDKxuPm4qbTZwB7AjeWCSgnzrJ9A+CQO8bVKaCAAgooMGoBrwFs1zWA0d8PBw4HdgZ+UkYEvzjkQojw98/AEcBHgFtL2IuQeR/gGuBcICal7AJcArwNeEeP/TAADrlzXJ0CCiiggAKjFjAAti8AjrrPY/2XAp8DXtdjY4cApwLbA/FkkliOBF5ZQml3EwNgjR5zGwoooIACCgxRwADYvgAYzwOOU69bd/VzjMgNY7k3EE8YOR14UrnNTMw+fjNwIbAG2BXYr7GxGDH8bLkv4a+6dsIAOIxecR0KKKCAAgpUFDAAtisAvhw4s0wA+XWjDu4AIhgOY9mh3GbmeuAZwFeBZwPnA/sChwFbAS9obGxFuUZwJ+A6A+AwusF1KKCAAgoosHgCBsB2BcDvA0eVkbhRVUWM2MXEjrcAxzY28nHgK8AWCxkBXLVqFUuXLt20upUrV276z0UBBRRQQAEF2iOwbt064r9YNm7cyNq1a+Of8dSxDe3Zy3p70qb7AEYwi6d/jHr5Vnm6SK8A+E3gtK5rAF9VJozExJTuxVPAo+4t16+AAgoooMCQBRwBbNcI4PuBdwKfHnI/d68uJnUcU67z+3p5/FycAo6nkFwFXF1mAZ9cJn5cXK4ZdBbwiDvG1SuggAIKKFBDwADYrgD4diBm4V7Q41q744ZcEEcDq8rQb4wIngB8rGxj93L7mZiMsh44G3jjLNt3BHDIHePqFFBAAQUUGLWAAbBdATCeBNJriUkgTxl1MSxw/QbABcLZTAEFFFBAgcUSMAC2KwAuVh1ktmsAzOjZVgEFFFBAgUUQMAAaALNlZwDMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJWcAzAraXgEFFFBAgcoCBkADYLbkDIBZQdsroIACCihQWcAAaADMlpwBMCtoewUUUEABBSoLGAANgNmSMwBmBW2vgAIKKKBAZQEDoAEwW3IGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkZALOCtldAAQUUUKCygAHQAJgtOQNgVtD2CiiggAIKVBYwABoAsyVnAMwK2l4BBRRQQIHKAgZAA2C25AyAWUHbK6CAAgooUFnAAGgAzJacATAraHsFFFBAAQUqCxgADYDZkjMAZgVtr4ACCiigQGUBA6ABMFtyBsCsoO0VUEABBRSoLGAANABmS84AmBW0vQIKKKCAApUFDIAGwGzJGQCzgrZXQAEFFFCgsoAB0ACYLTkDYFbQ9goooIACClQWMAAaALMlZwDMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJWcAzAraXgEFFFBAgcoCBkADYLbkDIBZQdsroIACCihQWcAAaADMlpwBMCtoewUUUEABBSoLGAANgNmSMwBmBW2vgAIKKKBAZQEDoAEwW3IGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkZALOCtldAAQUUUKCygAHQAJgtOQNgVtD2CiiggAIKVBYwABoAsyVnAMwK2l4BBRRQQIHKAgZAA2C25AyAWUHbK6CAAgooUFnAAGgAzJacATAraHsFFFBAAQUqCxgADYDZkjMAZgVtr4ACCiigQGUBA6AB8CLg2cBTgU+V+tsDOAPYE7gROAc4cZbaNABWftO6OQUUUEABBbICBsDpDoAHA38CPK38FwHwPsA1wLnAScAuwCXA24B39Cg4A2D2XWh7BRRQQAEFKgsYAKc3AO4IfBZ4PPCDxgjgIcCpwPbA7aUejwReCexsAKz8DnVzCiiggAIKjEDAADi9AXAd8CHg3SXodU4BrwF2BfZr1Ns+JSwuB37VVYeOAI7gjekqFVBAAQUUGKWAAXA6A+Dh5bq/laW4YqSvEwDfBWwFvKBReCuAK4GdgOsMgKN8S7puBRRQQAEFRi9gAJy+APiQMpq3N3BtjwC4oBHAVatWsXTp0k2rW7ly5ab/XBRQQAEFFFCgPQLr1q0j/otl48aNrF27Nv4ZZ/c2tGcv6+3JknqbasWW4hq/vymd3Tn2+wLrgQ8ClwNvBR7QuAbwVcARXgPYiv5zJxRQQAEFFEgLOAI4fSOA9wK27aqcHwIHAZ8AfgNcXWYBn1xC38XA6c4CTr/fXIECCiiggAKtEDAATl8A7FV4tzVuAxM/3x04q9wHMEYGzwbeOEvFOgmkFW9ld0IBBRRQQIH+BQyABsD+q6X3Kw2AWUHbK6CAAgooUFnAAGgAzJacATAraHsFFFBAAQUqCxgADYDZkjMAZgVtr4ACCiigQGUBA6ABMFtyBsCsoO0VUEABBRSoLGAANABmS84AmBW0vQIKKKCAApUFDIAGwGzJGQCzgrZXQAEFFFCgsoAB0ACYLTkDYFbQ9goooIACClQWMAAaALMlZwDMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJWcAzAraXgEFFFBAgcoCBkADYLbkDIBZQdsroIACCihQWcAAaADMlpwBMCtoewUUUEABBSoLGAANgNmSMwBmBW2vgAIKKKBAZQEDoAEwW3IGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkZALOCtldAAQUUUKCygAHQAJgtOQNgVtD2CiiggAIKVBYwABoAsyVnAMwK2l4BBRRQQIHKAgZAA2C25AyAWUHbK6CAAgooUFnAAGgAzJacATAraHsFFFBAAQUqCxgADYDZkjMAZgVtr4ACCiigQGUBA6ABMFtyBsCsoO0VUEABBRSoLGAANABmS84AmBW0vQIKKKCAApUFDIAGwGzJGQCzgrZXQAEFFFCgsoAB0ACYLTkDYFbQ9goooIACClQWMAAaALMlZwDMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJWcAzAraXgEFFFBAgcoCBkADYLbkDIBZQdsroIACCihQWcAAaADMlpwBMCtoewUUUEABBSoLGAANgNmSMwBmBW2vgAIKKKBAZQEDoAEwW3IGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkZALOCtldAAQUUUKCygAFwOgPgm4FnAA8CfgV8Gng98MNG/e0BnAHsCdwInAOc2KM+DYCV37RuTgEFFFBAgayAAXA6A+DJwEeAK4AtgbOBRwCPLgV1H+Aa4FzgJGAX4BLgbcA7uorOAJh9F9peAQUUUECBygIGwOkMgN1l9ijgy8C2wHrgEOBUYHvg9vLiI4FXAjsbACu/S92cAgoooIACQxYwABoAo6Ti9O/LgYeU+loD7Ars16i3fYDPAsvLaePOjxwBHPKb0tUpoIACCigwagEDoAHwqcBFwAHAJ0rBvQvYCnhBowBXAFcCOwHXNb5vABz1u9T1K6CAAgooMGQBA+B0B8BnAu8rp3w/2qitgUcAV61axdKlSzetYuXKlZv+c1FAAQUUUECB9gisW7eO+C+WjRs3snbt2vhnnNnb0J69rLcnS+ptqlVbeiFwJvBc4JNde3YwcFrXNYCvAo7wGsBW9aE7o4ACCiigwIIEHAGczhHACHIxu/dZwGU9KidmAV9dZgHHjOGY+HExcLqzgBf0PrORAgoooIACrRIwAE5nAIyZvbcCt5RqjFHQO8qkj04g3B04q9wHMGYGx61i3tijer0GsFVvaXdGAQUUUECB+QUMgNMZAOevjP5fYQDs38pXKqCAAgoo0AoBA6ABMFuIBsCsoO0VUEABBRSoLGAANABmS84AmBW0vQIKKKCAApUFDIAGwGzJGQCzgrZXQAEFFFCgsoAB0ACYLTkDYFbQ9goooIACClQWMAAaALMlZwDMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJWcAzAraXgEFFFBAgcoCBkADYLbkDIBZQdsroIACCihQWcAAaADMlpwBMCtoewUUUEABBSoLGAANgNmSMwBmBW2vgAIKKKBAZQEDoAEwW3IGwKyg7RVQQAEFFKgsYAA0AGZLzgCYFbS9AgoooIAClQUMgAbAbMkZALOCtldAAQUUUKCygAHQAJgtOQNgVtD2CiiggAIKVBYwABoAsyVnAMwK2l4BBRRQQIHKAgZAA2C25AyAWUHbK6CAAgooUFnAAGgAzJacATAraHsFFFBAAQUqCxgADYDZkjMAZgVtr4ACCiigQGUBA6ABMFtyBsCsoO0VUEABBRSoLGAANABmS84AmBW0vQIKKKCAApUFDIAGwGzJGQCzgrZXQAEFFFCgsoAB0ACYLTkDYFbQ9goooIACClQWMAAaALMlZwDMCtpeAQUUUECBygIGQANgtuQMgFlB2yuggAIKKFBZwABoAMyWnAEwK2h7BRRQQAEFKgsYAA2A2ZIzAGYFba+AAgoooEBlAQOgATBbcgbArKDtFVBAAQUUqCxgADQAZkvOAJgVtL0CCiiggAKVBQyABsBsyRkAs4K2V0ABBRRQoLKAAdAAmC05A2BW0PYKKKCAAgpUFjAAGgCzJWcAzAraXgEFFFBAgcoCBkADYLbkDIBZQdsroIACCihQWcAAaADMlpwBMCtoewUUUEABBSoLGAANgNmSMwBmBW2vgAIKKKBAZQEDoAEwW3KbAuAHP7iek06Kf8Jb3jKzymOOmfn/Aw6ACy+crH9v2ABLlsDWW0/m8U1qv2WOyz6frPfwbLVgP09HPzf7f1r7/LbbNnDVVcuDIv5nQzYMjGP7JeO40y3a500BcMWK9Vx11UwA3G23mb278sqZ/99iC7jllsn7d6cPJvX4PK6716x9Prnv52a928/T0c/2eWQ+A2CL8tTY7YoBcEIDrgHQADiJf7j1U9cGQAPgJA9e3PkeMAA6ApjLnJsC4Otfv56LL54ZATz1VPjc5+C002ZWvP/+cPnlcP31M6dMt9xy5vs33ww33QTbbQePexxcdNF4vD5O/d5ww8y+br/9eO37uNsv5v7fdttM7UYomq2ODz0U3vOeyan1xfRerM+O+Fy68UbYfHPYay/4/OfH43NpGvtqWL9H+nlvj9PvqH5q4YorwFPAXgOYi3/w3yOA3/zmTACMZffd7zwF3H1KeNy/bo4gjPuxuP8z9dq5XGE+j/n6/hvfmOzan89n0n4+X38PUjuTZjNpxzNtfR2fVU4CMQDOFQBPBA5jJuR9CVgVvyu7GmwKgDvuuJ7ly++cBBITQJrXAMYo3w9+cOcptXH+OkYA77hjMo4lPvRG3Rch1TyV2Pz6gQ+cGS3r/LztXzf7vte+xsjgt789+/GG91weg/68bV4RCmJUv/Nej/2LZVy/jlHAGB2KpQ3HEvU1qZ+rNT6Luj/rmu/F7s/1Yb5PB31f13j9wx4Gxx23gYMO8hrA7CjYJLZ/HXAEsB/wHeB44GBglzh72zjgTQFwhx3W86Mf3TkJJGYCP+c5d/4i3Gkn+OlPJ+fr5odFu45tHTvttHLRrZctmz8IxQdQzL679to7w/T979/er7t/Qcy3r7ffvo5bb1256eBG8YHeRr8LLoBXvOLOPoz3RiydPl7cr2feG4Psz2abwe23zxzD4u77zPbPOmtSPlfXASs3HVPt3wvdn03TFvyawXbFCmcBew1g7/j6XWANcGb58ebAdcBq4P3dAbB7FvCknwpr7+mC1ey2W3Rb/6c1R3Eqp5/+7+c1i30cze0P2uc/+clqfvGLmb5YLOPafu3u08HfG4P2+ai92+07SJ3Hr5E1I3lfzNcH3YZt6+P59n+YPzcAegq4V/yLobwbgX2AzzVeEH+2XQG8tjsAdu4DGKd/4vRo/JXVuf9f53vRJv7aiqVzCnVcv+5MYInjvNe9Zo5pIccSp3QOPBDOO2/mL+E4RRETCRb69dKlqznqqDULbp/dfrSPUYro38MPn/2Y+nlNxqHXccQ2m5OTjj56Zj87k5Xm+7o5IWC+18bP161bzVe/OhMA+3n9IPsS63vsY+c2Xgy/2vs03zE2+/w3v1nNX/zFmpH2+aB9OOjra/tm3zOz1f0pp6zmHvdYM5L3xXym3Ybdkzsuu2w4n8Wdz8LMZ84oPjeaPq9+9QZOPdVTwL3HwKb3uzvGZTvArsDVDYbzy80iX9odAK+99lqWLVvG3nvDVVfN/HTFiplfuN3fi581XzPOXy9dChs33nm8gx7LsH222eZYttvulEX17eeY+nnNoJbzvT67zUH7+vrrj+WGG0757/fCfPvX9p9n/eLzoPYx3nWfj2XFipn+6PfzZ9A+H2Tdi+Ex6P6Nrs+PBU7Z9Dti0H0a5euH3d+j81v475ym3847b+Bb39p0nYY3gp7evHe3Ix9kBHAH4IfaKaCAAgoooMBYCsSgz4/Gcs+TO+01gL0Be10D+GPgqK5rAMNve+CmZD/YXAEFFFBAAQXqCmxdru8v97aou/HF3poBsHcPxHV+MQv4GUCEweOAFwEP75oFvNj95/YVUEABBRRQQIGBBQyAs5OdALwMiL8QvjjLfQAHBreBAgoooIACCiiw2AIGwMXuAbevgAIKKKCAAgpUFjAA5sD7eVpIbgvT3frN5TT8g4BfAZ8GXt818WYP4Axgz3L7nnOA6JfmYj8Nv47i6dXPBp4KfKqs3r4YvvN8a4zbVb0JeCwQz+2IpxU93v6Yj23oP9+m3Dv26XHvc+CbwBuAS+2LoVs3V3hQOTv3KOA+wD2BcvvyTS8bxmdSP+sY6UGOauUGwIXL9vu0kIVvwZYnAx8p91/cEjgbeATw6EITb/hrgHOBk8qTWi4B3ga8o7zGfhp+HcVTcf4EeFr5LwKgfTF85/nWGOHvn8v1yvE+ubX8IfQF+2M+uqH//O/igSnA/sANZcJgXEYUDwT8jZ9TQ/furDA+g7YF4vfDu7oC4DA+k/pZx8gObtQrNgAuXLjfp4UsfAu27BaIv/K+XN7w64FDgFPLTOzOX31HAq8Edi6N7afh1lHcMuGzZZQp7pfZGQG0L4br3M/aYnQpblYff+R0L/ZHP4LDe008JOCd5WxErHWrcneIvcs9ZU/zc2p42D3WtG85E9EcARzGe6CfdYz0wEa5cgPgwnQHuVfgwrZgq14Ccfr35cBDyg/jURNxw+54ZnNniVGRCChxc8/NBniqi+L9CcQTcT4EvLucaukEQPuiP79hvereJWCcDjwJeCjwPSAum7iwnI70vTEs7fnXE5eZPAV4LvCL8tjQl5RTkNEn9sX8hplX9AqAw/hMmm8dcWnS2C4GwIV13SBPC1nYFmzVLRBBI647OwD4RPlhDPnHX9ovaLw47q8f10HF6ZgIgP0+1UXx+QUOL9f9rSwvjVHXTgC0L+b3G+Yr4ib01wLXl+tkv1r6Jp5YFL8MD/O9MUzuedcVD8WMz6d4b8Qp31+W08H/Xk5N+jk1L2HqBb0C4DA+k+Zbx3WpvV7kxgbAhXWAI4ALc1toq2cC7yunfD/aWMl8f505ArhQ8bu3i1HXGFmNU1oRPGJpBkD7YnjW/ayp8xn0FiCeLdZZPg58pUxEcNSpH8nhvOYz5Tq/1WVkNj6z3gs8ETjUEcDhIM+xFkcAF0BsAFwAWmnS79NCFr4FW4bAC4Ezy6mVT3aRxGSE7mtrXlUuip/rGsBeT3VRe26BuBbmb8rzsDufG/cF4lrMDwKXA28FHtCYhWdfjLaqvgV8eJYAGLNQfW+M1r+z9ngf/KxMTvtaY5NfAmJENkZp7YvR9kWvAJj5/fBq4O+BftYx2iMb4doNgAvH9WkhC7frt2U8jSVm9z4LuKxHo5ihdXWZBRwzhiP0XQzEdVGdWcD2U7/ac78uTnHFbLvmEs/BjtswxCn5OO1lXwzHut+1xISnY8o1sF8v75MIHDHqdJX90S/jUF4XAwLxB+pryi2r4ilSEc7/EIhZ2b43hsJ8t5XEWZ6Y+BEBMO4AEQ9uiNshbSyXQGTd+/kdM5ojq7BWA2AO2aeF5Pzmax2nGOPWFreUF0a9xjMbY9JHJxDuDpxVbn8Ro1Fxq5g3dq3YfppPemE/jw/auA1D5z6A9sXCHDOtji73QYtJTzEiGLX+sbJC+yMjO1jbuD1VjPLFJRIRSOIyibeXyVKxJvtiMM9+Xx1nJs4rvxeiTed3xJPLPRiH4d7POvrd31a9zgDYqu5wZxRQQAEFFFBAgdELGABHb+wWFFBAAQUUUECBVgkYAFvVHe6MAgoooIACCigwegED4OiN3YICCiiggAIKKNAqAQNgq7rDnVFAAQUUUEABBUYvYAAcvbFbUEABBRRQQAEFWiVgAGxVd7gzCiiggAIKKKDA6AUMgKM3dgsKKKCAAgoooECrBAyAreoOd0YBBRRQQAEFFBi9gAFw9MZuQQEF6gm8oTwW6un1NumWtLMV0AAABA1JREFUFFBAgfETMACOX5+5xwpMq8BNjUc+bQFsDtzcePxT8xGBNYwOA/4YeGaNjbkNBRRQYJgCBsBharouBRSoJRDPe/494Cm1NthjO/Hw+Q8C71nEfXDTCiigwIIEDIALYrORAgosssBsAfB44KnAE8r+/StwBXB/IEYIbwReB1wDnA3sBlwJHAx8q7TZDDgK+DNg+/L9o4FPNY55OfBjYAfgBuB7JQjuU4Lpj4CXl1HKNcBvA/8BvBj4aVnPEcCrgfuVkcwIlC9ZZFc3r4ACUyJgAJySjvYwFZgwgbkC4O8DT2wEwEcBzwIuA44ETgE+A7y0hLjzga1KQIxmJwB/BDwP+DbwbOD9wCNL0IvXvAg4tITN+DoCYCyxnQiUbwVeWEJfBMlbgH8Bvg68AngY8DXgd4BvAlsCjwE+O2H95OEooEBLBQyALe0Yd0sBBeYUGCQARoj787K2ZWUU8PnAh8r3DgDOAe5bvo5RwucCn2jsQYS3fyvhMb59UQl0MYrYCYCxjgiXsewBfAWIEcHPl++tLqFwzzIiGEExQuTHgbi+0UUBBRSoJmAArEbthhRQYIgCgwTAGO07rmw7Jo7cCjwJuLR8byXwT8DScqr4J8AG4Pby8/icvAfwXmBVGa2L07gxihev7QTA2Kdzy9cPLaeZHwz8oHzvZcBrgF3K1zFaeDjwu+W1fwXEaKSLAgooMHIBA+DIid2AAgqMQGBUATBCYFzTF6FwttOxzynX7j2+cVxxCnjQANhpHtccHljC38PLaecRkLlKBRRQ4E4BA6DVoIAC4ygwqgAYFqeXU7cxIeMq4N7lWr2Y9BGnkz9QTuvGiF1n6RUAY1JJTP7oNQIYo4AxShijkL8ugfNiYOfG9YTj2C/uswIKjImAAXBMOsrdVECBuwj0GwBj5m6M5DVPAW8EnjzLKeDYSHwuxgzdmCSyI/BfwJeB1wLfAX5WrvH7fmOPvgu8aYBTwLsDf11mIcf2IiS+uYRLu1oBBRQYuYABcOTEbkABBSZIIG76HLOEY/auiwIKKDC2AgbAse06d1wBBRZBIG4xExNJYlawiwIKKDC2AgbAse06d1wBBRRQQAEFFFiYgAFwYW62UkABBRRQQAEFxlbAADi2XeeOK6CAAgoooIACCxMwAC7MzVYKKKCAAgoooMDYChgAx7br3HEFFFBAAQUUUGBhAgbAhbnZSgEFFFBAAQUUGFsBA+DYdp07roACCiiggAIKLEzAALgwN1spoIACCiiggAJjK2AAHNuuc8cVUEABBRRQQIGFCRgAF+ZmKwUUUEABBRRQYGwFDIBj23XuuAIKKKCAAgoosDABA+DC3GylgAIKKKCAAgqMrYABcGy7zh1XQAEFFFBAAQUWJvD/AffMbvyXniDCAAAAAElFTkSuQmCC\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "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,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4Xu3dDbildVnv8e9AzBgwTNHRTiCYGi/qKKmVclTMtyYDM2cQMw3LSM1BTcQa6cSbb6DCyaOohUFHMzFlvLLx6KRH05Sy0F5gLsES01FQM3G2Qs7QMOe6Z/9X87DcL2vve63/PM96vuu6YPaeve7n5fO/19q/+T8vawU+FFBAAQUUUEABBXolsKJXe+vOKqCAAgoooIACCmAAtAkUUEABBRRQQIGeCRgAezbg7q4CCiiggAIKKGAAtAcUUEABBRRQQIGeCRgAezbg7q4CCiiggAIKKGAAtAcUUEABBRRQQIGeCRgAezbg7q4CCiiggAIKKGAAtAcUUEABBRRQQIGeCRgAezbg7q4CCiiggAIKKGAAtAcUUEABBRRQQIGeCRgAezbg7q4CCiiggAIKKGAAtAcUUEABBRRQQIGeCRgAezbg7q4CCiiggAIKKGAAtAcUUEABBRRQQIGeCRgAezbg7q4CCiiggAIKKGAAtAcUUEABBRRQQIGeCRgAezbg7q4CCiiggAIKKGAAtAcUUEABBRRQQIGeCRgAezbg7q4CCiiggAIKKGAAtAcUUEABBRRQQIGeCRgAezbg7q4CCiiggAIKKGAAtAcUUEABBRRQQIGeCRgAezbg7q4CCiiggAIKKGAAtAcUUEABBRRQQIGeCRgAezbg7q4CCiiggAIKKGAAtAcUUEABBRRQQIGeCRgAezbg7q4CCiiggAIKKGAAtAcUUEABBRRQQIGeCRgAezbg7q4CCiiggAIKKGAAtAcUUEABBRRQQIGeCRgAezbg7q4CLRG4EzgfuLBsT3x9LvB9QPys7Y/h7d/f23te8TtwkQ15APBC4KHAA4v3YjX7e99cvwIKTEDAADgBVBepgAKLCvwU8GXg5vLMQYA5yAC4qN1cTzgCuCfwt4tUnw68HLgW+GHgRMAAuCxyixTotoABsNvj59YrMC0CBsD6IxlB8BwDYH1416hAGwQMgG0YBbdBgW4IHAO8BvgfwGHA14G/AZ5eZu0eDXwUOBV4EvDkEi7+vBx2/GZjN4cPoc4VAH8WeDdwZamP8u8vh46fChwJfAV4K/BqYE9Z/iHARcDPl1muHcA/AmcCnxsT9fD2j2rzFODnilG8//4f4CXAw4BLgBOALwBnAX+xhG0dHEI/YAk1BsAlYPlUBaZNwAA4bSPq/igwOYF/Bv69hKv4MwJYhJlnA/8JDALgduDDwLuACEYRzv4OeNwSAmAcqry8hL2oj0ccqvxL4Phy7uD1wMPLuW9vBF5anhd1pwAvA/4F+CHgEcB7RjhEOqrecAAc1eZfgc3AB4CTgN8Ffh94THGKQ+Lxdw8G7gU0Q/NC2zbXOYB/BITjfKHQADjqaPs8BaZQwAA4hYPqLikwAYEIUf9WZtW2zLP8QQCMcHNy4zm/BLwdeHyZIYwfLTQDeHY5T+15ZfZvsKhfBiLURHD6ZGP5cRgzLiCJc+C+AVwHbAViOZN6NLd/KTZ/CPx6Y6M+Dfw48Ejgr8vfx8UZMWP5rOI2yj7MFQBjZjTMVs2zAAPgKLI+R4EpFTAATunAulsKTEAgZtO+C/xemYmL75uPQQCMGcEIaoNHXNjxHyWkvar85XwB8DLgDOA0YDho/nE5/Byzis3HQ4BPNcLpFeXrOKQah1H/foQLS+K9sDlTFts3OKQ8F+Xw9o9qE4eur24s8B3lcHkcUm967QT+JzDwWmw4R70KuLkcA+Biqv5cgSkWMABO8eC6awqMWeBHyyHZOOz738q5aq8F3lLWMwiAMfsXs4DNR8we/imwcYEAGCEmZvDiHLifLqGxuYwIczGLONcjwloEzzinLs4BjFnBCFv3BW4F3gb8zhzLHCwrzjOMGbd4xLJiObG8+R7DAXBUm9j+jzQWGuuNQ+NHD60olv+KEppHGUYD4ChKPkcBBf5LwABoMyigwHIEHlQuqojZuieWQ67jmAGMw6EfAm4sy729sXHvBH6yBLu53rvi/Lrhc+aOKhdcXFwusojzAud6RACLUDt4RBD90hICYPOpC9kYAJfTbdYooMDYBQyAYyd1gQr0RmA1EFfYxsUXcbh1EAA/WC4OGUA8o8zAjXoO4LFlliwOq0a4vK0sKGbo/qDcwHipV/PGuXZxgUVcnTyOx2I3gp7LJmb+nuAM4Dj4XYYCCmQFDIBZQesV6IdAXJjw+nJlbwSzuCL3V4H15UrcOM+ueRXw/wOuAo4rhzLjKuDm4dvFbgMT5/nFFb9xODhuB/Od8qkVMTsYP4vAGRdKrAR+rHHbmThH8RrgfeVikKiLw8lxPt2LgbhaeByP5vYvxabmDGBccBJXAcc5mINH3EYnDuHHI27fE7elifMt4xEzqBGUfSigQA8EDIA9GGR3UYExCNy93AMwPjkirraNoBVX276y3PIlVhEBMGa5NpSLMH6hBMUIYy8aOjy7u5xPGBcixCPOYYvbn0SgG3wUXJy/F/cVjNvKrCshMH6+CfhF4N5ldvDz5YKR2JaojdvG/AxwnxIabyozh3GBybgeze1fis1cM4CPLbd8aW5bLD/OAQyXUR4Dv/govcEjzi+Mq4Cbfxe3lolQPdcFLoud9zjKdvgcBRToiMC0BcB4448T0ONNLv7l/zHgt8pHTg2GJH5BxC+vuG9Z7H+8EcYvtW2NMbugXIkYV+bFv4jjxPXmzzsyvG6mAlUFBgFwOORU3QhXpoACCiiwuMC0BcCYAYibvcbMxMHAm4H7l5uqNgNgXHUXMwtzPeJ8pvjEgDj3KGYW4l/WcRglzktqnpC+uK7PUKBfAoNDwMOHOful4N4qoIACHRCYtgA4TB4fq/QZ4PBysnr8PGYAF/oFFYeLLm2cKxTnOsXJ4/HRTHHPLh8KKDC3gDOAk+mMeA9a6BGHi30ooIACSxKY9gAYh3/j0wTiXKDBIwLgV8uJ0V8s9zCLO+bHIw75fqscEo4byw4e8akCMas4yU8WWNLA+WQFFOiFwGBWdb6djVNY4lzIhW5Z0wsod1IBBZYmMM0BMGb53luuUowrBweP+MzNuEow/tUc5yrFrF7cGyw+jzNObo830vuV+5ANauJqxhngOUO84XcE8O2lsftsBRRQYCSBOJVl+JNPhgvjM5GdBRyJ0ycpcBeBuF1THOFb6FN/ppZsWgNgfBB8fPZo3DcsrkBc6BGfIRpXDMZncS51BvDIoQtMprZR3DEFFFBAAQWmUCAmfr4yhfu16C5NYwCMm87Gvb7iY6A+vKjA7IfIxy0mHlGeO9c5gLeUe4gNnwMYgXHH9u3bOeyw5kd5jrBWnzJ2gXPOOYdXvWrUj04d++pdYEPAsWhXOzge7RkPx6IdYzEzM8NRR8UHBbGmHOFrx4ZV3IppC4Bx9e6F5aawn5zD8cHl1i9xPl+cCxhXA8fHS0UIHNwjLM7zi+XE7WQiDMbPnlluaDt8FfDeALhjxw4DYMWmnW9VZ511FpdeGtfv+NjfAo7F/h6Bu67f8WjPeDgW7RiLCIBr1kT2MwC2Y0TyWxGh7g5gZ1nU4D5/cUuXCIRxaPg15Vy/uA9gXATyJuDyoVWfDzwXiPMDrl3gPoAGwPyYjW0JvrGOjTK9IMciTTjWBTgeY+VMLcyxSPGNrdgAOHsjZB/LFzAALt9u7JVbt25l3bo4mu9jfws4Fvt7BO66fsejPePhWLRjLAyABsBsJxoAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLbcogFwyxbYtGl2NevXw+bN3f96ZgZWrIDVq6dnn6ZpfCaxL45591+3o/SF49yPcW72Ql/HfPfuGW64wfsAZkNQn+sXDYBr18K2bbNEq1bBznKHwq5/PRj0ru+H2z96bzrm0/Uanq/3Hed+jHNz/Ps55jPlHtDeCLrPIS6z7wbAKQq1hsGFw2A/f0kYBqbtH6++zr/3dd7P17YB0BtBZ+IfLBoAm4eAN2yAq6+eXWGXv24eMujyfnR9HGpuv2Pe/dftKP3iOPdjnJu90Ncx9xCw5wDm4t8IATC7AusVUEABBRRQYLwCXgRiAMx21KIzgNkVWK+AAgoooIAC4xUwABoAsx1lAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy23YADcsgU2bZpdxfr1sHlz97+emYEVK2D16unZp2kan0nsi2Pe/dftUvvCMe/HmPd5nM89d4anPW1NDHT8byYbBrpYv6KLG92ibV4wAK5dC9u2zW7tqlWwc+d0fD3wn6Z9cl8W7k3HfLpew6P0u2PejzHv6zgff/wMN9xgAGxRnurcphgApyTUjvILsc/P6esvCce8HyHIce7fOBsAPQScTZwjHwLesAGuvnp2dV3+unnIoMv70fVxqLn9jnn3X7dL7RfHvB9j3udxPu+8GU47zRnAbAjqc70XgfR59N13BRRQQIFOCngRiDOA2cY1AGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAKcvAL4aOBm4F/Ad4GPAbwFfbvTWg4A3AA8FvgVcDlww1Hvx/RlABLxPAxuBbXP0pwGw8ovW1SmggAIKKJAVMABOXwB8JfAe4DrgYODNwP2BB5dmORT4HHAFcCFwLPAB4HXA68tzXgqcCTwR+DxwHnB6ee7tQ01nAMy+Cq1XQAEFFFCgsoABcPoC4HALnQB8Bjgc2AE8C7gYOAK4szz5hcALgGPK9zcBlwJvLN8fCNwMnAW8wwBY+VXq6hRQQAEFFBizgAFw+gNgHP59HnCf0jsR7O5XZvcG7XQi8AlgDXBAOSwcf/epRr9tLbOKZxsAx/wqdHEKKKCAAgpUFjAATncAfDzwXmA98KHSW28FDgGe3ui148v5fUeVAPilEhJvbDznKmAGeI4BsPKr1NUpoIACCigwZgED4PQGwFOAt5dDvu9r9I0zgGN+Ebk4BRRQQAEFuiZgAJzOAPiMcv7eU4EPDzVlXMzxmqFzAF9ULvpY6BzAW4AXz3cO4MaNG1m5cuXeVa1bt27vfz4UUEABBRRQoD0CW7duJf6Lx65du7jsssviyzj9K47w9e6xYsr2OK7ejat7nwR8co59i6uA49BuXAUcVwxH6Hs/cEnjKuA4zy+WE7eTiQtCzgWeCRwHeBXwlDWMu6OAAgoo0D8BZwCnbwYwruy9A9hZ2jkC7p5y0ccgEK4F3lTuAxhXBsetYl4+1P7nA88FVgPXeh/A/r05uMcKKKCAAtMrYACcvgBYu1u9D2BtcdengAIKKKBAUsAAaABMttDeTwrZsWPHDg47LL70oYACCiiggAJtFzAAGgCzPWoAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAATl8AfBqwETgBOBQ4CLiz0Vfx9XeB/2R23/cAJwLbGs+5ADgDiHD36bK85s+bbWoArPyidXUKKKCAAgpkBQyA0xcAnwAcDhwMvHWeAPg44KPzNM9LgTOBJwKfB84DTgeOBW6fo8YAmH0VWq+AAgoooEBlAQPg9AXAQQs9GvjIPAHw8eVnc7XbTcClwBvLDw8EbgbOAt5hAKz8CnV1CiiggAIKTEDAANjPAPjVEgy/CLylzBRGe8Vs3rfKIeFPNfptK3AdcLYBcAKvQhepgAIKKKBAZQEDYP8C4GOAa4DdQBwujlm9lwG/D9wT+BJwP+DGRi9eBcwAzzEAVn6FujoFFFBAAQUmIGAA7F8AHG6jc4GfAR7pDOAEXmEuUgEFFFBAgRYKGAANgBEA1wGPKP051zmAtwAvXugcwI0bN7Jy5cq9i1i3bt3e/3wooIACCiigQHsEtm7dSvwXj127dnHZZZfFl2vKUb72bGilLYlboUzT44Byfl9cBPIBYHU53LsL+PFy65c4ny9uBxNXA78TiBC4twvKeX5xFfDJQITB+NkzgeO8Cnia2sR9UUABBRTos4AzgNM3A/gs4Mpyf7/o7cG9/uLcv7jI4zXlXL+4D2BcBPIm4PKhF8H5wHNLeLzW+wD2+S3CfVdAAQUUmEYBA+D0BcDafep9AGuLuz4FFFBAAQWSAgZAA2CyhfbOKu7YsWMHhx0WX/pQQAEFFFBAgbYLGAANgNkeNQBmBa1XQAEFFFCgsoAB0ACYbTkDYFbQegUUUEABBSoLGAANgNmWMwBmBa1XQAEFFFCgsoABcP8HwC80rthdaPjvU7k3Rl3dggFwyxbYtGl2URddNPvn4Pv162Hz5tm/69rXMzOwYgWsXt29be+idxu22THv5mt1qb3jOPdjnJt90dcx3717hhtuiFsAeh/AUQPPuJ/3a40FHgU8H3gbEMHw3sAvl1u1XDDuFY9peQsGwLVrYdu22TU94AGzfw6+X7UKdu6c/bsufj3w6+K2u83L6zvHvJuv1aX2u+Pcj3Fu9kU/xzw+4dUAOKYslF7Mh4HfBf66saSHA68AHp9e+mQWMG8AjNm/U0/dF/KOPhq+9rV938cM2p49sxvVxa8HnF3cdrd5eX3nmHfztbrUfp/rrXKpy/D5y3uN7S+3fr62DYBt+iSQGI0fLJ/cMejHA4Fby02cJxPhckudNwA2Z//iX1r3uAds396tN4X53oz6+WbRj1/+jrnjPPyWuL9Cieut8/uir+/nK1fOsGuXM4C5CDS+6n8CXg/8YWORzy6fw/vA8a1mrEsaKQB6+Hd5hxyXeujK50/OuZ+HiPp3KNDDgv0b876+to8/3nMA2zQDuA74MyA+qzfOAfxRIILfU4APjjW2jW9hCx4Cfv7z4etfn71YIv41GyfbDr6+NeY1gYMPhttv3/d1nBfY5ufs3j176PqAA+DQQ2cPabd5e7to3LZtdsy78/rM9s6uXfDtb8/ub7zG430r3qN8nc9adPV9e76+WLkSbrtt9j39bneDO+7oR6/HEblLLpnhtNOcARxfHMovKS78eHr5vN4vA+8sYTC/5MksYeSLQOb6V9Z8mzTXv8KHn9um57RpW8LJ7Vm42cfhM45lOFb7zgnuwnuBYz7511XbXhPTOuZxVO6aa2ZYs8YAOJlo1I+lGgANXP91NXcXfomP6xfMtP5icL/mf+PWxgA4l0AX+8IAODuS+/sQcNzm5e2lqeJ8v/keV7Q0T450H8DmfZY2bICrr549HDw4ybn5ZxxObfNz/uM/Zg8RxXbGoaG2b28Xjdu2zY757CGywb0v2/z6zPbOlVfuO20lXt992e8uvyc75kvv04svhkc9yhnA/R0ArwfWlnAX5/3N9YibpXTyRtAtDa1ulgIKKKCAAr0W8JNA9v8MYNcb0I+C6/oIuv0KKKCAAr0TMAAaALNNbwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDYPsC4KOBnwJWD/XCuZV7Y9TVGQBHlfJ5CiiggAIKtETAANiuAHghsAn4B+C2Ro/EfQAf25KeGd4MA2BLB8bNUkABBRRQYD4BA2C7AuAtwFOAv+lQyxoAOzRYbqoCCiiggAIhYABsVwD8BnB3IGb8uvIwAHZlpNxOBRRQQAEFioABsF0B8A3Ax4F3d6hDDYAdGiw3VQEFFFBAAWcAZ3tgf38WcLMT3wn8AvBJ4OahFj29pS1rAGzpwLhZCiiggAIKzCfgDGC7AuCVC7Tqr7a0jQ2ALR0YN0sBBRRQQAED4Pw90KYZwC52qgGwi6PmNiuggAIK9FrAGcB2zQAODkn/JHA08EXg2pZfFGIA7PVbiDuvgAIKKNBFAQNguwLgjwDvAx4C3Ar8IPAZ4MlznBPYln4zALZlJNwOBRRQQAEFRhQwALYrAF5VLkr5DeCbwOHAZcABwNNGHNPaTzMA1hZ3fQoooIACCiQFDIDtCoBxI+jj4v6MjXFdA9wI/PfkWE+q3AA4KVmXq4ACCiigwIQEDIDtCoBfBX4M+E5jvOMzgf8F+OEJ9UB2sQbArKD1CvI3OGkAACAASURBVCiggAIKVBYwALYrAL4LuAN4fpkFjNm/NwKrgNMq98aoqzMAjirl8xRQQAEFFGiJgAGwXQHwCOD9wAPLOYBxEcj1wCnAV1rSM8ObMW8A3LIFNm2affr69bB583R9PTMDK1bA6tXTuX/TOm5L3S/Hebpet6OMv2PumE/b76vhvr/oIjjppBnWrIl5JuJ/zVPPWho3xr9ZbbsPYFzw8VPAUcB24G+BO8e/22Nb4pwBMMLfqafCzp2z6znoILgj5jaBAw+E3bun4+uB4jTtk/vyvb3pOE/X63aUHnfMHfNp+33V7Pujj4brrjMAtiUAriy3fPkJ4Ltji2eTX9CcAXDtWti2bd/KY6Zsz57Jb4xrUEABBRRQQIGFBVatgq9/3QDYlgAYoxWHee8N7OpQ8y4aAKPR4jDpN74xu1cHHAB3ljnNrn89GKeu74fbv3BvOs7T9bodpd8dc8d82n5fNfv+qKPg+usNgG0KgBcCcdD0lV0PgM3z/y6+eHb2b3A+4IYNcPXVs3vY5a+b5wl1eT+6Pg6T3n7Hufuv1aX2iGPumE/D76iF+j5+Lz/qUQbANgXAvwIeBnytfAxc89y/k1oaCr0KuKUD42YpoIACCigwn4BXAbfrKuDzFmjVC1raxgbAlg6Mm6WAAgoooIABcP4eaNMMYBc71QDYxVFzmxVQQAEFei3gDGC7ZgDjFjDzPdp6KxgDYK/fQtx5BRRQQIEuChgA2xUAI+TNd7OUA1vaYAbAlg6Mm6WAAgoooICHgLtxCPjRQ5t5JPAS4HLgLS1tYwNgSwfGzVJAAQUUUMAA2I0AONdW3gd4G/DIlraxAbClA+NmKaCAAgooYADsbgCMQ7+3AhG02vgwALZxVNwmBRRQQAEFFhDwHMB2nQMYs33NxyHAGcBjgAe1tJMNgC0dGDdLAQUUUEABZwC7MQM4fBFI3KLmJuB04JqWtrEBsKUD42YpoIACCihgAOxGALzX0GZ+G/hmy9vXANjyAXLzFFBAAQUUGBbwEHC7DgF3sUMNgF0cNbdZAQUUUKDXAgbAdgXAuBH0bwPPBu4BrAHWAUeXW8G0sVkNgG0cFbdJAQUUUECBBQQMgO0KgBcCTwJeC7wJ+AHgvsC7gYe0tJMNgC0dGDdLAQUUUECB+QQMgO0KgF8ATgK2l3P/DgdiVvAbQHzdxocBsI2j4jYpoIACCijgDOCCPRBX2rblEUHv7uXj4OLijwh9BwFfKYeE27Kdze0wALZxVNwmBRRQQAEFDICdCYAfLYd+45DvIACeCvx6ORewjc1sAGzjqLhNCiiggAIKGAA7EwAfBnwYeB/wFOBPgAiATwD+rqWdbABs6cC4WQoooIACCswn4DmA7ToHMMbpOOD5wDHAV8uM4LUtbmEDYIsHx01TQAEFFFBgLgEDYPsCYNc61QDYtRFzexVQQAEFei9gAGxfAIzPA34osHqoO69oabcaAFs6MG6WAgoooIACHgKevwfadBXw84A3lgtAbmts8h4ggmEbHwbANo6K26SAAgoooMACAs4AtmsG8IvAi4HNHepaA2CHBstNVUABBRRQIAQMgO0KgN8qn/7Rpe40AHZptNxWBRRQQAEFDIB7e6BNh4DfAfwB8LEOdacBsEOD5aYqoIACCijgDOBsD7QpAP4e8CzgauDmoRY9t6UtawBs6cC4WQoooIACCswn4CHgdgXA+CSQuR5xEchjW9rGBsCWDoybpYACCiiggAFw/h5o0wxgFzv1ewLgli2wadPsrlx00eyfg+/Xr4fN5RKXLn89MwMrVsDq1dDl/Yixcftne3QxB8d8NKdRLNv8HMe5H+Pc7MG+jvnu3TPccMOaoIj/zXQxgGS3edoC4NOAjcAJwKHAQcCdDaQHAW8o9xqMi04uBy4YQozvzwAi3H26LG/bPNB3CYAR/k49FXbunH32qlWzfw6+j9C0J+Yz49h7x78eeHR9P9z+hftxrr7XbDpew/ONo6/t7r8/L/U12s8xj8xnAMyGyDbVx+cGHw4cDLx1KABGIPwcEDeVvhA4FvgA8Drg9WUnXgqcCTwR+DxwHnB6ee7tc+zoXQLg2rWwbb6o2CYlt0UBBRRQQIFeCxgAp20GcNDOjwY+MhQA4wKTi4EjGrOCLwReUD57OGpvAi4tN6SO7w8sF6ScBcRVysOPvQHwXe/awWmnHYYBsNfvJu68AgoooEAHBGKG9MgjZ/jyl50B7MBwLXkT5wqAEezuV2b3Bgs8EfhEmQc+AIjDwvF3n2qscStwHXD2fAHw+ON38NnPHkYcAo5zqe64Y8nba4ECCiiggAIKVBK45z0NgH2aAYxDwocAT2/01/FAHLQ9CogA+KUSEm9sPOeqcoLocxYLgPHzu91t3zl/lfrY1SiggAIKKKDAEgRWrpxh1y5nAJdA1pmnVp0BPPnkjRx77Mq9OG9/+zq+8Y11e78+4AC4s1yCMk1fD7pgmvbJfZkd1fkcHPPpfT075v0b2+aY9+m1vWLFVvbsiYN6cMghu7jttsviS68C7ky0G21D5wqAcTHHa4bOAXxRuejjmLLYuc4BvKV8RvG85wDu2LGDww6L0wHZexh4cNuXDRvg6ritNTAtXzdvGTAt+zRN4zOJfXHMp+s1PEqPOOb9GPM+j/N5581w2mnOAI4WqbrxrDiMG7d+iQAYV/iuBnYDu8rh3zi0G1cBv7Jc+PF+4JLGVcBxnl9cBXxyuSAkPoHkmcBxwKJXAXeDyK1UQAEFFFCg3wJ+Eki7PglkHN0YV/peCZS77e39qLv4+jHAx4G1wJvKfQB3AG8GXj604vOB55bweO1S7gM4jh1wGQoooIACCigwWQED4PQFwMl2zPcu3Y+Cqy3u+hRQQAEFFEgKGAANgMkW2vtpITua5wBmF2i9AgoooIACCkxWwABoAMx2mAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwABoAMy2nAEwK2i9AgoooIAClQUMgAbAbMsZALOC1iuggAIKKFBZwADYzwB4HvC7wO3M7v8e4M+BZ5T+exDwBuChwLeAy4EL5ulNA2DlF62rU0ABBRRQICtgAOxvAHwccNIcDXQo8DngCuBC4FjgA8DrgNfP8XwDYPZVaL0CCiiggAKVBQyABsDhlnsWcDFwBHBn+eELgRcAxxgAK79CXZ0CCiiggAITEDAA9jcAnl0OAcdh4GuA3wH+FbgUuB/wxEa/nQh8AlgDfGeoD50BnMAL00UqoIACCigwSQEDYD8D4P2BbwPbgR8BXgs8DDgB+N/AIcDTG413PLANOAq42QA4yZeky1ZAAQUUUGDyAgbAfgbA4c5aCewAngT83HJmADdu3MjKlbEYWLdu3d7/fCiggAIKKKBAewS2bt1K/BePXbt2cdlll8WXcXRvpj1bWW9L4irYvj8GAfDny4zga4bOAXwRcKbnAPa9Tdx/BRRQQIFpEXAGsJ8zgE8FPgL8O/DD5RDwI4EHltvC3FiuAn5lCX3vBy7xKuBpedm7HwoooIACfRcwAPYzAP4Z8PByrt+twMfLfQFvKi+ItcCbyn0A49Dwm4GXz/Ni8SKQvr+LuP8KKKCAAp0TMAD2MwCOs1ENgOPUdFkKKKCAAgpUEDAAGgCzbWYAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLWcAzApar4ACCiigQGUBA6ABMNtyBsCsoPUKKKCAAgpUFjAAGgCzLbdgANyyBTZtml3FRRfN/tn8/pRToGvPGd6HadinLo5DzW12zKfrNTxK7zjmjvm0v7fv3j3DDTesid2M/81kw0AX61d0caNbtM0LBsC1a2HbttmtXbUK7nEP2L593/fvec9sIOzSc4b3YRr2qYvjUHObHfPpeg2P0juOuWM+/e/tkfkMgC3KU53blJED4CAE7ty5bx8f8IDZrwcBsAvPiSDb3Ifh77u4T27zwj3omO97zfal3x1zx3z6e90A6AxgLnMuegj41FP3Baajj4avfW3f9xE84tBwl54zvA/TsE9dHIea2+yYT9dreJTeccwd82l+b1+9GjwE7DmAufgHi14E0jzH7+KLYc+efecBxvcnn3zX8wC78JzhfZiGferiONTcZsd89q2iC6/PcfWFY+6YT+t7e7xGvAjEADjxAJhdgfUKKKCAAgooMF4BA6ABMNtRi84AZldgvQIKKKCAAgqMV8AAaADMdtS8AXD49i6j3Hqhi88JwGm71U0Xx2FS27x+PWzePPsyifNVJ7WeablN0jT4OOb96/c+jrnnABoAJxYAm7eAiZOur78ehv8uVj64Arirzxneh2nYp2kdq+XsV/NKwK726HL2e/jK8D7tu2M++2vBMb/r76fpe2/3KmCvAs5FwHlnAA2A/XoDndbAYBgwDBigpz0Izd6ndnB7r/4EXwOgAXBCAXD46t9xXZnXtuV4peBsA03r1aEbNsDVV+/bx7b1n9tz17sKjOOqTce8f/3exzH3ELCHgHPxb4TbwGRXYL0CCiiggAIKjFfAi0AMgNmO8irgrKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAF2q5C4AzgAh5nwY2AtuGCgyAlV+0rk4BBRRQQIGsgAHQADhfD70UOBN4IvB54DzgdOBY4PZGkQEw+yocY/3WrVtZt27dGJfoopYr4FgsV24ydY7HZFyXs1THYjlq468xABoA5+uqm4BLgTeWJxwI3AycBbzDADj+F+M4lnjWWWdx6aUxbD72t4Bjsb9H4K7rdzzaMx6ORTvGwgBoAJyrE2NW71vAicCnGk/YClwHnG0AbMcLeHgrfGNtz7g4Fu0Zi9gSx6M94+FYtGMsDIAGwLk68Z7Al4D7ATc2nnAVMAM8ZzgAbt++ncMOi9zoY38KnHPOObzqVa/an5vguouAY9GuVnA82jMejkU7xiIC4FFHHRUbs6b8bm/HhlXcihUV19WVVS1lBvBI4Mtd2TG3UwEFFFBAAQXuIhCTPl/po4kBcO5Rn+scwFuAFw+dAxh+RwDf7mPzuM8KKKCAAgp0WGB1Ob9/T4f3YdmbbgCcmy7O84urgE8GIgyeCzwTOG7oKuBlw1uogAIKKKCAAgrsLwED4Pzy5wPPBeJfCNfOcx/A/TVurlcBBRRQQAEFFFi2gAFw2XQWKqCAAgoooIAC3RQwAObGbZRPC8mtod/Vry6H4e8FfAf4GPBbQxfePAh4A/DQcvuey4EYl+bDcRp/H70XeDLweOAjZfGOxfidF1ti3K7qFcBPArvLpxU90vFYjG3sP//Bcu/YnwVWAZ8FXgZ83LEYu3VzgU8rR+dOAA4FDgLubDxhHO9Joyxjojs5qYUbAJcvO+qnhSx/DVa+EnhPuf/iwcCbgfsDDy408YL/HHAFcGH5pJYPAK8DXl+e4ziNv4/iU3F+CXhC+S8CoGMxfufFlhjh7/+W85XjdXJH+YfQ3zkei9GN/ed/DMQ9RZ4C3FouGIzTiI4G/tP3qbF7DxYY70GHA/H74a1DAXAc70mjLGNiOzfpBRsAly886qeFLH8NVg4LxL/yPlNe8DuAZwEXlyuxB//qeyHwAuCYUuw4jbeP4pYJnwBilinulzmYAXQsxus8ytJidiluVh//yBl+OB6jCI7vOfEhAX9QjkbEUg8pd4d4WLmn7Gt8nxof9hxLenQ5EtGcARzHa2CUZUx0xya5cAPg8nSXcq/A5a3BqrkE4vDv84D7lB/G577FDbvjM5sHj5gViYASN/c8YAmf6qL4aALxiTh/CvxhOdQyCICOxWh+43rW95eAcQnw08B9gS8AcdrE5nI40tfGuLQXX06cZvJY4KnAv5ePDX02EIcPY0wci8UNM8+YKwCO4z1psWXEqUmdfRgAlzd0S/m0kOWtwaphgQgacd7ZeuBD5Ycx5R//0n5648nHl/Og4nBMBMBRP9VF8cUFnl/O+1tXnhqzroMA6Fgs7jfOZ8RN6LcDXyvnyf5DGZv4xKL4ZXiGr41xci+6rLuV96d4bcQh32+Ww8F/XQ5N+j61KGHqCXMFwHG8Jy22jJtTW72fiw2AyxsAZwCX57bcqlOAt5dDvu9rLGSxf505A7hc8e+ti1nXmFmNQ1oRPOLRDICOxfisR1nS4D3oIuCcRsEHgb8vFyI46zSK5Hie81flPL+zysxsvGe9DTgJ+BVnAMeDvMBSnAFcBrEBcBlopWTUTwtZ/hqsDIFnAG8sh1Y+PEQSFyMMn1vzonJS/ELnAM71qS5qLywQ58L8fvnMzMH7xg8BcS7mu4BrgNcCP9K4Cs+xmGxX/TPw7nkCYFyF6mtjsv6Dpcfr4N/KxWn/2Fjlp4GYkY1ZWsdismMxVwDM/H74TeBPgFGWMdk9m+DSDYDLx/XTQpZvN2plfBpLXN37JOCTcxTFFVo3lquA44rhCH3vB+K8qMFVwI7TqNoLPy8OccXVds1HfA523IYhDsnHYS/HYjzWoy4lLnjaVM6B/afyOonAEbNONzgeozKO5XkxIRD/QH1JuWVVfIpUhPOfA+KqbF8bY2H+noXEUZ648CMCYNwBIj64IW6HtKucApF1H+V3zGT2rMJSDYA5ZD8tJOe3WHUcYoxbW+wsT4x+jc9sjIs+BoFwLfCmcvuLmI2KW8W8fGjBjtNi0sv7ebzRxm0YBvcBdCyW55ip+u1yH7S46ClmBKPXt5QFOh4Z2aXVxu2pYpYvTpGIQBKnSfxeuVgqluRYLM1z1GfHkYkry++FqBn8jnhMuQfjONxHWcao29uq5xkAWzUcbowCCiiggAIKKDB5AQPg5I1dgwIKKKCAAgoo0CoBA2CrhsONUUABBRRQQAEFJi9gAJy8sWtQQAEFFFBAAQVaJWAAbNVwuDEKKKCAAgoooMDkBQyAkzd2DQoooIACCiigQKsEDICtGg43RgEFFFBAAQUUmLyAAXDyxq5BAQUUUEABBRRolYABsFXD4cYooIACCiiggAKTFzAATt7YNSigQD2Bl5WPhfrZeqt0TQoooED3BAyA3Rszt1iBvgp8u/GRT6uAA4HbGx//1PyIwBpGZwC/AJxSY2WuQwEFFBingAFwnJouSwEFagnE5z0/AnhsrRXOsZ748Pl3AX+0H7fBVSuggALLEjAALovNIgUU2M8C8wXA84DHA48q2/dR4DrgHkDMEH4LeCnwOeDNwAOAbcDpwD+XmgOAFwO/BhxR/v63gY809nkNcAtwJHAr8IUSBE8swfQrwPPKLOWlwI8CfwP8MvD1spwzgd8E7l5mMiNQPns/u7p6BRToiYABsCcD7W4qMGUCCwXAxwEnNQLgCcCTgE8CLwReBfwV8JwS4q4CDikBMcrOB34eOA34F+DJwDuAB5agF895JvArJWzG9xEA4xHriUD5WuAZJfRFkNwJ/AXwT8BvAD8G/CPwE8BngYOBhwCfmLJxcncUUKClAgbAlg6Mm6WAAgsKLCUARoj79bK0w8os4C8Cf1r+bj1wOfBD5fuYJXwq8KHGFkR4+8sSHuOv31sCXcwiDgJgLCPCZTweBPw9EDOCf1v+7qwSCh9aZgQjKEaI/CAQ5zf6UEABBaoJGACrUbsiBRQYo8BSAmDM9p1b1h0XjtwB/DTw8fJ364A/B1aWQ8VfBWaAO8vP433y+4C3ARvLbF0cxo1ZvHjuIADGNl1Rvr9vOcx8b+BL5e+eC7wEOLZ8H7OFzwceXp77v4CYjfShgAIKTFzAADhxYleggAITEJhUAIwQGOf0RSic73DsqeXcvUc29isOAS81AA7K45zDDSX8HVcOO0+AzEUqoIAC+wQMgHaDAgp0UWBSATAsLimHbuOCjBuA7y/n6sVFH3E4+Z3lsG7M2A0ecwXAuKgkLv6YawYwZgFjljBmIW8rgfP9wDGN8wm7OC5uswIKdETAANiRgXIzFVDgLgKjBsC4cjdm8pqHgHcBj5nnEHCsJN4X4wrduEjknsB3gc8AZwOfB/6tnOP3xcYW3QS8YgmHgNcCbylXIcf6IiS+uoRLh1oBBRSYuIABcOLErkABBaZIIG76HFcJx9W7PhRQQIHOGg3WlgAAAQJJREFUChgAOzt0brgCCuwHgbjFTFxIElcF+1BAAQU6K2AA7OzQueEKKKCAAgoooMDyBAyAy3OzSgEFFFBAAQUU6KyAAbCzQ+eGK6CAAgoooIACyxMwAC7PzSoFFFBAAQUUUKCzAgbAzg6dG66AAgoooIACCixPwAC4PDerFFBAAQUUUECBzgoYADs7dG64AgoooIACCiiwPAED4PLcrFJAAQUUUEABBTorYADs7NC54QoooIACCiigwPIEDIDLc7NKAQUUUEABBRTorIABsLND54YroIACCiiggALLEzAALs/NKgUUUEABBRRQoLMCBsDODp0broACCiiggAIKLE/g/wNSh1R0WoiW/AAAAABJRU5ErkJggg==\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "No spikes on lsm_i: 2\n",
      "No spikes on lsm_i: 3\n",
      "No spikes on lsm_i: 4\n"
     ]
    }
   ],
   "source": [
    "for lsm_i in range(number_of_liquids):\n",
    "    if len(times[lsm_i]):\n",
    "        plt.figure()\n",
    "        plt.plot(times[lsm_i],spikes[lsm_i],'.')\n",
    "        plt.ylabel('neuron id')\n",
    "        plt.xlabel('Time/ms')\n",
    "        plt.title('spikes - lsm_i:'+str(lsm_i))\n",
    "        plt.show()\n",
    "    else:\n",
    "        print \"No spikes on lsm_i:\",lsm_i"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "hi, stdout\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "hi, stderr\n"
     ]
    }
   ],
   "source": [
    "#\n",
    "# All messages generated by C code will not be printed inside the notebook!\n",
    "#\n",
    "# Proof:\n",
    "\n",
    "import sys\n",
    "import ctypes\n",
    "libc = ctypes.CDLL(None)\n",
    "\n",
    "from __future__ import print_function\n",
    "print('hi, stdout', file=sys.stdout)\n",
    "print('hi, stderr', file=sys.stderr)\n",
    "ret_value = libc.puts(b\"hello\") # the C code prints only on the terminal..."
   ]
  },
  {
   "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
}