500K4L-rprop.ipynb

{
 "metadata": {
  "name": "",
  "signature": "sha256:aca40a0bb389e6f76f159d965a7c3f33b408c82da0f89b0cacdf79dda6167f83"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import theanets\n",
      "import logging\n",
      "from sklearn.datasets import make_classification\n",
      "from sklearn.metrics import confusion_matrix\n",
      "import matplotlib.pyplot as plt\n",
      "import time\n",
      "import numpy as np\n",
      "from theanets.trainer import SupervisedPretrainer\n",
      "import math\n",
      "from random import uniform\n",
      "from math import cos,tan,sqrt\n",
      "from numpy import float32,int32\n",
      "from sklearn.cluster import KMeans\n",
      "import pickle\n",
      "\n",
      "\n",
      "nlayers = 4\n",
      "module_size=20 #cm\n",
      "module_cell = 10 #per edge\n",
      "distance_between_layers = 5 #cm\n",
      "ncells=module_cell*module_cell*nlayers\n",
      "nhidden=ncells\n",
      "nhidden_layers = 10\n",
      "layers = (ncells,) + tuple([nhidden for i in range(nhidden_layers)])\n",
      "#net = theanets.Regressor((ncells,ncells,ncells,ncells,ncells,ncells,ncells,ncells,ncells,ncells,ncells,(ncells,'logistic')))\n",
      "trainer = theanets.Experiment(theanets.feedforward.Regressor(layers=(ncells,ncells,ncells,ncells,ncells,(ncells,'logistic'))))\n",
      "maximum_hight = nlayers * distance_between_layers\n",
      "theta_max = math.atan(0.5 * module_size / maximum_hight)\n",
      "theta_min = -1 * theta_max\n",
      "y = []\n",
      "y.append(0)\n",
      "y_prime = []\n",
      "y_prime.append(0)\n",
      "z = []\n",
      "z.append(0)\n",
      "x = []\n",
      "x.append(0)\n",
      "for i in range (nlayers) :\n",
      "    y.append((i+1) * distance_between_layers)\n",
      "    y_prime.append((i+1) * distance_between_layers)\n",
      "    z.append(0)\n",
      "    x.append(0)\n",
      "spent = []\n",
      "list_error = []\n",
      "color_list = ['bo','yo','ro','co','mo','go','ko','wo']\n",
      "min_ntracks = 4\n",
      "max_ntracks = 4\n",
      "nsamples = 500000\n",
      "#nsamples = 2 * ncells * ncells * (nhidden_layers + 1)\n",
      "#print nsamples\n",
      "training_ratio = 0.9\n",
      "ds_input = np.zeros((nsamples,ncells),dtype = float32)\n",
      "ds_output = np.zeros((nsamples,ncells),dtype = float32)\n",
      "cut = int(nsamples * training_ratio)\n",
      "train = ds_input[:cut], ds_output[:cut]\n",
      "valid = ds_input[cut:], ds_output[cut:]\n",
      "\n",
      "def addSample (sample_number = 1,sampleInput = np.zeros (ncells,dtype = float32),sampleOutput = np.zeros (ncells,dtype = float32)) :\n",
      "    for i in range (ncells) :\n",
      "        ds_input[sample_number][i] = sampleInput[i]\n",
      "        ds_output[sample_number][i] = sampleOutput[i]\n",
      "    \n",
      "def sample_creater(sample_number = 1,display = False) :                  # creates a sample with a random number of tracks\n",
      "    sampleInput = np.zeros (ncells,dtype = float32)\n",
      "    sampleOutput = np.zeros (ncells,dtype = float32)\n",
      "    ntracks = int( uniform(min_ntracks, max_ntracks ) ) #should be a better model than flat\n",
      "    sample_matrix = [[0 for t in range(nlayers + 1)] for t in range(ntracks)]\n",
      "# first argument of each row is track's radius and the other \"nlayers\" ones are hit indices of that track\n",
      "# then we would be able to sort these hits by their radius\n",
      "    \n",
      "# theta is the angle between the circle plane and xy plane. we randomly choose two points in the second and 4th layer and then fit a circle through these two points and also the coordinate origin\n",
      "# as y is fixed for each layer and z can be derived easily by y*tan(theta) , to specify a point in a layer we only need x of it\n",
      "    track = 0\n",
      "    while track < ntracks :\n",
      "        theta = uniform (theta_min,theta_max)\n",
      "        x[2] = uniform (-1*module_size/2.0,module_size/2.0)\n",
      "        x[4] = uniform (-1*module_size/2.0,module_size/2.0)\n",
      "        for i in range (4) :\n",
      "            y_prime[i+1] = y[i+1] / cos(theta)\n",
      "            z[i+1] = y[i+1] * tan(theta)\n",
      "        # center point of our circle\n",
      "        x_c = 0.5 * (((x[2]**2 + y_prime[2]**2) * y_prime[4]/y_prime[2] ) - x[4]**2 - y_prime[4]**2 )/ (x[2] * y_prime[4]/y_prime[2] - x[4])\n",
      "        y_c = 0.5 * (x[2]**2 + y_prime[2]**2 - 2 * x[2] * x_c) / y_prime[2]\n",
      "        R = sqrt(x_c**2 + y_c**2)\n",
      "# y_prime is the y component when we place our coordinate on the circle plane and as the layers are parallel, it would be y/cos(theta) for each layer\n",
      "        if x[2] > 0.5 * x[4] :\n",
      "            x[1] = x_c + sqrt(R**2 - (y_prime[1] - y_c)**2)\n",
      "            x[3] = x_c + sqrt(R**2 - (y_prime[3] - y_c)**2)\n",
      "        else :\n",
      "            x[1] = x_c - sqrt(R**2 - (y_prime[1] - y_c)**2)\n",
      "            x[3] = x_c - sqrt(R**2 - (y_prime[3] - y_c)**2)\n",
      "        sample_matrix [track][0] = R\n",
      "        x_edge = (-1) * module_size / 2.0\n",
      "        z_edge = (-1) * module_size / 2.0\n",
      "        \n",
      "        track_hits = np.zeros(nlayers,dtype = float32)\n",
      "        no_good = False\n",
      "        for layer in range (1,nlayers+1) :\n",
      "            # local indices on the module\n",
      "            index_hit1 = int( (x[layer] - x_edge ) / module_size * module_cell )\n",
      "            index_hit2 = int( (z[layer] - z_edge ) / module_size * module_cell )\n",
      "            # transform into the global index\n",
      "            hit = module_cell**2 * (layer - 1) + index_hit2 * module_cell + index_hit1\n",
      "            track_hits[layer-1] = hit\n",
      "            if sampleInput[hit] == 1:\n",
      "                    #discard this track\n",
      "                    no_good = True\n",
      "        if no_good:\n",
      "            #print \"discarding an overlaping track\"\n",
      "            continue\n",
      "        else:\n",
      "            for (layer,hit) in enumerate(track_hits):\n",
      "                sampleInput[int(hit)] = 1\n",
      "                sample_matrix [track][layer+1] = int(hit)\n",
      "        if display:\n",
      "            for i in range (0,5) :\n",
      "                plt.figure(1)\n",
      "                plt.plot(x[i],y_prime[i],color_list[track])\n",
      "                axis([-module_size/2. + 1 , module_size/2. + 1, 0, maximum_hight/cos(theta) + 1])\n",
      "                plt.suptitle('x - y_prime (perpendicular view of the circle plane)')\n",
      "                plt.xlabel('x')\n",
      "                plt.ylabel('y_prime')\n",
      "                plt.figure(2)\n",
      "                plt.plot(x[i],z[i],color_list[track])\n",
      "                axis([-module_size/2. + 1, module_size/2. + 1,-module_size/2. + 1, module_size/2. + 1])\n",
      "                plt.suptitle('x - z')\n",
      "                plt.xlabel('x')\n",
      "                plt.ylabel('z')\n",
      "                plt.figure(3)\n",
      "                plt.plot(z[i],y[i],color_list[track])\n",
      "                axis([-module_size/2. + 1, module_size/2. + 1, 0, maximum_hight/cos(theta) + 1])\n",
      "                plt.suptitle('z - y')\n",
      "                plt.xlabel('z')\n",
      "                plt.ylabel('y')\n",
      "        track += 1\n",
      "    sample_matrix.sort(key=lambda x: x[0])\n",
      "    for i in range (ntracks) :\n",
      "        for j in range(1,nlayers+1) :\n",
      "            sampleOutput[sample_matrix[i][j]] = (i+1)/float(ntracks)\n",
      "    if display == False :\n",
      "        addSample(sample_number,sampleInput,sampleOutput)\n",
      "        \n",
      "def new_samples() :\n",
      "    ds_input = np.zeros((nsamples,ncells),dtype = float32)\n",
      "    ds_output = np.zeros((nsamples,ncells),dtype = float32)\n",
      "    for i in range (nsamples) :\n",
      "        sample_creater(i,False)\n",
      "    train = ds_input[:cut], ds_output[:cut]\n",
      "    valid = ds_input[cut:], ds_output[cut:]"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Couldn't import dot_parser, loading of dot files will not be possible.\n"
       ]
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "def clusterer (in_array,revised_number = 0) :\n",
      "    n_lines = len(in_array)\n",
      "    n_gaps = n_lines - 1\n",
      "    gap_matrix = np.zeros((n_gaps,3),dtype=float32)\n",
      "    length_array = np.zeros(n_gaps,dtype=float32)\n",
      "    for i in range (n_gaps) :\n",
      "        gap_matrix[i][2] = i\n",
      "        gap_matrix[i][0] = in_array[i+1] - in_array[i]\n",
      "        length_array[i] = in_array[i+1] - in_array[i]\n",
      "        gap_matrix[i][1] = (in_array[i+1] + in_array[i])/2.0\n",
      "    length_array = sorted(length_array)\n",
      "    gap_matrix.view('f32,f32,f32').sort(order=['f0'], axis=0)\n",
      "    out_matrix = np.zeros((int(n_lines),2),dtype=int32)\n",
      "    if revised_number == 0 :\n",
      "        col = np.reshape(length_array,(len(length_array),1))\n",
      "        algo = KMeans(n_clusters=2)\n",
      "        gap_candidates = algo.fit_predict(col)\n",
      "        if gap_candidates[0] == 0:\n",
      "            accepted = list(gap_candidates).count(1)\n",
      "        else :\n",
      "            accepted = list(gap_candidates).count(0)\n",
      "    else :\n",
      "        accepted = revised_number-1\n",
      "    for i in range (n_gaps-accepted) :\n",
      "        gap_matrix[i][1] = in_array[n_gaps]\n",
      "    gap_matrix.view('f32,f32,f32').sort(order=['f1'], axis=0)\n",
      "    if accepted == 0 :\n",
      "        out_matrix [0] = [0,n_gaps]\n",
      "        return accepted+1,out_matrix\n",
      "    out_matrix [0] = [0,gap_matrix[0][2]]\n",
      "    out_matrix [accepted] = [gap_matrix[accepted-1][2]+1,n_gaps]\n",
      "    for bunch in range (1,accepted) :\n",
      "        out_matrix [bunch] = [int(gap_matrix[bunch-1][2]+1),int(gap_matrix[bunch][2])]\n",
      "    return accepted+1,out_matrix\n",
      "\n",
      "def clustering (in_array) :\n",
      "    central_tracks ,_ = clusterer (in_array)\n",
      "    col = np.reshape(in_array,(len(in_array),1))\n",
      "    min_n = max(1,int((0.75*central_tracks)))\n",
      "    max_n = int(1.25*central_tracks)\n",
      "    if max_n <= min_n + 1 :\n",
      "        return clusterer (in_array,max_n)\n",
      "    goc = np.zeros(1+max_n-min_n)\n",
      "    for nclu in range(min_n ,max_n +1):\n",
      "        algo = KMeans(n_clusters=nclu)\n",
      "        track_candidates = algo.fit_predict(col)\n",
      "        track_candidate_values = [np.array(in_array)[np.where(track_candidates==n)] for n in range(nclu)]\n",
      "        goc[nclu - min_n] = algo.score(col)\n",
      "    d_goc = []\n",
      "    for i in range (max_n-min_n):\n",
      "        d_goc.append(goc[i+1] - goc[i])\n",
      "    col = np.reshape(d_goc,(len(d_goc),1))\n",
      "    algo = KMeans(n_clusters=2)\n",
      "    gap_candidates = algo.fit_predict(col)\n",
      "    if gap_candidates[0] == 0:\n",
      "        accepted = list(gap_candidates).count(1)\n",
      "    else :\n",
      "        accepted = list(gap_candidates).count(0)\n",
      "    revised_number = max_n-accepted\n",
      "    return clusterer (in_array,revised_number)\n",
      "def similarity (exptrack,nntrack) :\n",
      "    sim_counter = 0\n",
      "    for exphit in range (len(exptrack)) :\n",
      "        for nnhit in range (len(nntrack)) :\n",
      "            if exptrack[exphit] == nntrack[nnhit] :\n",
      "                sim_counter += 1\n",
      "    return sim_counter\n",
      "def got (exp,nn,cut = 0.1,requested_efficiency = 1) : # Dear reader, get ready for awful variable names and no comment !\n",
      "    goodTrack = 0\n",
      "    # creating MidPut\n",
      "    for i in range (len(nn)) :\n",
      "        if nn[i] < cut :\n",
      "            nn[i] = 0\n",
      "    exp_Output = np.array([[i,exp[i]] for i in range (len(exp)) if exp[i] != 0],dtype = float32)\n",
      "    nn_Output = np.array([[i,nn[i]] for i in range (len(nn)) if nn[i] != 0],dtype = float32)\n",
      "    exp_Output.view('f32,f32').sort(order=['f1'], axis=0)\n",
      "    nn_Output.view('f32,f32').sort(order=['f1'], axis=0)\n",
      "    n_exp,exp_matrix = clustering (exp_Output[:,1])\n",
      "    n_nn,nn_matrix = clustering (nn_Output[:,1])\n",
      "    expM = np.array([[exp_Output[i][0] for i in range (exp_matrix[track][0],exp_matrix[track][1]+1)] for track in range (n_exp)])\n",
      "    nnM = np.array([[nn_Output[i][0] for i in range (nn_matrix[track][0],nn_matrix[track][1]+1)] for track in range (n_nn)])\n",
      "    similarities = np.zeros(n_exp,dtype = float32)\n",
      "    # doesn't work quite well for efficiencies less than 1,I know the reason,but don't have time to fix it !\n",
      "    remained_NNtracks = np.array([i for i in range (len(nnM))])\n",
      "    for exptrack in range (n_exp) :\n",
      "        for nntrack in remained_NNtracks :\n",
      "            if similarity(expM[exptrack],nnM[nntrack]) >= requested_efficiency * len(expM[1]) :\n",
      "                goodTrack +=1\n",
      "                remained_NNtracks = np.array([remained_NNtracks[i] for i in range (len(remained_NNtracks)) if i != nntrack]) #like !\n",
      "        #similarities[exptrack] = max([similarity(expM[exptrack],nnM[nntrack]) for nntrack in range (len(nnM))])\n",
      "    for exptrack in range (n_exp) :\n",
      "        if similarities[exptrack] >= requested_efficiency * len(expM[1]) :\n",
      "            goodTrack +=1\n",
      "    return float(goodTrack)/len(expM)\n",
      "def max_finder (listt) :\n",
      "    m = 0\n",
      "    for i in range (len(listt)) :\n",
      "        if listt[i] > listt[m] :\n",
      "            m = i\n",
      "    return m\n",
      "def compare (out,learn,n = 4,m = 4,efficiency = 1) : # n is nlayers , m is ntracks\n",
      "    counter = 0\n",
      "    output = np.zeros(len(out) , dtype=float32)\n",
      "    Otracks = np.zeros((m,n),dtype = float32)\n",
      "    Ltracks = np.zeros((m,n),dtype = float32)\n",
      "    for i in range (len(out)) :\n",
      "        output[i] = out[i]\n",
      "    learned = learn\n",
      "    for track in range (m) :\n",
      "        for layer in range (n) :\n",
      "            start = int(float(layer) * len(output) / n)\n",
      "            end = int(float(layer+1) * len(output) / n)\n",
      "            subOutput = output[start:end]\n",
      "            subLearned = learned[start:end]\n",
      "            Otracks[track][layer] = max_finder(subOutput)+start\n",
      "            Ltracks[track][layer] = max_finder(subLearned)+start\n",
      "            output[max_finder(subOutput)+start] = 0\n",
      "            learned[max_finder(subLearned)+start] = 0\n",
      "    for track in range (m) :\n",
      "        if max([similarity(Otracks[track],Ltracks[i]) for i in range(m)]) >= efficiency * n :\n",
      "            counter += 1\n",
      "    return counter"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 2
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "new_samples()"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 3
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "training_error = []\n",
      "validation_error = []\n",
      "e = []\n",
      "for i in range(10) :\n",
      "    e.append(0)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 4
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#trainer.load(\"network/500K4L-rprop\")"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 10,
       "text": [
        "<theanets.feedforward.Regressor at 0x7f3b38d17bd0>"
       ]
      }
     ],
     "prompt_number": 10
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# we try to create some training samples and train on them and then check the mean error in some testing samples\n",
      "# a function of iteration and train per iterations\n",
      "niter=1\n",
      "ntrain=2\n",
      "for i in range (0,niter+1) :\n",
      "    # training on training samples\n",
      "    for j in range (ntrain) :\n",
      "        start = time.mktime(time.gmtime())\n",
      "        print \"training started ... \"\n",
      "        e[i-1] = trainer.train(train,valid, algo='rprop',learning_rate= (1.0/(j+1)) * 10**(-1*i),batch_size=1000)\n",
      "        print e[i-1]\n",
      "        print \"\\ttraining ended\"\n",
      "        stop = time.mktime (time.gmtime())\n",
      "        spent.append (stop - start)\n",
      "        print \"training time (seconds) \\t\" , stop - start"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "training started ... \n",
        "(OrderedDict([('loss', 0.67500000000000571), ('err', 0.67500000000000571)]), OrderedDict([('loss', 0.67500000000000027), ('err', 0.67500000000000027)]))"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        "\ttraining ended\n",
        "training time (seconds) \t2020.0\n",
        "training started ... \n",
        "(OrderedDict([('loss', 0.67500000000000571), ('err', 0.67500000000000571)]), OrderedDict([('loss', 0.27440216401852935), ('err', 0.27440216401852935)]))"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        "\ttraining ended\n",
        "training time (seconds) \t360.0\n",
        "training started ... \n",
        "None"
       ]
      }
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "plt.plot(training_error,'b')\n",
      "print training_error[-1]"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "0.00637731631524\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEACAYAAABCl1qQAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XuYneO9//H3JwdxiIqkIpVEQjW6CUJKiGAkQQhCquyi\n3Q4ll9aulh9Bm5r0wO5G2erQltKdICIR6hRRatqK9BKEHESVTWRCjhqJopLM9/fHc0eWsWbWJJnk\nWTPzeV3XujzrXvf9rO/Syjf38VFEYGZmVp9WeQdgZmblz8nCzMxKcrIwM7OSnCzMzKwkJwszMyvJ\nycLMzEoqmSwkDZE0S9LLkkYW+bydpPGpzlRJPQo+uyy1myXpyILyDpImSHpJ0lxJB6by8ZJmpNcb\nkmak8p6SPiz47ObG+flmZtYQber7UFI74BZgALAImCbp8YiYUVDtfOCdiDhF0gnADcAwSX2B4cBe\nQBfgaUm9ImIVcCswKSLGSWoFtAeIiFMKvvsaYHnB97wWEftu5O81M7MNUKpn0Q+YExELImI1MB4Y\nWqvOMcDYdP0g0D8lgKHAPRGxJiIWAHOAfpI6AX0iYhxARNRExIrCG0oScDIwbiN+m5mZNZJSyaIb\nML/gfXUqK1onImqAZUBnoGuqX7vtbsASSfdKmi1pjKT2te55CLAoIl4vKOsp6UVJz0ga2IDfZmZm\njaRUsmjss0ACaA3sD1wdEb2Bd4FRtep9Hbi74P3bQNeI6AN8BxgrqUMjx2ZmZnWod86CrDfQveB9\ndz7d01hbZ2dgcRp+6gQsKdK2Wyp7C1gQEdNT+UQKkoWkNsCJwH5ryyLiY+DjdD1D0mzgy8BfCwOR\n5IOuzMw2QESovs9L9SymA70ldZXUlmweYXKtOo8Cp6frYcC0iFiTyk+R1EZSN6A38GxEVANLJfVK\nbQYDcwvuNxiYGxFvry2Q1DElIiT1TPd6rVjAEVH2ryuuuCL3GBynY3ScjnPtqyHq7VlExEeSzgOm\nkCWWsRHxgqTRwHMR8RBwI9mw0CxgJXBqavu8pPuBmUANMCKylVAAZwN3SdoamAecVvC1p/DZie3D\ngdGSagAB342IpQ36hWZmttFKDUMREZOp1ZuIiCsKrv9F1uMo1vZK4Moi5S+RzVsUa3NmkbL7gPtK\nxWpmZpuGd3DnoKKiIu8QGsRxNp6mECM4zsbWVOJsCDV0vKopkBTN6feYmW0OkoiNnOA2MzNzsjAz\ns9KcLMzMrCQnCzMzK8nJwszMSnKyMDOzkpwszMysJCcLMzMrycnCzMxKcrIwM7OSnCzMzKwkJwsz\nMyvJycLMzEpysjAzs5KcLMzMrCQnCzMzK8nJwszMSiqZLCQNkTRL0suSRhb5vJ2k8anOVEk9Cj67\nLLWbJenIgvIOkiZIeknSXEkHpvLxkmak1xuSZpS6l5mZbXpt6vtQUjvgFmAAsAiYJunxiJhRUO18\n4J2IOEXSCcANwDBJfYHhwF5AF+BpSb0iYhVwKzApIsZJagW0B4iIUwq++xpgeboudq/dI+Ljjf9X\nYGZmpdSbLIB+wJyIWADZ3/yBoUBhsjgGuCRdPwjcmhLAUOCeiFgDLJA0B+gnaS7QJyK+BhARNcCK\nwi+VJOBk4PBUVOxeBwBP1w64c2fo0AG22y57Fbuu6/PttoM2pf6NmJm1QKX+aOwGzC94Xw1U1FUn\nImokLQM6A12BP9Zq2w1YBSyRdC+wB/AC8O2IeL+g7iHAooh4Pb2v616fMXMmvPceLF+e/bP29Wuv\nFS9fvhxWrICttmpYYqkrCbVvD6r3sedmZk1PqWQRjfx9AbQG9gcuiIjpkq4HRgGF8yFfB+7ekC/o\n0iV7bVBwAe+/vy6J1JVw5s2rO+F89BFsu239iWWnnWCXXaBnT+jRI0tQZmblrFSyqAa6F7zvzqd7\nGmvr7AwsTsNPnYAlRdp2S2VvAQsiYnoqn0iWLACQ1AY4Edivnjhq93g+UVlZ+cl1RUUFFRUV9fy8\nT5OyP+i33Ra6dy9dv5jVq7MeSl0JZ/lymDED7r8f3ngD5s+H7bfPEsfaBFJ4vfPO0K7dhsViZlZM\nVVUVVVVV69VGEXV3HiRtCbwCHAwsBp4BRkTECwV1LgK6R8T3JJ0InBkRx6dJ6V8BB5EmpYEvRcQq\nSc8Bp0bEq5IqgQ4R8b10vyHAyIg4vOA76rxXrXijvt9Tjmpq4J13ssTx5pvr/rn2esEC2GGHupNJ\n9+7Qtm2OP8DMmjxJRES9A+j1Jot0k6OBq8mW2Y6NiKskjQaei4iH0oqpscC/ASvJksCbqe3lwOlA\nDXBRRExJ5fsAtwFbA/OA0yLiH+mzO4BpEfGbWnEUvVetOk0uWZSyejW8/fZnk8jafy5cmA271ZVM\nunb1pL2Z1a9RkkVT0hyTRSmrVkF1dfFk8uabsHhxljDqSiZf+AK0bp1f/GaWPycL41//yuZF6uqZ\n/OMf2VBW7SSy666w776eLzFrCZwsrKQPP4S33vpsMvn737NlxoccAkcdlb169fKyYLPmyMnCNsq7\n78KTT8KUKdmrdWs48sgscQwalC0DNrOmz8nCGk0EzJ27LnFMnQp7772u1/GVr3juw6ypcrKwTebD\nD+Evf1mXPN55BwYPzhLHkUdCt6L7682sHDlZ2GZTXQ1/+EOWOJ54IlvOuzZxHHqod6mblTMnC8vF\nmjXw/PPreh0vvQT9+68bstpjD0+Um5UTJwsrC8uXwx//uC55rFmzbqJ88GDo2DHvCM1aNicLKzsR\n8Oqr6xLHX/6S9TTW9joOOMA7zs02NycLK3v/+le2smpt8njrLRg4cF3Po0eP0vcws43jZGFNzsKF\n8PjjWeL4wx+yIaq1vY7DDoNttsk7QrPmx8nCmrSamuw497W9jhdegH791iWPvfbyRLlZY3CysGZl\n5Up46ql1yeODD+Doo+HCC2HPPfOOzqzpcrKwZu3112HCBLjuumxV1Y9+BLvvnndUZk1PQ5JFq80V\njFlj++IX4dJLswMP99gDBgyAM87IkoiZNS4nC2vytt0WfvCDLGn07JnNa5xzTvasdDNrHE4W1mxs\ntx1UVmb7OHbYAfbbD77znezRtGa2cZwsrNnp2BGuvBJeeQW23jpbNfX972fLcs1swzhZWLO1ww5w\n9dUwZ062DHePPeCSS2Dp0rwjM2t6SiYLSUMkzZL0sqSRRT5vJ2l8qjNVUo+Czy5L7WZJOrKgvIOk\nCZJekjRX0kEFn/1nKp8l6epU1lPSh5JmpNfNG//TraX4whfgf/4HZs7Mlt/uvjv88IfZI2XNrGHq\nTRaS2gG3AEOAvYGTJO1bq9r5wDsRsRdwNXBDatsXGA7sldr/WlLb1OZWYFJE7APsCcxJbYYCRwF9\n0/3+q+B7XouIfdPr2xv6g63l6tYNbrklOxF34UL40pdg9Gh47728IzMrf6V6Fv2AORGxICJWA+OB\nobXqHAOMTdcPAv0ltUr17omINRGxgCwh9JPUCegTEeMAIqImIlak9t8Cfp6+i4hYtpG/z+wzevaE\n226Dv/41W2b7pS/BVVfB++/nHZlZ+SqVLLoB8wveV6eyonUiogZYBnQGuqb6tdvuBiyRdK+k2ZLG\nSGqf6nwZOErSi5KmSepf0L5nKn9G0sD1+I1mRe22G4wZA3/+c/bMjd12g2uvzXaGm9mnlToMurG3\nQwfQGtgfuCAipku6HhgFjCRLXttGRB9J+wP3pTmQt4GuEbEiDYM9LGnPiFhe+wsqKys/ua6oqKCi\noqKRf4I1N1/+MtxzD8yalS29vfbabLPfuefCllvmHZ1Z46uqqqKqqmq92tR73IekQ4CREXFsen8x\nsEVE/KygzpOpznNp+GkR0AW4HPgwIq5J9R4GrgLmAU9HRM+C7/hhRBwl6QngJxHxp/TZa8CAiPjU\nokdJU4ArIuKvtcp93IdttBkzsqNDXnwx2+x31lmwxRZ5R2W26TTGcR/Tgd6SuqbJ6ZOBybXqPAqc\nnq6HAdMiYk0qP0VSG0ndgN7AsxFRDSyV1Cu1GQTMTdePAANT8L2ArYHFkjqmRISknuler5WI3WyD\n7LsvPPQQ3HcfPPAA9OoFv/0trFqVd2Rm+Sl5kKCko8lWObUCxkbEVZJGA89FxENpxdRY4N+AlcCp\nEfFmans5WSKpAS6KiCmpfB/gNrJkMA84LSL+kRLS7cDaFVcXRsTjkr4KjE73EVAZEfcVidU9C2t0\nU6dmPY158+CKK+DUU6F167yjMms8PnXWrBFVVcGoUdmmvspK+NrXoJW3tVoz4GRh1sgisif4jRqV\nrZoaPRpOPNEPYbKmzcnCbBOJgEceyYanAH78Yxg61EnDmiYnC7NNLCKbBP/Rj7Lng//4x3DEEU4a\n1rQ4WZhtJjU12VP7Kivh85/Pksbhh+cdlVnDOFmYbWZr1sC4cdlcRvfuWdIYMCDvqMzq58eqmm1m\nrVvD6afD3LnZP7/xDTjqKJg+Pe/IzDaOk4XZJtCmTbbz+29/g+HD4fjjs3mN1avzjsxsw3gYymwz\nWLgQvvlN+Oc/4e67oUeP0m3MNhcPQ5mViS5d4LHH4IQTYP/9YeLEvCMyWz/uWZhtZs8+C1//Ogwe\nDNddlz0n3CxP7lmYlaEDDshOtn3/ffjKV7LHvZqVOycLsxx87nNw553ZczMGDYKbbso2+JmVKw9D\nmeXs1VezYanu3bOj0Dt1yjsia2k8DGXWBPTqBc88A1/8IvTpA3/6U94RmX2WexZmZWTy5Gx/xjnn\nZPsy2pR68LFZI/BxH2ZNkPdk2ObmYSizJsh7MqwcuWdhVsa8J8M2B/cszJo478mwclEyWUgaImmW\npJcljSzyeTtJ41OdqZJ6FHx2WWo3S9KRBeUdJE2Q9JKkuZIOKvjsP1P5LElXl7qXWXPnPRlWDuod\nhpLUDngFGAAsAqYB50bEjII6FwHdI+J7kk4AzoyIYZL6Ar8CDgS6AE8DvSJilaQJwKSIGCepFdA+\nIlZIGgqcB5wQEasldYqIZXXca/eI+LhWvB6GsmbNezJsU2iMYah+wJyIWBARq4HxwNBadY4Bxqbr\nB4H+KQEMBe6JiDURsQCYA/ST1AnoExHjACKiJiJWpPbfAn6evouIWJbKi93rgBKxmzU73pNheSmV\nLLoB8wveV6eyonUiogZYBnQGuqb6tdvuBiyRdK+k2ZLGSGqf6nwZOErSi5KmSeqfyuu6l1mL064d\nXHst/OY38O//7udk2OZRastPY4/pBNAa2B+4ICKmS7oeGAWMJEte20ZEH0n7A/dJ6rk+X1BZWfnJ\ndUVFBRUVFY0SuFm5OfrobPL7m9+Eww7zngxruKqqKqqqqtarTalkUQ10L3jfnU/3NNbW2RlYnIaf\nOgFLirTtlsreAhZExNoHTd4H/DBdzwcmAaRE8jGwYx33qh0H8OlkYdbcrd2Tce212Z6Mm2+Gk07K\nOyord7X/Ij169OiSbUoNQ00HekvqKqktcDIwuVadR4HT0/UwYFpErEnlp0hqI6kb0Bt4NiKqgaWS\neqU2g4C56foRYCBA+nxrson1ovcq+evMWoBWreDii+Hhh2HkSBgxAj74IO+orLmpN1lExEdkq5Om\nAC+RrWB6QdJoScelajcCO0maBVwMfDe1fR64H5gJPAaMiIhVqc3ZwF2S1k5Ujy64166SZpP1MM5I\nE+D13cvM8J4M27S8g9usmYmAsWPhoougshK+/W1QvYsiraXzQYJmLZj3ZFhD+bgPsxbMezKsMbln\nYdYC+DkZVh8PQ5nZJxYuhG98I1sp5T0ZVsjDUGb2iS5dYMoUGDbMz8mw9eeehVkL5OdkWCH3LMys\nqLV7Mlau9J4MaxgnC7MW6nOfg7vuynZ9Dxzo52RY/TwMZWa8+mp2gu2uu8KYMR6Wamk8DGVmDdKr\nF0ybBlttBYcfDosW5R2RlRsnCzMDsudkjBkDQ4bAQQfB3Lml21jL4WEoM/uM3/0um8sYPx78SJjm\nz8NQZrZBzjgj27h38slw5515R2PlwD0LM6vT7Nlw7LFw9tnwwx/69Nrmysd9mNlGe+edLGHssw/8\n+tfQtm3eEVlj8zCUmW20L3whO7F2yRI45hh47728I7I8OFmYWUnt28MDD8Duu8PBB8O8eXlHZJub\nk4WZNUjr1vDLX2bzF/37w/PP5x2RbU6eszCz9TZpEowYAXfckc1nWNPWKHMWkoZImiXpZUkji3ze\nTtL4VGeqpB4Fn12W2s2SdGRBeQdJEyS9JGmupANTeaWkakkz0mtIKu8p6cOC8pvX51+EmTWu4cPh\noYeyhynd7P8aW4R6n5clqR1wCzAAWARMk/R4RMwoqHY+8E5EnCLpBOAGYJikvsBwYC+gC/C0pF4R\nsQq4FZgUEeMktQLap3sF8IuI+EWRcF6LiH03/KeaWWM68ECYOjWb9P6//4P//m9o5YHtZqvU/7T9\ngDkRsSAiVgPjgaG16hwDjE3XDwL9UwIYCtwTEWsiYgEwB+gnqRPQJyLGAURETUSsKLifV3KbNRG7\n7po953v6dPja17Kn8FnzVCpZdAPmF7yvTmVF60REDbAM6Ax0TfVrt90NWCLpXkmzJY2R1L6g3nfS\n0NSdkjoWlPeU9KKkZyQNbOgPNLNNq2NHePxx2HLL7KjzxYvzjsg2hVKPbW/s2eIAWgP7AxdExHRJ\n1wOjgJHATcCPU91KsiGt04G3ga4RsULSvsDDkvaMiOW1v6CysvKT64qKCip8sI3ZJteuXXYsyI9+\nlB1C+Oij2TJbK09VVVVUVVWtV5t6V0NJOgQYGRHHpvcXA1tExM8K6jyZ6jyXhp8Wkc1RXA58GBHX\npHoPA1cB84CnI6JnKh8AjIqIo2p9907AUxHxmf/LSZoCXBERf61V7tVQZjm7/Xa47DKYMAEOPTTv\naKwhGmM11HSgt6SuktoCJwOTa9V5lOxv/wDDgGkRsSaVnyKpjaRuQG/g2YioBpZK6pXaDAbmpoA7\nF9z3q2TzHEjqmBIRknqme71WInYzy8FZZ2W9jJNOyg4jtOah3mGoiPhI0nnAFLLEMjYiXpA0Gngu\nIh4CbgTGSpoFrAROTW2fl3Q/MBOoAUaklVAAZwN3SdqarKdxWiq/VtLewBap/OxUfjgwWlIN2QT4\ndyNiaSP8fjPbBI44Ap58MtuD8cYbcPnlPoSwqfOmPDPbZN5+O0sY++4Lv/qVDyEsVz5I0MxytdNO\n8Oc/w8KFMHSoDyFsypwszGyTat8efv972G03OOQQmD+/dBsrP04WZrbJtWkDN90E//Ef2dLaGTNK\nt7Hy4jkLM9usJk6E886D//3f7KgQy5/nLMys7Jx0Ejz4YHbU+S235B2NNZR7FmaWi9dfz3oWxx8P\nP/+5DyHMk5/BbWZlbdkyOOEE6NIFxoyBrbbKO6KWycNQZlbWOnWCP/whmwAfNCh7zreVJycLM8vV\nllvCXXfB4YdnK6VefTXviKyYUqfOmpltcq1awc9+Brvskh0+OHEiDBiQd1RWyD0LMysb3/pWtqT2\nxBPhnnvyjsYKeYLbzMrOzJnZmVLf/jaMHOlDCDc1r4YysyZrwYIsYey/f7b724cQbjpeDWVmTVbX\nrtkhhNXVcNxxsGJF3hG1bE4WZla2tt022+3ds2d2CGF1dd4RtVxOFmZW1tq0yY4FOe20bGntiy/m\nHVHL5DkLM2sy7r0XvvMdGDsWhgzJO5rmw3MWZtasnHwyPPBAdtT5nXfmHU3L4k15ZtakHHwwPPVU\n1rNYtAguuijviFqGkj0LSUMkzZL0sqSRRT5vJ2l8qjNVUo+Czy5L7WZJOrKgvIOkCZJekjRX0oGp\nvFJStaQZ6XV0qXuZWcuzxx4wdSr89rdw8cVQU5N3RM1fvXMWktoBrwADgEXANODciJhRUOcioHtE\nfE/SCcCZETFMUl/gV8CBQBfgaaBXRKySNAGYFBHjJLUC2kfECklXACsj4he14ih2r90j4uNa9Txn\nYdaCLFuW7cXo1Qtuu817MTZUY8xZ9APmRMSCiFgNjAeG1qpzDDA2XT8I9E8JYChwT0SsiYgFwByg\nn6ROQJ+IGAcQETURUbiCuljAxe51QInYzayZ69QJnnhi3VHn//xn3hE1X6WSRTeg8PHq1amsaJ2I\nqAGWAZ2Brql+7ba7AUsk3StptqQxktoX1PtOGpq6U1LHVFbXvcyshdtmG7j/fthhh+yY82XL8o6o\neSo1wd3YYzoBtAb2By6IiOmSrgdGASOBm4Afp7qVwA3A6evzBZWVlZ9cV1RUUFFRsbExm1mZa9sW\n7rgDLr00O612yhTYeee8oypfVVVVVFVVrVebUsmiGuhe8L47n+5prK2zM7A4DT91ApYUadstlb0F\nLIiI6al8IlmyICKWrq0s6dfAU3XEUbvH84nCZGFmLYeUPZ51xx2zhDF5Muy5Z95Rlafaf5EePXp0\nyTalhqGmA70ldZXUFjgZmFyrzqOs+9v/MGBaRKxJ5adIaiOpG9AbeDYiqoGlknqlNoOBuQCSOhfc\n96tkcxPUda+Sv87MWpwLL4SrroKBA7MVU9Y46u1ZRMRHks4DppAllrER8YKk0cBzEfEQcCMwVtIs\nYCVwamr7vKT7gZlADTAiIlalW58N3CVpa2AecFoqv1bS3sAWqfzsBtzLzOxTTjsNPv/5bNL79tuz\ngwht4/i4DzNrtp59Fo4/Hq68Es46K+9oyldDls56B7eZNVsHHAB/+hMcdVS22/vSS/0gpQ3lnoWZ\nNXsLFmTHgwwaBL/4RfbMb1vHT8ozM0uWL8+GpLp1g9/9DrbYIu+IyodPnTUzSzp0yPZffPBBdkTI\nypV5R9S0OFmYWYux1VYwcSL06JEtrV28OO+Img4nCzNrUdq0gd/8JpvDGDAA3ngj74iaBq+GMrMW\nR4Kf/CTb7X3IIfDII7DPPnlHVd6cLMysxTr/fOjcGY44AiZMgMMOyzui8uVhKDNr0U4+GcaNg5NO\ngkmT8o6mfLlnYWYt3qBB2UqpY4+FJUtgxIi8Iyo/ThZmZsB++8Gf/7xut/eoUd7tXcib8szMCixc\nCEcfDQcdBL/8JbRunXdEm553cJuZbYAVK7ITazt2hDvvhC23zDuiTcs7uM3MNsDnPgePPpoNQx19\nNLz3Xt4R5c/JwsysiC23hHvugT32gIqKbHiqJXOyMDOrQ+vWcOONMHw4HHwwvPZa3hHlx6uhzMzq\nIWUro3bcEQ49FB56CPr2zTuqzc89CzOzBjj33KyXMWQIPPFE3tFsfk4WZmYNNHx4dmrtqafCvffm\nHc3mVTJZSBoiaZaklyWNLPJ5O0njU52pknoUfHZZajdL0pEF5R0kTZD0kqS5kg6qdc+LJNVI6pje\nV0h6T9KM9Prhxv1sM7MNc9hhWc/iwguznkZLUe+chaR2wC3AAGARME3S4xExo6Da+cA7EXGKpBOA\nG4BhkvoCw4G9gC7A05J6RcQq4FZgUkSMk9QKaF/wnd2BI4B5tcL5U0QcvzE/1sysMey9N/zlL9lu\n74ULsxNsm/tu71I9i37AnIhYEBGrgfHA0Fp1jgHGpusHgf4pAQwF7omINRGxAJgD9JPUCegTEeMA\nIqImIlYU3O8XwCVFYmnm/1OYWVOyyy4wdWp2ptQ558Dq1XlHtGmVShbdgPkF76tTWdE6EVEDLAM6\nA11T/dptdwOWSLpX0mxJYyS1B5A0DKiOiJm1viOAg9Jw1pOSfPK8meVuhx3gj3+Et97KTq398MO8\nI9p0Si2dbeyzMwJoDewPXBAR0yVdD4ySVAlcTjYEtdba3sTzQLeI+CjNfTwA7FLsCyorKz+5rqio\noKKiopF/gpnZOttuCw8/DGecAUceCQ8+CNtvn3dU9auqqqKqqmq92tR7NpSkQ4CREXFsen8xsEVE\n/KygzpOpznNp+GkR2RzF5cCHEXFNqvcwcBXZXMTTEdEzlQ8ARgEXAU8CH6RbdwMWAAdExKeelCvp\nb8BhEbGwVrnPhjKzXNTUZJPeTz4Jjz0GXbvmHVHDNcbZUNOB3pK6SmoLnAxMrlXnUeD0dD0MmBYR\na1L5KZLaSOoG9AaejYhqYKmkXqnNYGBuRMyOiB0jYpeI2IVs2Gq/iFgs6fMFP6ovsA3gR62bWdlo\n1Qquuw5OOy3b7f3KK3lH1LjqHYZKwz7nAVPIEsvYiHhB0mjguYh4CLgRGCtpFrASODW1fV7S/cBM\noAYYkVZCAZwN3CVpa7Kexmkl4vy6pHPT9cfAqWl+xMysbEhw6aXZbu+KCvj976Ffv7yjahw+otzM\nbBN4+GE480y47TYYNizvaOrnI8rNzHJy7LEweTKcd172EKWmzj0LM7NN6I034JhjstfVV2dzG+XG\nT8ozMysD774LJ56Y7csYOxa22irviD7Nw1BmZmWgY0d4/HFo2xYGD4alS/OOaP05WZiZbQbt2sFd\nd2XPxOjfH15/Pe+I1o8ffmRmtpm0agVXXQU9e8KAAXD//XDggXlH1TDuWZiZbWYjRmRLao87LksY\nTYGThZlZDoYOzZbWnn8+3HBD3tGU5tVQZmY5evPNbFntkCFwzTX5LK310lkzsybgH//IltZ26gR3\n3rn5l9Z66ayZWROw/fbZQ5S23BIGDYIlS/KO6LOcLMzMykC7dtmGvYqKbGnta6/lHdGneemsmVmZ\naNUKrrwSevRYt7T2oIPyjirjnoWZWZkZMQJuvx2OPx4mTco7mox7FmZmZeiYY7J5jOOOy57x/b3v\n5RuPV0OZmZWxefOyxHHEEXDttdC6deN/h5fOmpk1A8uXZ0trt98+O1+qsZfWeumsmVkz0KEDPPYY\nbL01DByYz9JaJwszsyZg7dLaQYOypbV///vm/f6SyULSEEmzJL0saWSRz9tJGp/qTJXUo+Czy1K7\nWZKOLCjvIGmCpJckzZV0UK17XiSpRlLHgrIbJM2R9IKkfTf8J5uZNU0S/PSncMklcMgh8Mwzm++7\n600WktoBtwBDgL2Bk4r8QX0+8E5E7AVcDdyQ2vYFhgN7pfa/ltQ2tbkVmBQR+wB7AnMKvrM7cAQw\nr6Dsq8DOEbEncDZwxwb9WjOzZuCcc+COO+CEE+C++zbPd5bqWfQD5kTEgohYDYwHhtaqcwwwNl0/\nCPSX1CpCA2OzAAAI10lEQVTVuyci1kTEArKE0E9SJ6BPRIwDiIiaiFhRcL9fAJfU9R0RMQNoI6nb\n+vxQM7Pm5Oijs6W1F1wA11+/6b+vVLLoBswveF+dyorWiYgaYBnQGeia6tduuxuwRNK9kmZLGiOp\nPYCkYUB1RMzcgDjMzFqUfffNhqJuuy1LGmvWbLrvKrUpr7HXoQbQGtgfuCAipku6HhglqRK4nGwI\nai3VcV1nbJWVlZ9cV1RUUFFRsdFBm5mVq513hqefhuHD4aSTsqW1W29df5uqqiqqqqrW63vq3Wch\n6RBgZEQcm95fDGwRET8rqPNkqvNcGn5aBHQh+4P/w4i4JtV7GLiKbC7i6YjomcoHAKOAi4AngQ/S\nrbsBC8iGwq4EJkfExNRmNnBUGt4qjNf7LMysRfr4Yzj77GyV1IMPQufODW/bGPsspgO9JXVNk9Mn\nA5Nr1XkUOD1dDwOmRcSaVH6KpLXzC72BZyOiGlgqqVdqMxiYGxGzI2LHiNglInYhG2raLyIWpXud\nln7UfsCa2onCzKwl22ILGDMm2+ndvz+8+mrj3r/eYaiI+EjSecAUssQyNiJekDQaeC4iHgJuBMZK\nmgWsBE5NbZ+XdD8wE6gBRkTEqnTrs4G7JG1N1tM4rdjXF8Rxn6TDJc0B/gWcueE/2cyseZLgJz+B\nnj3h0EOzQwj792+kezenYRsPQ5mZZaZMgW98A26+OZvLqI/PhjIza8FefBGOPRYuvBC+//2s51GM\nk4WZWQv31lvZqbUDB8J11xU/tdbJwszMWL4cvvpV2HZbuPvuzy6t9amzZmZGhw4weTJstx0cfjgs\nXrz+93CyMDNrAbbYAn73OxgyJHuu99/+tn7t/VhVM7MWQoLRo6FHDzjsMJg4EQYMaFhbJwszsxbm\nrLOga9fsiJCbbmpYG09wm5m1UC++CMcdB9XVXg1lZmb1ePdd6NTJycLMzErw0lkzM2sUThZmZlaS\nk4WZmZXkZGFmZiU5WZiZWUlOFmZmVpKThZmZleRkYWZmJZVMFpKGSJol6WVJI4t83k7S+FRnqqQe\nBZ9dltrNknRkQXkHSRMkvSRprqQDU/lPU9lsSX+WtGsqr5D0nqQZ6fXDxvn5ZmbWEPUmC0ntgFuA\nIcDewEmS9q1V7XzgnYjYC7gauCG17QsMB/ZK7X8tqW1qcyswKSL2AfYEXk7l/xUR+0REb2ACcEXB\n9/wpIvZNr59u2M8tD1VVVXmH0CCOs/E0hRjBcTa2phJnQ5TqWfQD5kTEgohYDYwHhtaqcwwwNl0/\nCPSX1CrVuyci1kTEAmAO0E9SJ6BPRIwDiIiaiFiRrt8vuG974J2C9/VuRW9Kmsr/gRxn42kKMYLj\nbGxNJc6GKJUsugHzC95Xp7KidSKiBlgGdAa6pvq12+4GLJF0bxpuGiOp/dpKkn4m6S3gDOCqVBzA\nQWk460lJ+6zHbzQzs41UKlk09ql8AbQG9geuTsNN7wKjPqkQ8YOI2Bm4A7g+FT8PdEtDXT8HHmjk\nuMzMrD4RUecLOAR4uOD9xcAPatV5EvhKum4FLCFLCKOA/1dQ72HgYLLexZsF5QOAKUW+e2fglTri\n+hvQpUh5+OWXX375tf6v+nJBRJR8Ut50oLekrsBi4GRgRK06jwKnA88Bw4BpEbFG0qPAryRdD3QB\negPPRsQqSUsl9YqIV4HBwFwASbtExBvpvsOAWan88xGxNF33BbZJ8XxKqSN2zcxsw9SbLCLiI0nn\nAVPIeg1jI+IFSaOB5yLiIeBGYKykWcBK4NTU9nlJ9wMzgRpgRESsSrc+G7hL0tbAPOC0VH6tpC8C\nbYE3gG+l8q9LOjddfwycmuZHzMxsM2hWDz8yM7NNo1ns4C61cbBcSLpd0qLUCytLkrqnDZGzJP1N\n0iV5x1SMpC0lTU+bNF+VdF3eMdVHUusU60N5x1IXSW9KmpnifDbveOpSZFPvQXnHVJuk3Qs2Ec9I\nm4q/m3dctUkanf77eUXSxDTaU7xuU+9ZpI2Dr5BNlC8CpgHnRsSMXAMrQtIhwPvAmLSyq+xI2hHY\nISJmpyXNLwBfi4iXcg7tMyRtFREfSmoDPA1cFhFP5R1XMZIuBPoC20bE8XnHU4ykN4C+EfFu3rHU\nR9IEsk2949KervZr92qVoxTjAuCAiJhfqv7mImk34HHgyxHxsaTxwOMR8dti9ZtDz6IhGwfLQkT8\nBfhH3nHUJyIWRcTsdP0+2ZzTTvlGVVxEfJgutyBbgbcox3DqJKkb2ebV2yj/zaVlHV99m3rL2GDg\n9XJKFMm7wCpgm/QXrrVzyEU1h2TRkI2DtgEk9STbE/N0vpEUJ6mVpBfJksRTEfFyqTY5uY5s2Xm5\nL8oI4A9pKOr8vIOpw5eoZ1Nvmfp34O68g6gt9SCvBd4C3gaWR8QTddVvDsmiaY+jlan0H+AE4IKI\nWJl3PMWkv1X2IfvLwaGSKnIO6TMkHQssTsOiZf23duDAiNgPGAScKWlw3gEV0Yp6NvWWG0lbAMeR\n/bdUVtLK0+8BPclGD9pLOq2u+s0hWVQD3Qved+fTPQ1bT+nAx/uAuyOi7HfLR8R7wCPAgXnHUkR/\n4Pg0HzAOGChpTM4xFRURi9M/lwATyf5QLjfzgQURMT29nwj0yTGeUo4Gnk//TsvNAcAzEbEsDeFP\nIpv7Lao5JItPNg6mP+ROBibnHFOTJUnAb4GXI6JsVxhJ6iRp23S9FXAEaRNnOYmIyyOie0TsQjYc\n8ceI+GbecdUmaeu1K2EkbUN2UvScfKP6rDTuv1RSr1T0yabeMvV1sr8klKPXgAMlbZX+ux+cyooq\ntYO77NW1cTDnsIqSNA44DOgkaT7wo4i4I+ewajuYbEf+TElrV5RdFhGP5RhTMTsBY9L/ybck6wU9\nknNMDVGuw6Y7Ag9ICrKJznsi4sGcY6pLXZt6y0pKuoOBc/KOpZiImC5pIus2Ts8AbqqrfpNfOmtm\nZptecxiGMjOzTczJwszMSnKyMDOzkpwszMysJCcLMzMrycnCzMxKcrIwM7OSnCzMzKyk/w9qRxeN\nOw+4qwAAAABJRU5ErkJggg==\n",
       "text": [
        "<matplotlib.figure.Figure at 0x7f3b04fc0210>"
       ]
      }
     ],
     "prompt_number": 12
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "plt.plot(validation_error,'b')\n",
      "print validation_error[-1]"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "0.0064090880801\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEACAYAAABCl1qQAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xu41GW99/H3h7OKqVFoAoqlWAoqogKitlRUBAUFFQXT\nyp7c7szqseLCHYm70y46mFlWmhZknBTPIuZhVaAmCMrR7aO5VdAANQ+VboT1ff64f8SwWKxZwKz1\nm5n1eV3XXMzcc9+/+Q6X+J37+FNEYGZm1pg2eQdgZmblz8nCzMyKcrIwM7OinCzMzKwoJwszMyvK\nycLMzIoqmiwkDZG0RNJySeMaeL+jpOlZnXmS9i14b3zWbomkkwvKd5c0U9JTklZIGpCVT5e0KHs8\nL2lRsWuZmVnza9fYm5I6AtcBxwCrgUcl3R8RiwqqXQq8EhGjJZ0BXAOMkNQPGAn0AfYC5krqFRHv\nAdcDsyJiqqQ2QGeAiBhd8NnfB97Injd0rQMjYt2O/xWYmVkxxXoW/YFlEbEqItYD04Fh9eoMBaZk\nz+8Ejs4SwDBgWkRsiIhVwDKgv6QuwGERMRUgIuoi4q3CC0oScA4wNStq6FpHbcf3NTOz7VAsWXQH\nXip4vTIra7BORNQBrwFdgW5Z/fpt9wfWSpohaamkyZI617vmscDqiHgue721a5mZWQsolixKfRZI\nAG2BI4FJEdEbeB2YUK/eecDvSvzZZma2nRqdsyD9gu9R8LoHm/c0NtbZB1iTDT91AdY20LZ7VvYi\nsCoi5mflt1CQLCS1A84EDm8kjvo9no1tfdCVmdl2iAg19n6xnsV8oLekbpLak+YRZtercy9wfvZ8\nBPBoRGzIykdLaiepO9AbeDwiVgKvSuqVtRkMrCi43mBgRUS8XO8ztrhWQwFHRNk/rrzyytxjcJyO\n0XE6zo2Ppmi0ZxER70q6BJhDSixTImKhpKuABRFxF3AtMEXSEuBtYEzW9glJtwGLgTrg4kgroQAu\nAm6WtDPwAjC24GNHs2limyZcy8zMmlmxYSgiYjb1ehMRcWXB8/8l9Tgaavtt4NsNlD9FmrdoqM2n\ntuVaZmbW/LyDOwc1NTV5h9AkjrN0KiFGcJylVilxNoWaOl5VCSRFNX0fM7OWIInYwQluMzOz4nMW\nlWanneB970uP3Xbb9HxrZVt73b593t/EzKx8VN0w1D//Gbz1Frz5Jrz11qZH4evG3tv4un37bUsu\nW3vtpGNm5a4pw1BVlyxK8X0i4J13tj/RFL7u0GHLZNK1K+y996ZHt26bnu+8cwn+IszMtoGTRc4i\n4J//3DyRvPkmrFkDL7+cHqtWbXr+8svQqdOWCaR+Utlrr5SEzMxKwcmiwkTA3/62ZQKpn1hWr4Y9\n9mg4qRS+/uAHoW3bvL+VmZU7J4sqtWEDrF3bcM+ksOyNN9KQV2O9lL33TolHjf5nYmbVzMmilVu3\nDv761+JJ5d13t0wg3brB4YdD//6wyy55fxMza05OFtYk//gHvPLK5knkxRdh/nx46in42Mdg0KBN\nj27d8o7YzErJycJ22LvvwoIFMG/epseuu25KHMccAwcf7LkRs0rmZGElV1cH//3fmyePNWtgwIBN\nCcRDV2aVxcnCWsSaNfDII5uSh4euzCqLk4Xlov7Q1SOPbD50NWiQh67MyomThZWFiE1DV3PneujK\nrNw4WVjZamzo6phj0p977513lGatg5OFVQwPXZnlx8nCKlbh0NXGx+rVHroyaw5OFlZVGhq66tMH\nLr0Uzj0X2lXd3VnMWoaThVW1d9+Fhx6C730v7TgfNw4uvDCd3GtmTedkYa3GvHnw7W/Dk0/C5ZfD\nZz8LnTvnHZVZZfA9uK3VGDQI7rkH7r4bHnsMPvxh+MY30pHvZrbjnCysqvTtCzNmwJ/+BM8/D/vv\nn4anVq/OOzKzyuZkYVXpwAPhxhth0aJ0t8KPfSxNhL/wQt6RmVUmJwuravvsAz/5CaxYkeYwDj8c\nPvlJePrpvCMzqyxOFtYq7Lkn/Nd/wbPPpqGp446Ds8+GhQvzjsysMjhZWKuyxx7wta/BX/4CRx8N\nw4fDqaemM6vMbOuKJgtJQyQtkbRc0rgG3u8oaXpWZ56kfQveG5+1WyLp5ILy3SXNlPSUpBWSBha8\n9/msfImkSVlZT0nvSFqUPX6241/dWrPOneFLX4LnnoORI9P+jOOOgzlz0u5xM9tco/ssJHUEngaO\nAVYDjwKfjYhFBXUuB3pExBclnQF8KiJGSOoH/BwYAOwFzAV6RcR7kmYCsyJiqqQ2QOeIeEvSMOAS\n4IyIWC+pS0S8JqkncFdE9Gn0y3ifhW2n9evTKqpvfxs6doQrroAzz4Q27ntbK1CKfRb9gWURsSoi\n1gPTgWH16gwFpmTP7wSOzhLAMGBaRGyIiFXAMqC/pC7AYRExFSAi6iLiraz9Z4DvZp9FRLzWpG9q\ntoPatYMxY2DxYvj61+G7300HF06eDO+9l3d0Zvkrliy6Ay8VvF6ZlTVYJyLqgNeArkC3rH79tvsD\nayXNkLRU0mRJG/fafhQ4RdKTkh6VdHRB+55Z+SOSTtiG72jWZG3awIgR8Oc/w7XXwq9/DQccAD/7\nWTpexKy1KpYsSj2mE0Bb4EhgUkT0Bl4HJhTEs2tEHAZcBkyTJOBloFtW/jlgiqTdSxyb2b9IcOKJ\n6eypqVPhvvvSrvBJk+Dtt/OOzqzlFTuncyXQo+B1DzbvaWyssw+wJht+6gKsbaBt96zsRWBVRMzP\nym8FvpY9fwmYBRAR8yWtA/aMiL8C67LyRZKWknohj9UPeOLEif96XlNTQ01NTZGvaNa4gQPhzjvT\nENV3vpOSxuc+B5//PHTpknd0ZtuutraW2trabWpTbIK7E2mCexCwBngEuDgiFhbUKZzgPpM0wT28\nYIJ7IJsmuA/IJrgXAGMi4hlJE4Hds/Zfyp5fKakXUEtKMrsDb0REXTbZPQ84NCJerRevJ7it2f2/\n/5fmNGbNgk9/Oh1c+KEP5R2V2fbb4QnuiHiXtDppDvAUaQXTQklXSTo9q3YtsLekJcBXSMNHRMQT\nwG3AYuA+UpLZOFV4EXCzpGXAUcBVBdf6cNZzmAV8MpsHOR5YLGkxcBdwWf1EYdZSDjgAbrgh3U9j\n/fo0EX7JJeksKrNq5SPKzXbQmjXw4x/Dz38OQ4fC+PFw0EF5R2XWdD6i3KwFdO0K3/pW2hV+0EFw\n/PFpo9+CBXlHZlY6ThZmJbLbbqlX8fzzUFOTNvWdcgr84Q/eFW6Vz8NQZs1k3TqYMiUdYNi9e9rg\n16NH8XZmLc3DUGY56tABLrooHYd+yilw1FHwwAN5R2W2fdyzMGshDz8MY8emPRrjx/vcKSsfTelZ\nOFmYtaBVq9J9ND7wAfjNb9KR6WZ58zCUWZnp1g1qa9Mu8COOgCefzDsis6ZxsjBrYR06wNVXp+W2\nJ50EN92Ud0RmxXkYyixHy5fDqFFw7LFwzTXQqVPeEVlr5GEoszJ30EHw+OPwxhswaJCPDLHy5WRh\nlrNdd4Xp0+ETn4ABA+Dee/OOyGxLHoYyKyNz58K556bTbK+8Etq2zTsiaw28dNasAv31rylhdOwI\nN9+cltmaNSfPWZhVoL32Sju9DzsM+vVLcxpmeXOyMCtD7dqlGyxdfTWcdlo6/tydZsuTh6HMytwz\nz6TltX37pqSx8855R2TVxsNQZlWgVy947LHUsxgwIN3W1aylOVmYVYBddklHnF9ySdqPcfvteUdk\nrY2HocwqzOOPp8MIzz03HRnSrl3eEVml89JZsyr16qswZgy89x5MmwZ77pl3RFbJPGdhVqU+8AGY\nPTudKdWvH8ybl3dEVu3cszCrcPfck3Z8X3EFXHYZqNHfh2Zb8jCUWSvx/PNpee0BB8ANN6Tzpsya\nysNQZq3EfvvBI4+kJNG/P6xYkXdEVm2cLMyqRKdOqVdx+eVw3HEwY0beEVk18TCUWRVatCgNSw0f\nDt/7Xro7n9nWeBjKrJXq2xeeeAKefRaOPx5Wrco7Iqt0ThZmVWqPPeDOO2HoUDjySHj44bwjskrm\nYSizVuD3v0934vvSl+CrX/XyWttcSYahJA2RtETScknjGni/o6TpWZ15kvYteG981m6JpJMLyneX\nNFPSU5JWSBpY8N7ns/IlkiYVu5aZFXfSSTB/PsyaBSNHwptv5h2RVZpGk4WkjsB1wBDgEOAsSX3r\nVbsUeCUi+gCTgGuytv2AkUCfrP0vJLXP2lwPzIqIQ4GDgWVZm2HAKUC/7Hr/1ci1PGVntg169IA/\n/hG6dYMjjoDFi/OOyCpJsZ5Ff2BZRKyKiPXAdGBYvTpDgSnZ8zuBoyW1yepNi4gNEbGKlBD6S+oC\nHBYRUwEioi4i3srafwb4bvZZRMRrWXlD1zpqO7+zWavVsSNcey1MnAgnnghTphRtYgYUTxbdgZcK\nXq/MyhqsExF1wGtAV6BbVr9+2/2BtZJmSFoqabKkzlmdjwKnSHpS0qOSjs7Kt3YtM9sOY8fCQw/B\nN76Rjj3/3//NOyIrd8UONy71bHEAbYEjgS9ExHxJVwMTgHGk5LVrRBwm6UjgVkk9t+UDJk6c+K/n\nNTU11NTUlCRws2rTp0+ax/jUp9KBhLfcAvvsk3dU1hJqa2upra3dpjbFksVKoEfB6x5s3tPYWGcf\nYE02/NQFWNtA2+5Z2YvAqoiYn5XfCnwte/4SMAsgSyTrgD23cq36cQCbJwsza9xuu8Gtt8IPfgBH\nHQW//S0MHpx3VNbc6v+Qvuqqq4q2KTYMNR/oLalbNjl9DjC7Xp17gfOz5yOARyNiQ1Y+WlI7Sd2B\n3sDjEbESeFVSr6zNicDGk2zuAU4AyN7fGVi9tWsV/XZmVpQEX/5yui/G2LHpFFuz+hrtWUTEu5Iu\nAeaQEsuUiFgo6SpgQUTcBVwLTJG0BHgbGJO1fULSbcBioA64OCLeyy59EXCzpJ2BF4CxWfm1wI2S\nlmavP5nNgzR2LTMrgZqatInv9NPh179Om/nMNvKmPDPbzGOPpTOlfvMbOPXUvKOxluCzocxsmw0Y\nkHoYF14I992XdzRWLpwszGwLAwbAHXfABRfAnDl5R2PlwMnCzBo0cCDcfns6U+r++/OOxvLmZGFm\nW3X00XDbbXD++ekwQmu9nCzMrFGDBqUDCMeMgQceyDsay4uThZkVdcwxmxLGgw/mHY3lwcnCzJpk\n45Eg552XzpWy1sXJwsya7LjjYOZMOPdc33mvtXGyMLNt8vGPw4wZMHo0bONZdFbBnCzMbJvV1KSE\ncc458Ic/5B2NtQQnCzPbLjU16fDBs89Od+Cz6uZkYWbb7YQTYOpUOOss+NOf8o7GmpOThZntkBNP\nhN/9DkaOhLlz847GmouThZntsMGDNyWMefPyjsaag5OFmZXESSelO+2deSY88kje0VipOVmYWcmc\nfDJMmQJnnAGPPpp3NFZKThZmVlKnnAKTJ8OIEU4Y1cTJwsxKbsiQTQnjscfyjsZKwcnCzJrFkCHp\nXt7Dh8Of/5x3NLajnCzMrNkMHQo33QSnnw6PP553NLYjnCzMrFkNGwY33pgSxvz5eUdj28vJwsya\n3WmnwQ03pD8XLMg7GtseThZm1iJOPx2uvz71NJ54Iu9obFu1yzsAM2s9hg+HiJQw7r0XDj8874is\nqZwszKxFjRiREsapp8Ls2U4YlcLJwsxa3BlnpIQxdGhKGH375h2RFeNkYWa5OPPMTT2M++6Dww7L\nOyJrjJOFmeVm5MiUMIYMgTlz4NBD847ItqboaihJQyQtkbRc0rgG3u8oaXpWZ56kfQveG5+1WyLp\n5ILy3SXNlPSUpBWSBmTlEyWtlLQoewzJyntKeqeg/Gel+fpmlrdRo+Daa1PCWLw472hsaxrtWUjq\nCFwHHAOsBh6VdH9ELCqodinwSkSMlnQGcA0wQlI/YCTQB9gLmCupV0S8B1wPzIqIqZLaAJ2zawXw\nw4j4YQPhPBsRHtk0q0JnnZV6GKecAvffD3365B2R1VesZ9EfWBYRqyJiPTAdGFavzlBgSvb8TuDo\nLAEMA6ZFxIaIWAUsA/pL6gIcFhFTASKiLiLeKrieduwrmVklOvtsuPrqdMz50qV5R2P1FUsW3YGX\nCl6vzMoarBMRdcBrQFegW1a/ftv9gbWSZkhaKmmypM4F9T6XDU39VtL7C8p7SnpS0iOSTmjqFzSz\nyjF6tBNGuSo2wR0l/rwA2gJHAl+IiPmSrgYmAOOAnwL/mdWdSBrSOh94GegWEW9J6gvcLengiHij\n/gdMnDjxX89ramqoqakp8Vcws+Y0ejTU1aWE8fvfw8EH5x1R9amtraW2tnab2ihi6/lA0rHAuIg4\nLXv9FaBDRHyroM6DWZ0F2fDTatIcxRXAOxHx/aze3cB3gBeAuRHRMys/BpgQEafU++y9gYcj4sAG\n4poDXBkRj9Urj8a+j5lVjt/9Dr78ZXjgATjooLyjqW6SiIhGpwCKDUPNB3pL6iapPXAOMLtenXtJ\nv/4BRgCPRsSGrHy0pHaSugO9gccjYiXwqqReWZvBwIos4K4F1x1FmudA0vuzRISkntm1ni0Su5lV\nsDFjYNKkdG/vFSvyjsYaHYaKiHclXQLMISWWKRGxUNJVwIKIuAu4FpgiaQnwNjAma/uEpNuAxUAd\ncHG2EgrgIuBmSTuTehpjs/IfSDoE6JCVX5SVHw9cJamONAF+WUS8WoLvb2ZlbOzYtEpq8GB48EH4\n6Efzjqj1anQYqtJ4GMqsOk2ZAuPHpyEpJ4zSa8owlHdwm1nZ+8QnNu9hHLjFTKY1NycLM6sIF1yQ\nEsaJJ8JDD0GvXsXbWOk4WZhZxbjwQli/Ph0N8thj0LVr8TZWGp6zMLOK8/Wvpz0YDz0EO+2UdzSV\nrylzFk4WZlZxItJKqQ0bYOpUaOMbRO+QUuyzMDMrOxLceCOsWgUTJuQdTevgZGFmFalTJ7jtNpg+\nHW66Ke9oqp8nuM2sYn3wg3D33fDxj0PPnnD88XlHVL3cszCzivbRj8K0aXDuufD003lHU72cLMys\n4h1/PHz3u3DaabB2bd7RVCevhjKzqvEf/wG1tWmXd6dOeUdTObx01sxalbq6dFptmzZw881p1ZQV\n56WzZtaqtGmTVkY9/zxceWXe0VQXJwszqyo77QR33AG//S1Mnpx3NNXDS2fNrOp07Qr33AM1NbDv\nvmlpre0Y9yzMrCp97GPp1qyjR8Mzz+QdTeVzsjCzqnXiifDNb8KwYfCq7625Q7waysyq3rhx8Mgj\n6U57HTvmHU358dJZMzPSktpzzkmT35Mne0ltfV46a2ZGWlI7eXKau/jP/8w7msrk1VBm1irsvHNa\nUjtgAOy/f7ofhjWdk4WZtRp77ZVOqT3hBNhnHzj22LwjqhwehjKzVqV377Rh7+yz4dln846mcjhZ\nmFmrc/LJcNVVaUnt66/nHU1l8GooM2u1vvxlWLAA7r8fOnTIO5r8eOmsmVkjNmyAs86C970Pfv3r\n1ruk1ktnzcwa0bZtmr9Ytgy+9a28oylvXg1lZq3aLrvAXXelJbUf+Qicd17eEZWnoj0LSUMkLZG0\nXNK4Bt7vKGl6VmeepH0L3huftVsi6eSC8t0lzZT0lKQVkgZk5RMlrZS0KHucWuxaZmY76kMfSgnj\nC19Ix4LYlhpNFpI6AtcBQ4BDgLMk9a1X7VLglYjoA0wCrsna9gNGAn2y9r+Q1D5rcz0wKyIOBQ4G\nlmflAfwwIvpmj9mNXKsVT0eZWakdcgj85jcwahQ891ze0ZSfYj2L/sCyiFgVEeuB6cCwenWGAlOy\n53cCR0tqk9WbFhEbImIVsAzoL6kLcFhETAWIiLqIeKvgeg1NsjR0raOa/jXNzIo79VSYMCEtqf3b\n3/KOprwUSxbdgZcKXq/MyhqsExF1wGtAV6BbVr9+2/2BtZJmSFoqabKkzgX1PpcNTf1W0vuzsq1d\ny8yspP7931PSGDUK1q3LO5ryUWyCu9TrUANoCxwJfCEi5ku6GpgAjAN+Cmw85msiaUjr/G35gIkT\nJ/7reU1NDTU1NTsas5m1Mt//PowcCf/2b/CrX1Xfktra2lpqa2u3qU2j+ywkHQuMi4jTstdfATpE\nxLcK6jyY1VmQDT+tBvYCrgDeiYjvZ/XuBr4DvADMjYieWfkxwISIOKXeZ+8NPBwRB0qa0NC1ImJe\nvTbeZ2FmJfH3v8Nxx6VjQcaPzzua5lWKfRbzgd6SumWT0+cAs+vVuZdNv/5HAI9GxIasfLSkdpK6\nA72BxyNiJfCqpF5Zm8HAiizgrgXXHUWam2Br1yoSu5nZduvcOa2Quu46mDEj72jy1+gwVES8K+kS\nYA4psUyJiIWSrgIWRMRdwLXAFElLgLeBMVnbJyTdBiwG6oCLI+K97NIXATdL2pnU09h4WPAPJB0C\ndMjKL2rCtczMmkW3bnDnnXDSSemU2gED8o4oPz7uw8ysiHvugc98Ju3B2G+/vKMpPR/3YWZWAsOG\nwRVXpD/feCPvaPLhnoWZWRNddhksXw6zZ0P79sXrVwqfOmtmVkIbNsCIEel4kF/+snqW1HoYysys\nhNq2halTYf58mDQp72halk+dNTPbBrvumu7jPXBgOqV21Ki8I2oZHoYyM9sOCxfCKaeklVJHVfhJ\ndR6GMjNrJocfno4COeMMeOGFvKNpfh6GMjPbTsOHw1/+kpbUzpsHu+2Wd0TNx8NQZmY7IAIuvRSe\nfTbNZVTikloPQ5mZNTMJfvzjtFLq859PyaMaOVmYme2gdu1g2rR0HMgPf5h3NM3DcxZmZiXwvvdt\nvqT2jDPyjqi0PGdhZlZCCxakO+3Nng1HHJF3NE3jOQszsxZ2xBFw/fXpWJAXX8w7mtLxMJSZWYmd\ncQY89xycdhrMnZuGqCqdh6HMzJpBBFxySdqwd9ddaRK8XHkYyswsJxL85CdQV5eONq/037FOFmZm\nzaR9+3T/7j/+Ea6+Ou9odkwZd4zMzCrfbrulwwYHDoQPfzhNfFciz1mYmbWAxx9PZ0jddx/065d3\nNJvznIWZWZk46ij4+c9Tz+Kll/KOZtt5GMrMrIWMGrX5ktpdd807oqbzMJSZWQuKgM9+Fl5+Ge64\nozyW1HoYysyszEjws5/BunXwpS/lHU3TOVmYmbWw9u1h5kx46CG45pq8o2maMugAmZm1Prvvnk6p\nHTQoLak97bS8I2qc5yzMzHL05z+nRHH//dC3bz4xeM7CzKzM9e8P112X7ue9alXe0Wxd0WQhaYik\nJZKWSxrXwPsdJU3P6syTtG/Be+OzdksknVxQvrukmZKekrRC0sB617xcUp2k92evayS9KWlR9vja\njn1tM7PycdZZ8LnPwemnw9//nnc0DWs0WUjqCFwHDAEOAc6SVL+jdCnwSkT0ASYB12Rt+wEjgT5Z\n+19I2ngr8+uBWRFxKHAwsKzgM3sAJwEv1PucP0RE3+zxzW3+pmZmZWzcODj8cBgzBjZsyDuaLRXr\nWfQHlkXEqohYD0wHhtWrMxSYkj2/EzhaUpus3rSI2BARq0gJob+kLsBhETEVICLqIuKtguv9EPhq\nA7E0Op5mZlbJpDQc9Y9/wOWX5x3Nlooli+5A4cb0lVlZg3Uiog54DegKdMvq12+7P7BW0gxJSyVN\nltQZQNIIYGVELK73GQEMzIazHpR0aJO/oZlZhWjfHm65BebMgZ/+NO9oNlcsWZR6aVEAbYEjgUkR\n0Rt4HZggaSfgCuDKgvobexNPAN2zoa7vAreXOC4zs7Kwxx7plNpvfhPuvTfvaDYpts9iJdCj4HUP\nNu9pbKyzD7AmG37qAqxtoG33rOxFYFVEzM/KbwEmAB8BegJPSdpY/wlJR0XEmo0XiYj7Ja2TtFdE\n/LV+wBMnTvzX85qaGmpqaop8RTOz8vLhD8OsWWmF1AMPwKElHkupra2ltrZ2m9o0us9CUifgaWAQ\nsAZ4BLg4IhYW1Lkc6BERX5R0JvCpiBieTXD/HBgI7AXMBQ6IiPckLQDGRMQzkiYCu0fEF+t99vNA\nv4h4XdIHIuLVrLwfcAewTzbsVdjG+yzMrGpMnw5f+Qo89hjsvXfzfU5T9lk02rOIiHclXQLMIQ1Z\nTYmIhZKuAhZExF3AtcAUSUuAt4ExWdsnJN0GLAbqSEnmvezSFwE3S9qZtOppbJHvcp6kz2bP15ES\nTV1jDczMKt3o0fDss2lJ7R//CLvskl8s3sFtZlbGIuDTn4bXX09DU23blv4zvIPbzKzCSfCLX8Db\nb6chqbw4WZiZlbkOHeDWW9PqqOuuyycGnzprZlYBNi6pPeYY2G8/GDKkZT/fPQszswrxkY+kTXsX\nXABLlrTsZztZmJlVkEGD4Mc/Tseav/JKy32uk4WZWYU57zz4zGfSpr1//KNlPtNLZ83MKlAEfPKT\naZXULbdAmx346e+ls2ZmVUqCX/4SXnstHW/e3JwszMwqVMeOaaPeHXekxNGcvHTWzKyCdemSltQe\neyz07Aknn1y0yXZxz8LMrMIdcADMnAnnnw9LlzbPZzhZmJlVgWOPhR/9KB06+Nctbt6w45wszMyq\nxNixaYXU8OHwz3+W9tpeOmtmVkUi0g7vd96BGTOatqTWS2fNzFoZCW64AdasgfHjS3ddJwszsyrT\nsSPcdltaVnv99aW5ppfOmplVocIltfvtB4MH79j13LMwM6tSvXqlJbVjxsDy5Tt2LScLM7Mqdtxx\n8P3vw7BhsHr19l/HycLMrMpdcAF84hMwYkRaJbU9vHTWzKwViEj7MNavh2nTNl9S66WzZmYGpCW1\nN94IL78MX/vatrd3sjAzayU6dYLbb0+b9W66advaeumsmVkr8oEPpCW1xx0H++4LJ5zQtHbuWZiZ\ntTIHHgjTp6fbsz79dNPaOFmYmbVCNTXwve+lJbVN4WRhZtZKXXghfPGLTavrpbNmZq1cSZbOShoi\naYmk5ZK2uC24pI6Spmd15knat+C98Vm7JZJOLijfXdJMSU9JWiFpYL1rXi6pTtL7C8qukbRM0kJJ\nfYvFbWZmpdNospDUEbgOGAIcApzVwP+oLwVeiYg+wCTgmqxtP2Ak0Cdr/wtJ7bM21wOzIuJQ4GBg\nWcFn9gDKmQjHAAAFVUlEQVROAl4oKBsF7BMRBwMXAdu46Ku81NbW5h1CkzjO0qmEGMFxllqlxNkU\nxXoW/YFlEbEqItYD04H60yFDgSnZ8zuBoyW1yepNi4gNEbGKlBD6S+oCHBYRUwEioi4i3iq43g+B\nr27tMyJiEdBOUvdt+aLlpFL+A3KcpVMJMYLjLLVKibMpiiWL7sBLBa9XZmUN1omIOuA1oCvQLatf\nv+3+wFpJMyQtlTRZUmcASSOAlRGxeDviMDOzZlIsWZR6tjiAtsCRwKSI6A28DkyQtBNwBXBlQX1t\n5XlzxGZmZlsTEVt9AMcCdxe8/grwH/XqPAgckT1vA6wlJYQJwJcL6t0NDCL1CP6noPwYYA7QG1gN\nPJ893gP+B9gT+BVwVkGbpUC3BuINP/zwww8/tv3RWC6IiKLHfcwHekvqBqwBzgEurlfnXuB8YAEw\nAng0IjZIuhf4uaSrgb1IyeDxiHhP0quSekXEM8BgYEVELCUlBgAkPQ/0i4jXs2udD9wi6XBg4zzI\nZoot/TIzs+3TaLKIiHclXUL65d8GmBIRCyVdBSyIiLuAa4EpkpYAbwNjsrZPSLoNWAzUARdHxHvZ\npS8Cbpa0M2nV09iGPr4gjlslHS9pGfC/wKe2/yubmdm2qqpNeWZm1jyq4riPYhsHy4WkGyWtznph\nZUlSD0l/zP4+/1tS/WXMZUFSJ0nzJS2S9IykH+UdU2Mktc1ivSvvWLZG0v9IWpzF+Xje8WxNsU29\n5UDSgdnf48bHm5Iuyzuu+iRdlf37eVrSLdloT8N1K71nkW0cfJo0Ub4aeBT4bLYfo6xIOhb4OzA5\n28RYdiTtCXwwIpZmS5oXAmdHxFM5h7YFSTtFxDuS2gFzgfER8XDecTVE0v8F+gG7RsTwvONpSOE8\nYd6xNEbSTNKm3qnZnq7O9fZqlZUsxlXAURHxUrH6LUXS/sD9wEcjYp2k6cD9EfGrhupXQ8+iKRsH\ny0JE/An4W95xNCYiVmeLDYiIv5PmnPbON6qGRcTGuwl3IK3A24Hb0TefbAPpUOAGtlwCXm7KOr4m\nbOotR4OB58opUWReJ6063SX7wbVxDrlB1ZAsvGGvmUjqSdoTMzffSBomqY2kJ0lJ4uGIWJ53TFvx\nI9Ky87q8AykigN9nQ1GX5h3MVhzAVjb1lrFzgd/lHUR9WQ/yB8CLwMvAGxHxwNbqV0OyqOxxtDKV\n/QOcCXwhIt7OO56GZL8qDyP9ODhOUk3OIW1B0mnAmmxYtKx/tQMDIuJw4ETgU5IG5x1QA9rQwKbe\nfEPaOkkdgNNJ/5bKiqSPAF8EepJGDzpLamhlKlAdyWIl0KPgdQ8272nYNsoOfLwV+F1E3J53PMVE\nxJvAPcCAvGNpwNHA8Gw+YCpwgqTJOcfUoIhYk/25FriF9D/lcvMSsCoi5mevbwEOyzGeYk4Fnsj+\nTsvNUcAjEfFaNoQ/izT326BqSBb/2jiY/U/uHGB2zjFVLEki7ZhfHhFlu8JIUhdJu2bPdyKdVFx2\nq8wi4oqI6BER+5GGIx6KiAvyjqs+STtvXAkjaRfSSdHLGm/V8rJx/1cl9cqKBgMrcgypmPNIPxLK\n0bPAAEk7Zf/uB2dlDSq2g7vsbW3jYM5hNUjSVODjQBdJLwFfj4ibcg6rvkGk3fKLJW1cUTY+Iu7L\nMaaG7A1Mzv4j70TqBd2Tc0xNUa7DpnsCt0sK0kTntIi4M+eYtqYpm3pzlyXdwcD/yTuWhkTEfEm3\nsGnj9CLgp1urX/FLZ83MrPlVwzCUmZk1MycLMzMrysnCzMyKcrIwM7OinCzMzKwoJwszMyvKycLM\nzIpysjAzs6L+P0wgHdIrtQ0bAAAAAElFTkSuQmCC\n",
       "text": [
        "<matplotlib.figure.Figure at 0x7f3b07812790>"
       ]
      }
     ],
     "prompt_number": 13
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "e"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 11,
       "text": [
        "[(OrderedDict([('loss', 0.0085986440083445168), ('err', 0.0085986440083445168)]),\n",
        "  OrderedDict([('loss', 0.0086501840236289668), ('err', 0.0086501840236289668)])),\n",
        " 0,\n",
        " 0,\n",
        " 0,\n",
        " 0,\n",
        " 0,\n",
        " 0,\n",
        " 0,\n",
        " 0,\n",
        " 0]"
       ]
      }
     ],
     "prompt_number": 11
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "spent"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 33,
       "text": [
        "[18.0, 15.0, 38.0, 201.0, 210.0, 1113.0, 2674.0]"
       ]
      }
     ],
     "prompt_number": 33
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "trainer.save(\"network/500K4L-rprop\")"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 15
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "Predict = trainer.network.predict(valid[0])"
     ],
     "language": "python",
     "metadata": {},
     "outputs": []
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "gots = np.zeros(5 , dtype=float32)\n",
      "start = time.mktime(time.gmtime())\n",
      "for i in range (2000) :\n",
      "    n_goodTracks = int(4*got(ds_output[cut+i],Predict[i],0.1,1))\n",
      "    if n_goodTracks == 4 :\n",
      "        print i\n",
      "    gots[n_goodTracks] += 1\n",
      "stop = time.mktime (time.gmtime())\n",
      "spent.append (stop - start)\n",
      "print \"\\t\" , stop - start\n",
      "print \"gots :\\t\",gots\n",
      "print sum([i*gots[i] for i in range (5)])/(2000*0.04)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": []
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "counter = np.zeros(5 , dtype=float32)\n",
      "start = time.mktime(time.gmtime())\n",
      "for i in range (2000) :\n",
      "    n_good_tracks = compare(ds_output[cut+i],Predict[i],4,4,1)\n",
      "    if n_good_tracks == 4 :\n",
      "        print i\n",
      "    counter[n_good_tracks] += 1\n",
      "stop = time.mktime (time.gmtime())\n",
      "spent.append (stop - start)\n",
      "print \"\\t\" , stop - start\n",
      "print \"A-error: \",counter\n",
      "print sum([i*counter[i] for i in range (5)])/(2000*0.04)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": []
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "plt.plot(trainer.network.predict(train[0])[4])\n",
      "plt.xlabel('hit index')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 10,
       "text": [
        "<matplotlib.text.Text at 0x7f7551933050>"
       ]
      },
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAEPCAYAAABIut/fAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXuYXWV97z+/yYTMJCTkIohNIlADNZBSwaOoFR1btbRc\nelooaq3PqY9ai1Vr9UE87TmQtPUR2x59qlSs4qVHqeIRrxWKx9OOFBEKJOWSYGwVaBIil8kNkplJ\nZuY9f7x7uddes9bea++19prMer+f59nPvq3Z88477/vd3/V9L8uccwghhKgfA3NdACGEEP1BAi+E\nEDVFAi+EEDVFAi+EEDVFAi+EEDVFAi+EEDWlo8Cb2afN7DEzu7/NMR8xs61mttnMziq3iEIIIXoh\nj4P/DHBe1ptmdjHwbOfcGcCbGscLIYSYYzoKvHPuX4C9bQ75NeBzjWO3AINmtqac4gkhhOiVMjL4\nNcCO2POdjdeEEELMIWUNslriufY/EEKIOWawhM/YCawF7mw8X9N4rQUzk+gLIUQPOOeSJjoXZTj4\nm4DXA5jZ2cC0c25X2oHOOd2c46qrrprzMhwtN9WF6kJ10f5WhI4O3sy+ALwceIaZ7QCuAhY2BPtv\nnXM3mtkrzGwrMAm8sVCJhBBClEJHgXfOvS7HMW8vpzhCCCHKQitZ54CRkZG5LsJRg+qiieqiieqi\nHKxoxpP7F5m5qn6XEELUBTPDzeEgqxBCiKMQCbwQQtQUCbwQQtQUCbwQQtQUCbwQQtQUCbwQQtQU\nCTzw0Y/Chz8816UQQohyKWOzsXnP7t3w6KNzXQohhCgXOXjg8GF4+OG5LoUQQpSLBB44cgQeeWSu\nSyGEEOUigccL/M6dMDU11yURQojykMDjI5qpKeXwQoh6IYHHO3hQTCNEVTz9NBw6NNelqD8SeLzA\nL1gggReiKq6+Gq69dq5LUX80TRIf0ZxyimbSCFEVhw7BokVzXYr6IwePd/CnnioHL0RVTE1pUkMV\nSODxDn7dOgm8EFUhga8GCTxNB6+IRohqkMBXgwQeL/Dr1sGOHTAzM9elEaL+HDnSnL0m+ocEHh/R\nLFsGS5fC44/PdWmEqD9y8NUggcc7iWOOgZNPVkwjRBVI4KtBAo8X+IUL4aSTNNAqRBVI4KtBAo+P\naI45xgu8HLwQ/UcZfDVI4Gk6+JNPloMXogrk4KtBAo938IpoRL+ZmICDB+e6FEcHEvhqkMDTHGRV\nRCP6hXPwm78Jf/VXc12SowMJfDVI4Jk9yOrcXJfo6GBsDD74wbkuRT34+tfh5pth7965LsnRgQS+\nGiTwNCOa447z93v2zHWJjg5+/GO4/vq5LsX859AheNe74NJLFdFEaJC1GiTwNCMaUEwTR52wHK6+\nGl7wArjoIgl8hBx8NUjgaUY0oJk0cdQJi/OjH8Hf/A186EOwZIkuchGhtlUNwQv8zAxMT8NgY2d8\nzaRpIgdfnPe+Fy6/HNau9QIvB++RwFdD8AIfuXcz/1wRTRN1wuLcdx/8xm/4xxL4JjIP1SCBj8Uz\noIgmjjphceLtSwLfROahGoIX+GibgghFNE3UCYszNSWBT0Ntqxo6CryZnWdm95vZNjO7IuX955rZ\nnWb2QOOYX+9PUftD0sGvWqW5yhFy8MWRg09HAl8NbQXezBYB1wLnAWcCl5jZWYnD/gfwaefcBuBi\n4Jp+FLRfRHPgIxYvVieMUCcszpEjzQF8ta0mMg/V0MnBnwNsdc7tcs5NATcA5yeO2QEc13i8HJhX\nAUd8DjxoKlscdcLipDl4rZSWeaiKwQ7vr8ELeMROYCRxzAeA75vZO4AlwC+XVroKSEY0w8MwPu6n\nTw4EPkIRdULnmrOMRHfEM/jBQX+bnIShobkt11wjga+GTgKfx2t8CLjOOfdhM3sR8HngjLQDN27c\n+NPHIyMjjIyM5CtlH0lGNAMDvvONj3vHFTKRe4+vExDdkTQQkYuXwEvgsxgdHWV0dLSUz+rUbXcC\na2PP19Lq6AFeClwF4Jy7w8yGzOwE59ysq5vGBf5oIRnRQDOmCV3gow44NSWB7wXnZn85RgK/atXc\nletoQAKfTdL8btq0qefP6hRC3AVsMLPVZrYQuBS4OXHMj4BXApjZenxMM9ZziSom6bBAg2ERkYNX\nDt8bU1OwYEFrvKWZNB6N71RDW4F3zk0AlwG3APcCX3HObTazTWZ2YeOwdwO/b2ZbgRuBNzvnpvtZ\n6DJJzoMHdcKIuIMX3ZNmHtS2mmc2alf9p+OJt3PuZhKu3Tl3VezxduDF5RetGrI6oWbSyMEXRQKf\nznTD/kng+0/g80RmD7KCIpoIOfhixOfAR0jg1a6qJHiBbzfIGjpy8MWIT5GMkMD79rRokdpVFUjg\ndRqdiZxWMdS20pma8tNEp6e16KvfBC/wimjgfe9L/3vl4IuhiCad6MxmwYJmHi/6Q/ACr4gGPvtZ\neHzWqgU5+KLIwacTrasYHFTb6jcSeM2Dz5yTLAdfDGXw6UjgqyN4gc+aBx+Sg88SeDn4YsjBpxNF\nV4ODMg/9JniBVyeUg+8XyuDTic5sFi6Ueeg3EnhFNHLwfULmIR1FNNURvMCHHtFEy8bl4MtHGXw6\nEvjqCF7gQ3dZ7URcDr4YobetLJTBV0fwAh/6PPh2Ai8HXwxl8Okog6+O4AU+9HnwnRz8wIA6Ya9k\nRTShtK0sFNFUhwQ+8EHWTg5+eFgOvldCb1tZSOCrI3iBD32QtZODHx5WJ+wVZfDpKIOvjuAFPvRO\nKAffP5TBpyMHXx0S+MBPo+Xg+4emSaajQdbqCF7gFdG03iffk4PvnTTzsGiRLjgtB18dwQt8Wicc\nHvYCH8Je1YcP+/t2Aq9O2BtpEY2ZXLwEvjqCF/i0efALFninNTExN2WqkjwRjRx8b6SZB5DAa5C1\nOoIX+LR58BBOJ4w6WOTkk+9J4HsnLYOHcNpWFsrgq0MCn+GyQhlo1SBr/5CDT0cRTXUEL/Bpg6wQ\nzkBrp0HWxYvl4HslLYMHCbwEvjqCF/jQXZYcfP8IvW1loQy+OiTwimha7pPvKYPvHWXw6UQOXhl8\n/wle4BXRtN7HkYMvhhx8OtEXnyKa/hO8wMvBt94n35OD7x1l8Okog68OCXwblyUHLwdfBEU06Ujg\nqyN4gW8X0YTQCeXg+4cimnQ0yFodwQu8Ihq/cjfZ0aJrtQ4NyWX1igQ+HS10qo7gBT5tqwIIK6JJ\nm+s+NeWFf+FCuaxeUQafjiKa6ghe4LVVQbbAy2UVQxl8OhL46pDAt4loQnbwykmLo4gmHbWt6ghe\n4DXIKgffLyTw6ahtVUfwAq9BVjn4fpGVwYfStrJQRFMdHQXezM4zs/vNbJuZXZFxzKVmtsXM7jOz\nvy+/mP1hetrfL1gw+z0NssplFUUZfDoS+OpI8RdNzGwRcC3wUuAx4Ptm9m3n3JbYMb8AvAd4qXPu\noJmt7GeByyQrnoFwOuHhw17gn3669XU5+OIooklHbas6Ojn4c4Ctzrldzrkp4Abg/MQxbwSucc4d\nBHDO7Sm/mP0hqwNCOKfRcvD9Q9Mk05GDr45OAr8G2BF7vrPxWpyfA55nZneb2T1mdlGZBewnWXPg\nQRGNXFZxQt8GIwuZh+poG9EAeS47PQCcjHf7a4Hbzey2+eDks+bAQzgu68gR/7fKwZePMvh05OCr\no5PA78SLdsRaWh09jee3OeemgYfNbBtwGnBH8sM2btz408cjIyOMjIx0X+IS6RTRhOCyjhyBFSvk\n4PtBuxla4+MwMwMDAc5jk8C3Z3R0lNHR0VI+q5PA3wVsMLPVwOPApcBbE8d8C/h14LNm9gxgPfCj\ntA+LC/zRgAZZmw4+edFtOfjiZGXwAwN+j5/xcV/3odEP8/CZz8COHXDlleV83lySNL+bNm3q+bPa\n+gfn3ARwGXALcC/wFefcZjPbZGYXNo75KjBmZluB24D3Oeee6LlEFZLHwbs8IdU8Rhl8/8iKaCAc\nA5FGP8zDT37ib6KVTg4e59zNwM2J165KPH8PfqrkvKKdwEcbbU1OerdVVzSLpn+0a1+hC3zZEc3h\nw7PPQkXgK1nbRTQQRieUg+8fEvh0+iHwk5P+JloJWuDbdUAIY6BVDr5/ZGXwELbA98M8HD4sgU9D\nAt9G4EPohHLw/UMZfDqRgy/TPMjBpxO0wCuikYPvJ4po0onaljL4/hO0wCuiab3uanzGkBx8cSTw\n6SiDr46gBV4O3ovQokV+bna0uybIwZeBMvh0+jWLRgI/m6AFPo+Dr3snjOogee1VOfhiRBctl4Of\nTT/a1uSkIpo0JPAdBllDiGjSBD6ek05P13/BV9lEFy03S38/ZIHvx9mhHHw6QQu8IprODt7MC5Vi\nmu7QDK1slMFXR9ACr0HW9gIf1Y1y+O5pl7+DBF4rWatBAh+4y4r2xM+KaEA5fC+0mwMPYbStLPqV\nwcvBzyZoge8U0WiQ1T9Ovic6I/OQjTL46gha4DXI2nmQFbRvdy90imhCMA9Z9CuDV0Qzm6AFvt0l\n+yAMlxVd1UoOvlwU0WSjefDVEbTAt7tkH4Q9yCoHXwxFNNn0K4M/csRfJUs0CV7gQ++EyuD7g9pW\nOjMzfk3FwED5EQ2onSYJWuBDnwcfbU0QXdwky8FrmmT3aJpkOlG7Mit/kBUU0yQJWuBDnwefnOue\n5eA1TbJ7OmXwQ0NhilGUv0P5Dn7RojDrtB0S+IBPo+N//zHHtM5CkIMvRqe2NTQEExPVledooV8C\nf/gwLF2qmTRJghb40OfBy8H3Dwl8OvF2NTDgM/miA6MzMz5uXLJEDj5J0AIf+jz4dgIvB1+MThl8\nyBFN1K7KyuEjo6aIZjYSeEU0gBx82XTK4CNhi+/BHwLxiAbKiWkmJ5sCr4imlaAFPk9Ec+hQfbfK\nnQsH/4MflPM5RzudzINZmC6+HwJ/+LAXdzn42QQt8J064eCgn0JYV1dQtYPftQte8YrinzMf6NS2\nIMwcPhldldG2ohk0xxwjgU8StMB3cvBQ75imagd/6BCMjxf/nPlApwwevCiFJvBJB192Bl9XM9Yr\nQQt8HpdV57nwVTv4iYlwBK1TBg9hOvhkvXQT0Tz6KLzwhbNfjxy8IprZSOA7dEI5+PIc/MSE74B1\nHdOIkzeiCU2QimTwTz7pRT6JZtFkE7TA54lo6jwXfi4cPITRCZXBp1NE4Ccm0iO+eAaviKaVoAU+\nTyccHq5vJ4xvl1yVg4cwcvikkKURosAXGWTNivjk4LORwAfssubKwde1PuOE3raySGbw3ZiH8XFf\nX8mITxl8NkELfJ6IZni4vo5zLjL4+H2dkcCnUzSimZmZbTaifqyIZjZBC3zonTB+wRM5+HLRNMl0\nigp8/D5CDj6boAU+r4Ovayes2sFHZ0J1rc84miaZTpEMPqv9KIPPJmiBz+vgQ4xo5OCLoWmS6RTJ\n4Ds5eEU0s5HAB+yylMH3j9DbVhZlRDRJwxXfbCy0L8xOBC3wGmSVg+8XeTJ4CXx3At8uolEGn05H\ngTez88zsfjPbZmZXtDnuYjObMbOzyy1i/wjdZcnB9w9l8OkUnQcfv4+IHLwimtm0FXgzWwRcC5wH\nnAlcYmZnpRy3FPhD4I5+FLJfhL7QSQ6+f4RuHrLoR0QjB59NJwd/DrDVObfLOTcF3ACcn3LcnwFX\nA5OAlVvE/uCcBlnl4PtHnrYV6jTJIgudINvBS+Bn00ng1wA7Ys93Nl77KY1IZrVz7qbGS/NiK6mp\nKb/X+0CHGqizy5oLB79wYX3rM07eDD40QerHPPjIwSuimU2HJtherM1sAPgQ8N/iLxctVBXEF/m0\nI5RB1mTn6JeDX7EiDIFXBp9OPwZZJydh2TI5+DQ6CfxOYG3s+VpaHf1S4Axg1MwATgS+YWYXOuc2\nJz9s48aNP308MjLCyMhIT4Uug/hGW+2ocyfs5OCj98p08MuX1/cLM44y+HSKDrIuXlz/DH50dJTR\n0dFSPquTwN8FbDCz1cDjwKXAW6M3nXP7geOj52b2z8B70sQdWgV+rsnTAaHenfDIEb/fPbQK/MyM\nv0XxVZkOfvny+tZnHAl8OkUXOqW1n7rNokma302bNvX8WW0TaOfcBHAZcAtwL/AV59xmM9tkZhf2\n/FuPAvLMgYdwIpq4wEed0BphW9kOPgRR03bB6RSNaNIivro5+DLp5OBxzt0M3Jx47aqMY+fNJZXl\n4LMFPnkaXaaDP/74+tZnHDn4dIoOsq5Ykb6SVQKfTrArWfMKfIjz4NOumykH3x2aJplO0Qw+rf1o\nu+BsghX4vBFNnefBZ13RKc3BS+C7Q5uNpdOPiEYOPptgBV4RTXcOXtMku0MZfDpFB1mzMngtdEpH\nAt+BEAdZ5eCLoww+naIZfNo0W20XnE2wAt9NRFPXTph1Rad+OngJfJM6t60s0gS+mwt+yMF3R7AC\nr0FWOfh+ou2C00kbZC0a0SiDzyZYgdcga7UZvHNhCby2Kkin1ww+3n7SVrJGDl4RTSvBCnxeBz84\n6Fd1lhFRHG1U6eCj37V4cRiilneaZGiOs9cM/vBhf+ySJe0v2Tc56b8MhEcC3wGz+sY0eR18GQud\nJia8Yw3FtXYzDz4kQepV4Nu1n2irgmh32DqasV4JVuDzRjRQX1HK6+DLWOgU76B1jbzi5MngBwfD\nE6Re29b4uDdaaWYr2qoAFNMkCVbg8zp4CEfgp6a8m5SDL06eDB7CqY+IpIPP27baGYTIwUMzphEe\nCXwO6joXPl4HZv4Ud2qqPw5+fDwsgddCunR6HcBvZxCSDl4C3yRYgVdEM1uEophGDr44Evh0es3g\no4imXQYPimiSBCvwcvDZAt/PDH7hQpiern/unGerAghP4HttW1H7SfbF6NrKimjSCVrg5eCrd/DR\nrKS6d8K8BiK0SKHsWTTRhnnRtQtCq89OBCvweS/ZB+EJfD8dPNS3PiPSBqqzqHtdJOnVPCQjmmhq\naTx/Bwl8kmAFXhFNusAfPtxfBw/1F7WpKT9gbTkuP1/3ukhS1MEPDLRO6Y3n76ANx5IEK/AaZJ0t\n8MccIwdfBnnzd6h/XSQpIvDDw/5xfKqkHHx7ghV4Ofi5yeCh/qKmNRbZFFnolNZ+kg5eAt+KBD4H\nde2EyXGIrAx+YMDvxzMz0/vvksCnU/e6SNKrechqP0kHr4imlWAFXhFNfgdvVtzFhybw3UQ0ITnO\nMiKa+Bm1HHx7ghX40COarC0J0hw8FM/hQxL4vDNoILwLbxdZ6JTHwUvgWwla4EN28FNTPnoZiLWA\nLAcfvVemg6/bF2YcRTTZFF3oBLMzeEU02QQr8N3Mg6/jdsFpX3By8OUggc+mrIgm7uAV0WQTrMB3\n2wnr5jjT/v4qHXydRa3bDL7OdZGkyEKntDPApIOXwLcSrMCHPsjaTuDl4IvRTQZf97pIUtY8+CwH\nr4imlWAFPvRB1rl08HWMvOIooskmTeCLzoOXg89GAp+DOnbCuXDwaQ6sjnTbtkISpLRB1m7nwWua\nZH6CFXhFNMrg+0U3WxWEOE2y14VOWRFN0sEromkSrMAromnv4JPvKYPPT+hnh+0oex582mZjcvBN\nJPA5qGMn7OTgkw5UDj4/obetLKLFdQsWNF/rNaLRQqd8BCvw3UQ0cvBy8N0ggU9nenr24rpuFjql\n7SaZlsEromkSrMCH3gm7HWRNOvgvfhGuvz7/7wtJ4LVdcDpFzgzzblWgiKaVoAVeg6ytr7UbZE06\nrbvugrvvzv/7QhJ4zaJJJ6tdFdlNUrNo2hOswB861Dzl60SIEU0np7V3r7/lJcuB1ZGiZ4c33gjX\nXlt+ueaaNAe/YIF/PboEXxZZu0lqFk17cp5I1o+9e2HFinzH1lGQijr4ffu62x9eDj6dtGmSW7bA\n7t3ll2uuSTMOUSY/Pd0+1tIsmt7I5eDN7Dwzu9/MtpnZFSnvX25mW83sATO71cxOKb+o5TE9DQcP\nwrJl+Y6Xg2+9DiZ07+BD2k2yaAY/NtZd3c4XsuolT0zTzTz4bgT+wQfzHzsf6SjwZrYIuBY4DzgT\nuMTMzkocdgdwtnNuA/D3wIfKLmiZ7NvnxX0gZ0AVNZpOp5HzibTdNLMuug2zO+G+ff6WFzn4dNLq\nYs+e+gp8Wr10GmidmWkdM4tPkywyi2ZsDM4+u179OkkeiTsH2Oqc2+WcmwJuAM6PH+Cc+xfnXPS9\n+T1gdbnFLJdu4hloXsm9TtlemghlXXQbijl451qdlvaiaSIH39nBT0z4tmPmn7e76HY3Ec2TT/rP\nPngw3/HzkTwCvwbYEXu+s/FaFm8Fvl6kUP2mW4GH+sU0vWTwvTr4yGXFO2jdBV4RzWyKCHx8QkRZ\ns2jGxlrv60ieZpj7BMbMXg+cDbw87f2NGzf+9PHIyAgjIyN5P7pUehH4uolSkQx+ehqefto/zpM3\nx+MZqF9dJul2u+CkINU1osn64uu02Ck+wAqdV7J2E9FE9yedlO9nqmB0dJTR0dFSPiuPwO8E1sae\nr6XV0QNgZq8E/gR4mXMu9d8VF/i5pFeBr5uDz7qiUycHf+AALF3qHfn+/bBqVfvfFZrAdxPRRPuX\nO9c8wxkb87FBN58zH+g1g09rP1krWbuJaPbs8fdHm4NPmt9Nmzb1/Fl5Ipq7gA1mttrMFgKXAjfH\nD2gMun4cuNA592TPpamIXiOaOolSEQe/dy8sX+5veZymBD4bs1ZRmpz0gr9yZXeD2POBfkQ0RWbR\nRMIeCX0d6SjwzrkJ4DLgFuBe4CvOuc1mtsnMLmgc9hfAEuDLZrbFzL7WtxKXgCKaYhn8vn2+/las\nyCdCoQl8N9MkobU+9uzx4r5iRf1iml4FPhnRdLrgRy8RTV3J1QydczeTcO3Ouatij19Vcrn6yt69\nnWOFJCENsuZ18Ga9CXy0uCceS9SJI0fg2GPzHx8XrLExL/BLltRP4HvN4JPtp91Ft7uJaMbG/Bdp\nnQU+yK0K5OCLXfAj7uB7iWgGB/3U0yLbDx/NdJudJx38qlVy8HHSIpoyLro9NgannVZvgQ9yqwIN\nsha7ZF9RBw9NUavTIGJELwIfidLYmBf4oaF6Cnwvg6zJiCa+8DA5yDo46BdGTU+37jufxp49cOqp\ngWfwdUSDrNU7+OTGbnU7I4pTJIOPIho5+CbJ9jMw0Jx9lBxkjQat8+TwY2Ne4Ovs4CXwOambIJXh\n4JcvL+7g64gimnTaCXw3GTw0z6iTDh7yxzQhRDQS+JyENMja7wwe6hd5xelW4OM7SkYRTR0Fvt0g\nazcRDTS/FJMOHvLPpJHA15Q9e+Tg59rB1y3yilPUwdc5oul1oVNWxJfm4PPMpBkf9zn9s59db4EP\nbpA1WmZ/3HHd/VxIAp/XwZsVc/B1qs84RTP4aApv3Qb/imTwWQYhy8F3Evh4FHbgQL5B2flIcALf\n7VbBESFFNFXOoqkjRRx8NMjqXD0dfJkRzfj47HnwkC+iiep5wQK/7ca+fd2vjZkPBCfwveTv4BvU\noUPll2eukIPvH0WmSUbOcnq6fgJfZKFT8uI8Q0Pw1FNeoJNmLY+Dj58prVrV+rxOSOBzMjRUr6wu\nS+CjjpHsNHLw+Snq4Fet8p9RN4HPcvB5MvgTTmh9bXjYb3SXdO+QL4OPvkihKfB1RAKfk7oNCmYJ\n/KFDnQfCIgcPvj47bTkwMTG7zuss8L1m8M41o4O6Cnxa28oT0aQNsu7fPzt/h/wRTQgCH9wsmiIO\nvk6ClHXJvvHx9qfR0TjE0JC/LVjQeWxCDr490TTJQ4d8fQ4P+1x4fLx9dDHfKHOQNRL4NAefN6JZ\nudI/XrWqfgPaERL4nIQyyNrJwe/b5+OZiDzT+drNY64jvUY0cdExyz8Ndb5Q1gU/wD8/cCDdweeJ\naOIOfuVKOfjaIAfvyRL4mZn2nTBZf3lEKDQH32tEE8+FoX5z4cvaqgDaZ/CKaJpI4HNSt5WXWVd0\nit8n38ty8BL4VnqdRZOcyVFHgS/jik7QOYPXIKtHAp+TEAZZI8Hv1sF3EiEJfHvSIhrw9VynbLjs\nefBFZtEog68pimg8aSIUreSTgy9GrwKviMZz991+HUBEVkSTlcEromkigc9JCIOsZv61rKlsZTr4\nutVnnF4z+DQHXyeBzzvIetFFcMcdzef9mkWjQdYaIgfvyXKZCxe2X4wiB9+ZXqdJJjP4lSvrJfB5\nFjodOAC7d8ODDzbfbxfR9DKLJtoGIh7RSOBrggZZPe0EvgoHL4FvElJE02mh0/bt/j4u8Fm7Sfa6\n0OnAAf95kfuXwNcIDbJ65OD7R5mDrHUT+E4Z/PbtfqfXbdua72dFfL1GNMl6PvZY/z+rY3sMSuCn\np/0GRd1uFQz1EyQ5+P7RSwY/OVl/B58ng9++Hc4/f3ZEk+bge13olIzCzLzg13EmTVACv3+/XwLe\ny77PIUU0/XDwIV2TVQ4+nbwO/ld/FR57DA4e9K9lGQTnelvolLZzZF1jmqAEvtd4BpqrPNvN151P\nyMH3jyICn3TwdXKVeRY6bd8Op58O69b5x9HryS+GyDD0stApeaYEEvhaUETgzeolSmU5eG1V0Ipz\nvUU0hw61zuyA8Bz8zAz8+7/766SefrqPadLiGWi2p14WOiXPlEACXwuKCDzUa6C1LAffSYRmZtJ/\nV10Ffnra76XfzRXDFi2Cxx+HxYtb6ykkgT9yBP7zP73wHnssrF/vB1rTzAE0X+tlFk1WRFOns6UI\nCXwX1EmUenHwk5N+kDp+dZ1ly/w1buMrD+NMTvp6S+4XX6e6jNNtPAO+Lg4enC06S5f6OqrLlsHt\nBlmnpnwk83M/519bv947+E4C3+ssmmRd13WxkwS+C+oy0OqcF+QsIc9y8Pv2zR6kHhjwIn/gQPrv\natdBJfCeqH6SsYFZvVx8pww+TeCzIpp2GXyniEYZfE1RROOJnFTaVZjaOXhozd8j2g20hibwWSLW\njqh+0q4JWjeBz+vgTzsNHnrInzH2w8Erg68himg8aVdzimjn4CG9/tpNlcwS+Lp8WSbJiiHaMTjo\nz4ok8E2BHxqCtWvhgQeUwRdBAt8Fddkgq12MUKWDr0NdJuklogFfH0lXCfUS+E4LnbZvh+c+t/n6\n+vWweXMlPUj9AAAMxElEQVT7iKasDF4OvgbIwXs6CXzaewMDzUw4SS8OPnJZzuUv93ygiMCH6uAX\nLvTtZ2wMnv3s5uuRwHfr4HvJ4DXIWgM0yOrpxcFH75Xl4M3yXZhhvtHtHPiIRYuyBb4u0UG7zca2\nbfOLm+LTS9evh3/7t3QHHwl7tytZp6b8rK/kdiVy8DVAg6yetMv1RWQ5ePAdMU3ge3HwUJ8zojiK\naLJpl8H/5CfN/D3i9NO9GGcZhEWLul/JumePb8PJdQpRBl+3M8p5KfAf/Sh8/evd/5wiGk8RB59W\nf704eKhPfcbZu1cRTRbtMniYLfBRHp/VfuJb/sZpd2aYlr9HPzM05Gft1ImOAm9m55nZ/Wa2zcyu\nSHl/kZnd0Djme2Z2Un+K6vnXf4U//VN485v9suZuUETj6SWDBzn4Ttx6K1xyCbztbd3/7HHHwYkn\nzn69TgLfzsHDbIFftgxWr06PaMC3nzQHPzzs2+P+/bPfyxJ4qGcO31bgzWwRcC1wHnAmcImZnZU4\n7O3AbufczwN/CXykHwUFLwa/+7vwkY/Apk1w6aX5BWJmxi/GSROovJQV0YyOjhb/kAJU6eDTrsYT\nMTQEt9462rG884GPfxx+67fg7/4Ofv/3u//5m26C/ftHZ72evKrT7t3NXRaL4pzfHqAq09JuoRO0\nCnzUR9avb99+0hz8iSfC7/wOnHsu7NzZ+l7aAGvEqlV+y4jbboM//mPYuHH+jxF1cvDnAFudc7uc\nc1PADcD5iWN+Dfhc4/E3gJeYpS2hKc7Gjf607bWvhcsug+c8By6/PN/P7t/v97joZavgiLIc/NEs\n8Mcc0178y3bwt98+2qG01TEzA/fcA4884h/n4ZFH4HWv86bje9+DV7+6t9+9fDl897ujs16PHLxz\n8Nd/DWecAc96lt8z/W//FkZH4XOfg/e/H664wj9PbhsxNQU//CF8+9vwiU/Ae98Lr3oVPOMZ8MIX\nepf8e7/ny++c/1/+8Id+BkveekgyPe3r5Fd+xZdvYqI7B59H4IeH0x28mf/db3gDvOQlcP/9vs1/\n97vw+c/7vzuN44+HkRF4+9v9Z9x3n6+frVtz/9mZ7NoFX/0q/NM/wZYtsGNHNXl/p/H+NcCO2POd\nwEjWMc65GTMbA04AHiupjADceSd85jO+0qOvj+uug7PPhrPOgosvbn8hj6LxDPgG9dRT8K1vwfXX\nw+23wy/+Ipx3nu/Yz3xmsc+vAuf8aWgvDn5wsPwMfts2uPJK+NGP/IDaypX+dsIJfpBtwwY46aTu\nNu/qlt27fdu67jr/9x886P+edet8+zr3XH9bt67Z9vbvhw98AD75SfiDP/D3xx5bftlWrPBicOGF\n8MQTcNddXqD+8R/hG9/w7XDtWj+9cHgY/uiP/F7ql1ziRfaee7zAnXginHIKnHyyv3/Xu+D5z/ev\n79zphS+KPZcs8fU/MOBF+W1vgze+cfYg8Pi4n7v+4x/7zzzjDG8Q7r0X3vIWv3nam97ky/jud/v/\nb5bAn3hiev995zuzhTDLwYP/P11+OaxZA694ha+L5zzH99W3vCX9Zz7+cf//X7PGP3cOPv1pL/qX\nXw7Pe57/m4aG/N98553+9thjzXr92Z/1j08+2f9P7rwTPvUprxUvfrFvW/v2+Ta3ZAlccIH/377g\nBf7vL9sadxL4Ur9jXvayxoc6f4u7A7PmLXoef337du9g4iK6fDl86Uv+lPid7/Qd7NRT/bf6zIy/\nRZ31ySfhzDOLlX/JEh8NvfjF8PrXw5/8if/Hfe1r/lt/etr/k5Yt8w0hKnu0u+DAgD+DeOQR/02e\nJF4vMzOtPx+vm3j9dOMCxsbg4Ye9EPz2b6cfs3Jl9hfhsmXpX2LHH+9dyYYNzdcGB31neeIJLzZp\nXHCB//+Z+Y63dKk/hd6zx8+q+NjHvDjt3+/r88gRf1uwwJfluOP8z0Tin6yjqD6TmPkvnqee8rdD\nh3y88qUvedEz82K0fbsX1O98B666ypfLzP9vpqd9G7jvPu+A+8Xxx/v28prX+LGn6Iv5Na/xtyQb\nN/o9XG680bfX177WC9PSpdm/Y80aeN/7/BnA4cNNV+ycH/O65hovjmvXNuv4qafg0Uf966ec4gXv\noYd8/9u9G66+2n8pmHkn/R//4ftJWl0997l+4kQa69Zll/tjH/NfwO143eu8Cz/22M4G7JRTWp+b\n+S+oc8/19fqd7/i2Mj7u6+Kcc7werF7t/0c//rG/bd7s+9nDD/uzkje9CW64wf8/Ipzzbfsf/sG3\nrQce8G3qZ37GfyF94hPty5oXc20UwszOBa5wzl3QeH45cIxz7v2xY/5f45i7zWwA79yf6ZybSXxW\nzSYgCSFENTjnevL2nRz8XcAGM1sNPA5cCrw1ccxNwO8AdwO/Dnw/Ke5FCiiEEKI32gq8c27CzC4D\nbsEPyH7OObfZzDYBdzvnvglcA3zOzO4HngIyTv6FEEJUSduIRgghxPyl7ytZOy2UqjNmttbMbm38\n/dvN7L2N11ea2f81s/vM7BYzKzA7f35hZgvMbIuZfbPx/BQz+36jjr5oZj2sA51/mNlyM/s/Znav\nmT1oZi8KtV2Y2SYz+6GZ/cDMvmxmi0NpF2b2aTN7rJGARK9ltgMz+4iZbTWzzSlrkmbRV4HPuVCq\nzhwG3tZYBPZ84M1m9gvAJuBbzrkzgZsbz0PhD4FtNGdofQT4YKOOfoJfOBcCnwS+4pz7BeAMfJ0E\n1y7MbB3wBmCDc+65wDTwOsJpF5/B62Oc1HZgZhcDz3bOnQG8qfGzbem3g8+zUKq2OOcec8490Hj8\nNHAfsJrWxWGfJ5A6MbM1+L/9Ov/UFgAvcs59rXFIEHVhZquA5znnvgB+/Yhz7gBhtos9wBFgiZkN\nAouB/ySQduGc+xcguYokqx2cH73unNsCDDb6VCb9Fvi0hVJtC1RXzOxk4AXAbcDxzrkxAOfck/iF\nYSHwYeByIJpldQLwZOz9XYTRPk4FnjCzL5nZA2b2v81sKQG2C+fcHuB/4UX9UWAf8ABhtouIrHaw\nmi71tN8CrxFcwMyOBb4M/GHDqQWHmV0APN5wHtGU2VCnzg7gv+z/0jm3Ae9i/+fcFmluMLPnAO8C\nTgZ+BjgWeNVclukoJ9ln2mpsvwV+J7A29nwtrd9AtacxOHQjcH3slPMJM3tG4/3j8WsM6s5LgIvM\n7CHgC8AvAR8E4juDrMG3mbqzA9jlnLur8fzLwPOAxwNsFy8EbnfOjTVi3K8ALyPMdhGRpQ9JPe1Y\nL/0W+J8ulGoI3aX4QYMgaGy69ilgm3Puw7G3osVhNO5vqrpsVeOc+2Pn3Frn3CnAa4F/cs69AbjD\nzP5r47BQ6mIH8KSZndZ46ZXAg/i+EVS7AP4DeJGZDTf6yyuBHxBgu4iRpQ83Aa8HMLOzgWnn3K52\nH9T3efBm9qv4bYSjhVIf6OsvPIows5cCt+IHV6OK/u/Av+IHnJ+JnyFwqXMuYz/G+mFmLwfe45y7\nyMxOAf4ef2q+FXiDc+7InBawAhqzqa7DDyo+gu+4RoDtwsw24v/+GWAL8LvAswigXZjZF4CX489Y\nHgOuBL5ORjsws2uAVwCTwJudc5vbfr4WOgkhRD2Zl5fsE0II0RkJvBBC1BQJvBBC1BQJvBBC1BQJ\nvBBC1BQJvBBC1BQJvJg3mNnJ8W1VE+9tMrNfajx+l5kNZxz3STNb3+Xvfbr70gox92gevJg3NDZs\n+2ZjC9l2xz0E/Jdow6YSfu9Tzrk2l60W4uhEDl7MNxaY2ccbuzCOmtkSADP7rJldbGbvwG9a9c/m\nLwjfQuNnzm48ftrM/rxxAZItZvasxuunNZ7fY2Z/nvj5KxsXYnjQzD7QeO03zOw7jcfPMn9xl9rv\nBCmOfiTwYr5xKnBNYxfGx4DfarzuAOec+yh+29kR59wvp/x8/JR1MXCbc+4s4Ns0Lyj/UeAvnHPP\nx29jC4CZXQSsblyI4Qzg583slc65rwK7zeztwCeAK51zIWwUJo5yJPBivvFQdBEV4B5ad9frlsPO\nuX9M+awXA19qPP5i7PhXA682sy2N40/Db3ML8A78PkMTzrkbCpRJiNIYnOsCCNElk7HH0xQzKfHN\nq2YSn5U1OPVnzrlPp7y+tlGeZ5qZOQ1uiaMAOXhRR8aBJQV+/nb81tYAr4m9fgvwRjMbAjCzZ5rZ\nMxqXmvsUfhvkHwDvLvC7hSgNOXgx30g64zSn/Cn8IOsjGTl82s+62PN3ADeY2eXAt6LXnXPfNLPT\ngc1mdhiYAH4buAy41Tl3u5ndB9xlZv/gnNve7R8nRJlomqQQQtQURTRCCFFTJPBCCFFTJPBCCFFT\nJPBCCFFTJPBCCFFTJPBCCFFTJPBCCFFTJPBCCFFT/j8SD9kze/I2swAAAABJRU5ErkJggg==\n",
       "text": [
        "<matplotlib.figure.Figure at 0x7f7553b5bb10>"
       ]
      }
     ],
     "prompt_number": 10
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "plt.plot(ds_output[80004])\n",
      "plt.xlabel('hit index')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 11,
       "text": [
        "<matplotlib.text.Text at 0x7f10d6b30310>"
       ]
      },
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAEPCAYAAABIut/fAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnX2wJFd53n/vvdKuWEloASFwpEVQYGIHmS/HEbgIXGzs\nUgzCiUWU2A6VgKkCUlAYU4QkELRbdsoxroCNqRKJAbuCA5aDZRCxFEGCrzFflkAC9AGiAgQkxZH4\n2qtdSXv3fpz80dO6Pb3dPad7zuk+PfP8qqb2zkzfub3vOf3Me57zntPmnEMIIcTisTL0CQghhIiD\nBF4IIRYUCbwQQiwoEnghhFhQJPBCCLGgSOCFEGJBmSnwZvY+M7vHzG5pOOadZnabmd1kZs8Ie4pC\nCCG64JPB/wFwSd2bZnYZ8Djn3FOAX5kcL4QQYmBmCrxz7q+AHzQc8nPA+yfH3gycZmYXhDk9IYQQ\nXQnhwV8A3Fl4ftfkNSGEEAMSapLVSs+1/4EQQgzMaQE+4y7gEPDXk+cXTF6bwswk+kII0QHnXDmJ\n9iJEBn8t8MsAZvZMYMc5d3fVgc65hX38xE84Pv5xv2OvuOKKwc83lYdPLJ7/fMfVVw9/rinEIoXH\nNdc4nvtcxSLG4xOfcKytTb82DzMzeDP7IPA84FwzuxO4Ajh9Itj/yTn3p2b2fDO7DdgEXjbXGY2U\n7e3sIcKj2KaF2iMe29twWghfZcLMj3LO/aLHMa8JczrjRZ0+HoptWqg94hFa4LWSNRBtOv3a2lrU\ncxkTPrFYFkEZS7/ooz3GEovQSOATRQLfDQn8HmPpFxL4eEjgE2VZRGgIFNu0UHvEY3sbVlfDfZ4E\nPhA7O9lDhEexTQu1Rzx2dpTBJ4mymngotmmh9oiHLJpEUaePh2KbFmqPeEjgE0WdPh6KbVqoPeIh\ngU8Udfp4KLZpofaIhwQ+UdTp46HYpoXaIx6qokkUVRbEQ7FNC7VHPFRFkyjKauKh2KaF2iMesmgS\nxLnsm1edPg4SlLRQe8RDAp8g+XBVnT4OEpS0UHvEQwKfIHlnV6ePgwQlLfL2mHOrclGBBD5BJPDx\nkP2VHnlb7O4Oex6LiKpoEiS3aFRZEJ5cRBTbdFB/j4eqaBJEGXw8FNv0UJvEQxZNgqjDx0OxTQ+1\nSTwk8AmiDh8PxTY91CbxkMAniDp8PBTb9FCbxEMCnyDq8PFQbNNDbRIPVdEkiBY6xUMVG+khgY+H\nqmgSJO/oEqHwSEzSQ1+68ZBFkyASoXgotumhNomHBD5B1OHjodimh9okHhL4BFGHj4dimx5qk3hI\n4BNEHT4eim16qE3ioSqaBFEVTTw0oZceEvh4qIomQVRFEw+JSXroSzcesmgSJG8UiVB4FNv0UJvE\nQwKfINvbsH+/OnwMFNv0UJvEQwKfIOrw8VBs00NtEg8JfIJsb8MZZ6jDx0CxTQ+1STxURZMgOzvq\n8LHIY6sJvXSQwMdDVTQJknd4iVB4JCbpoS/deMiiSRCJUDwU2/RQm8Sjd4E3s0vM7BYzu93M3lTx\n/o+Y2V+b2a2TY34+3OmNA006xUOxTQ+1STx6FXgz2w9cCVwCPBV4iZk9o3TYW4D3OecuAi4D3hXu\n9MaBOnw8FNv0UJvEo+8M/mLgNufc3c65beAq4IWlY+4Ezpn8fBD4VrjTGwcassZDsU0PtUk8Qgv8\nrI+6gEzAc+4C1krH/CbwWTN7LXAm8NPBzm4kqIomHprQSw8JfDx2dsKWSc4SeOfxGW8H3uOce4eZ\nPQv4I+ApVQcePnz4oZ/X1tZYW1vzO8vEUYePh2KbHmqTeGxvww03rPO+960H+bxZAn8XcKjw/BDT\nGT3Ac4ArAJxznzOzM8zsPOfcveUPKwr8IqEyyXhITNJDo6p4bG/Dc5+7xktesvbQa0eOHOn8ebM8\n+BuBi8zsfDM7HbgcuK50zNeBFwCY2Y+S2TTf63xGI0STTvFQbNNDbRKPXidZnXMngFcD1wNfAq52\nzt1kZkfM7NLJYb8GvMrMbgP+FHiFc26pvtvV4eOxvQ2nnw7Owe7u0GcjQP09Jn1PsuKcu45S1u6c\nu6Lw8x3As8Od0viQjRCPXODz7Wn37Rv6jJab3d3sy3bfPvX3GGgla4JI4OORd3jtP54GeXucfrra\nIwbabCxBdnayjCbPbkQ48s2XTjtNk3opoC/cuGizsQQp2ggSobBIUNJCX7jx2N3NHisBVVkCHwCJ\nUDwU27RQe8Qj//I0C/eZEvgA5L7Z6qo6fWgU27RQe8Qj9AQrLLHAb22F+yxlNd3Y2po9Z6HY+hOy\nT9eh9vBne7tdaa8EPiAXXwzf/naYz1Kn78aLXwyf/3zzMYqtH/ffD096Uvy/o/bw5/Wvh6uu8j8+\ndAUNLLHAf//7sLER5rOKE0/q9P784Adw9GjzMZrU8+OBB+Cee+L/HQm8Pz79u0joChpYYoHf2go3\npC12eomQPz5tIEHxI2R/bkJfuP60bRNZNAGJJfASIX8k8OHY2tors4uJ2sMfCfyAhBZ4VRa0x1fg\nFdvZ5HGMncWrPfyRwA+IMvjhUQYfjj4FXu3hhwR+QCTwwyOBD4cEPj26CLyqaAIRUuBVRdMNnzbQ\npJ4fEvj0aKsxqqIJxM5OtsBGGfywKIMPhwQ+PWTRDMTJk9P/zovKJLtx8uTsNpCg+BG6T9ehEZU/\nPv27iAQ+EKGzHVUWdENVNOFQFU16KIMfiBgCryyzPbJowiGLJj0k8AMhgR8e57JYSeDDIIFPD1XR\nDEToi0FVNO3J46QqmjBI4NNDVTQDoQx+eHzbQLH1QwKfHrJoBiKmwCvL9EMCH5a+BF4jKn8k8AOh\nKprhaSPwiu1sVEWTHhL4gYhZB69O70cee2XwYfCN57yoPfw5eVICPwjy4Icnj70WOoXBN57zovbw\nZ2ur/UInVdEEQFU0w+PbBvJ8/dAka1r4lgEXURVNIJTBD48mWcMigU8L3zLg8u9I4AOgKprhkcCH\nRVU0adGlPSTwgVAGPzwS+LAog08LCfyAqExyeFQmGRaVSaaFBH5AlMEPjzL4sCiDT4uuAq8qmgBs\nbYFZuJIyVdG05+TJrA1URRMG33jOiwTejy7toSqaQJw8CQcOKIMfkq2trA1UBx8G33jOi9rDjy7t\nIYsmEHnwQwj87m72WFlRp2+DbxtIUPwI2aebUHv40aU9JPCBCHkx5MMqM9kIbZDAh6UvgZdl5sdo\nBN7MLjGzW8zsdjN7U80xl5vZzWb2ZTP7QNhTDE/Ii6E4MaLKAn/aCLyqNmbTZwa/upqNWPOb14tT\nSUXgGz/OzPYDVwLPAe4BPmtmH3PO3Vw45mnAG4DnOOfuN7NHhj3F8IQW+LxRTjsNHnhg/s9cBvI2\nOHas/hjnMvtrdVUZ/Cz6tmjMsnaJMTG4CHQV+DPOCHseszL4i4HbnHN3O+e2gauAF5aOeRnwLufc\n/QDOue+HPcXwxLBoQCLUBp822NnJRET212z6FnhQf29iaysT63yOzoe8v4dklsBfANxZeH7X5LUi\nfxt4upl93sy+YGYvDnmCMYiZwavD++HTBoqtPxL4tNjagtNPzx6+bTKEB+/jsK0AjyfL9i8D3p26\nTROypEwdvhs+paqKrT+hS3/rUJv4cfIk7Ns3vMDP+ri7gEOF54eYzuiZPP+Uc24H+D9mdjvwZOBz\n5Q87fPjwQz+vra2xtrbW/owDcPIknHlmnAxeNoIfW1uz20Bi4o9PPENQtiTV36vJM/h9+9oL/Pr6\nOuvr60HOY5bA3whcZGbnA/cClwOvLB3z58DPA39oZucCPwp8verDigI/JCEvBlXRdCNvg6ZRlGLr\nTx7P++6L+3fUJn4ULRpfpyAX+HLye+TIkc7n0WjROOdOAK8Grge+BFztnLvJzI6Y2aWTY/4M+J6Z\n3QZ8CvjXzrnvdD6jHpAHPzzy4MMiDz4tUvHgZ36cc+464LrSa1eUnr+BrFRyFGxtwTnnSOCHRAIf\nFgl8WnQVeG02FgCVSQ6PTxmZ/F5/JPBp0UXgtdlYIGTRDI/PBaDY+iOBT4tULBoJ/JyoiqYbEviw\n9L0XDai/NyGBH5DQdfCqKmjPyZN+Aq/Y+tFnHbzaZDap1MEvpcAXL4Z5N0tSltmNra3sAmiqE1Zs\n/emrDl5t4sc8dfAhWUqBzyf48h3x5kEdvhs+dcKKrT8+6wpCoDbxo2sdvKpoAtDFH6tDVTTd8GkD\n+b3+aJI1LVRFMyAhBV4dvhuaZA2LBD4tNMk6IBL44ZHAh0UCnxYS+AGJKfCyEfyQwIfDuSw2D3uY\nyiRTQQI/IKEFXmVj7fEVeMV2Nnmc2lRszPu3QG3ShAR+QGTRDI8y+HCE7M+zUJv40bUOXlU0ASjW\nYM9bVqYqmm7kF0BT1ik7wA+fWIZCAu9Hlzp4VdEEwmcVpS/q8N1QBh8OZfDpIYtmQGTRDI8WOoWj\ny6KarqhN/Jjnhh8hkcCrimYQlMGHI4/l6mrz9sshkG3mhzL4AVEVzfCoiiYceSzN4ts0ahM/JPAD\n0hR85+BrX/P/LGWZzfzgB3Dvvae+rgy+G1//+qnxymMJ/Qj82NrkG9+IPzdRRnd0GpCm4N9xB1x6\nqf9nqYqmmfe+F972tlNf11403Xj5y+HTn55+TQLfzMtfDp/6VL9/U3vRDEhT8B94IHv4MsYO3yd1\n8VQG342qeErgm2l7TYcgL8WWRdMzeRa4ulo9w725mT18GWOH75O6ePre8COPbe73zrt//9ipimce\nS5DAV9H2mg5Bl1JsCXwA8kUhUL0I4cSJ7OGLqmiaqYtn2xt+rKxkj5gVImOgKp55LCHuYifnsv6d\n+8Rj6e9tr+kQ6IYfAzFrONtF4FVVUE+TwLepogHFF+oFvo8Mfnc3q9RZmajGWNpjSIFXBt8zsy6G\nzc3sNd9McYxD1j7xEXifhU6g+EIWy7LdUO7TsRY7jbU9qmIWG93RaSB8Mnjw7xDFme88s1l2G6FI\nlf+Zb2/bpooGxmMJxGRzc7gMfqwCXxWz2KiKZiB8Bd63Q4y10/dFVQafZyqzFuYotqcypEUz1vYY\ng0WTJ4UrgRVZAl9h0YAEPhTzCJJiO83ubjbcn2XRSOD3cG4cGXwM/x0k8HNbNFWdftlthCJVFo0E\nvhu5lztUBj9GyyyP2RAefJs6eAl8IGZNSHWxaFTpUU9VBu9bt63YTlPXN/uqgx9je7S9nkPRtg5e\nAh8Inzr44r+zUJbZzDx124rtNHV9s686+DG2x1AC37YOPkYFDSyhwPt68F2qaGAcnb5P5rFoxmgJ\nxKSub8qDr6ft9RyKth58jAoakMArg4+MzySr6uD9aMrgVQdfzRAZfLkM2Kc9ZNEEQgLfL6qiCYev\nwCuD32MIgfctAy7/jgQ+AKEtmjF2+j6pWhksge+GLJr2DGHRdGkPCXwgYmTw5cqCZfaJy1SVnbYR\n+LFVbcRk6Ay+uNEYjKOvD5HBS+AHxEfgzzhDFk0oquKpDL4bdX1TGXw9ba/nEBSrmpIXeDO7xMxu\nMbPbzexNDcddZma7ZvbMsKcYllkTUpubcPCgqmhC4FwW73POmY6nb922qmimqeubfdbBj62vt72e\nQ9ClPcqjo1A0CryZ7QeuBC4Bngq8xMyeUXHc2cDrgM+FP8WwzKoZPnEiEyRl8POzuZnFuCqDVx18\ne+r6purg62l7PYegmES2qYMfIoO/GLjNOXe3c24buAp4YcVxvw78B2ATsLCnGJZZ364S+HCcOAH7\n98uiCUWTwCuDr2ZogU/dorkAuLPw/K7Jaw8xsWTOd85dO3kp6Zuq+VTRtBnSjbHT98XmZibuZ5zR\nfZJVsd2jrm9K4OsZwqJJSeBnfWSjWJvZCvB24J8XX573pGLiM8l68OB8Gfwy+8RF8gmuWRm8Fjr5\nUdc3+1roNMY5kbbXcwi6tMdQAn8XcKjw/BDTGf3ZwFOAdTMDeCxwjZld6py7qfxhhw8ffujntbU1\n1tbWOp30PPgIfFuLRqV81eQWzf793S0axXYPX4smlpiNsT3GaNGsr6+zvr4e5FxmCfyNwEVmdj5w\nL3A58Mr8TefcBvDo/LmZ/QXwhipxh2mBHwofi+acc+DBB/0+T1lmPb4ZvCwaP3KxkkXjT13MYtJV\n4PMvz3Lye+TIkc7n0ujBO+dOAK8Grge+BFztnLvJzI6Y2aWd/+qA+Fo0KpOcn3k9+DFaAjHZ3ISz\nz876VzEOEvh6NjfhrLOyldR9nWuxqml1Nfvbs27jGWuzsZkf6Zy7Driu9NoVNcc+P9B5RcNnP3hV\n0YShzqJpsx+8YrvHiROZwO/fnwnXgQPZ66qDr+fECXjUo/ZiFkNEyxTbo7gfzf799b+jlayB8LVo\nJPDzI4smLMV4dhkRzcsY26MuZjEptgf4tYkEPhBNN/zY3c2eP/zhumVfCJosGi10ak8xnl0Wjs3L\nGC2zupjFpCzwPm0igQ9EU7ZTt/KyiTFWFvSFqmjCMm8852WM7VEXs5gogx+QWQLf9tteWWY9qoMP\ni69Foxt+7JGKRTOrTXTLvkA0CXyXzqAqmnpURROWJotGHnw1KVg0Pm2iW/YFYpbAtx3OjbHT94WP\npdBURqbYTpOCRTO29hjKosnnREAWTa/IoukPH4um6bZmiu00qqJpTyoWjQS+J0JbNGPs9H3hY9GA\nBN4XWTTtGcKiKdbBgwS+V5ompLpaNGO7jVlf5F+YTQudoFngx1a1EZN5F47Ny1hv2acqmiWi6nZa\nbrJnpiyasDTtB1/0KOvqhBXbaZosGt3wo5oULBrfOnhV0QSgmO2srsLKyl4WoiqasMxr0aiKZhpZ\nNO1RFc2S0RT88nDOzbh1SV75sVKI4hg6fV/UWTTy4LuhKpr2yKJZMpqCn3/br6z4L04oN8oYOn1f\n+NyyD+pjPUZBiYkWOrUnBYumq5aEQAJfyuDPOCP72WdIJ4FvxkeQQBm8L/Nu3jYv5fbIR66ztsId\nkrqYxUR18APiY9HA3vaiTVRNjIyhsqAv8hFRF4tmdzezyIr21zJX0TiXxbNsNziXxWSIKhpIu7/X\nxSw2smgGxMeiAWXwIWhj0ZQvAMV2mq2tTExXV6dHRHmSYZM7IfeZwUPabbK9ncXltNPGUQevKpoA\nNPljbS2aqpnvlDt839QNj30ugLrYppotxqaub3YRk66MTeDLMUt5JauqaAIR2qIZU4fvm6JF09aD\nV2ynya0GmLYbuohJV8bWJnUxi40smgEp3vADphchyKIJyzwLnRTbaeqyUd9FYyEYW5u0HZGHQjf8\nGBBV0fSHT9UHKIP3oa5vKoOvZ0wWjQQ+EH1YNMvqE5dpY9GU64THJiax8bVo8qqWGKWLY5sXScmi\nUR18T4SuoqkqG1tWESpTXMm6ubm3Mtg3g1ds92iyaIqxbNp+eV7G1iZDWjRd6uBVRROAkBaNqmia\nyUdE5ZXBPgI/tmwxNr4WDcQV+DH19zFZNKqiCUQbgVcVzXyUR0R1Wac8+NnUjS4l8PW0HZGHQh78\ngDT5Y209u7F1+D5x7tQ5jTa124rtNHXzQ+VYggQ+p67/xUY3/BiIfHhf9LpURROH4spLaF/5odhO\nI4umPWOyaCTwASjXwMN0jaosmnAUh8fQvnZbsZ2myaJp6tMhGVubpGLRqA6+J2ZlO10smjFtvtQn\nxS9LaL/6cmwVG7Gps2j6zODHttnYUBaNMvie+OIX4Wtf23s+62JQFU0zH/6w/zC3eHFBe4sm9Sqa\nq6/ut53HaNEcOwbXXhv+PHwZk0VT9eUZgoUW+Pe8B666au95W4GXRTPN614HX/mK37F1Fk15e1sY\n50KnV74SvvGN/v5e2yqaGDf9aNsmN9wAb3lL+PPwZUiLplwHP9RCpwgfmQ4bG9Mi4yPwqqKpZ2Mj\ne/hQZ9GUt7eF8U2yOtcuFiGo65spZ/B9x6iMLJoFF/ijR6dtAh8PXlU01ezuwn33ZTH1oc6iqWuD\ncqxTju2DD2b/D99YhKD4hbl/f5YROpe2wB892m+MypRjtrWVxayYXMQgJYFfaIumnEHMGs6qTLKe\nY8f2Mlcfyhl8btH41m2nHNs8Bn1n8Hk8iyuDU66Dz6+/WTevj0UxZmaZbdKHD686+J7wFfh5Nhsb\nU1XBPLQVtbIHnw+RfTPOlKtohhD4YoUXtI9nCNpW0WxsZO898ED4c/GhLmaxUQbfE+UhYlPNcPH+\njaAqmjJ5HENYND512ylX0bSNRQjqRkQp18EPEacidTGLTZc6eFXRdKCcwTcNZ7e2ss5atfKyjpRt\nhNC0zVqrLq62GXyqsR3aooH28QxBF4um+G/f1MUsNqPL4M3sEjO7xcxuN7M3Vbz/RjO7zcxuNbNP\nmtkTwp9qO/JJQV+Lppxx6pZ90+Rx9M3GqiyaPOMcu8APkZmmYNGMLYMfwqKpKwNOVuDNbD9wJXAJ\n8FTgJWb2jNJhnwOe6Zy7CPgA8PbQJ9qWY8eygJUtmrqLoapuWxn8HkePZrFSBp/FoE0sQtBk0aQq\n8EPEqcgQFk0eo1llwHW/FxqfDP5i4Dbn3N3OuW3gKuCFxQOcc3/lnMtD92ng/LCn2Z6NDXjMY7JM\n3mdZd5fhXMoiFJqNDXjc49oJfJsyyTEtdGobixC0/cJMYaHTEHEqMoRF07U9hhT4C4A7C8/vmrxW\nxyuBj8xzUiE4ehTOOSd75B0shkWzTFU0j3tcvxZNqlU0R4+2i0UIUrBo2lbRDBGnIkNYNF3bY8iF\nTt5VrGb2y8AzgedVvX/48OGHfl5bW2Ntbc33o1uzsQEHD2aNfPQonHde87erLJpmjh6FCy+Ez3zG\n7/i67GlR6uAvvDDb66gvmtYVHDgwfWxKFs2FF6aVwce2aLquSygmNOvr66yvrwc5Hx+Bvws4VHh+\niOmMHgAzewHwZuC5zrnK/05R4GOzsZFl7ydO+GfwbQV+mcok84v1uuv8jm9r0YypTHJjA574RPjL\nv+zvb46tiiZfFCeLxn+zsTy25eT3yJEjnc/Hx6K5EbjIzM43s9OBy4Gpy3wy6fpu4FLn3Hc7n01A\nihZNPkRUFU13hrBoyrHN7YChVkbmLKtF06a/P/BA9t6jHy2LJulJVufcCeDVwPXAl4CrnXM3mdkR\nM3vR5LC3AWcCHzKzm83sw+FPtR25RXPw4F4G0XTDD1k0zRw9Co997PSkdRNNGafPwpyq2JqlMceR\nf9ndd19/XzZDL3Ta3c0eZQ++rr9XXX99M4RFUyXwQ97ww+sjnXPXUcranXNXFH7+mcDnNTd9WDTL\nJPD5BZtPWp93XvPx81o0dR0+j2+Mi8GXjQ0499zs/3f//XDWWfH/5tAWTT7BWt6oq0ngy0UOfTMm\ni0ZbFbRkXosmD3aTWDdV0QxtI4Qmv2APHvQbcseoooE0Kmmq+lZshrZo2rbHEDEqM5RFUx5RSeAj\nUDVEbLPQCWZ/49fZCCsr2XB2kTh6dDqDn8W8dduzMvghGcJ+aLvQKXQdfNWkN9RPfMui2WNWezin\nvWhaUzVEbGPRwGyBb+r0Q4tQaNoOuWNU0cDwlTS7u9kq6Yc/vF/7YWiLpu0XriyaPWa1x+5ulhTG\n2Kd+YQV+XosGZg/pUs4yQ5KPcM4809+iqbqj07xVNDB8bI8dy+rOV1f7sx+2tzMRKMZjCIumTXsM\nbdHs7Jy6J0wfFk1VHfzq6t4kdRUx55QWVuDbWDRVC51g9pAuVREKTZ6NmflnZHVVSfMsdILhY5v3\nK+jPfshjWczw2t5AZV66ZPAHD2Yjnfvv73/UVRezITJ4s+Y2kcB3oMtWBSE8eBhehEKTCzz4Z2Sx\nq2iGIu9X0F922tQ3U8/gV1ayKqP77gt7PrOoi9kQHjxI4INT9ADzi7CpDl4CX09R1Hyz1hh18DB8\nbMtfdn1k8LMEvo86+K4ePAzjw3e5nkNQJ/BNbSKB70CXKpoqD36WRVNXOjb0YpyQFG2JrhbNopRJ\nDmXR1PXNvuvgy9T19SHiVKQuZkMJ/KwMPkYFDSyowBcnBVVFMz/FbKzrJOuiVNHIoplmlkUDw0y0\nDmnRlEdU0NwmdX09BAsp8FWTgs7Jg+9KXgMP3csku9x0O8XYDpGZjlHgh87gU7No5MEHpJhxnn76\n3pLyGBZNiiIUmi5+atmiye9y8+CD1WVkOzvTZWSpxnYIb7mLRRN6odPYPPgULZq6NpHAt6Q4PIS9\nDqYMvhshLBrInh875ldGlmpsZdFMI4tmmqqyVVAGH5Ti8BD2RKkp25HA1xPCooF6gYfxCHxqFk3q\ndfAgiwYk8EEpZpzgl8GriqaetsPt7e1szqPcaffvbyfwqVbRjMGiGbKKZncXjh+Hs8/OnsuiURVN\nUKosmqYMfmsr84ZVRVNNlag17ZZZtYoQ2mXwqqLZY2wWzX33ZRVsuWgtk0XTReBVRdOSKotmY6N6\nodPqarba7oEHZNHUUbRoipPWdVTZM5DF8/jx6jKy8kKQVGNb7Ftnn531m9jnM7aFTnXXX5+kZtFo\noVNA2kyyQvbasWOqoqmjbHnNmmit2tcH2ls0Kca22LdWVjKRj70MPwWLpk171I2g+2RIi6ZtHbwE\nviVVHnydRQN7Aq8Mvpq6OY06qrInWJxJ1jaxCMHYLJohYlRmTBaNBL4ldUPEpuAfPy6Br6No0YCf\nwDdZND51wqnGdgj7oYvAD1kHL4vm1NdVBx+QPi2aFCs9QuLccBZNarEtboGR04f90GQ3lPc8h3gZ\nvG97pG7RxLydZtc6eFXRtKCLRbOzE7aKZlHKJE+cyLzm4sUyj0WzszPeKpqNjWx/82J10FAWTb4y\nuOpG2LHKJH3bI1WLZnU1foKgKpoe6GLRgCyaKsr2DMwn8DBeD77cr2A4iway19qKSVfmsWgOHMiO\nie1/F2mKWUybRh58D3SxaEC37KuinI3BbIumzoPPXxurwJf7FfRjPzTFM0WBL8epzZ3AQtEUMwn8\nyKmzaKrq4CF7bWXl1CDrln3VAj/rYq3z4PPXxloH3yUWIWiKp08sQzBPFQ30L/BNMYs5klAdfGTq\nJgVnZfB1Ky+XPYOvsmhm2RKyaMIydosG+q+kGdKiUR18RB58MMvGi43rY9F0Gc6lWOkRmrpsrA+L\nJrXYjsV1MPiLAAAKSElEQVSiWV3N9oMJWS0yTxUN9F9JI4smY+EEvip7OPvsbGn9iRPNGXyZWcO5\nFCs9QlMVz3ktmjFX0QyRmTbFsyqWZlmcQmbxbatohs7gU7NoZlXRqEzSk6rsIV9S3lSD3WU4l6KN\nEJqqePpMsnYR+NQXOg2ZwbcReAi/2GmeSVYYJoNPrYpGC50CUGUpwPQdnsrMY9GkJkKh6TJhtqhV\nNENNHra1aCC8Dz+2SdahLBrd8CMyVcNDmN4NsUxXiyZFEQrNUBZNirEdi0UDwwu8LJppJPCBqBoe\nwt5NuKu8Llk09fRp0aQu8GOzaIYQ+JMns8eBA9PHyaKRwAehyaKpKl+C7HVZNNVUxfOss7J90Osm\nO2dZNKqDb0dTPJv69BACn8eoXHK8LBaN6uAj02TRNGU7bYdzzrW7jdlYqYrnykq2J0vdPuihqmhS\ni21VLB72sKwkMaZopGDR+LZH0/W3LBZN2zp4VdG0oMmiaSvw+bdulaDs7mZCV85UYLEy+KZ41g25\nl8mi6WMZ/pgsmi79JQayaDIWTuCbLJq2Am9W/43f1CiLJPBN8awTtab94GG8At8lFiFIWeDzhU75\noqqhYlRGAp8xU+DN7BIzu8XMbjezN1W8v9/Mrpoc82kzuzDOqfrR1aKpEiSo3xN+WQS+aqsCaB5y\nN+0HX7W9LUxfAM7Vr5wcKrZVW2DkxLYfqvY2hzTKJFdWssfubvY8JYumLmaqoplgZvuBK4FLgKcC\nLzGzZ5QOew3wN865HwN+G3hnjBP1JaRFA/Xf+PMI/Pr6ev2bCbG7my0Oe/jDT32vq0VTboM8FsWF\nILn9tVLRO4cS+KotMHJC2Q9V/WJ3N4tL3Yho6IVOMN0moSyaea4R5+oFPnYG31QHn+JCp4uB25xz\ndzvntoGrgBeWjvk54P2Tn68BftKsKkfrh5AWDSy3wB8/npW7VWXSsyyaLgKfZzgpjo7q+hWEsx+q\n+kUuVFVXVAoWDUy3ySyLxnd/nHmukc3NvR1iy8iimeYC4M7C87smr1Ue45zbBb4HnBfqBNvSZCl0\ntWjqBL5u5ntRNhuriyU018LPu/Iyxdh2jcW81MUS+rdofNqkLk779mWP++8Pd051zIpZigIfq4pm\n1vdG0LsXXnppyE+r5utfr+5gj3xkc5Zex4ED8KpXnWpT1G1cBllHvuYa+OY3q9+/4w74whfq/2Yq\nHD9eL2qPeARceSV8/OOnvnfrrVn5YJkDB+pjvX8/fOQj8I1vVN9nNGffPvjYx/rpS0WaBP4Rj4Df\n+Z3s/Oehql9sblbHEmbH881vht/93fnOKeeLX4RnP7v6vX374PLLs39vvx1e85rq4x7xCPiFX6gX\n3yLzXCOzYvbhD2c6EYP7768uk9y/H264obrf3norvP71cc7HXMOYycz+PvAm59yLJs/fCOxzzv37\nwjH/a3LM581sBbgHeMwkmy9+VsRb3QohxOLinOtke8/K4G8ELjKz84F7gcuBV5aOuRb4Z8DngZ8H\nPlsW93lOUAghRDcaBd45d8LMXg1cT+bXv985d5OZHQE+75z7KPAu4P1mdgtwDPil2CcthBBiNo0W\njRBCiPESfSXrrIVSi4yZHTKzT07+/3eY2b+avP5IM/u4mX3ZzK43s5rpu8XDzFbN7GYz++jk+RPM\n7LOTGP2xmdVMry4WZnbQzP6bmX3JzL5iZs9a1n5hZkfM7Gtm9lUz+5CZHViWfmFm7zOzeyYOSP5a\nbT8ws3ea2W1mdlPFmqRTiCrwngulFpmTwL+cLAL7ceAVZvY04Ajw5865pwLXTZ4vC68DbmevQuud\nwG9NYvT/yBbOLQO/D1ztnHsa8BSymCxdvzCzJwEvBS5yzv0IsAP8IsvTL/6ATB+LVPYDM7sMeJxz\n7inAr0x+t5HYGbzPQqmFxTl3j3Pu1snPx4EvA+czvTjsj1iSmJjZBWT/9/dkT20VeJZz7sOTQ5Yi\nFmb2KODpzrkPQrZ+xDl3H8vZL74PbAFnmtlpwAHg2yxJv3DO/RXwg9LLdf3ghfnrzrmbgdMm11Qt\nsQXeZ6HUUmBmjwd+AvgU8Gjn3PcAnHPfZcCFYT3zDuCNQF5ldR7w3cL7d7Mc/eOHge+Y2Z+Y2a1m\n9l/M7GyWsF84574P/EcyUf+/wFHgVpazX+TU9YPzaamnsQVeM7iAmZ0FfAh43SRTWzrM7EXAvZPM\nIy+ZXdbS2RWyL/vfds5dRJbF/rthT2kYzOyJwK8Cjwf+FnAW8DNDnlPilK+ZRo2NLfB3AYcKzw8x\n/Q208Ewmh/4U+K+FIed3zOzcyfuPJltjsOj8JPBiM/sm8EHgp4DfAs4tHHMBWZ9ZdO4E7nbO3Th5\n/iHg6cC9S9gv/h7wGefc9yY27tXAc1nOfpFTpw9lPZ0Zl9gC/9BCqYnQXU42abAUTDZdey9wu3Pu\nHYW38sVhTP69tu9z6xvn3L91zh1yzj0B+KfAJ5xzLwU+Z2b/cHLYssTiTuC7ZvbkyUsvAL5Cdm0s\nVb8A/jfwLDN72OR6eQHwVZawXxSo04drgV8GMLNnAjvOububPih6HbyZ/QOybYTzhVK/GfUPJoSZ\nPQf4JNnkah7ofwPcQDbh/BiyCoHLnXM93u9mWMzsecAbnHMvNrMnAB8gG5rfBrzUORdwm6w0mVRT\nvYdsUvFbZBeusYT9wswOk/3/d4GbgX8B/BBL0C/M7IPA88hGLPcAbwU+Qk0/MLN3Ac8HNoFXOOdu\navx8LXQSQojFZOFu2SeEECJDAi+EEAuKBF4IIRYUCbwQQiwoEnghhFhQJPBCCLGgSODFaDCzxxe3\nVS29d8TMfmry86+aWeVdOc3s983sR1v+3ePtz1aI4VEdvBgNkw3bPjrZQrbpuG8CfzffsCnA3z3m\nnDs7xGcJ0SfK4MXYWDWzd092YVw3szMBzOwPzewyM3st2aZVf2HZDeGnmPzOMyc/Hzez35jcgORm\nM/uhyetPnjz/gpn9Run33zq5EcNXzOw3J6/9IzP7n5Off8iym7ss/E6QIn0k8GJs/DDwrskujPcA\n/3jyugOcc+73yLadXXPO/XTF7xeHrAeATznnngF8jL0byv8e8Dbn3I+TbWMLgJm9GDh/ciOGpwA/\nZmYvcM79GfA3ZvYa4D8Db3XOLcNGYSJxJPBibHwzv4kK8AWmd9dry0nn3P+o+KxnA38y+fmPC8f/\nLPCzZnbz5Pgnk21zC/Basn2GTjjnrprjnIQIxmlDn4AQLdks/LzDfElKcfOq3dJn1U1O/bpz7n0V\nrx+anM9jzMycJrdEAiiDF4vIg8CZc/z+Z8i2tgb4J4XXrwdeZmZnAJjZY8zs3Mmt5t5Ltg3yV4Ff\nm+NvCxEMZfBibJQz46pM+b1kk6zfqvHhq37XFZ6/FrjKzN4I/Hn+unPuo2b2d4CbzOwkcAL4JeDV\nwCedc58xsy8DN5rZf3fO3dH2PydESFQmKYQQC4osGiGEWFAk8EIIsaBI4IUQYkGRwAshxIIigRdC\niAVFAi+EEAuKBF4IIRYUCbwQQiwo/x+/npQTl0T9/AAAAABJRU5ErkJggg==\n",
       "text": [
        "<matplotlib.figure.Figure at 0x7f10dbb882d0>"
       ]
      }
     ],
     "prompt_number": 11
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "plt.hist([val for val in trainer.network.predict(train[0])[4] if val > 0.1 ], bins = 50)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 12,
       "text": [
        "(array([ 1.,  0.,  1.,  0.,  0.,  0.,  0.,  1.,  0.,  0.,  1.,  0.,  0.,\n",
        "         0.,  1.,  0.,  0.,  0.,  0.,  0.,  0.,  1.,  0.,  0.,  1.,  0.,\n",
        "         1.,  0.,  1.,  0.,  0.,  0.,  1.,  0.,  0.,  0.,  0.,  1.,  0.,\n",
        "         1.,  0.,  0.,  0.,  1.,  0.,  0.,  0.,  0.,  0.,  3.]),\n",
        " array([ 0.11243238,  0.13007094,  0.14770949,  0.16534804,  0.1829866 ,\n",
        "         0.20062515,  0.2182637 ,  0.23590226,  0.25354081,  0.27117936,\n",
        "         0.28881792,  0.30645647,  0.32409502,  0.34173358,  0.35937213,\n",
        "         0.37701068,  0.39464924,  0.41228779,  0.42992634,  0.4475649 ,\n",
        "         0.46520345,  0.482842  ,  0.50048056,  0.51811911,  0.53575766,\n",
        "         0.55339622,  0.57103477,  0.58867332,  0.60631188,  0.62395043,\n",
        "         0.64158898,  0.65922754,  0.67686609,  0.69450464,  0.7121432 ,\n",
        "         0.72978175,  0.7474203 ,  0.76505886,  0.78269741,  0.80033596,\n",
        "         0.81797452,  0.83561307,  0.85325162,  0.87089018,  0.88852873,\n",
        "         0.90616728,  0.92380584,  0.94144439,  0.95908294,  0.9767215 ,\n",
        "         0.99436005]),\n",
        " <a list of 50 Patch objects>)"
       ]
      },
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAXYAAAEACAYAAACnJV25AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEOtJREFUeJzt3X2MZXV9x/H3BxdBQGkVtZHlwbapWkp1MURF1CFpFVEw\nirEVNT5sI9Wgpg8+1mSHNG1NTbFZSRVT8GFVQtWIUKDYWJcHpQbcFWaXSGOiZqG6oKUCjbAg3/4x\nZ2H27sx9GO6dufz2/Upucs49v3vOZ87MfObM78zdTVUhSWrHfqsdQJI0Xha7JDXGYpekxljsktQY\ni12SGmOxS1Jj+hZ7kgOTXJ9ka5L/SvKxRcYckOSiJHNJvpXkqMnFlSQN0rfYq+pe4MVVtQ74XeAF\nSU7qGXYW8JOqOhb4KLBxIkklSUMZOBVTVb/sFh8LPAbY2TPkFGBTt3wJcEKSjC2hJGkkA4s9yX5J\nvsd8oX+zqm7uGbIW2AFQVQ8CPweeMu6gkqThDHPF/mBVPYf5An9xkpmJp5IkLduaYQdW1S+SXAY8\nH9i8YNOtwJHA7Un2A54E3NH7+iT+ozSStAxVNdL09qC/inlSksd3y48D/hCY6xl2OfDGbvlVwHXd\nlMxi4abqsWHDhlXPYKa2cpnJTKM8umZc5EHPmNEMumJ/GvC57mbogcAXq+qyJGcDN1TVpcC5wKYk\nc8DdwBnLSiJJGou+xV5Vc8C6RZ7fsGD5PuB1448mSVqOffqdpzMzM6sdYS9mGt405jLTcMw0WVnu\nHM7IB0pqpY4lSY8G87Pci/ViHppfT0KN8+apJOnRx2KXpMZY7JLUGItdkhpjsUtSYyx2SWqMxS5J\njbHYJakxFrskNcZil6TGWOyS1BiLXZIaY7FLUmMsdklqjMUuSY2x2CWpMRa7JDXGYpekxljsktQY\ni12SGmOxS1JjLHZJaozFLkmNsdglqTEWuyQ1xmKXpMb0LfYkRyS5OslckluSvG+RMTNJfpFka/f4\n8OTiSpIGWTNg+y7gnVW1LckhwJYkV1bVjT3jrqqq0yYTUZI0ir5X7FW1s6q2dcv3ADcBT1tkaCaQ\nTZK0DEPPsSc5GjgeuLZnUwEv6KZrvpHk2eOLJ0ka1aCpGAC6aZgvAe+pqrt7Nn8XWFtV9yZ5KXAx\n8PTxxpQkDWtgsSfZH/gK8MWqurh3ezdFs3v560l2JfmNqvpp79jZ2dmHlmdmZpiZmVlmbElq1WZg\nz74cVapq6Y1JgM8CP6+qP1tizGFV9bNu+bnA14Ajq+rBnnHV71iStK+Zr9jFejHs7sskVNVI9zEH\nXbG/EHgjcFOSrd1zHwKOBKiq84DXJ3l7t20XcEZvqUuSVk7fK/axHsgrdknaw6Su2H3nqSQ1xmKX\npMZY7JLUGItdkhpjsUtSYyx2SWqMxS5JjbHYJakxFrskNcZil6TGWOyS1BiLXZIaY7FLUmMsdklq\njMUuSY2x2CWpMRa7JDXGYpekxljsktQYi12SGmOxS1JjLHZJaozFLkmNsdglqTEWuyQ1xmKXpMZY\n7JLUGItdkhpjsUtSY/oWe5IjklydZC7JLUnet8S4jUm2J9mSZN1kokqShrFmwPZdwDuraluSQ4At\nSa6sqht3D0hyOnBkVR3TlfqngedMLrIkqZ++V+xVtbOqtnXL9wA3AU/rGXYKsKkbsxVYk2TtBLJK\nkoYw9Bx7kqOB44FrezatBXYsWL+1e06StAoGTcUA0E3DfAl4T1XdvdiQnvVabD+zs7MPLc/MzDAz\nMzNUSEnad2wG9uzLUaVq0Q5+eECyP/CvwL9V1ccW2X4+cEVVfblb3wa8rKpu6xlXg44lSfuSJCx+\nHRx292USqqr34rmvQX8VE+B84ObFSr1zOfCGbvxxwK96S12StHL6XrEnORG4mvmbprsHfgg4EqCq\nzuvGnQucBNwH/ElVbVlkX16xS9ICk7piHzgVMy4WuyTtaVWmYiRJjz4WuyQ1xmKXpMZY7JLUGItd\nkhpjsUtSYyx2SWqMxS5JjbHYJakxFrskNcZil6TGWOyS1BiLXZIaY7FLUmMsdklqjMUuSY2x2CWp\nMRa7JDXGYpekxljsktQYi12SGmOxS1JjLHZJaozFLkmNsdglqTEWuyQ1xmKXpMZY7JLUmIHFnuSC\nJDuTzC2xfSbJL5Js7R4fHn9MSdKw1gwx5tPAx4HP9RlzVVWdNp5IkqRHYuAVe1VdA9w5YFjGE0eS\n9EiNY469gBckmUvyjSTPHsM+JUnLNMxUzCDfBdZW1b1JXgpcDDx9sYGzs7MPLc/MzDAzMzOGw0tS\nSzYDe/blqFJVgwclRwOXVtWxQ4y9BXhJVf205/ka5liStK9Iwvykx15b2N2XSaiqkaa7H/FUTJLD\nFiw/FzgYuP2R7leStDwDp2KSXAi8BDgsyQ5gA7A/QFWdB7w+ydu74buAM6rqwQnllSQNMNRUzFgO\n5FSMJO1haqdiJEnTxWKXpMZY7JLUGItdkhpjsUtSYyx2SWqMxS5JjbHYJakxFrskNcZil6TGWOyS\n1BiLXZIaY7FLUmMsdklqjMUuSY2x2CWpMRa7JDXGYpekxljsktQYi12SGmOxS1JjLHZJaozFLkmN\nsdglqTEWuyQ1xmKXpMZY7JLUGItdkhozsNiTXJBkZ5K5PmM2JtmeZEuSdeONKEkaxTBX7J8GTl5q\nY5LTgSOr6hhgfTdekrRKBhZ7VV0D3NlnyCnApm7sVmBNkrXjiSdJGtWaMexjLbBjwfqt3XO3LvWC\n+++/nzvuuGPJHR566KEcfPDBY4gmSfuecRQ7QHrWa7FBs7OzAOzcuZPzzvsUj3vcU/cas2vX//LJ\nT36c9evX732Q9B6m56C16GHHatwZVupjmoZztxwreb6X2te4X9PvdaOahs/rcjLsS+dndJuBh/ty\nOTLMB5bkaODSqjp2kW3nA1dU1Ze79W3Ay6rqtp5xtftYc3NznHjiGdx11973Yw86aD0bN57Qp9iX\nypsV/CIeX4aV+pim4dwtx8qd76X3Nd7X9H/dqKbh87qcDPvS+elnmK+tJFRV/59QPcbx546XA2/o\nAhwH/Kq31CVJK2fgVEySC4GXAIcl2QFsAPYHqKrzquorSU5Ksh24D3jrJANLkvobWOxV9fohxpw1\nnjiSpEfKd55KUmMsdklqjMUuSY2x2CWpMRa7JDXGYpekxljsktQYi12SGmOxS1JjLHZJaozFLkmN\nsdglqTEWuyQ1xmKXpMZY7JLUGItdkhpjsUtSYyx2SWqMxS5JjbHYJakxFrskNcZil6TGWOyS1BiL\nXZIaY7FLUmMsdklqjMUuSY2x2CWpMQOLPcnJSeaS3Jzk/Ytsf0uSO5Js7R5vm0xUSdIw1vTbmOQA\n4BPAicBO4LokX6+qrQuGFXBhVb17cjElScMadMX+PGB7Vd1WVQ8AFwGv6BmT7iFJmgKDin0tsGPB\n+q3dcwsV8Jok25NckuSocQaUJI1mULHXEPu4BDiqqo4BvgZ84RGnkiQtW985duav0I9YsH4Ee17B\nU1V3Llg+P8k/LrWz2dlZAHbu3MkDD9wzalZJ2gdsBh7uy+VI1dIX5UkOBL4PvBC4Hfg2cGZVbVkw\n5slVdUe3fCpwdlUdt8i+avex5ubmOPHEM7jrrrm9jnnQQevZuPEE1q9fv1gelv4lIvT7WMZl3BlW\n6mOahnO3HCt3vpfe13hf0/91o5qGz+tyMuxL56efYb62klBVI93H7HvFXlX3JnkHcCXz0zabqmpL\nkrOBG6rqUuAvkpwCPAa4E3jTKAEkSeM1aCqGqroCuKLnuQ0Llj8AfGD80SRJy+E7TyWpMRa7JDXG\nYpekxljsktQYi12SGmOxS1JjLHZJaozFLkmNsdglqTEWuyQ1xmKXpMZY7JLUGItdkhpjsUtSYyx2\nSWqMxS5JjbHYJakxFrskNcZil6TGWOyS1BiLXZIaY7FLUmMsdklqjMUuSY2x2CWpMRa7JDXGYpek\nxljsktSYgcWe5OQkc0luTvL+RbYfkOSibsy3khw1maiSpGH0LfYkBwCfAE4Gfh94bZJ1PcPOAn5S\nVccCHwU2TiLoJGzevHm1Iyxi82oH2Mt0nqfpNI3nahozTaOWztOgK/bnAdur6raqegC4CHhFz5hT\ngE3d8iXACUky3piTMZ2fyM2rHWAv03meptM0nqtpzDSNWjpPg4p9LbBjwfqt3XOLjqmqB4GfA08Z\nV0BJ0mjWDNhekzrwL3/5I57whFP3en7Xru8BJ0zqsJLUvFQt3d1JXgS8v6pe2a2/F3hsVf3NgjHf\n6MbckGQ/YCfw1O7qfeG+JvZDQpJaVlUjTW8PumK/Hvi9JIcDtwOvA87sGXM58EbgBuBVwHW9pb6c\nYJKk5elb7FV1b5J3AFcyPx+/qaq2JDkbuKGqLgXOBTYlmQPuBs6YdGhJ0tL6TsVIkh59xv7O0yHe\n0PTiJFuS3J/k9HEff5mZ3ptke5JtSa5O8vQpyPTOJDcmuSnJDUmeu9qZFow7PcmDSY5b7UxJ3pLk\njiRbu8fbVjtTN+Z1XZ6bknxx0pmGyZXknAXn6ZYkd05Bpmcm+U73vXdzkldNQabfTHJtl+mb3VT0\nJPNckGRnN+ux1JiNXUdtWeS9RHurqrE9gAOAHwKHMz/Ncz2wrmfMUcCxwGeB08d5/EeQ6UXAAd3y\nnwJfnYJMhyxYPhW4arUzdeMeD1wNfBs4brUzAW8GNk7662jETM8GvgMc3K0/cRpy9Yw/C/jn1c4E\nfB44s1t+FrBjCjJdCrypWz4J+PKEM70IWAfMLbH9dODibnkd8L1B+xz3FfvANzRV1Y+rag7Y6wbr\nhAyT6Zqquq9b/Rbzn/TVznTPgtVDgJ+sdqbOXwMfAe4DJn1DfJhMWYEco2Z6K3BuVf0fQFX9z5Tk\nWugM4MIpyLQDOLRb/jXgx1OQ6RnAf3TLm4GXTfJNl1V1DdDvt6eH3gRaVVuBNUl630+0h3EX+zBv\naFppo2Y6E/jaRBMNmambjvkBcA7wwdXO1E29HF5Vl3dPTfoGzTDnqYDXdL+mXrIC/1bRMJmeATyn\nm0L7bpLTJpxp2FwAdOfoaB4ur9XM9HfAm5PsAC4D3jUFmeaYv0oGeDVwMKv7psuRe3XcxT6Nd2KH\nzpTkDcBxwN9PLg4wZKaq+qeq+m3gz4ELJhupf6buPQrnAH+58OmJJhruPF0CHFVVxzD/A/kLk400\nVKb9mC/O5zFfEJ9M8sRJhmK0770/Br5U3e/2EzTM/s9hfkroCOavTD8/2UhDZXo38NIk24GXAz8a\n8nWT1Pu91jfPuIv9VuCIBetHsOdPml4rcbKGypTkD4C/Ak6rqvunIdMCFwHHTzTR4EyPB44BNif5\nIfB84JIJ30AdeJ6q6s7uV2qq6nzm57cnaZjP3Q7g0qr6VVX9CLgZ+J0pyLXbHzH5aRgYLtOJwL8A\nVNV/AgcmmeTV8TBfU7dV1Su7i4V3AQdW1e0TzDRIb+a13XNLG/NNgAOZ/+l2OLA/8zcmFr3BBnyG\nlbl5OjAT8zckfgD81qTzjJDp6AXLpwLXr3amnvHf7Ld9Bc/Tk3vO05YpyPRq4DPd8mHAbQtzrubn\nD3gm8MNJZhnxXF0GvLlbfhbwU+Axq5zp13n4T8E/DHx0Bc7V0fS/efrVbvk44MaB+5tAwJcD25i/\nSvlg99zZwKnd8vHM/4S8B/jZUh/MCmV6Zbf878zfnNzaPS5exUy7z9PHgZu6MVcBz1rtTD1jJ17s\nQ56nj3TnaTtwLXDMamfq1v+hy/R9ur+wmJJcG4C/XYk8Q37+ngFc152rm3d/T65yptcCt3RfV58C\n9p9wnguB/wZ2dd34Nubv9Z25YMy53TnaMsz3nW9QkqTG+F/jSVJjLHZJaozFLkmNsdglqTEWuyQ1\nxmKXpMZY7JLUGItdkhrz/1n1F+f8VmbNAAAAAElFTkSuQmCC\n",
       "text": [
        "<matplotlib.figure.Figure at 0x7f3de003e150>"
       ]
      }
     ],
     "prompt_number": 12
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "plt.plot(trainer.network.predict(train[0])[4],'bo')\n",
      "plt.plot(ds_output[4],'y')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 9,
       "text": [
        "[<matplotlib.lines.Line2D at 0x7fa6d0d602d0>]"
       ]
      },
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAEACAYAAAC57G0KAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXuYHFWZ/79vd0/PhJAQruGSxMRRVpeIIa4bXCEMApng\nmFVXDb8VFW8o/jQTHnwQJWGTrER/eN1M5MH1xqIuyMriijYhg7hDB7kbbgkQZILkwh1yIZe+1vn9\ncbq6q6uruqu7zqk+VXU+z5Mn0901VTWnTr/11vd8z3uIMQaNRqPRRI9Et09Ao9FoNHLQAV6j0Wgi\nig7wGo1GE1F0gNdoNJqIogO8RqPRRBQd4DUajSaitAzwRPQzInqRiB5rss0IEW0moo1EdIrYU9Ro\nNBpNJ3jJ4K8FsNDtQyL6EIAZjLGTAHymsr1Go9FoukzLAM8Y2wBgV5NN3gvgF5VtHwKQIqJpYk5P\no9FoNJ0iQoOfBmC75fWOynsajUaj6SKiBlnJ9lrXP9BoNJoukxKwjx0ApgO4r/J6WuW9OohIB32N\nRqPpAMaYPYn2hIgM/lYA5wMAEc0FUGaM7XTakDEWi3/btn0PTzzxadfPV6xY0fVzVOVfu22xffta\nbNr0ka6ftwptocq/O++ciFxup26LDv89+eSFePbZb7t+7oeWGTwR3QDgDABHEdF2ACsA9FQC9r8z\nxv6biM4kos0A8gA+5euMIgBjJTBW6vZpRBLDyOm2VQzd3/0hs/1aBnjG2D972OZLYk4nGugOLw8d\n4NVD93d/yGw/PZNVAq0u2MDAQHAnozjttgVj+cgGkzD2Cy4hlIVfkzC2RafoAB8ydID3TrttEeUM\nPoz9grFy5X8d4DtFB/iQoR9Z5RHlAB9GzGuhr0nn6AAfMhgrVzMbjVgMI6/bVinMDF5fk07hbSen\n/XSAl4DO4OWhM3i10Bm8f3QGHzJ0gJeHDvBqoQO8f3SADxGZTBa33HIXHn74aQwOLkcmk+32KUUK\nLtHoYKIKOsD7Rwf4kJDJZLF06Xrs3Pl3eP31aRgdvRJLl67XQV4gOoNXCx3g/aMDfEgYGRnF+Phq\nJJMlJJP8go2Pr8batbd3+cyigw7waqEDvH90gA8J+TyfGJxIlJFI1EbFc7lkt04pcmgXjVrUfPD6\nmnSKTNedDvAC6e3ld2FrBg8AfX2684tCZ/BqoTN4/+gMPiQMDy9Af/+yugDf3385liw5p8tnFh10\ngFcLHeD909ViYxrvDA3NBwBs334RJk16GYODV2DJkoXV9zX+iXItmjCiA7x/dIAPEUND87Fp099g\n//4ibrvt690+ncjBM3i9dowq6ADvHx3gQ4ae6CQPw8gB0IPWqqADvH+0Bh8ydC0aeRhGHrLqdmja\np9bP9TXpHPHllk10Bi8BncHLgTEGw8iBSGfwqqAzeP/weKFtkqFBB3g5MFYEwHTbKsQ99/wZAPCd\n7/xGl+boEC3RhAwd4OVgGHkQpXXbKkImk8U119wPANi69VxdmqNDdIAPGTrAy8EwckgmD620r3bS\ndJuRkVG8+OLHAUCX5vCBDvAhQwd4OfAAfwgAAmB0+3RiTz6fqgb2RKLW33VpjvbQAT5kcBeNDvCi\nMYwcEok+EKW0S0kBentL1ZpLujRH58iMFzrAS0DmqHicGRv7E7Zt24NCgWHRoiu01ttlhocX4IQT\nfg4ASCZ5f9elOdpHu2hCBr8bG2BMywiiyGSy+O53s3jttRkoFvswNrZMD+h1maGh+fjUp94GAOjv\nvx2Dg1dgzRpdmqMdMpksDKOIbdtekOJC0gFeAjVvsM7iRTEyMoqXXvo0CoU+lMtc+9UDet1nzpw3\nAwA++cnTcdttX9fBvQ0ymSwuvvg2EAH5/BQpLiQd4CWgJ3+IJ59PIZ3Oo1jsrQZ4QA/odRvd1ztn\nZGQUf/3rSgDyXEg6wEtAd3rx9PaWkE7n6jJ4QA/odRvd1zvH6kKyDlKLTFp0gJdAbZUb3elFMTy8\nANOn/6wa4BOJsh7QUwLd1zslCBeSrkUjAa3Bi2doaD6SydsxPv4bpNOv48wzr8InPrFYa75dRvf1\nzhkeXoDnn18BANVAz5OWhcKOoQO8BBgrVbzaOqsRydvfPgszZ/49du/O4dprl+KQQ97c7VOKPbqv\ndw5PTnajXF6D3t59UhYI0gFeArzT9+pOL5j6iU66bVVA93V/nH32O3HffROQTpOUBYK0Bi8BxkpI\nJHSnF41h5JFI9OoArxC6r/tDdvvpAC8BftH6dKcXjM7g1UP3dX/Ibj8d4KVQ1p1eAroWjXowpvu6\nH3j79UqrkKoDvGAYY9XHLr2MmVh4PfheECV1QFGEWgaq+3on8DGMNHgoFl/aRAd44RgAEnphCglo\niUY9tETjD9OFJOuptGWAJ6KFRPQYET1ORJc5fP4WIrqPiDZVtnm/8LMMEfyCJXWWKQEd4NVDD7L6\nQ3a8aBrgiagXwDUAFgI4GcCHiegU22bLAfyMMTYbwIcA/ED4WYaI+juy7vQiYUy7aFRDB3h/yI4X\nrTL4eQA2M8Z2Mn70GwEM2bbZDuCwys9TADwr9hTDhQ7w8tAZvHpoicYfsuNFq4lO08ADuMkOAAO2\nbb4J4B4iWgJgIoCzhJ1dCGGsrAO8JKwBXg9gq4HpoimVXu/2qYQUufGiVYD34tv5HoCfMMa+T0Sn\nAvglgJOcNly5cmX154GBAQwMDHg7yxChM3h56IlO6lHL4Hd1+1RCiVO8GBsbw9jYmJD9twrwOwBM\nt7yejvqMHgBOA7CCnyy7l4j6iOgYxthL9p1ZA3xUkT0qHmdqGbwewFYFM8DrJ6rOcIoX9uR31apV\nHe+/lQb/AIDZRHQCEfUAWAxgnW2bcQBnAwARvRVcpnm14zMKOdpFIw+twauHrkXjj1q86MIgK2Ms\nB+ALANYDeATAzYyxjUS0iogWVTa7BMBFRLQZwH8D+CyLceqqJRp58ACvJRqV0C4af9TihZyEsGU1\nScbYOtiydsbYCsvPWwC8S/iZhRQd4OXBNXidwauEdtH4o9s2SU2b6AAvDy3RqIcO8P7QAT5kaJuk\nPAwjV6lFowewVUEXG/OH7HihA7xgtItGHjWJRg9gq4IuNuaPrtei0bSHlmjkoSUa9dASjT+0RBMy\n+EXSNkkZaBeNemgXjT9MmyTQhWJjUcYwilL2G9cMXlZ71vZfAsBi0baMGaGRPMIQ4A2jKGUxDRHo\nDF4SGzfOQy63Tfh+4xjgc7lt2LhxntRj8EqSfSCiyLftc8/9O555ZkXrDRUgDBLNpk0fwOuv39/t\n03BEB3hJFIuvoFTaI3y/cXPRZDJZXHjh17Fjx1YMDi5HJpOVchxTngEQ+QHsUuk1lEq7u30aHlHf\nRVMq7UKx+Fq3T8MR7aKRhGHkwJh4WSFOLppMJoulS9fjz3/+LIrFiRgdvRJLl66XEuRNBw0n2uMb\nsvqmDMLgomGsCMPIdfs0HNEZvCSCCfDRDUIAMDIyivHx1Uinc0ileFuOj6/G2rW3Cz+W6aABEPm2\nDWeAV/d6GEY4AryMgm0xDvB5KQODskfFVSKf55Uuenry1QAPALlcUvixGiWa6LatYeRDFeBVLzbG\nM/h8t0/DEdmuu1gGeO5SKOgM3ie9vfzvs2bwANDXJz4TsUo0UW9bw8hJdyWJIgwumjhLNC2LjUUR\n824uO8Cr2qlEMTy8AOPjy5BOvx3JJG/L/v7LsWTJQmHHyGSyGBkZxdFHP49zz30ZL76YxcknRz/A\nhymDTyTSAAwwZoBIvZxRB/iYITPAy16CSyWGhuYDAP74x++ip6eAwcErsGTJwur7fjEHccfHV2PO\nnP/F3/3dOFasWI+rr96PWbPUHdTzS7gkmjKIeqqmAhUDPPfBqyrRlKXOzo5pgM9V/tcSjV+Ghubj\nlFO24KmnbsG6dSsr4w9iMAdxAa7zFwp9GB9fjbGx92LmzBnCjqMaYcvg6/t7T7dPqYE4Z/Dq3W4D\nwLzYjBWE71v2qLiK1G6YYtvTHMQFuM5fKPRV3o/2ADbX4MX3TRmEwRbMWCEkAV67aIRgPq7J0+Dj\nVYtGluRlDuICPMAXi9xFk0wi0m0bLolGfdcYt0mqKtHIjRexDPBaohFL7YlIbHsODy9Af/8yAEA6\nzSWa/v7LcfbZJ0e6bcMt0ahHnCWaWGvw2iYpBlk3THOwdu3aK3DSSQ/hmGOew5o1n8GcOTvxyivj\nQo+lEnxhk3B8NcPQ33WAjxlybZLxcdGYyJS8hobmY2hoPnbsWIODB7fizW+ej5de+jWiPL7B21PN\n6od2VO/vXNdmyko0sl13MQ3wOoMXicz2rB0jH6OZrDkkEtTt0/CE6v3dfKrUGXyMCE6Dj26WaSWY\nAG+tRaPugJ4IwirRqPhUZfbJuAb4mA6yaheNSMz2lDm9Pl7FxuTUSRINYwYAA0BC2f5eC/BqSjS1\neKFtksIIygevYoeXgcz2tB4jThJNGFw0PCAllV6ExZxPoH4Gr22SwtA2SbEEp8FHP4NnjFUCvPoT\nncy+Dqh7TbREE0O0i0Ys3ZFo1NN7RcD7DAuJRFMOTYBXuRaNDvCC4Xdz0i4aQchsT+sxiKIv0QTR\nlqIITwZPOoOPE4aRQzJ5aAABPppZph2Z7Vk7Rj4WLpog2lIUjQFevf5uGEUkk4fqAB8nGMsjmZyk\nNXhBGAZvz+BsktFt2z/84U7s3s3AWBGDg8ukLWIugrBk8Py7rqpEo4uNCUd2Bi9zCS4VMdtTvgYf\nbYkmk8niyivvwN69U1EuJ/CHP6yStoi5CGqFxtR9quIBXvUMXl6xthgHeDkZZzwzeHntacJY9F00\nIyOjeP75L6JQ6EO53INUqihtEXMRhCuDVznAa4lGKIaRQyo1SUrN7Xi6aHh7yvfBR9tFk8+nqmWR\nS6We6jq3MhYxF4P6LhrDKCCZnAjGipWJWWqhXTQS4JpxEIOs6nV4GfAxDS3R+KW3t1Qti8wDPL9h\nyljEXARhyODvu28jHn54O0qlJN73vq8pJ3fpWjQSMCWFcnlf3fvmAs/5fAq9vSUMDy9oe33ROAb4\nICQaq4tG5cUl/DA8vADXXHO1JcAXhS9iLhLVA3wmk8VPfnI/5s17C44++jVks1/Fli3fAQBh6wb7\nRQd4CRhGDun0cSiVdlXfsy7wbDI+zhebaKczxM0myWdeaheNCIaG5iOZvAc7dvwEyWQOZ531TVxw\nwYeVCUZ2VLdJjoyMYv/+f0a5fC0KhT6k07nKmMYVyrSp7GJtLSUaIlpIRI8R0eNEdJnLNouJ6CEi\nepSIrhd+loIxA5JVUrAu8GzSyQBXGJYwEwljRRAlkUj0BRDgoy3RAMA73vFmnHrq2zB16pG47rql\nygQiJ1R30eTzKaRSRRSLaRQKvejp4VZJlcY0ZLvummbwxKcOXgPgNAAvAriHiEYZYw9Ztnk7gC8D\nOI0xtp+IjhB+loKpDQrWApJ1gWcr7XaGuEk0ZmZN1CNZg4++iwYIrj1FoLpE09tbQrlcRLncU83g\nAbXGNLrtopkHYDNjbCfjR78RwJBtm08B+AFjbD8/Yfaa8LMUjJNmbF3g2Ur7nSFeLhozICUSPQFL\nNOp8SUViDfCqz2ZVvRbN8PACHH/8DSiVelAs8gDPxzTO6fapVem2i2YagO2W1zsq71n5GwBziOhB\nIvozEf2jyBOUgZNmbF3g2aSTzhC/DD4Pol6pAYkxA4wVQJQGoGYwEYW5cpXsG6YIVM/gh4bm4/zz\n34pjj30UfX2v4PTTR7BmzUKlZK9uD7J6WRgyAWAmeLY/HcDdRHSXypl8beZlzbdtXeA5l0uir6+M\nJUva7wzxC/BWSUHOdHDDKFRuInwZuyi3bbgyeLUDPADMnj0Ts2adiv37N2No6AJMmRJMcPfqyOt2\ngN8BHrRNpqM+o0fl9V2MPzP/lYgeB3AigHvtO1u5cmX154GBAQwMDLR/xgJws/WZCzz7QfUlzERT\nH5D2tf4FH8cwUXFATxTh1uDV6+/cBNCDRKIvsNms7TjynAL82NgYxsbGhJxLqwD/AIDZRHQCgJcA\nLAbweds2GQDvB/AfRHQUgLcCGHfamTXAdxPG8g2DrOL2Ha8l+4KQFKwOGgAV50YZjLFqVh8V6iUv\ntRf9sLpoVHWNGYYZ4HsDKzjm7shrtGc6LdlnT35XrVrV8bk01eAZYzkAXwCwHsAjAG5mjG0kolVE\ntKiyzW8AvEpEmwHcBeCrjLGXOz6jAJBZHCvOEo08DT5vy+ATABJKZox+CWrQWgRhkGi6kcG348iT\nvWRfy4lOjLF1ANbZ3lthe/1lcKtkKNDFxsQRhKRgl2gAa0CJ1ly9WnumQyjRqNffGSsikQg2wLfj\nyOu2TTJyyJ55WbM98TsyY17GqcMLz657kUikA5NoAER2jMNsz3AMsqptkwTMDD4NouAkmvYceXJt\nktFKfzwge+Zl7Y7MZQTAAKDOzDnRBCHRNM/go4WWaMRiGEWkUhMDzeC9OvIYYzYNXgd438iWFJyc\nBbWBqOgRTIDPOwR4NQf1/BJum6R6qyZ1Q6IBvDryDAAJECW0RCMK06XAMyQZ9eDVrs8hGtNFw2+Y\nclwfbhJNFNu2vj3VD/Dm06mqfZ1PkAvWReOV+lihl+wTgj1DEq2Rh+GxVSRBSApxlGjCmcGrdz1q\nNslgM3gv1LefXrJPCLUvkByrXRg6vUi6J9FEs23rb5hh8MGr3de7YZP0ShDtF8MAn7eUnRUflMLg\nLBCJdWKOTJskL2xaQ9WZk34Jl0Sjfl+vafC9yo0RBNF+MQzwtcd9GbJCGLIakWgXjVisPngt0fhH\nZ/Axo77srPgvURg6vUhqkoJMH3z8XDTaJikGrsGnYxvgY2eTtE57l/EYHIYCTCIxg6/8DD5OLpqw\nDrKq19dNiQboU9RFI7f9YhfgrcFClkSjunVMJGZ7dq9UQbQwxxtkll8WRRgswTWJBopm8HKLtcVc\notEavF+CcH3Yi40B0W1bbZMUC19LIAwafK1CqkhiH+BFTs7hF8ewTV5Qr9OLJLhB1ri4aMKjwZt1\nVAB1+7rag6zW9iPwLF5sn45hgJdnk+QXJ4k4rDxkYi2OpSUa/4Sr2Jj6GbypwQdZbMwr1vYD5LRh\nTDV4OTbJIC6YauiJTuIwjBIYM8CL1ckr/SCKsAR4oh7w8suqZfDy40UMM3h5GrzzBYuejGClW6UK\nVF1ByA/mWAMRSbWdiiIMLpqwlCoAdIAXgjnzEhBvk6wfFQeiGITsWGdeapukP2TPshZNuFw0ako0\n1lLi5kCrSGIu0TRmSYwxHDz4FxxyyIlt7zveEo28FYjiI9HYDQBi2zOTyWJkZBSHHrob+/ZNwfDw\nAl+LzIdFokkkujPRKZ9/DsnkJKRSk1zOTWvwwjGMHFKpwwE4Z0kHDmzBpk3vx7x5WzrYezm2AV62\nRBOHWjQy5cNMJoulS9fj2WdXYN26iTjnnALGx5cDQMdBPiy1aLol0TzzzHIcdti7cdxxn3E5N/nx\nIpYSTbPHYMM4AMM40NG+45nB14qNySi/zI8RDxeNtW+KvmGOjIxifHw10ukcUqkSenoKGB9fjbVr\nb+94n2HI4GsafPDFxsrlAyiX3WOJ1uAlYBg5PP74dgwOLsfdd49j5cprkclkLZ/nO9bq4hngreWX\nSUpWHVeJRmSAz+d5v+zpydf9n8t1vtpYGAI8Y4VqNUnDyElJQNyP3TyWaIlGAtu2PYvrr09idPT3\neNe7tmDTpo9g3br1APijqmHkOn6Uc7pgUVwY2opzUBLbrZwzeDUH9fwgU4Pv7eVtlU7nqv8fODAZ\nfX2d988wuGhqNslkpX8WGuQ+WbSKJTqDl8Bf/rINO3d+FABQLvcglSrWPar6D/D1o+JRC0J2ZMoK\ntWPEx0VjdXiJbMvh4QXo719WF+D7+y/HkiXndLzPIGqp+MWUaAAE7qTxFuCt8UL8TTJ2GTyRgUKB\nZ0mlUg+SSf4lMh9VDSNf0ZKNiuzgnThLNIA8a19catE0OrzETXQyB1JvuOE7AID3vOcqfOxj58XC\nRVML8OZA6+RAjs0loXYkGvE3ydgF+HS6UBfgUykekMxHVfOOaxh5JJMT2tp3EKPiKmEYJQDM8iWX\nU67AbZA1avKXbJvk0NB8nH76Idi48Uf48Y8vwqGHvs3X/sLiouHlghG4k4aP5zU7nnbRCGfmzMNw\n1FG/BFAL8NZH1VqAb78jxC2Dr9VNMWvvaInGD2Z7AnLb0vq/H1TP4E25ozYZS0WJRg+yCuXww/vw\nxS++G6XSFZg6dSMmTtyKCy9cVn1UNR+pdIBvjT2zljW9Pj42SbnLSZrHsP7vB9UDvLmak0nwGbwO\n8IFjGDmcdtqpOPfcd+Dppw8inT4WM2bMr/uc/9/+nT5utWjsgVdOfX1WN5BbO5aag3p+CGo8gx/L\nfyZrD/CqSWZWeQbolkTTXZtk7CQaq6faaZEK/xJNfFw0dn+6jAqIfJAsaavxo2bG6Bdre8ocz7D+\n7wfVVy/jlkh7gFdNorG7aHSA94V12rvzTFatwXvFro3LkBWcK0lGs22t7ak1eP9YLZKAaZNUWaLR\nC374ppVTwbzDdzKtOX4uGvkSjdMsVn6s6MlfwWjwnfdvO6q7aKwWSUBFiUa7aIRjfwzWGXznOEs0\n4jN4p5mHUWzb+r6ZFi538WPEJ4N30uCDqkfDLcTlrg+yxjDA1x6DnbIkHeC9Yw++cjJ4LdGIPob1\nfz+EIcB3S6Lx0s46wEuglVPBfKQS56JRq9OLpNEmGaREo96gnl+ClGhkuGhUux6NGnxwEo0Xt5K2\nSQqmceZl4yIVfjN4+wotUdOJgdrCEbNmPY25c5/Bjh1ZDA3NB5F4H7zTJCdAzYDiFyf5kDFWnUgm\n5hiiM3iri0atvm4u9mESpIvGewYv10UTqwBvn3npJtF0eqd3uiN3WlteVcyFI8bHV+PMM3+FqVMZ\nvvUtXo1zxgw5GnycJJqawysBINEwECfiGKIy2foMPlF5r/0aTrLotkTTqp3t8YIXbAvYRUNEC4no\nMSJ6nIgua7Ldh4jIIKK5Qs9QIM6uD7sPPo9UakqHEk30XTTmwhEAr0hYKPRVq3E6zSvwi1OhMSD6\nLhpAniup0/7diNr93TCcfPBBBXjezs2LjXXZRUM8nbgGwEIAJwP4MBGd4rDdJABLAdwr9OwE4/QF\ncpJoksnD9CCrC+bCEUAtwAO8Gqe8Qda4SDRBjGl03r/tqN7fu2mT9NLOKgyyzgOwmTG2k/Ej3whg\nyGG7rwP4fwDyAMQJhoKxT3l3s0mmUjrAu2EuHAEA6XQexSJvz76+sjSbZFwkGmuxMUCeK6nT/m1H\n9f5ut0kGWWys1s55uK0ipUKAnwZgu+X1jsp7lpOiuQBOYIzdWnkruDWx2sRLhsSYH4lG7Q4vAnPh\nCADo6eEZvFmNM/iJTtFqWy9PmH7x078b96V2f+9uBp9HMnlIpU2cr6EKLpqmwZr4aMr3AFxgfdvv\nScnCi8bJ77xTBGbw0dKJzaqba9degRNP3IB0uow1a1ZjaGg+nnrqhsAkGlVXEPKDc3VOsWMafvq3\nHdX7ezdtkua1NI9pdfOYqBDgdwCYbnk9HfUZ/SQAJwEYqzhTjgVwCxEtYoxttO9s5cqV1Z8HBgYw\nMDDQ0Ul3inVJNMBdg/cn0ViLYkUvCAE8yA8NzcfWrZcjmZyIN7yBB3090ckfXiRE/8fIobf3aBQK\nL/nel+rF9ZxdNMFJNES9FudO4ypSjbZqXpFzbGwMY2NjQs6jVYB/AMBsIjoBwEsAFgP4fO0E2R4A\nR9dOkP4XwJedgjtQH+C7gbcMPo9U6jCUywfb3r/qj6yiMYwcenqOrL6Wo8FriUbsMXj/zue3t964\nBar3926WC67P4J1vKm5L9tmT31WrVnV8Hk01eMZYDsAXAKwH8AiAmxljG4loFREt6vioXcLLAhXm\nI2xnNSvUto2JJogFP5pn8OrIASIIxiYpUqJRu7/zDL47C36YiUmzYwbRfi1nUDDG1gFYZ3tvhcu2\nZwo6Lyk4PQLbCzppm6R3nCQvGZpxvIqNyS+/HBebZKMPvrfDxK2TY+cqkyrdJ1ep4KKJFO1INDrA\ntyYoSSHOEo2ciU6H+daiGTMA1Gaw8p/Vuibhk2h0gPdFK5skYwYYKyKVmizMJqnaMmYiCUpSiGOx\nMUBe+WUREk3jNHv1ZLNu2iTNGdjNJRod4IVin/ZuD0hcckj7rEWjrqtANPbsWt6KTtGXaBhjDZqx\nyhOdGh1jgGquMWebZLASTbP6N3rJPsHYg4U9Q/IyMNKMeEo0cm19zWvRRKdtzeTCWjlSjg9elETj\nlsGrc01UKDbWbJERvWSfYFpJCl50s2ao7ioQTVCSgluAj5L85fR3qi3RNFa5VK2/d1eD95Is6iX7\nhNJKUvAy8t2MuGXwrSQvEcRHoml8UpHVniJcNOHM4FWUaHSAF0YrScHLwEgz4hbg7e0Z7IpO0Wpb\npxuZ6PY0F7xJJifGQqKxa/CdJm6dHVu7aAKnUVKoX9HJv0QTxwBvb0/xPvg4uGjcJBqRAb6WwPRU\njtl5+4UhwDNWQOOKTupINDrAC8ZZoilYPm/9WNUMZxdNdHRiO/b21BJN5zg9qYjW4K1t6XfA0V5H\nBUCl76vT37sv0fRpF02Q2GdFWte95J9riaYdtEQjDqcbmegbprUtm7k7vBCGDL7RJpkCwHw9uXg/\ntllsrN1aNNpF0zGNkkJt3Uvr5512fu2iCbqapDrZol+c/k7xGnyuLsD7y+DD4aKxBnggOJlGlVo0\nMQvwzZ0KdonGbSUWN+KXwTdKNDJsfXGoRePcN8UWbxMt0YQhwFttkoD/JxevmDdTXYsmQFo5Fcwv\nGVECnRTOilOA52UdCgHMvIy3RCNy0Nou0fjRo8MS4Bsz+GCcNLVkUbtoAqOVU8HvI2ycArxZSbJ+\n5qVe8KNTwifRqB/g7Ro8EKREk/Mg0egALxS3ySSmrGCVAzpZ/aVxhZboumhayV3ijuPmoomaTVJ+\nezb2b78BvtFFo1J/d9fg5Us0NUuqdtEEhrNEk4ZdouHv6wy+Gc6SQpA++Gi1rdNYgxybZHwyeLsP\nHuiWRONBxJ31AAAfZElEQVQtg+fF2rSLpmPcJZpCw+edBfj4uGiCkRT4zEt7IAHi4aKRbZP0l8mq\n76JRRaJxLzamXTRCaVXQyZqV8tHv9iWaOAd4OTMv63X+2rGi1bbBafDW/h31DF4Fm6R20QSGfUk0\noD4oWTV6LdE0xwy+VmRKCnai1rZu7akHWTvHySbZSeLWCZ1INDrA+6RVlqRdNN4JWlJoJAHAgLl0\nXNhxbs+00BtmYwITR5uk/AzeXBmutniQtkkGgjebpD8XTeMFi45ObCVoScEOl23Ucm34wb09Rfrg\ng5jopM716JYGX5tPQ12XaBpHryJMq4JOYlw0VutYtKx8VoKy9bln8NYvRI/rNmEh6PYUI9E42STV\n6e/2DD6TyWJ8/Als3ZrDE09sxvDwAgwNzRd+XK/t3Gir1gG+Y5xmXgJiJzo1rtDCq+sxxhwHCsOM\n+8zLoCQa9TJGPwTdnv6LjanvorFq8JlMFkuXrsf7338qtm49BaOjF2F8fBkACA/y9fMNmkk0zvFC\nJLGRaJxmXgLNNfh2vwCNj1zRkhGsBKPBO9ehqR1PrYDSKZlMFqOj9+Kqq9ZhcHA5MpksgDBk8OoH\neDODHxkZxfj4ahQKfUin+d89Pr4aa9feLvy4XpQAXueqLH2iU2wyePeaJrXJOVaXTScaZfNOH62m\ndmpPUzMW9cTiXaIJL2Zmef75b8Wjj34IY2OLq5nlvHniA3xPz5EAeP8ul1/veF9ufd0wDvg6R5EY\nRu2JPZ/n51oo9KGnp5a45XJJx9/1d9xav3Wzo/KkL1n3PdGDrB2SyWRx3nkrsWdPqS5DAuTWouH7\nD38QcsJ55mUS1vLL/o/hRaIJd9uamWU6nUOxyNvTzCxVnugUhr5ulWh6e/l5FYu91QweAPr6xD9d\n1w9mO7dzUO0X+QBvZkh3330x9u07EqOjV2Lp0vXVIN/cJulPogHU6/SicC8hIC4otc7g1RrU6wQz\ns0yncygUan9rLpfUpQp8YpVohocXoL9/WZ1E099/OZYsOUfCcb1INDrAC6GWIeWrXyCr9lafwefh\npxiTm7NApWXMROGWXYu0SjazSQLqBZROMDPLnp58XYDv6yvX1UkSgexiY6qNN1ltkkND87FmzSBm\nzfpfvPGNYxgcvAJr1iwMwEXjvIqUc6zQGnzbWDMk8xEYqGlvjaUKxNWi4fsPfxBywr3Ko7gA71T9\ns/5YKYT95jk8vADj48sqGTxvT55ZLqz0HVUlGvX7ut0mOTQ0H3PmfAD792/C1772dWnHbVzK0jRs\nWNvLqf3E3yAjH+DNDGnChH3I5SZW3ze1t2YafLH4WlvHCsNjqyjK5X1IpSY3vC9SVojDIKuZQR44\n8GPMnn0dZszIYMkSnlnu3fsAxC74ET+JpnFFp4kol/dLPa796dZs62SyFn+Car/IB3gzQ5o48d3Y\nt28KgFqGBDhp8J0XYwpDpxdFubwHfX0zGt4Xr8FHW6IBeJC/555e/Pznl6Ovb3r1fTk2yXgUGzMz\nYbsMkkpNQam0R+qx7YmJkxymA7wgzAzpjju+j4kTt2Nw8IpqhgTILTbG969OpxdJqbQbqdRhDe+L\n1eCj76IxcWpP0aUf4lSLxqlMAQCkUoehVNot+diNEo29rXWAFwjX3jZj376p+NKX6rU367qX2ibp\nnVJpD1KpKQ3vi1xHNA4uGoBnm+XyQSSTh9a9ryc6dQ6fte4U4LuRwTe2tXbRCIYHJKeM07qik4hi\nY2ovYyaKUmkPksnG9iQS5/yIi0RTKu1FKjUJRPVfRzk2SXEuGmsdFUAtxxjX39MN76dSh6Fclh3g\n7Rq8m0Rjb79E5TNxFVJjFOCdJQUzS2KM2Waytpfh8IvCYG/SqAQhO+pINGoEFD+USrtdbpZyJzpF\nuRZNtyUa6yRAZ4mmsf0A8U+lsQnw5XIzSaFYsVSlqnfV9gN8ufL79VP0Ver0Imkm0QQ30SkabevW\nljJ88PGRaJwDfDJ5GEql1yu1YOTQqUQDiE9aPAV4IlpIRI8R0eNEdJnD55cS0WYi2kREWSKaJewM\nBeEm0ZgBqXFgpLetDKf5BVOj04uE3zCd21OkTTIOxcaataUsF428evBqXA8niyTAJx4lEn0ol/dJ\nO7b9ydPJsRRUvGgZ4Il/w64BsBDAyQA+TESn2Da7F8BcxthsANcD+J6wMxSE22OwKSm4eVe9EqcA\nbxgllMsHGgYFAdEZfDxcNM3kQ5nlgqPsonHL4AH5Mo1zBt/aRQN0IcADmAdgM2NsJ+NHvhHAkHUD\nxtgGVkt3/wTgBGFnKIhWkoKXx6pmxCnAl8t7kUpNbhgUBGSUKohDgHeTaGSs6BQPicZNgwfkO2mc\nbZKKZvAApgHYbnm9o/KeG58H8Fs/JyWDZhKNYRQdBkbac9E4uQr4/qPnonFz0ABBT3SKhk2ylXwo\nQi+2L3gjoxaNSn29VQYv00ljL7Hh1UUDiA/wXnzwnnsXEZ0PYC6AM5w+X7lyZfXngYEBDAwMeN21\nb5q7aApaomkDt7YExPrgvdSiiULburtoEjDLLzv1rXbgtdFrC95wO2sJjBmOT2KtUD2DZ6zgqMED\n3ZFo7ON5zeLFnXduwN13Py7kXLz0mh0AplteT0d9Rl85MTobwDIA85lLCmcN8EHj5qIxnQr+JRo3\n25M6nV4UbpICINb54U2iUSNj9EOptAfp9LGOn9WeiPwG+Pq2rC0InUcyOaGDPYbBJtnogwdUkWjc\nbtpJnH7632PBgsXVd1atWtXxuXi5dT8AYDYRnUD8mWcxgHXWDSqDrj8EsIgx9krHZyOJVoOCpkRj\nvSh85Fu7aJxwc30AuhZNJzR7IhI1puG85mvnMo36Gby7RMOtkjIDvFeJRgENnjGWA/AFAOsBPALg\nZsbYRiJaRUTvq2z2LQATAdxERA8R0f8IO0MBlMt7kUw2zhQErDqnlmi80lqi0S6adgjihukkd/lx\n0oQhwKsk0XTLRePpuY8xtg62rJ0xtsLys/hlUQTSXFLQLpp2adaeeqJT+zRvz7SQG6ZTW/px0oQh\nwKvkorH77pXJ4KNAq4zTeaITb3z7SixuuI2K81Vu1Oj0omjmoglyRaeotG23JBo/Thp3F40a16O5\nTVKui8arROPsutMBvm1aZZxcg3d7hPX2BQhq6rEK8IAURAYfD4mm9RORf1eSW//utB6NewavRl9v\nncHHQ6KJRYD3onH6f4SNm4smCA0+Pi6aZvMK1JRo1HbRtNbg5Uo0jcXGvLloRFfkjEWA9/II7PcR\nNk4afOsbprh68HGoReNFQvSLHIlG3QDPff/dctHYLakKu2iiQGuJpnGiE9DeI2ycAnwziUaUZmzO\nvIy6TbJczgFgrk8qYjV4PcgKBCHRtHbk6QAvELeZgkBtgQq/X4B4BfhmGbwo10f9zEvnY4W/bc2n\nIbe/U1wGH0ebpNtEpyAy+PaX7AN0gO+IZgHJzSbJP9MB3okgatG0dtCo5drolGZ9E1BZg1c/wHfP\nReO3HrwO8G3hxbdtGPkGvbedgmPuxYPUWcZMFEFINK0GWAG1AkqnNOubgLjSD07jGaJtkty2qkZf\nb2aTTCYPRbl8wLMFuh0YYw3SYjeLjcUkwLeeeek/g9cuGkDuzMvGY6UQ9ptnM/kQUFmiUd9F4xbg\niRJIpSajXN4r4bh8cNc6a955LE8v2SeMIGyScZFo+KCg4Rp8ZUoKjccKf9t6kWj0IGv7NLNJAqaT\nRvxAq7NbSfEl+8KOl1IFPGN0KsakXTRWzKqc8gcFvWjw4W9btyqnJqLKL+tiY/XIKlfg/KSkbZJS\nae6i0Rl8O7TKOEWtQtRqFisQjbZtJh8CIm2nQbho+HiTzAWtvdLMBw/Ic9K4xxHtopFG65mXBR3g\nPeJFM9YSjXeiJdEQVBlobZ3BqyzR6ADfFl5dNM4z/fy7aFTo8KJo5frQNsn28NKe4m6YYl00bktU\nqjDw3UqDT6WmSLFKtifRaBeNb1oNCtZ88Ae1i8YDrSUaUbY+rxJN94OJH7rtoum82Jja/b3Zik6A\nChKNdtEIofWgYBJAAuXyAS3ReKCZBx6QKyk0Hiv8bRukD162RAOoc028DbLKkWjs8w3M2fKMGZbz\n0y4aIbQaxAJ4UCqXX3d8hNW1aOppZjkFxEoKzQqN8WOFv21b9c8wFRsD1Lkm3mySwUg0RFQp4VGL\nJVqDF0SzafUmiYQZ4HUG34rgBgXj4aIJ7oYp30UDqHNNvA2yBiPRAI1trdSSfWGm1SMwYGbw+xy/\nAMXiy56Oo3qHF0WptBsTJvS7fq5LFbRHa4lGlO003BJNJpPFyMgo8vkUentLGB5egKGh+a7bNytV\nAMiVaJzMAfa21gFeEF4lmmLxVZdaHd4lGjdXQdiDkJVWT0QyJ+Y0Hiv8bRukRCO/Fo2ca5LJZLF0\n6XqMj6+uvjc+vgwAXIO8WTLADVkFx9xKbNjbWrtoBNFKUgBQ6QhliS6acDs9rAQl0XitRRPmtmWM\noVTai2Rysus2qko0zVcwE3tNRkZG64I7AIyPr8batbe7/k5rm2TwEo11PC8oF1LkA3yrqeAAqh1B\na/Ct0S4acZTL+5FI9DUNRCpOdOIVE4PJQAEgn3cWGnI5p0XuOd2VaNw0+NYSjeiJYpEP8N4kGu6X\n9VOrw+mCZTJZXHnlr3HnnY9hcHA5MplsG2euJq0GBcXa+sLtoslkshgcXI6BgZWO199L3xQ7piGq\nFo0BIAFrxcTaPsVfk95e5/319bkHQj7I6u6Dl+eicdPgnSQarcH7plTag76+Nzbdht/pEw0N3s4j\nLGOluju3qRsee+zHMXXqzzA6emVL3TAMBLdARbhdNF5041Y3S6DmofaLyFo07tmnnGsyPLwA4+PL\n6tqyv/9yLFmysMk5dkuicdPgtYtGCl5cNIlEDxKJvobJUH4kGlM3PProDJJJfsG4bnhFyAO8lmi8\n4K4b166/V4eXais6BR3gzfZau/YK5HJJ9PWVsWTJwqbfo1Y2Sd4WDOVyDslk837WDn4lGh3g28Sr\ni8bLY1Uz7JqkqRuWy8lqgAea64aq42VQUKak0Ii6LhovunG3JZpOXTRu+jsgz9k0NDS/rcSolQZP\nRFUnjegA39NzZMP7ThKNs+tOB/i28Oqi8fJY1Zz6UXFTNyyXU0gkalphM91QdfigYG9Ag4LhdtF4\n0Y29TMJTsRaNmwMEUOeatMrggZpMk05P7egYTt78E0/0KtEE46KJQYBvLikArQJ8ZxKNqRuWy+dU\nM/hWuqHqeG1LXS7Ym27svT39zStgjHkuY+ttf8FKNJ3QSoMH/Dlp3MZYrrnmGZx00t82bO9dohH7\nBBT5AO9lIItr8CIkmlpzmo+Tv/nNT3HEEVsxOHhFS91QdbwNCoqbeRnmWjRedGOvT5d+M3iezSYb\nZBVzIh9jrGH8qfn+uh/gW81sbbXgB+DPSeM2xvLkkydj9uzGfmt3LDXT4EWWW458gPc6kOVlckIz\nnC7Y0NB8nHZaGk8/vQW33fZ17yetKF4HrHnlvPaChp2wu2iA1rqx1zkafgO8W1sSJSpt2NxSaKfb\nAd6LQ8mbRNP5sn1uYyxEzqXJu+WiibQPng8K7mk6KAjIkWhq+1Y7CLWDtwHrJADyrcOGXaLxglcD\ngP8A796Wncg03Q7wXma2etfgO5No3MZY0umiUi6a0AX4VpNHrJgFxFppcYlE2lEO8CPRmKgWhNpp\nPzteBgUBMd7tKEx0aoU3iSbte0yjmdzViZOm2wHei0OJa/DNn0r81KMZHl6A/v5lde/191+ON75x\nSscTnTKZLK666mb84Q8PC5sYGSqJpt2iQ14kBUCMi4aPiqu7hBnQWdEmK17bsyYrTOj4XL3VovHX\ntnYd913vOh733POc54qFfmlH8vJDM7mr3clOmUwW//mf1+ODH9yFf/mX5Q5tJG6qvZvO7sWh1Mom\nCfiTaNzGWI466h6Pcm99vDC/m9OmfRSHHfZfdRMj/RCqAO9l8ogVL4/AgHuAr63EUsatt/6p6aBO\nGDL4dtvPjpf2zGSySCRKWLToShhGuuMgKVuiabzZZfHHP16PUumH1W1kzzxutVwfoJZEY7aZYXwR\nZ5/9J8fZ2W7XpN1yv82SES8OJa8SzcGD4x7+cmecxlg2bmxs60wmi/vuG0Nv70FkswzDwwtwzDHO\nEyOPP/7mhomRfghVgG+36JAX1wfgHuBvvXUD0ukU5s79HLZs6cHBg+5f/jAE+E6KNllp1Z7ml/Jb\n35qM++67BLt2HdtxkJQd4BtvdqN1wR2QP/O4lUSTyWTx61//B84551l87WtO2bL771mD6cUXz8DR\nRzcP8F4CsNlmb3rTQyiXeV+yt5HTNWnnydE8jwceeBq7dv2q7jPzWKZhwcmhZP7+V79awNDQvzZ1\nrsmoR2OXFs2//eSTP4Bp056q3hR/9KO9dfGiNjGyft6M34mRLQM8ES0E8G3waVfXMcausn3eC+Dn\nAP4WwF4AH2WMPevrrFxofDTLAhjFo49ux+Bg4xegPUmhXjczL8y//dtEPP30FBw8+N26z+0dOwwB\n3vnRNotNm57AwMDKlplVqbQbvb3TXfdvBoBy+TqkUjzrtGYh7S3YkJeqwTfe7Pzd/DqhmYvG7H/A\nME47baPnWkZOwZToE7j00kYZJpPJYu/eV7BixQqMjR3TNIEBam2WTJaqAR6obyPzmlhvGJs2PYFX\nX72x7thON8/6c1/p+PeZx3LKns3ff+aZf8Vll30D69d/A08/7d5mbj74dp82rNjlMPM78da3/gjp\ndL76t+/efY3rxEjrzHe/EyObBvhK8L4GwGkAXgRwDxGNMsYesmz2JQDPM8bOI6IPABgB8H5fZ2XD\nbPCdO1/GhAkXVTpiFsB6AKuxaxcwOtrYKb08AmcyWTzxxIMADFx8ce0mYV6YQuGnSDeM1fAby733\n7sDg4HK8613HY/r0x3HvvTdi27b76zqEUxAaGxvDwMCAz1ZpjlMnbXy0zSKVuh6vvnoj7ryTv9Ms\ns/rrXzfgySdfxFNPPefY6c0AUCr1oKenALOdNmzYhGz25boA8uijn8Fxx/0KhnEAiUQKQBqTJx9T\nPddJk+Rm8I03O/fkYe7czzacn/m3dxIMzN+55JKX8JGPjOCiixZhaGh+Xb8w+9/MmZuRSvF5BV6e\nKJxkuJdfvgDbt29oOIelS9dj6dJ+bN8+qe7a8GMN4oILrsbs2X/E3r07AKTx17++BqAxwJtBKJPJ\nYnz8MTzyyA7ccMNbLPtc6XiuuVyyyY3AWTpxCnj2rD+dPohSqafyd7i3mVPBMb/jVPYnT/M7USj0\nIZ2uSWFEzGVi5BkNEyPXr7+y5XHdaJXBzwOwmTG2k58U3QhgCIA1wL8XwFcqP98C4MdERIwxZt+Z\n/YtiHdQyO9Hkycc0/Pz885Pxwgvfq+wliwkTzkMyWca+fTfV7d/eKf/hH57CzJn7ccklzo+35sVc\nsGA+kslSXZZUf2Gs2U/txrJnDzA6yrXbb3zjjXj44Qvw4IML6jqENQiZHXHLlrtwxBFvcvx7m7WL\nnzbbsOFq9Pcfh8mT92Hu3C9i0qSjHTKrLMbHCR//+E8xc+bPG/a3dOkb8OCDn8a9977PsdObQbNU\n6kEyeTeAJwGsxoEDywFcWXecF144Fi+8sBrApwBMBVC76WzY8AP8/vcHMTR0pesjdiaTxQ9/+Dt8\n6UsHqjfZdtqsUChbkgUAWIBU6qKKTGNNHrLYtWu95fxqN6dCoYytW8myD7Odb8bxxx/qeCPYufNl\nbN1KyOd/gK985ZvIZK7CE0/wJ5wHHqgFeOvN0nwaArK4//6/YGBgpeu1f+SR7baWyiKd/hFeeaVY\n95RbS2DORjpt/6ryv5/f+Gttwd9fhmRyYTXAm0HI/C598IPvxAsvPG+7YTjfPB988DF85CMvuNwI\nUgCW1bW700xwp6w/lSqiXK7p705PYZlMFjfe+Euce+4WfPWrtf7TKA3VvhPvfOdo04H4TCYLw3gZ\nixdfjUJhAoaHF1S/E/UBPotyOY/zzhtBPj+xbh+33PIjHHXUX8RNjOTF+53/AfgogGssr/8PgB/a\nttkC4BjL6ycBTHXYFwMuZwCr/LuTpVKfr/5c++xO23bLLD/X/h1++Cds7zXu46Mf/Qb73Oe+wgDG\n+vsvZ7///Z3MyoIFfN8XXngZ+/znL63ua3BwefWz6647kc2YcZ1l3/bz4a+/+90z2dy5f6jbB2OM\n5fMvsrvuOor9/vd3sv5+cx+fbPL3urVLp21m/6zWFmecsaJp+9n3t3btP7DZszc0/I0m5t947bV/\ny2bNutDy+yuYU5vxf6c3nEMqlWe3355yvW7mcXp797PbbuvrsM346wkTFrPZs5eywcHlbMWKq9ng\n4HJ2+OHnNel/1n14a+f6a89/Z9KkV9ktt0ypa8sVK1Y09M1jj93KbrjhDU2uT+trP3/+TWzVqg86\nXvtvfvO97NRTFze5Po1//5w5H2Nr105lg4PLq9fGPN9LLvkcW7RoqEmbubWf/fWKyrbL2eGHf6Lu\nWE7fYevvT578SkPbOvWf444bZ9dfP9PWR9y+E/a+5HyNf/vbw9nkya9UP1ux4mrW3385e/e7/4dd\neeWi6j5uuulYduSROxv6+K5dd7KNG0+vO18ept3jdLN/rTJ45u/2Uc/q1Y8CWFR59SSAt1Reu/0M\nAE8BeKRhXz09G1EsLrK807iPGTOexLp1nwbg/KhmzdKt5HJJXHrpezA+vgz5/CG45JKfYv/+IpLJ\n4wEkUC5bz4efX3//I9XHQnMfAHfilEq78dxzn8GnP83P7447/oSzznq3y9/brC06aTP7Z5znn/81\nzj+fYcGCP3s4Lt/fzJmbsG9fTTO2Z0Zm277++s/w5S//Fnv2PG/Z958tW9bO7447NuOss+r/vmTy\n/cjnucXS6bqZ2WcyWUQqVcDq1Z/poM1qHHPMrZg3728APIwPfQi4++6n8Oqrbv3Pug9v7cwYq157\n83d6ew82bUvzkX337i/iiCNeaPI32o/7KoA3V97jnx111Dps2/bWuvbs7WWV4x6CT37yESxaZP6O\n/W9u/P5NmfIyensn1s3ONr9LxWIaH/jA/Tj11Pp2AF5FOn0CGAOKxbkO7Vd/7nfcsQVnnfVnTJz4\nJE466WhMnfoagIfx2GPfrtvrP/3TUzjjjEfqfr+nZxaKRT5+45T1m/3n0EN34cgjn3NoW6fvhNNr\njvUaT5y4txpPxsdX4957r8CaNYP43e/+HXPm3IFvf/tBFArvwKRJrznKSKLH7MhBSal9SHQ6gMsY\nY++rvL4UQJoxttqyzR2VbR4kvsTLi+AZvGHbl9CbhUaj0cQFxlhHdT9aZfAPAJhNRCcAeAnAYgCf\nt21zK4CPAXgQfHD1Hntw93OCGo1Go+mMpgGeMZYjoi+Aj7IkAPyCMbaRiFYBeJAx9jsAPwDwCyJ6\nDMDr4Lq9RqPRaLpMU4lGo9FoNOFFerExIlpIRI8R0eNEdJns46kEEU0nomzl799CRF+pvH8EEd1O\nRI8S0Xoiaj0bKyIQUZKIHiKi31VezyKieypt9CtqNb88IhDRFCL6NRE9QkRPENGpce0XRLSKiJ4i\noieJ6CYiOiQu/YKIfkZEL1YUEPM9135ARCNEtJmINhLRKa32LzXAWyZKLQRwMoAPezmpCFEA8H8Z\nY28D8A4AnyWitwNYBSDDGDsZwLrK67iwFMDjqDm0RgBcVWmjF8AnzsWBHwO4mTH2dgAngbdJ7PoF\nEb0JwMcBzGaMvQW8etw/Iz794lrw+GjFsR8Q0YcAzGCMnQTgM5XfbYrsDL46UYpx7485USoWMMZe\nZIxtqvy8D8CjAE4Anxz2i8pmv0RM2oSIpoH/7T/hLykJ4FTG2P9UNolFWxDRkQDmMMZuAADGmMEY\n24t49ovXABQBTCQ+tfMQANsQk37BGNsAYJftbbd+MGS+z3g1gVTlO+WK7AA/DYB1at2Oynuxg4hm\nAngngLsAHM0YexUAGGOvADime2cWKN8HcCkA02V1DIBXLJ/vRDz6x5sBvExE/0VEm4jo50Q0CTHs\nF4yx1wB8FzyoPwdgN4BNiGe/MHHrByegzXgqO8DrEVwARHQogJsALK1karGDiN4H4KVK5mFaZuNq\nnU2A3+y/zRibDZ7F+qsLG1KIqB/AxQBmAjgewKEAzunmOSmO/TvTNMbKDvA7AFjLD05H/R0o8lQG\nh/4bwH9aHjlfJqKjKp8fDT7HIOr8A4B/JKJnANwA4D0ArgJwlGWbaeB9JupsB7CTMfZA5fVNAOYA\neCmG/eLvAdzNGHu1IuPeDGA+4tkvTNzigz2etmwX2QG+OlGqEugWgw8axALiq07/FMDjjLHvWz4y\nJ4eh8v+tQZ9b0DDGLmeMTWeMzQKvafRHxtjHAdxbqUIKxKcttgN4hYhOrLx1NoAnwL8bseoXAJ4G\ncCoRTah8X84GrwkQu35hwS0+3ArgfAAgorkAyqxSCNIN6T54IjoXvJ68OVHqm1IPqBBEdBp46bxH\nUXuU+hqA+8EHnKeCOwQWM8Y6W/03hBDRGQC+zBj7RyKaBeB68EfzzQA+zvwuYRQCKm6qn4APKj4L\n/sUlxLBfENFK8L/fAK9U+0kAxyEG/YKIbgBwBvgTy4sA/gXAb+HSD4joBwDOBJAH8FnG2Mam+9cT\nnTQajSaaSJ/opNFoNJruoAO8RqPRRBQd4DUajSai6ACv0Wg0EUUHeI1Go4koOsBrNBpNRNEBXqPR\naCKKDvAajUYTUf4/riwi5B3MQAQAAAAASUVORK5CYII=\n",
       "text": [
        "<matplotlib.figure.Figure at 0x7fa6d0e74a90>"
       ]
      }
     ],
     "prompt_number": 9
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "def max_finder (listt) :\n",
      "    m = 0\n",
      "    for i in range (len(listt)) :\n",
      "        if listt[i] > listt[m] :\n",
      "            m = i\n",
      "    return m\n",
      "def compare (out,learn,n = nlayers,m = 4) : # m is ntracks\n",
      "    counter = 0\n",
      "    output = np.zeros(len(out) , dtype=float32)\n",
      "    for i in range (len(out)) :\n",
      "        output[i] = out[i]\n",
      "    learned = learn\n",
      "    for track in range (m) :\n",
      "        subcounter = 0\n",
      "        for layer in range (n) :\n",
      "            start = int(float(layer) * len(output) / n)\n",
      "            end = int(float(layer+1) * len(output) / n)\n",
      "            subOutput = output[start:end]\n",
      "            subLearned = learned[start:end]\n",
      "            if max_finder(subOutput) == max_finder(subLearned) :\n",
      "                subcounter += 1\n",
      "            output[max_finder(subOutput)+start] = 0\n",
      "            learned[max_finder(subLearned)+start] = 0\n",
      "        if subcounter == n :\n",
      "            counter += 1\n",
      "    return counter"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 63
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "counter = np.zeros(5 , dtype=float32)\n",
      "start = time.mktime(time.gmtime())\n",
      "for i in range (200) :\n",
      "    counter[compare(ds_output[cut+2000+i],trainer.network.predict(valid[0])[2000+i],4,4)] += 1\n",
      "stop = time.mktime (time.gmtime())\n",
      "spent.append (stop - start)\n",
      "print \"\\t\" , stop - start\n",
      "print counter"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "ename": "KeyboardInterrupt",
       "evalue": "",
       "output_type": "pyerr",
       "traceback": [
        "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[1;31mKeyboardInterrupt\u001b[0m                         Traceback (most recent call last)",
        "\u001b[1;32m<ipython-input-74-ef500169a963>\u001b[0m in \u001b[0;36m<module>\u001b[1;34m()\u001b[0m\n\u001b[0;32m      2\u001b[0m \u001b[0mstart\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mtime\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmktime\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mtime\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mgmtime\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m      3\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mi\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mrange\u001b[0m \u001b[1;33m(\u001b[0m\u001b[1;36m200\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m----> 4\u001b[1;33m     \u001b[0mcounter\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mcompare\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mds_output\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mi\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0mtrainer\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mnetwork\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mpredict\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mtrain\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mi\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m4\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m4\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;33m+=\u001b[0m \u001b[1;36m1\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m      5\u001b[0m \u001b[0mstop\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mtime\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmktime\u001b[0m \u001b[1;33m(\u001b[0m\u001b[0mtime\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mgmtime\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m      6\u001b[0m \u001b[0mspent\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mappend\u001b[0m \u001b[1;33m(\u001b[0m\u001b[0mstop\u001b[0m \u001b[1;33m-\u001b[0m \u001b[0mstart\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
        "\u001b[1;32m/tmp/theanets/theanets/graph.pyc\u001b[0m in \u001b[0;36mpredict\u001b[1;34m(self, x)\u001b[0m\n\u001b[0;32m    549\u001b[0m             \u001b[0mvariables\u001b[0m\u001b[1;33m.\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    550\u001b[0m         '''\n\u001b[1;32m--> 551\u001b[1;33m         \u001b[1;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mfeed_forward\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mx\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0moutput_name\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m    552\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    553\u001b[0m     \u001b[1;32mdef\u001b[0m \u001b[0mscore\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mx\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0my\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mw\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mNone\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
        "\u001b[1;32m/tmp/theanets/theanets/graph.pyc\u001b[0m in \u001b[0;36mfeed_forward\u001b[1;34m(self, x, **kwargs)\u001b[0m\n\u001b[0;32m    530\u001b[0m             )\n\u001b[0;32m    531\u001b[0m         \u001b[0mlabels\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mf\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_functions\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mkey\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 532\u001b[1;33m         \u001b[1;32mreturn\u001b[0m \u001b[0mdict\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mzip\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mlabels\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mf\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mx\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m    533\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    534\u001b[0m     \u001b[1;32mdef\u001b[0m \u001b[0mpredict\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mx\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
        "\u001b[1;32m/usr/local/lib/python2.7/site-packages/theano/compile/function_module.pyc\u001b[0m in \u001b[0;36m__call__\u001b[1;34m(self, *args, **kwargs)\u001b[0m\n\u001b[0;32m    593\u001b[0m         \u001b[0mt0_fn\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mtime\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mtime\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    594\u001b[0m         \u001b[1;32mtry\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 595\u001b[1;33m             \u001b[0moutputs\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mfn\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m    596\u001b[0m         \u001b[1;32mexcept\u001b[0m \u001b[0mException\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    597\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mhasattr\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mfn\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m'position_of_error'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
        "\u001b[1;31mKeyboardInterrupt\u001b[0m: "
       ]
      }
     ],
     "prompt_number": 74
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "gots = np.zeros(5 , dtype=float32)\n",
      "start = time.mktime(time.gmtime())\n",
      "for i in range (200) :\n",
      "    gots[int(4*got(ds_output[cut+i],trainer.network.predict(valid[0])[i],0))] += 1\n",
      "stop = time.mktime (time.gmtime())\n",
      "spent.append (stop - start)\n",
      "print \"\\t\" , stop - start\n",
      "print \"gots :\\t\",gots"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\t831.0\n",
        "gots :\t[ 10.  42.  50.  54.  44.]\n"
       ]
      }
     ],
     "prompt_number": 71
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "print counter"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "[ 38.  33.   3.   1.   0.]\n"
       ]
      }
     ],
     "prompt_number": 75
    }
   ],
   "metadata": {}
  }
 ]
}