RunCanny.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "import pickle\n",
    "from canny_cf import *\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline\n",
    "%load_ext autoreload\n",
    "%autoreload 2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Ycol shape: (28066, 126280)\n",
      "Ytrain before removing empty: (28066, 126280)\n",
      "Ytrain after removing empty: (27892, 124574)\n",
      "Ytest before removing empty: (5439, 23550)\n",
      "Ytest after removing empty: (3566, 9824)\n",
      "n=27892\n",
      "m=124574\n",
      "k=10\n",
      "Starting Iterations\n",
      "========================\n",
      "evaluating user: 0\n",
      "evaluating user: 1000\n",
      "evaluating user: 2000\n",
      "evaluating user: 3000\n",
      "evaluating user: 4000\n",
      "evaluating user: 5000\n",
      "evaluating user: 6000\n",
      "evaluating user: 7000\n",
      "evaluating user: 8000\n",
      "evaluating user: 9000\n",
      "evaluating user: 10000\n",
      "evaluating user: 11000\n",
      "evaluating user: 12000\n",
      "evaluating user: 13000\n",
      "evaluating user: 14000\n",
      "evaluating user: 15000\n",
      "evaluating user: 16000\n",
      "evaluating user: 17000\n",
      "evaluating user: 18000\n",
      "evaluating user: 19000\n",
      "evaluating user: 20000\n",
      "evaluating user: 21000\n",
      "evaluating user: 22000\n",
      "evaluating user: 23000\n",
      "evaluating user: 24000\n",
      "evaluating user: 25000\n",
      "evaluating user: 26000\n",
      "evaluating user: 27000\n",
      "evaluating user: 28000\n",
      "evaluating user: 29000\n",
      "evaluating user: 30000\n",
      "evaluating user: 31000\n",
      "evaluating user: 32000\n",
      "evaluating user: 33000\n",
      "evaluating user: 34000\n",
      "evaluating user: 35000\n",
      "evaluating user: 36000\n",
      "evaluating user: 37000\n",
      "evaluating user: 38000\n",
      "evaluating user: 39000\n",
      "evaluating user: 40000\n",
      "evaluating user: 41000\n",
      "evaluating user: 42000\n",
      "evaluating user: 43000\n",
      "evaluating user: 44000\n",
      "evaluating user: 45000\n",
      "evaluating user: 46000\n",
      "evaluating user: 47000\n",
      "evaluating user: 48000\n",
      "evaluating user: 49000\n",
      "evaluating user: 50000\n",
      "evaluating user: 51000\n",
      "evaluating user: 52000\n",
      "evaluating user: 53000\n",
      "evaluating user: 54000\n",
      "evaluating user: 55000\n",
      "evaluating user: 56000\n",
      "evaluating user: 57000\n",
      "evaluating user: 58000\n",
      "evaluating user: 59000\n",
      "evaluating user: 60000\n",
      "evaluating user: 61000\n",
      "evaluating user: 62000\n",
      "evaluating user: 63000\n",
      "evaluating user: 64000\n",
      "evaluating user: 65000\n",
      "evaluating user: 66000\n",
      "evaluating user: 67000\n",
      "evaluating user: 68000\n",
      "evaluating user: 69000\n",
      "evaluating user: 70000\n",
      "evaluating user: 71000\n",
      "evaluating user: 72000\n",
      "evaluating user: 73000\n",
      "evaluating user: 74000\n",
      "evaluating user: 75000\n",
      "evaluating user: 76000\n",
      "evaluating user: 77000\n",
      "evaluating user: 78000\n",
      "evaluating user: 79000\n",
      "evaluating user: 80000\n",
      "evaluating user: 81000\n",
      "evaluating user: 82000\n",
      "evaluating user: 83000\n",
      "evaluating user: 84000\n",
      "evaluating user: 85000\n",
      "evaluating user: 86000\n",
      "evaluating user: 87000\n",
      "evaluating user: 88000\n",
      "evaluating user: 89000\n",
      "evaluating user: 90000\n",
      "evaluating user: 91000\n",
      "evaluating user: 92000\n",
      "evaluating user: 93000\n",
      "evaluating user: 94000\n",
      "evaluating user: 95000\n",
      "evaluating user: 96000\n",
      "evaluating user: 97000\n",
      "evaluating user: 98000\n",
      "evaluating user: 99000\n",
      "evaluating user: 100000\n",
      "evaluating user: 101000\n",
      "evaluating user: 102000\n",
      "evaluating user: 103000\n",
      "evaluating user: 104000\n",
      "evaluating user: 105000\n",
      "evaluating user: 106000\n",
      "evaluating user: 107000\n",
      "evaluating user: 108000\n",
      "evaluating user: 109000\n",
      "evaluating user: 110000\n",
      "evaluating user: 111000\n",
      "evaluating user: 112000\n",
      "evaluating user: 113000\n",
      "evaluating user: 114000\n",
      "evaluating user: 115000\n",
      "evaluating user: 116000\n",
      "evaluating user: 117000\n",
      "evaluating user: 118000\n",
      "evaluating user: 119000\n",
      "evaluating user: 120000\n",
      "evaluating user: 121000\n",
      "evaluating user: 122000\n",
      "evaluating user: 123000\n",
      "evaluating user: 124000\n",
      "iter: 0 \tlam_diff:  1317996.5464\tpsi_diff:  11.1685\n",
      "evaluating user: 0\n",
      "evaluating user: 1000\n",
      "evaluating user: 2000\n",
      "evaluating user: 3000\n",
      "evaluating user: 4000\n",
      "evaluating user: 5000\n",
      "evaluating user: 6000\n",
      "evaluating user: 7000\n",
      "evaluating user: 8000\n",
      "evaluating user: 9000\n",
      "evaluating user: 10000\n",
      "evaluating user: 11000\n",
      "evaluating user: 12000\n",
      "evaluating user: 13000\n",
      "evaluating user: 14000\n",
      "evaluating user: 15000\n",
      "evaluating user: 16000\n",
      "evaluating user: 17000\n",
      "evaluating user: 18000\n",
      "evaluating user: 19000\n",
      "evaluating user: 20000\n",
      "evaluating user: 21000\n",
      "evaluating user: 22000\n",
      "evaluating user: 23000\n",
      "evaluating user: 24000\n",
      "evaluating user: 25000\n",
      "evaluating user: 26000\n",
      "evaluating user: 27000\n",
      "evaluating user: 28000\n",
      "evaluating user: 29000\n",
      "evaluating user: 30000\n",
      "evaluating user: 31000\n",
      "evaluating user: 32000\n",
      "evaluating user: 33000\n",
      "evaluating user: 34000\n",
      "evaluating user: 35000\n",
      "evaluating user: 36000\n",
      "evaluating user: 37000\n",
      "evaluating user: 38000\n",
      "evaluating user: 39000\n",
      "evaluating user: 40000\n",
      "evaluating user: 41000\n",
      "evaluating user: 42000\n",
      "evaluating user: 43000\n",
      "evaluating user: 44000\n",
      "evaluating user: 45000\n",
      "evaluating user: 46000\n",
      "evaluating user: 47000\n",
      "evaluating user: 48000\n",
      "evaluating user: 49000\n",
      "evaluating user: 50000\n",
      "evaluating user: 51000\n",
      "evaluating user: 52000\n",
      "evaluating user: 53000\n",
      "evaluating user: 54000\n",
      "evaluating user: 55000\n",
      "evaluating user: 56000\n",
      "evaluating user: 57000\n",
      "evaluating user: 58000\n",
      "evaluating user: 59000\n",
      "evaluating user: 60000\n",
      "evaluating user: 61000\n",
      "evaluating user: 62000\n",
      "evaluating user: 63000\n",
      "evaluating user: 64000\n",
      "evaluating user: 65000\n",
      "evaluating user: 66000\n",
      "evaluating user: 67000\n",
      "evaluating user: 68000\n",
      "evaluating user: 69000\n",
      "evaluating user: 70000\n",
      "evaluating user: 71000\n",
      "evaluating user: 72000\n",
      "evaluating user: 73000\n",
      "evaluating user: 74000\n",
      "evaluating user: 75000\n",
      "evaluating user: 76000\n",
      "evaluating user: 77000\n",
      "evaluating user: 78000\n",
      "evaluating user: 79000\n",
      "evaluating user: 80000\n",
      "evaluating user: 81000\n",
      "evaluating user: 82000\n",
      "evaluating user: 83000\n",
      "evaluating user: 84000\n",
      "evaluating user: 85000\n",
      "evaluating user: 86000\n",
      "evaluating user: 87000\n",
      "evaluating user: 88000\n",
      "evaluating user: 89000\n",
      "evaluating user: 90000\n",
      "evaluating user: 91000\n",
      "evaluating user: 92000\n",
      "evaluating user: 93000\n",
      "evaluating user: 94000\n",
      "evaluating user: 95000\n",
      "evaluating user: 96000\n",
      "evaluating user: 97000\n",
      "evaluating user: 98000\n",
      "evaluating user: 99000\n",
      "evaluating user: 100000\n",
      "evaluating user: 101000\n",
      "evaluating user: 102000\n",
      "evaluating user: 103000\n",
      "evaluating user: 104000\n",
      "evaluating user: 105000\n",
      "evaluating user: 106000\n",
      "evaluating user: 107000\n",
      "evaluating user: 108000\n",
      "evaluating user: 109000\n",
      "evaluating user: 110000\n",
      "evaluating user: 111000\n",
      "evaluating user: 112000\n",
      "evaluating user: 113000\n",
      "evaluating user: 114000\n",
      "evaluating user: 115000\n",
      "evaluating user: 116000\n",
      "evaluating user: 117000\n",
      "evaluating user: 118000\n",
      "evaluating user: 119000\n",
      "evaluating user: 120000\n",
      "evaluating user: 121000\n",
      "evaluating user: 122000\n",
      "evaluating user: 123000\n",
      "evaluating user: 124000\n",
      "iter: 1 \tlam_diff:  920300.4925\tpsi_diff:  2.5167\n",
      "evaluating user: 0\n",
      "evaluating user: 1000\n",
      "evaluating user: 2000\n",
      "evaluating user: 3000\n",
      "evaluating user: 4000\n",
      "evaluating user: 5000\n",
      "evaluating user: 6000\n",
      "evaluating user: 7000\n",
      "evaluating user: 8000\n",
      "evaluating user: 9000\n",
      "evaluating user: 10000\n",
      "evaluating user: 11000\n",
      "evaluating user: 12000\n",
      "evaluating user: 13000\n",
      "evaluating user: 14000\n",
      "evaluating user: 15000\n",
      "evaluating user: 16000\n",
      "evaluating user: 17000\n",
      "evaluating user: 18000\n",
      "evaluating user: 19000\n",
      "evaluating user: 20000\n",
      "evaluating user: 21000\n",
      "evaluating user: 22000\n",
      "evaluating user: 23000\n",
      "evaluating user: 24000\n",
      "evaluating user: 25000\n",
      "evaluating user: 26000\n",
      "evaluating user: 27000\n",
      "evaluating user: 28000\n",
      "evaluating user: 29000\n",
      "evaluating user: 30000\n",
      "evaluating user: 31000\n",
      "evaluating user: 32000\n",
      "evaluating user: 33000\n",
      "evaluating user: 34000\n",
      "evaluating user: 35000\n",
      "evaluating user: 36000\n",
      "evaluating user: 37000\n",
      "evaluating user: 38000\n",
      "evaluating user: 39000\n",
      "evaluating user: 40000\n",
      "evaluating user: 41000\n",
      "evaluating user: 42000\n",
      "evaluating user: 43000\n",
      "evaluating user: 44000\n",
      "evaluating user: 45000\n",
      "evaluating user: 46000\n",
      "evaluating user: 47000\n",
      "evaluating user: 48000\n",
      "evaluating user: 49000\n",
      "evaluating user: 50000\n",
      "evaluating user: 51000\n",
      "evaluating user: 52000\n",
      "evaluating user: 53000\n",
      "evaluating user: 54000\n",
      "evaluating user: 55000\n",
      "evaluating user: 56000\n",
      "evaluating user: 57000\n",
      "evaluating user: 58000\n",
      "evaluating user: 59000\n",
      "evaluating user: 60000\n",
      "evaluating user: 61000\n",
      "evaluating user: 62000\n",
      "evaluating user: 63000\n",
      "evaluating user: 64000\n",
      "evaluating user: 65000\n",
      "evaluating user: 66000\n",
      "evaluating user: 67000\n",
      "evaluating user: 68000\n",
      "evaluating user: 69000\n",
      "evaluating user: 70000\n",
      "evaluating user: 71000\n",
      "evaluating user: 72000\n",
      "evaluating user: 73000\n",
      "evaluating user: 74000\n",
      "evaluating user: 75000\n",
      "evaluating user: 76000\n",
      "evaluating user: 77000\n",
      "evaluating user: 78000\n",
      "evaluating user: 79000\n",
      "evaluating user: 80000\n",
      "evaluating user: 81000\n",
      "evaluating user: 82000\n",
      "evaluating user: 83000\n",
      "evaluating user: 84000\n",
      "evaluating user: 85000\n",
      "evaluating user: 86000\n",
      "evaluating user: 87000\n",
      "evaluating user: 88000\n",
      "evaluating user: 89000\n",
      "evaluating user: 90000\n",
      "evaluating user: 91000\n",
      "evaluating user: 92000\n",
      "evaluating user: 93000\n",
      "evaluating user: 94000\n",
      "evaluating user: 95000\n",
      "evaluating user: 96000\n",
      "evaluating user: 97000\n",
      "evaluating user: 98000\n",
      "evaluating user: 99000\n",
      "evaluating user: 100000\n",
      "evaluating user: 101000\n",
      "evaluating user: 102000\n",
      "evaluating user: 103000\n",
      "evaluating user: 104000\n",
      "evaluating user: 105000\n",
      "evaluating user: 106000\n",
      "evaluating user: 107000\n",
      "evaluating user: 108000\n",
      "evaluating user: 109000\n",
      "evaluating user: 110000\n",
      "evaluating user: 111000\n",
      "evaluating user: 112000\n",
      "evaluating user: 113000\n",
      "evaluating user: 114000\n",
      "evaluating user: 115000\n",
      "evaluating user: 116000\n",
      "evaluating user: 117000\n",
      "evaluating user: 118000\n",
      "evaluating user: 119000\n",
      "evaluating user: 120000\n",
      "evaluating user: 121000\n",
      "evaluating user: 122000\n",
      "evaluating user: 123000\n",
      "evaluating user: 124000\n",
      "iter: 2 \tlam_diff:  499980.9452\tpsi_diff:  0.3594\n",
      "evaluating user: 0\n",
      "evaluating user: 1000\n",
      "evaluating user: 2000\n",
      "evaluating user: 3000\n",
      "evaluating user: 4000\n",
      "evaluating user: 5000\n",
      "evaluating user: 6000\n",
      "evaluating user: 7000\n",
      "evaluating user: 8000\n",
      "evaluating user: 9000\n",
      "evaluating user: 10000\n",
      "evaluating user: 11000\n",
      "evaluating user: 12000\n",
      "evaluating user: 13000\n",
      "evaluating user: 14000\n",
      "evaluating user: 15000\n",
      "evaluating user: 16000\n",
      "evaluating user: 17000\n",
      "evaluating user: 18000\n",
      "evaluating user: 19000\n",
      "evaluating user: 20000\n",
      "evaluating user: 21000\n",
      "evaluating user: 22000\n",
      "evaluating user: 23000\n",
      "evaluating user: 24000\n",
      "evaluating user: 25000\n",
      "evaluating user: 26000\n",
      "evaluating user: 27000\n",
      "evaluating user: 28000\n",
      "evaluating user: 29000\n",
      "evaluating user: 30000\n",
      "evaluating user: 31000\n",
      "evaluating user: 32000\n",
      "evaluating user: 33000\n",
      "evaluating user: 34000\n",
      "evaluating user: 35000\n",
      "evaluating user: 36000\n",
      "evaluating user: 37000\n",
      "evaluating user: 38000\n",
      "evaluating user: 39000\n",
      "evaluating user: 40000\n",
      "evaluating user: 41000\n",
      "evaluating user: 42000\n",
      "evaluating user: 43000\n",
      "evaluating user: 44000\n",
      "evaluating user: 45000\n",
      "evaluating user: 46000\n",
      "evaluating user: 47000\n",
      "evaluating user: 48000\n",
      "evaluating user: 49000\n",
      "evaluating user: 50000\n",
      "evaluating user: 51000\n",
      "evaluating user: 52000\n",
      "evaluating user: 53000\n",
      "evaluating user: 54000\n",
      "evaluating user: 55000\n",
      "evaluating user: 56000\n",
      "evaluating user: 57000\n",
      "evaluating user: 58000\n",
      "evaluating user: 59000\n",
      "evaluating user: 60000\n",
      "evaluating user: 61000\n",
      "evaluating user: 62000\n",
      "evaluating user: 63000\n",
      "evaluating user: 64000\n",
      "evaluating user: 65000\n",
      "evaluating user: 66000\n",
      "evaluating user: 67000\n",
      "evaluating user: 68000\n",
      "evaluating user: 69000\n",
      "evaluating user: 70000\n",
      "evaluating user: 71000\n",
      "evaluating user: 72000\n",
      "evaluating user: 73000\n",
      "evaluating user: 74000\n",
      "evaluating user: 75000\n",
      "evaluating user: 76000\n",
      "evaluating user: 77000\n",
      "evaluating user: 78000\n",
      "evaluating user: 79000\n",
      "evaluating user: 80000\n",
      "evaluating user: 81000\n",
      "evaluating user: 82000\n",
      "evaluating user: 83000\n",
      "evaluating user: 84000\n",
      "evaluating user: 85000\n",
      "evaluating user: 86000\n",
      "evaluating user: 87000\n",
      "evaluating user: 88000\n",
      "evaluating user: 89000\n",
      "evaluating user: 90000\n",
      "evaluating user: 91000\n",
      "evaluating user: 92000\n",
      "evaluating user: 93000\n",
      "evaluating user: 94000\n",
      "evaluating user: 95000\n",
      "evaluating user: 96000\n",
      "evaluating user: 97000\n",
      "evaluating user: 98000\n",
      "evaluating user: 99000\n",
      "evaluating user: 100000\n",
      "evaluating user: 101000\n",
      "evaluating user: 102000\n",
      "evaluating user: 103000\n",
      "evaluating user: 104000\n",
      "evaluating user: 105000\n",
      "evaluating user: 106000\n",
      "evaluating user: 107000\n",
      "evaluating user: 108000\n",
      "evaluating user: 109000\n",
      "evaluating user: 110000\n",
      "evaluating user: 111000\n",
      "evaluating user: 112000\n",
      "evaluating user: 113000\n",
      "evaluating user: 114000\n",
      "evaluating user: 115000\n",
      "evaluating user: 116000\n",
      "evaluating user: 117000\n",
      "evaluating user: 118000\n",
      "evaluating user: 119000\n",
      "evaluating user: 120000\n",
      "evaluating user: 121000\n",
      "evaluating user: 122000\n",
      "evaluating user: 123000\n",
      "evaluating user: 124000\n",
      "iter: 3 \tlam_diff:  170580.4056\tpsi_diff:  2.0028\n",
      "evaluating user: 0\n",
      "evaluating user: 1000\n",
      "evaluating user: 2000\n",
      "evaluating user: 3000\n",
      "evaluating user: 4000\n",
      "evaluating user: 5000\n",
      "evaluating user: 6000\n",
      "evaluating user: 7000\n",
      "evaluating user: 8000\n",
      "evaluating user: 9000\n",
      "evaluating user: 10000\n",
      "evaluating user: 11000\n",
      "evaluating user: 12000\n",
      "evaluating user: 13000\n",
      "evaluating user: 14000\n",
      "evaluating user: 15000\n",
      "evaluating user: 16000\n",
      "evaluating user: 17000\n",
      "evaluating user: 18000\n",
      "evaluating user: 19000\n",
      "evaluating user: 20000\n",
      "evaluating user: 21000\n",
      "evaluating user: 22000\n",
      "evaluating user: 23000\n",
      "evaluating user: 24000\n",
      "evaluating user: 25000\n",
      "evaluating user: 26000\n",
      "evaluating user: 27000\n",
      "evaluating user: 28000\n",
      "evaluating user: 29000\n",
      "evaluating user: 30000\n",
      "evaluating user: 31000\n",
      "evaluating user: 32000\n",
      "evaluating user: 33000\n",
      "evaluating user: 34000\n",
      "evaluating user: 35000\n",
      "evaluating user: 36000\n",
      "evaluating user: 37000\n",
      "evaluating user: 38000\n",
      "evaluating user: 39000\n",
      "evaluating user: 40000\n",
      "evaluating user: 41000\n",
      "evaluating user: 42000\n",
      "evaluating user: 43000\n",
      "evaluating user: 44000\n",
      "evaluating user: 45000\n",
      "evaluating user: 46000\n",
      "evaluating user: 47000\n",
      "evaluating user: 48000\n",
      "evaluating user: 49000\n",
      "evaluating user: 50000\n",
      "evaluating user: 51000\n",
      "evaluating user: 52000\n",
      "evaluating user: 53000\n",
      "evaluating user: 54000\n",
      "evaluating user: 55000\n",
      "evaluating user: 56000\n",
      "evaluating user: 57000\n",
      "evaluating user: 58000\n",
      "evaluating user: 59000\n",
      "evaluating user: 60000\n",
      "evaluating user: 61000\n",
      "evaluating user: 62000\n",
      "evaluating user: 63000\n",
      "evaluating user: 64000\n",
      "evaluating user: 65000\n",
      "evaluating user: 66000\n",
      "evaluating user: 67000\n",
      "evaluating user: 68000\n",
      "evaluating user: 69000\n",
      "evaluating user: 70000\n",
      "evaluating user: 71000\n",
      "evaluating user: 72000\n",
      "evaluating user: 73000\n",
      "evaluating user: 74000\n",
      "evaluating user: 75000\n",
      "evaluating user: 76000\n",
      "evaluating user: 77000\n",
      "evaluating user: 78000\n",
      "evaluating user: 79000\n",
      "evaluating user: 80000\n",
      "evaluating user: 81000\n",
      "evaluating user: 82000\n",
      "evaluating user: 83000\n",
      "evaluating user: 84000\n",
      "evaluating user: 85000\n",
      "evaluating user: 86000\n",
      "evaluating user: 87000\n",
      "evaluating user: 88000\n",
      "evaluating user: 89000\n",
      "evaluating user: 90000\n",
      "evaluating user: 91000\n",
      "evaluating user: 92000\n",
      "evaluating user: 93000\n",
      "evaluating user: 94000\n",
      "evaluating user: 95000\n",
      "evaluating user: 96000\n",
      "evaluating user: 97000\n",
      "evaluating user: 98000\n",
      "evaluating user: 99000\n",
      "evaluating user: 100000\n",
      "evaluating user: 101000\n",
      "evaluating user: 102000\n",
      "evaluating user: 103000\n",
      "evaluating user: 104000\n",
      "evaluating user: 105000\n",
      "evaluating user: 106000\n",
      "evaluating user: 107000\n",
      "evaluating user: 108000\n",
      "evaluating user: 109000\n",
      "evaluating user: 110000\n",
      "evaluating user: 111000\n",
      "evaluating user: 112000\n",
      "evaluating user: 113000\n",
      "evaluating user: 114000\n",
      "evaluating user: 115000\n",
      "evaluating user: 116000\n",
      "evaluating user: 117000\n",
      "evaluating user: 118000\n",
      "evaluating user: 119000\n",
      "evaluating user: 120000\n",
      "evaluating user: 121000\n",
      "evaluating user: 122000\n",
      "evaluating user: 123000\n",
      "evaluating user: 124000\n",
      "iter: 4 \tlam_diff:  21527.8918\tpsi_diff:  2.5420\n",
      "evaluating user: 0\n",
      "evaluating user: 1000\n",
      "evaluating user: 2000\n",
      "evaluating user: 3000\n",
      "evaluating user: 4000\n",
      "evaluating user: 5000\n",
      "evaluating user: 6000\n",
      "evaluating user: 7000\n",
      "evaluating user: 8000\n",
      "evaluating user: 9000\n",
      "evaluating user: 10000\n",
      "evaluating user: 11000\n",
      "evaluating user: 12000\n",
      "evaluating user: 13000\n",
      "evaluating user: 14000\n",
      "evaluating user: 15000\n",
      "evaluating user: 16000\n",
      "evaluating user: 17000\n",
      "evaluating user: 18000\n",
      "evaluating user: 19000\n",
      "evaluating user: 20000\n",
      "evaluating user: 21000\n",
      "evaluating user: 22000\n",
      "evaluating user: 23000\n",
      "evaluating user: 24000\n",
      "evaluating user: 25000\n",
      "evaluating user: 26000\n",
      "evaluating user: 27000\n",
      "evaluating user: 28000\n",
      "evaluating user: 29000\n",
      "evaluating user: 30000\n",
      "evaluating user: 31000\n",
      "evaluating user: 32000\n",
      "evaluating user: 33000\n",
      "evaluating user: 34000\n",
      "evaluating user: 35000\n",
      "evaluating user: 36000\n",
      "evaluating user: 37000\n",
      "evaluating user: 38000\n",
      "evaluating user: 39000\n",
      "evaluating user: 40000\n",
      "evaluating user: 41000\n",
      "evaluating user: 42000\n",
      "evaluating user: 43000\n",
      "evaluating user: 44000\n",
      "evaluating user: 45000\n",
      "evaluating user: 46000\n",
      "evaluating user: 47000\n",
      "evaluating user: 48000\n",
      "evaluating user: 49000\n",
      "evaluating user: 50000\n",
      "evaluating user: 51000\n",
      "evaluating user: 52000\n",
      "evaluating user: 53000\n",
      "evaluating user: 54000\n",
      "evaluating user: 55000\n",
      "evaluating user: 56000\n",
      "evaluating user: 57000\n",
      "evaluating user: 58000\n",
      "evaluating user: 59000\n",
      "evaluating user: 60000\n",
      "evaluating user: 61000\n",
      "evaluating user: 62000\n",
      "evaluating user: 63000\n",
      "evaluating user: 64000\n",
      "evaluating user: 65000\n",
      "evaluating user: 66000\n",
      "evaluating user: 67000\n",
      "evaluating user: 68000\n",
      "evaluating user: 69000\n",
      "evaluating user: 70000\n",
      "evaluating user: 71000\n",
      "evaluating user: 72000\n",
      "evaluating user: 73000\n",
      "evaluating user: 74000\n",
      "evaluating user: 75000\n",
      "evaluating user: 76000\n",
      "evaluating user: 77000\n",
      "evaluating user: 78000\n",
      "evaluating user: 79000\n",
      "evaluating user: 80000\n",
      "evaluating user: 81000\n",
      "evaluating user: 82000\n",
      "evaluating user: 83000\n",
      "evaluating user: 84000\n",
      "evaluating user: 85000\n",
      "evaluating user: 86000\n",
      "evaluating user: 87000\n",
      "evaluating user: 88000\n",
      "evaluating user: 89000\n",
      "evaluating user: 90000\n",
      "evaluating user: 91000\n",
      "evaluating user: 92000\n",
      "evaluating user: 93000\n",
      "evaluating user: 94000\n",
      "evaluating user: 95000\n",
      "evaluating user: 96000\n",
      "evaluating user: 97000\n",
      "evaluating user: 98000\n",
      "evaluating user: 99000\n",
      "evaluating user: 100000\n",
      "evaluating user: 101000\n",
      "evaluating user: 102000\n",
      "evaluating user: 103000\n",
      "evaluating user: 104000\n",
      "evaluating user: 105000\n",
      "evaluating user: 106000\n",
      "evaluating user: 107000\n",
      "evaluating user: 108000\n",
      "evaluating user: 109000\n",
      "evaluating user: 110000\n",
      "evaluating user: 111000\n",
      "evaluating user: 112000\n",
      "evaluating user: 113000\n",
      "evaluating user: 114000\n",
      "evaluating user: 115000\n",
      "evaluating user: 116000\n",
      "evaluating user: 117000\n",
      "evaluating user: 118000\n",
      "evaluating user: 119000\n",
      "evaluating user: 120000\n",
      "evaluating user: 121000\n",
      "evaluating user: 122000\n",
      "evaluating user: 123000\n",
      "evaluating user: 124000\n",
      "iter: 5 \tlam_diff:  2747.7937\tpsi_diff:  1.9919\n",
      "evaluating user: 0\n",
      "evaluating user: 1000\n",
      "evaluating user: 2000\n",
      "evaluating user: 3000\n",
      "evaluating user: 4000\n",
      "evaluating user: 5000\n",
      "evaluating user: 6000\n",
      "evaluating user: 7000\n",
      "evaluating user: 8000\n",
      "evaluating user: 9000\n",
      "evaluating user: 10000\n",
      "evaluating user: 11000\n",
      "evaluating user: 12000\n",
      "evaluating user: 13000\n",
      "evaluating user: 14000\n",
      "evaluating user: 15000\n",
      "evaluating user: 16000\n",
      "evaluating user: 17000\n",
      "evaluating user: 18000\n",
      "evaluating user: 19000\n",
      "evaluating user: 20000\n",
      "evaluating user: 21000\n",
      "evaluating user: 22000\n",
      "evaluating user: 23000\n",
      "evaluating user: 24000\n",
      "evaluating user: 25000\n",
      "evaluating user: 26000\n",
      "evaluating user: 27000\n",
      "evaluating user: 28000\n",
      "evaluating user: 29000\n",
      "evaluating user: 30000\n",
      "evaluating user: 31000\n",
      "evaluating user: 32000\n",
      "evaluating user: 33000\n",
      "evaluating user: 34000\n",
      "evaluating user: 35000\n",
      "evaluating user: 36000\n",
      "evaluating user: 37000\n",
      "evaluating user: 38000\n",
      "evaluating user: 39000\n",
      "evaluating user: 40000\n",
      "evaluating user: 41000\n",
      "evaluating user: 42000\n",
      "evaluating user: 43000\n",
      "evaluating user: 44000\n",
      "evaluating user: 45000\n",
      "evaluating user: 46000\n",
      "evaluating user: 47000\n",
      "evaluating user: 48000\n",
      "evaluating user: 49000\n",
      "evaluating user: 50000\n",
      "evaluating user: 51000\n",
      "evaluating user: 52000\n",
      "evaluating user: 53000\n",
      "evaluating user: 54000\n",
      "evaluating user: 55000\n",
      "evaluating user: 56000\n",
      "evaluating user: 57000\n",
      "evaluating user: 58000\n",
      "evaluating user: 59000\n",
      "evaluating user: 60000\n",
      "evaluating user: 61000\n",
      "evaluating user: 62000\n",
      "evaluating user: 63000\n",
      "evaluating user: 64000\n",
      "evaluating user: 65000\n",
      "evaluating user: 66000\n",
      "evaluating user: 67000\n",
      "evaluating user: 68000\n",
      "evaluating user: 69000\n",
      "evaluating user: 70000\n",
      "evaluating user: 71000\n",
      "evaluating user: 72000\n",
      "evaluating user: 73000\n",
      "evaluating user: 74000\n",
      "evaluating user: 75000\n",
      "evaluating user: 76000\n",
      "evaluating user: 77000\n",
      "evaluating user: 78000\n",
      "evaluating user: 79000\n",
      "evaluating user: 80000\n",
      "evaluating user: 81000\n",
      "evaluating user: 82000\n",
      "evaluating user: 83000\n",
      "evaluating user: 84000\n",
      "evaluating user: 85000\n",
      "evaluating user: 86000\n",
      "evaluating user: 87000\n",
      "evaluating user: 88000\n",
      "evaluating user: 89000\n",
      "evaluating user: 90000\n",
      "evaluating user: 91000\n",
      "evaluating user: 92000\n",
      "evaluating user: 93000\n",
      "evaluating user: 94000\n",
      "evaluating user: 95000\n",
      "evaluating user: 96000\n",
      "evaluating user: 97000\n",
      "evaluating user: 98000\n",
      "evaluating user: 99000\n",
      "evaluating user: 100000\n",
      "evaluating user: 101000\n",
      "evaluating user: 102000\n",
      "evaluating user: 103000\n",
      "evaluating user: 104000\n",
      "evaluating user: 105000\n",
      "evaluating user: 106000\n",
      "evaluating user: 107000\n",
      "evaluating user: 108000\n",
      "evaluating user: 109000\n",
      "evaluating user: 110000\n",
      "evaluating user: 111000\n",
      "evaluating user: 112000\n",
      "evaluating user: 113000\n",
      "evaluating user: 114000\n",
      "evaluating user: 115000\n",
      "evaluating user: 116000\n",
      "evaluating user: 117000\n",
      "evaluating user: 118000\n",
      "evaluating user: 119000\n",
      "evaluating user: 120000\n",
      "evaluating user: 121000\n",
      "evaluating user: 122000\n",
      "evaluating user: 123000\n",
      "evaluating user: 124000\n",
      "iter: 6 \tlam_diff:  518.5323\tpsi_diff:  1.4241\n",
      "evaluating user: 0\n",
      "evaluating user: 1000\n",
      "evaluating user: 2000\n",
      "evaluating user: 3000\n",
      "evaluating user: 4000\n",
      "evaluating user: 5000\n",
      "evaluating user: 6000\n",
      "evaluating user: 7000\n",
      "evaluating user: 8000\n",
      "evaluating user: 9000\n",
      "evaluating user: 10000\n",
      "evaluating user: 11000\n",
      "evaluating user: 12000\n",
      "evaluating user: 13000\n",
      "evaluating user: 14000\n",
      "evaluating user: 15000\n",
      "evaluating user: 16000\n",
      "evaluating user: 17000\n",
      "evaluating user: 18000\n",
      "evaluating user: 19000\n",
      "evaluating user: 20000\n",
      "evaluating user: 21000\n",
      "evaluating user: 22000\n",
      "evaluating user: 23000\n",
      "evaluating user: 24000\n",
      "evaluating user: 25000\n",
      "evaluating user: 26000\n",
      "evaluating user: 27000\n",
      "evaluating user: 28000\n",
      "evaluating user: 29000\n",
      "evaluating user: 30000\n",
      "evaluating user: 31000\n",
      "evaluating user: 32000\n",
      "evaluating user: 33000\n",
      "evaluating user: 34000\n",
      "evaluating user: 35000\n",
      "evaluating user: 36000\n",
      "evaluating user: 37000\n",
      "evaluating user: 38000\n",
      "evaluating user: 39000\n",
      "evaluating user: 40000\n",
      "evaluating user: 41000\n",
      "evaluating user: 42000\n",
      "evaluating user: 43000\n",
      "evaluating user: 44000\n",
      "evaluating user: 45000\n",
      "evaluating user: 46000\n",
      "evaluating user: 47000\n",
      "evaluating user: 48000\n",
      "evaluating user: 49000\n",
      "evaluating user: 50000\n",
      "evaluating user: 51000\n",
      "evaluating user: 52000\n",
      "evaluating user: 53000\n",
      "evaluating user: 54000\n",
      "evaluating user: 55000\n",
      "evaluating user: 56000\n",
      "evaluating user: 57000\n",
      "evaluating user: 58000\n",
      "evaluating user: 59000\n",
      "evaluating user: 60000\n",
      "evaluating user: 61000\n",
      "evaluating user: 62000\n",
      "evaluating user: 63000\n",
      "evaluating user: 64000\n",
      "evaluating user: 65000\n",
      "evaluating user: 66000\n",
      "evaluating user: 67000\n",
      "evaluating user: 68000\n",
      "evaluating user: 69000\n",
      "evaluating user: 70000\n",
      "evaluating user: 71000\n",
      "evaluating user: 72000\n",
      "evaluating user: 73000\n",
      "evaluating user: 74000\n",
      "evaluating user: 75000\n",
      "evaluating user: 76000\n",
      "evaluating user: 77000\n",
      "evaluating user: 78000\n",
      "evaluating user: 79000\n",
      "evaluating user: 80000\n",
      "evaluating user: 81000\n",
      "evaluating user: 82000\n",
      "evaluating user: 83000\n",
      "evaluating user: 84000\n",
      "evaluating user: 85000\n",
      "evaluating user: 86000\n",
      "evaluating user: 87000\n",
      "evaluating user: 88000\n",
      "evaluating user: 89000\n",
      "evaluating user: 90000\n",
      "evaluating user: 91000\n",
      "evaluating user: 92000\n",
      "evaluating user: 93000\n",
      "evaluating user: 94000\n",
      "evaluating user: 95000\n",
      "evaluating user: 96000\n",
      "evaluating user: 97000\n",
      "evaluating user: 98000\n",
      "evaluating user: 99000\n",
      "evaluating user: 100000\n",
      "evaluating user: 101000\n",
      "evaluating user: 102000\n",
      "evaluating user: 103000\n",
      "evaluating user: 104000\n",
      "evaluating user: 105000\n",
      "evaluating user: 106000\n",
      "evaluating user: 107000\n",
      "evaluating user: 108000\n",
      "evaluating user: 109000\n",
      "evaluating user: 110000\n",
      "evaluating user: 111000\n",
      "evaluating user: 112000\n",
      "evaluating user: 113000\n",
      "evaluating user: 114000\n",
      "evaluating user: 115000\n",
      "evaluating user: 116000\n",
      "evaluating user: 117000\n",
      "evaluating user: 118000\n",
      "evaluating user: 119000\n",
      "evaluating user: 120000\n",
      "evaluating user: 121000\n",
      "evaluating user: 122000\n",
      "evaluating user: 123000\n",
      "evaluating user: 124000\n",
      "iter: 7 \tlam_diff:  174.2610\tpsi_diff:  0.9921\n",
      "evaluating user: 0\n",
      "evaluating user: 1000\n",
      "evaluating user: 2000\n",
      "evaluating user: 3000\n",
      "evaluating user: 4000\n",
      "evaluating user: 5000\n",
      "evaluating user: 6000\n",
      "evaluating user: 7000\n",
      "evaluating user: 8000\n",
      "evaluating user: 9000\n",
      "evaluating user: 10000\n",
      "evaluating user: 11000\n",
      "evaluating user: 12000\n",
      "evaluating user: 13000\n",
      "evaluating user: 14000\n",
      "evaluating user: 15000\n",
      "evaluating user: 16000\n",
      "evaluating user: 17000\n",
      "evaluating user: 18000\n",
      "evaluating user: 19000\n",
      "evaluating user: 20000\n",
      "evaluating user: 21000\n",
      "evaluating user: 22000\n",
      "evaluating user: 23000\n",
      "evaluating user: 24000\n",
      "evaluating user: 25000\n",
      "evaluating user: 26000\n",
      "evaluating user: 27000\n",
      "evaluating user: 28000\n",
      "evaluating user: 29000\n",
      "evaluating user: 30000\n",
      "evaluating user: 31000\n",
      "evaluating user: 32000\n",
      "evaluating user: 33000\n",
      "evaluating user: 34000\n",
      "evaluating user: 35000\n",
      "evaluating user: 36000\n",
      "evaluating user: 37000\n",
      "evaluating user: 38000\n",
      "evaluating user: 39000\n",
      "evaluating user: 40000\n",
      "evaluating user: 41000\n",
      "evaluating user: 42000\n",
      "evaluating user: 43000\n",
      "evaluating user: 44000\n",
      "evaluating user: 45000\n",
      "evaluating user: 46000\n",
      "evaluating user: 47000\n",
      "evaluating user: 48000\n",
      "evaluating user: 49000\n",
      "evaluating user: 50000\n",
      "evaluating user: 51000\n",
      "evaluating user: 52000\n",
      "evaluating user: 53000\n",
      "evaluating user: 54000\n",
      "evaluating user: 55000\n",
      "evaluating user: 56000\n",
      "evaluating user: 57000\n",
      "evaluating user: 58000\n",
      "evaluating user: 59000\n",
      "evaluating user: 60000\n",
      "evaluating user: 61000\n",
      "evaluating user: 62000\n",
      "evaluating user: 63000\n",
      "evaluating user: 64000\n",
      "evaluating user: 65000\n",
      "evaluating user: 66000\n",
      "evaluating user: 67000\n",
      "evaluating user: 68000\n",
      "evaluating user: 69000\n",
      "evaluating user: 70000\n",
      "evaluating user: 71000\n",
      "evaluating user: 72000\n",
      "evaluating user: 73000\n",
      "evaluating user: 74000\n",
      "evaluating user: 75000\n",
      "evaluating user: 76000\n",
      "evaluating user: 77000\n",
      "evaluating user: 78000\n",
      "evaluating user: 79000\n",
      "evaluating user: 80000\n",
      "evaluating user: 81000\n",
      "evaluating user: 82000\n",
      "evaluating user: 83000\n",
      "evaluating user: 84000\n",
      "evaluating user: 85000\n",
      "evaluating user: 86000\n",
      "evaluating user: 87000\n",
      "evaluating user: 88000\n",
      "evaluating user: 89000\n",
      "evaluating user: 90000\n",
      "evaluating user: 91000\n",
      "evaluating user: 92000\n",
      "evaluating user: 93000\n",
      "evaluating user: 94000\n",
      "evaluating user: 95000\n",
      "evaluating user: 96000\n",
      "evaluating user: 97000\n",
      "evaluating user: 98000\n",
      "evaluating user: 99000\n",
      "evaluating user: 100000\n",
      "evaluating user: 101000\n",
      "evaluating user: 102000\n",
      "evaluating user: 103000\n",
      "evaluating user: 104000\n",
      "evaluating user: 105000\n",
      "evaluating user: 106000\n",
      "evaluating user: 107000\n",
      "evaluating user: 108000\n",
      "evaluating user: 109000\n",
      "evaluating user: 110000\n",
      "evaluating user: 111000\n",
      "evaluating user: 112000\n",
      "evaluating user: 113000\n",
      "evaluating user: 114000\n",
      "evaluating user: 115000\n",
      "evaluating user: 116000\n",
      "evaluating user: 117000\n",
      "evaluating user: 118000\n",
      "evaluating user: 119000\n",
      "evaluating user: 120000\n",
      "evaluating user: 121000\n",
      "evaluating user: 122000\n",
      "evaluating user: 123000\n",
      "evaluating user: 124000\n",
      "iter: 8 \tlam_diff:  84.4125\tpsi_diff:  0.6773\n",
      "evaluating user: 0\n",
      "evaluating user: 1000\n",
      "evaluating user: 2000\n",
      "evaluating user: 3000\n",
      "evaluating user: 4000\n",
      "evaluating user: 5000\n",
      "evaluating user: 6000\n",
      "evaluating user: 7000\n",
      "evaluating user: 8000\n",
      "evaluating user: 9000\n",
      "evaluating user: 10000\n",
      "evaluating user: 11000\n",
      "evaluating user: 12000\n",
      "evaluating user: 13000\n",
      "evaluating user: 14000\n",
      "evaluating user: 15000\n",
      "evaluating user: 16000\n",
      "evaluating user: 17000\n",
      "evaluating user: 18000\n",
      "evaluating user: 19000\n",
      "evaluating user: 20000\n",
      "evaluating user: 21000\n",
      "evaluating user: 22000\n",
      "evaluating user: 23000\n",
      "evaluating user: 24000\n",
      "evaluating user: 25000\n",
      "evaluating user: 26000\n",
      "evaluating user: 27000\n",
      "evaluating user: 28000\n",
      "evaluating user: 29000\n",
      "evaluating user: 30000\n",
      "evaluating user: 31000\n",
      "evaluating user: 32000\n",
      "evaluating user: 33000\n",
      "evaluating user: 34000\n",
      "evaluating user: 35000\n",
      "evaluating user: 36000\n",
      "evaluating user: 37000\n",
      "evaluating user: 38000\n",
      "evaluating user: 39000\n",
      "evaluating user: 40000\n",
      "evaluating user: 41000\n",
      "evaluating user: 42000\n",
      "evaluating user: 43000\n",
      "evaluating user: 44000\n",
      "evaluating user: 45000\n",
      "evaluating user: 46000\n",
      "evaluating user: 47000\n",
      "evaluating user: 48000\n",
      "evaluating user: 49000\n",
      "evaluating user: 50000\n",
      "evaluating user: 51000\n",
      "evaluating user: 52000\n",
      "evaluating user: 53000\n",
      "evaluating user: 54000\n",
      "evaluating user: 55000\n",
      "evaluating user: 56000\n",
      "evaluating user: 57000\n",
      "evaluating user: 58000\n",
      "evaluating user: 59000\n",
      "evaluating user: 60000\n",
      "evaluating user: 61000\n",
      "evaluating user: 62000\n",
      "evaluating user: 63000\n",
      "evaluating user: 64000\n",
      "evaluating user: 65000\n",
      "evaluating user: 66000\n",
      "evaluating user: 67000\n",
      "evaluating user: 68000\n",
      "evaluating user: 69000\n",
      "evaluating user: 70000\n",
      "evaluating user: 71000\n",
      "evaluating user: 72000\n",
      "evaluating user: 73000\n",
      "evaluating user: 74000\n",
      "evaluating user: 75000\n",
      "evaluating user: 76000\n",
      "evaluating user: 77000\n",
      "evaluating user: 78000\n",
      "evaluating user: 79000\n",
      "evaluating user: 80000\n",
      "evaluating user: 81000\n",
      "evaluating user: 82000\n",
      "evaluating user: 83000\n",
      "evaluating user: 84000\n",
      "evaluating user: 85000\n",
      "evaluating user: 86000\n",
      "evaluating user: 87000\n",
      "evaluating user: 88000\n",
      "evaluating user: 89000\n",
      "evaluating user: 90000\n",
      "evaluating user: 91000\n",
      "evaluating user: 92000\n",
      "evaluating user: 93000\n",
      "evaluating user: 94000\n",
      "evaluating user: 95000\n",
      "evaluating user: 96000\n",
      "evaluating user: 97000\n",
      "evaluating user: 98000\n",
      "evaluating user: 99000\n",
      "evaluating user: 100000\n",
      "evaluating user: 101000\n",
      "evaluating user: 102000\n",
      "evaluating user: 103000\n",
      "evaluating user: 104000\n",
      "evaluating user: 105000\n",
      "evaluating user: 106000\n",
      "evaluating user: 107000\n",
      "evaluating user: 108000\n",
      "evaluating user: 109000\n",
      "evaluating user: 110000\n",
      "evaluating user: 111000\n",
      "evaluating user: 112000\n",
      "evaluating user: 113000\n",
      "evaluating user: 114000\n",
      "evaluating user: 115000\n",
      "evaluating user: 116000\n",
      "evaluating user: 117000\n",
      "evaluating user: 118000\n",
      "evaluating user: 119000\n",
      "evaluating user: 120000\n",
      "evaluating user: 121000\n",
      "evaluating user: 122000\n",
      "evaluating user: 123000\n",
      "evaluating user: 124000\n",
      "iter: 9 \tlam_diff:  48.1745\tpsi_diff:  0.4550\n",
      "evaluating user: 0\n",
      "evaluating user: 1000\n",
      "evaluating user: 2000\n",
      "evaluating user: 3000\n",
      "evaluating user: 4000\n",
      "evaluating user: 5000\n",
      "evaluating user: 6000\n",
      "evaluating user: 7000\n",
      "evaluating user: 8000\n",
      "evaluating user: 9000\n",
      "evaluating user: 10000\n",
      "evaluating user: 11000\n",
      "evaluating user: 12000\n",
      "evaluating user: 13000\n",
      "evaluating user: 14000\n",
      "evaluating user: 15000\n",
      "evaluating user: 16000\n",
      "evaluating user: 17000\n",
      "evaluating user: 18000\n",
      "evaluating user: 19000\n",
      "evaluating user: 20000\n",
      "evaluating user: 21000\n",
      "evaluating user: 22000\n",
      "evaluating user: 23000\n",
      "evaluating user: 24000\n",
      "evaluating user: 25000\n",
      "evaluating user: 26000\n",
      "evaluating user: 27000\n",
      "evaluating user: 28000\n",
      "evaluating user: 29000\n",
      "evaluating user: 30000\n",
      "evaluating user: 31000\n",
      "evaluating user: 32000\n",
      "evaluating user: 33000\n",
      "evaluating user: 34000\n",
      "evaluating user: 35000\n",
      "evaluating user: 36000\n",
      "evaluating user: 37000\n",
      "evaluating user: 38000\n",
      "evaluating user: 39000\n",
      "evaluating user: 40000\n",
      "evaluating user: 41000\n",
      "evaluating user: 42000\n",
      "evaluating user: 43000\n",
      "evaluating user: 44000\n",
      "evaluating user: 45000\n",
      "evaluating user: 46000\n",
      "evaluating user: 47000\n",
      "evaluating user: 48000\n",
      "evaluating user: 49000\n",
      "evaluating user: 50000\n",
      "evaluating user: 51000\n",
      "evaluating user: 52000\n",
      "evaluating user: 53000\n",
      "evaluating user: 54000\n",
      "evaluating user: 55000\n",
      "evaluating user: 56000\n",
      "evaluating user: 57000\n",
      "evaluating user: 58000\n",
      "evaluating user: 59000\n",
      "evaluating user: 60000\n",
      "evaluating user: 61000\n",
      "evaluating user: 62000\n",
      "evaluating user: 63000\n",
      "evaluating user: 64000\n",
      "evaluating user: 65000\n",
      "evaluating user: 66000\n",
      "evaluating user: 67000\n",
      "evaluating user: 68000\n",
      "evaluating user: 69000\n",
      "evaluating user: 70000\n",
      "evaluating user: 71000\n",
      "evaluating user: 72000\n",
      "evaluating user: 73000\n",
      "evaluating user: 74000\n",
      "evaluating user: 75000\n",
      "evaluating user: 76000\n",
      "evaluating user: 77000\n",
      "evaluating user: 78000\n",
      "evaluating user: 79000\n",
      "evaluating user: 80000\n",
      "evaluating user: 81000\n",
      "evaluating user: 82000\n",
      "evaluating user: 83000\n",
      "evaluating user: 84000\n",
      "evaluating user: 85000\n",
      "evaluating user: 86000\n",
      "evaluating user: 87000\n",
      "evaluating user: 88000\n",
      "evaluating user: 89000\n",
      "evaluating user: 90000\n",
      "evaluating user: 91000\n",
      "evaluating user: 92000\n",
      "evaluating user: 93000\n",
      "evaluating user: 94000\n",
      "evaluating user: 95000\n",
      "evaluating user: 96000\n",
      "evaluating user: 97000\n",
      "evaluating user: 98000\n",
      "evaluating user: 99000\n",
      "evaluating user: 100000\n",
      "evaluating user: 101000\n",
      "evaluating user: 102000\n",
      "evaluating user: 103000\n",
      "evaluating user: 104000\n",
      "evaluating user: 105000\n",
      "evaluating user: 106000\n",
      "evaluating user: 107000\n",
      "evaluating user: 108000\n",
      "evaluating user: 109000\n",
      "evaluating user: 110000\n",
      "evaluating user: 111000\n",
      "evaluating user: 112000\n",
      "evaluating user: 113000\n",
      "evaluating user: 114000\n",
      "evaluating user: 115000\n",
      "evaluating user: 116000\n",
      "evaluating user: 117000\n",
      "evaluating user: 118000\n",
      "evaluating user: 119000\n",
      "evaluating user: 120000\n",
      "evaluating user: 121000\n",
      "evaluating user: 122000\n",
      "evaluating user: 123000\n",
      "evaluating user: 124000\n",
      "iter: 10 \tlam_diff:  31.6967\tpsi_diff:  0.3030\n",
      "evaluating user: 0\n",
      "evaluating user: 1000\n",
      "evaluating user: 2000\n",
      "evaluating user: 3000\n",
      "evaluating user: 4000\n",
      "evaluating user: 5000\n",
      "evaluating user: 6000\n",
      "evaluating user: 7000\n",
      "evaluating user: 8000\n",
      "evaluating user: 9000\n",
      "evaluating user: 10000\n",
      "evaluating user: 11000\n",
      "evaluating user: 12000\n",
      "evaluating user: 13000\n",
      "evaluating user: 14000\n",
      "evaluating user: 15000\n",
      "evaluating user: 16000\n",
      "evaluating user: 17000\n",
      "evaluating user: 18000\n",
      "evaluating user: 19000\n",
      "evaluating user: 20000\n",
      "evaluating user: 21000\n",
      "evaluating user: 22000\n",
      "evaluating user: 23000\n",
      "evaluating user: 24000\n",
      "evaluating user: 25000\n",
      "evaluating user: 26000\n",
      "evaluating user: 27000\n",
      "evaluating user: 28000\n",
      "evaluating user: 29000\n",
      "evaluating user: 30000\n",
      "evaluating user: 31000\n",
      "evaluating user: 32000\n",
      "evaluating user: 33000\n",
      "evaluating user: 34000\n",
      "evaluating user: 35000\n",
      "evaluating user: 36000\n",
      "evaluating user: 37000\n",
      "evaluating user: 38000\n",
      "evaluating user: 39000\n",
      "evaluating user: 40000\n",
      "evaluating user: 41000\n",
      "evaluating user: 42000\n",
      "evaluating user: 43000\n",
      "evaluating user: 44000\n",
      "evaluating user: 45000\n",
      "evaluating user: 46000\n",
      "evaluating user: 47000\n",
      "evaluating user: 48000\n",
      "evaluating user: 49000\n",
      "evaluating user: 50000\n",
      "evaluating user: 51000\n",
      "evaluating user: 52000\n",
      "evaluating user: 53000\n",
      "evaluating user: 54000\n",
      "evaluating user: 55000\n",
      "evaluating user: 56000\n",
      "evaluating user: 57000\n",
      "evaluating user: 58000\n",
      "evaluating user: 59000\n",
      "evaluating user: 60000\n",
      "evaluating user: 61000\n",
      "evaluating user: 62000\n",
      "evaluating user: 63000\n",
      "evaluating user: 64000\n",
      "evaluating user: 65000\n",
      "evaluating user: 66000\n",
      "evaluating user: 67000\n",
      "evaluating user: 68000\n",
      "evaluating user: 69000\n",
      "evaluating user: 70000\n",
      "evaluating user: 71000\n",
      "evaluating user: 72000\n",
      "evaluating user: 73000\n",
      "evaluating user: 74000\n",
      "evaluating user: 75000\n",
      "evaluating user: 76000\n",
      "evaluating user: 77000\n",
      "evaluating user: 78000\n",
      "evaluating user: 79000\n",
      "evaluating user: 80000\n",
      "evaluating user: 81000\n",
      "evaluating user: 82000\n",
      "evaluating user: 83000\n",
      "evaluating user: 84000\n",
      "evaluating user: 85000\n",
      "evaluating user: 86000\n",
      "evaluating user: 87000\n",
      "evaluating user: 88000\n",
      "evaluating user: 89000\n",
      "evaluating user: 90000\n",
      "evaluating user: 91000\n",
      "evaluating user: 92000\n",
      "evaluating user: 93000\n",
      "evaluating user: 94000\n",
      "evaluating user: 95000\n",
      "evaluating user: 96000\n",
      "evaluating user: 97000\n",
      "evaluating user: 98000\n",
      "evaluating user: 99000\n",
      "evaluating user: 100000\n",
      "evaluating user: 101000\n",
      "evaluating user: 102000\n",
      "evaluating user: 103000\n",
      "evaluating user: 104000\n",
      "evaluating user: 105000\n",
      "evaluating user: 106000\n",
      "evaluating user: 107000\n",
      "evaluating user: 108000\n",
      "evaluating user: 109000\n",
      "evaluating user: 110000\n",
      "evaluating user: 111000\n",
      "evaluating user: 112000\n",
      "evaluating user: 113000\n",
      "evaluating user: 114000\n",
      "evaluating user: 115000\n",
      "evaluating user: 116000\n",
      "evaluating user: 117000\n",
      "evaluating user: 118000\n",
      "evaluating user: 119000\n",
      "evaluating user: 120000\n",
      "evaluating user: 121000\n",
      "evaluating user: 122000\n",
      "evaluating user: 123000\n",
      "evaluating user: 124000\n",
      "iter: 11 \tlam_diff:  23.7216\tpsi_diff:  0.2021\n",
      "evaluating user: 0\n",
      "evaluating user: 1000\n",
      "evaluating user: 2000\n",
      "evaluating user: 3000\n",
      "evaluating user: 4000\n",
      "evaluating user: 5000\n",
      "evaluating user: 6000\n",
      "evaluating user: 7000\n",
      "evaluating user: 8000\n",
      "evaluating user: 9000\n",
      "evaluating user: 10000\n",
      "evaluating user: 11000\n",
      "evaluating user: 12000\n",
      "evaluating user: 13000\n",
      "evaluating user: 14000\n",
      "evaluating user: 15000\n",
      "evaluating user: 16000\n",
      "evaluating user: 17000\n",
      "evaluating user: 18000\n",
      "evaluating user: 19000\n",
      "evaluating user: 20000\n",
      "evaluating user: 21000\n",
      "evaluating user: 22000\n",
      "evaluating user: 23000\n",
      "evaluating user: 24000\n",
      "evaluating user: 25000\n",
      "evaluating user: 26000\n",
      "evaluating user: 27000\n",
      "evaluating user: 28000\n",
      "evaluating user: 29000\n",
      "evaluating user: 30000\n",
      "evaluating user: 31000\n",
      "evaluating user: 32000\n",
      "evaluating user: 33000\n",
      "evaluating user: 34000\n",
      "evaluating user: 35000\n",
      "evaluating user: 36000\n",
      "evaluating user: 37000\n",
      "evaluating user: 38000\n",
      "evaluating user: 39000\n",
      "evaluating user: 40000\n",
      "evaluating user: 41000\n",
      "evaluating user: 42000\n",
      "evaluating user: 43000\n",
      "evaluating user: 44000\n",
      "evaluating user: 45000\n",
      "evaluating user: 46000\n",
      "evaluating user: 47000\n",
      "evaluating user: 48000\n",
      "evaluating user: 49000\n",
      "evaluating user: 50000\n",
      "evaluating user: 51000\n",
      "evaluating user: 52000\n",
      "evaluating user: 53000\n",
      "evaluating user: 54000\n",
      "evaluating user: 55000\n",
      "evaluating user: 56000\n",
      "evaluating user: 57000\n",
      "evaluating user: 58000\n",
      "evaluating user: 59000\n",
      "evaluating user: 60000\n",
      "evaluating user: 61000\n",
      "evaluating user: 62000\n",
      "evaluating user: 63000\n",
      "evaluating user: 64000\n",
      "evaluating user: 65000\n",
      "evaluating user: 66000\n",
      "evaluating user: 67000\n",
      "evaluating user: 68000\n",
      "evaluating user: 69000\n",
      "evaluating user: 70000\n",
      "evaluating user: 71000\n",
      "evaluating user: 72000\n",
      "evaluating user: 73000\n",
      "evaluating user: 74000\n",
      "evaluating user: 75000\n",
      "evaluating user: 76000\n",
      "evaluating user: 77000\n",
      "evaluating user: 78000\n",
      "evaluating user: 79000\n",
      "evaluating user: 80000\n",
      "evaluating user: 81000\n",
      "evaluating user: 82000\n",
      "evaluating user: 83000\n",
      "evaluating user: 84000\n",
      "evaluating user: 85000\n",
      "evaluating user: 86000\n",
      "evaluating user: 87000\n",
      "evaluating user: 88000\n",
      "evaluating user: 89000\n",
      "evaluating user: 90000\n",
      "evaluating user: 91000\n",
      "evaluating user: 92000\n",
      "evaluating user: 93000\n",
      "evaluating user: 94000\n",
      "evaluating user: 95000\n",
      "evaluating user: 96000\n",
      "evaluating user: 97000\n",
      "evaluating user: 98000\n",
      "evaluating user: 99000\n",
      "evaluating user: 100000\n",
      "evaluating user: 101000\n",
      "evaluating user: 102000\n",
      "evaluating user: 103000\n",
      "evaluating user: 104000\n",
      "evaluating user: 105000\n",
      "evaluating user: 106000\n",
      "evaluating user: 107000\n",
      "evaluating user: 108000\n",
      "evaluating user: 109000\n",
      "evaluating user: 110000\n",
      "evaluating user: 111000\n",
      "evaluating user: 112000\n",
      "evaluating user: 113000\n",
      "evaluating user: 114000\n",
      "evaluating user: 115000\n",
      "evaluating user: 116000\n",
      "evaluating user: 117000\n",
      "evaluating user: 118000\n",
      "evaluating user: 119000\n",
      "evaluating user: 120000\n",
      "evaluating user: 121000\n",
      "evaluating user: 122000\n",
      "evaluating user: 123000\n",
      "evaluating user: 124000\n",
      "iter: 12 \tlam_diff:  19.3482\tpsi_diff:  0.1364\n",
      "evaluating user: 0\n",
      "evaluating user: 1000\n",
      "evaluating user: 2000\n",
      "evaluating user: 3000\n",
      "evaluating user: 4000\n",
      "evaluating user: 5000\n",
      "evaluating user: 6000\n",
      "evaluating user: 7000\n",
      "evaluating user: 8000\n",
      "evaluating user: 9000\n",
      "evaluating user: 10000\n",
      "evaluating user: 11000\n",
      "evaluating user: 12000\n",
      "evaluating user: 13000\n",
      "evaluating user: 14000\n",
      "evaluating user: 15000\n",
      "evaluating user: 16000\n",
      "evaluating user: 17000\n",
      "evaluating user: 18000\n",
      "evaluating user: 19000\n",
      "evaluating user: 20000\n",
      "evaluating user: 21000\n",
      "evaluating user: 22000\n",
      "evaluating user: 23000\n",
      "evaluating user: 24000\n",
      "evaluating user: 25000\n",
      "evaluating user: 26000\n",
      "evaluating user: 27000\n",
      "evaluating user: 28000\n",
      "evaluating user: 29000\n",
      "evaluating user: 30000\n",
      "evaluating user: 31000\n",
      "evaluating user: 32000\n",
      "evaluating user: 33000\n",
      "evaluating user: 34000\n",
      "evaluating user: 35000\n",
      "evaluating user: 36000\n",
      "evaluating user: 37000\n",
      "evaluating user: 38000\n",
      "evaluating user: 39000\n",
      "evaluating user: 40000\n",
      "evaluating user: 41000\n",
      "evaluating user: 42000\n",
      "evaluating user: 43000\n",
      "evaluating user: 44000\n",
      "evaluating user: 45000\n",
      "evaluating user: 46000\n",
      "evaluating user: 47000\n",
      "evaluating user: 48000\n",
      "evaluating user: 49000\n",
      "evaluating user: 50000\n",
      "evaluating user: 51000\n",
      "evaluating user: 52000\n",
      "evaluating user: 53000\n",
      "evaluating user: 54000\n",
      "evaluating user: 55000\n",
      "evaluating user: 56000\n",
      "evaluating user: 57000\n",
      "evaluating user: 58000\n",
      "evaluating user: 59000\n",
      "evaluating user: 60000\n",
      "evaluating user: 61000\n",
      "evaluating user: 62000\n",
      "evaluating user: 63000\n",
      "evaluating user: 64000\n",
      "evaluating user: 65000\n",
      "evaluating user: 66000\n",
      "evaluating user: 67000\n",
      "evaluating user: 68000\n",
      "evaluating user: 69000\n",
      "evaluating user: 70000\n",
      "evaluating user: 71000\n",
      "evaluating user: 72000\n",
      "evaluating user: 73000\n",
      "evaluating user: 74000\n",
      "evaluating user: 75000\n",
      "evaluating user: 76000\n",
      "evaluating user: 77000\n",
      "evaluating user: 78000\n",
      "evaluating user: 79000\n",
      "evaluating user: 80000\n",
      "evaluating user: 81000\n",
      "evaluating user: 82000\n",
      "evaluating user: 83000\n",
      "evaluating user: 84000\n",
      "evaluating user: 85000\n",
      "evaluating user: 86000\n",
      "evaluating user: 87000\n",
      "evaluating user: 88000\n",
      "evaluating user: 89000\n",
      "evaluating user: 90000\n",
      "evaluating user: 91000\n",
      "evaluating user: 92000\n",
      "evaluating user: 93000\n",
      "evaluating user: 94000\n",
      "evaluating user: 95000\n",
      "evaluating user: 96000\n",
      "evaluating user: 97000\n",
      "evaluating user: 98000\n",
      "evaluating user: 99000\n",
      "evaluating user: 100000\n",
      "evaluating user: 101000\n",
      "evaluating user: 102000\n",
      "evaluating user: 103000\n",
      "evaluating user: 104000\n",
      "evaluating user: 105000\n",
      "evaluating user: 106000\n",
      "evaluating user: 107000\n",
      "evaluating user: 108000\n",
      "evaluating user: 109000\n",
      "evaluating user: 110000\n",
      "evaluating user: 111000\n",
      "evaluating user: 112000\n",
      "evaluating user: 113000\n",
      "evaluating user: 114000\n",
      "evaluating user: 115000\n",
      "evaluating user: 116000\n",
      "evaluating user: 117000\n",
      "evaluating user: 118000\n",
      "evaluating user: 119000\n",
      "evaluating user: 120000\n",
      "evaluating user: 121000\n",
      "evaluating user: 122000\n",
      "evaluating user: 123000\n",
      "evaluating user: 124000\n",
      "iter: 13 \tlam_diff:  16.5451\tpsi_diff:  0.0943\n",
      "evaluating user: 0\n",
      "evaluating user: 1000\n",
      "evaluating user: 2000\n",
      "evaluating user: 3000\n",
      "evaluating user: 4000\n",
      "evaluating user: 5000\n",
      "evaluating user: 6000\n",
      "evaluating user: 7000\n",
      "evaluating user: 8000\n",
      "evaluating user: 9000\n",
      "evaluating user: 10000\n",
      "evaluating user: 11000\n",
      "evaluating user: 12000\n",
      "evaluating user: 13000\n",
      "evaluating user: 14000\n",
      "evaluating user: 15000\n",
      "evaluating user: 16000\n",
      "evaluating user: 17000\n",
      "evaluating user: 18000\n",
      "evaluating user: 19000\n",
      "evaluating user: 20000\n",
      "evaluating user: 21000\n",
      "evaluating user: 22000\n",
      "evaluating user: 23000\n",
      "evaluating user: 24000\n",
      "evaluating user: 25000\n",
      "evaluating user: 26000\n",
      "evaluating user: 27000\n",
      "evaluating user: 28000\n",
      "evaluating user: 29000\n",
      "evaluating user: 30000\n",
      "evaluating user: 31000\n",
      "evaluating user: 32000\n",
      "evaluating user: 33000\n",
      "evaluating user: 34000\n",
      "evaluating user: 35000\n",
      "evaluating user: 36000\n",
      "evaluating user: 37000\n",
      "evaluating user: 38000\n",
      "evaluating user: 39000\n",
      "evaluating user: 40000\n",
      "evaluating user: 41000\n",
      "evaluating user: 42000\n",
      "evaluating user: 43000\n",
      "evaluating user: 44000\n",
      "evaluating user: 45000\n",
      "evaluating user: 46000\n",
      "evaluating user: 47000\n",
      "evaluating user: 48000\n",
      "evaluating user: 49000\n",
      "evaluating user: 50000\n",
      "evaluating user: 51000\n",
      "evaluating user: 52000\n",
      "evaluating user: 53000\n",
      "evaluating user: 54000\n",
      "evaluating user: 55000\n",
      "evaluating user: 56000\n",
      "evaluating user: 57000\n",
      "evaluating user: 58000\n",
      "evaluating user: 59000\n",
      "evaluating user: 60000\n",
      "evaluating user: 61000\n",
      "evaluating user: 62000\n",
      "evaluating user: 63000\n",
      "evaluating user: 64000\n",
      "evaluating user: 65000\n",
      "evaluating user: 66000\n",
      "evaluating user: 67000\n",
      "evaluating user: 68000\n",
      "evaluating user: 69000\n",
      "evaluating user: 70000\n",
      "evaluating user: 71000\n",
      "evaluating user: 72000\n",
      "evaluating user: 73000\n",
      "evaluating user: 74000\n",
      "evaluating user: 75000\n",
      "evaluating user: 76000\n",
      "evaluating user: 77000\n",
      "evaluating user: 78000\n",
      "evaluating user: 79000\n",
      "evaluating user: 80000\n",
      "evaluating user: 81000\n",
      "evaluating user: 82000\n",
      "evaluating user: 83000\n"
     ]
    }
   ],
   "source": [
    "k = 10 # number of latent factors\n",
    "FNAME = 'data/LasVegas_local.pck'\n",
    "PROPORTION = .2 # Train to test percentage\n",
    "ITERS = 20 # Number of times to run EM\n",
    "\n",
    "with open(FNAME, 'rb') as pickle_file:\n",
    "    Ycol = pickle.load(pickle_file)\n",
    "print('Ycol shape: {}'.format(Ycol.shape))\n",
    "\n",
    "# Split the dataset\n",
    "Ytest, Ytrain, test_row_ind, test_col_ind  = split_Y(Ycol, PROPORTION)\n",
    "\n",
    "# Train \n",
    "lam_diff, psi_diff, x, Lam = train(Ytrain, k, ITERS)\n",
    "\n",
    "# Get the correct rows and columns of Lambda and x\n",
    "x_test, Lam_test = split_others(Ycol, x, Lam, PROPORTION, test_row_ind, test_col_ind)\n",
    "\n",
    "# Test the result\n",
    "res = test(Ytest, x_test, Lam_test)\n",
    "pickle.dump(res, open('data/prediction_difference.pck', 'wb'))\n",
    "print('\\nResults\\n===========')\n",
    "print('Mean of test array: {}'.format(np.mean(abs(res))))\n",
    "print('Std of test array: {}'.format(np.std(abs(res))))\n",
    "\n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 162,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEWCAYAAABrDZDcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzt3XeYVNX5wPHvO9tmgd2lL72IiA1B\nXAv2YAGNsSUmGmtiYjTRaDRGk/iLxphEjTWJJVhiiS1qjL0i2AIKKAqKCChlkV53Wba/vz/Omd1h\nmN2dLTN3d+f9PM88c/t95+7sfeeee885oqoYY4xJX6GgAzDGGBMsSwTGGJPmLBEYY0yas0RgjDFp\nzhKBMcakOUsExhiT5iwRBERE7haR/4sav0BEVotIqYj0EpGDRGShHz8xyFhbSkQeEJHrgo7DxCci\nvxGRe4OOIxlE5FMROTzoODoKsXoEbU9ElgCFQDVQA3wGPARMVtXaOMtnAVuAA1T1Yz9tCvCcqt6e\nqrjbmog8ABSr6lWpXDdVRGQP4FagCPejajHwf6r6UqCBBUxEpgH/UtV7/cn4X6o6KIn7e4B2/l1p\n7+yKIHm+pap5wFDgeuAK4L4Gli0EwsCnUdOGxownTEQyW7KeabbngdeBfkBf4Oe4hG7aiH2XU0RV\n7dXGL2AJcGTMtP2AWmBPP/4AcB2wC7AVUKAUeBP3y7IW2Oan5QAFuESyEljh183w2zoHeA/363Q9\ncJ2f/kNgPrAReBUYGhWPAucDC4FNwB34K0Q//8d+3RLcFc04P30A8DSwFvgK+Hkjx+EB4G7cybIE\neCsmhl39vA3AAuC7fvp5QBVQ6T//88APgOej1l0IPBk1vhwY29h2/bwc4CZgGbDax5fr5x0OFAOX\nAWv8sf5BA5+ttz+G3RuYfw7wbsw0BXb2w7nAzcBSYDPwblQcBwP/83+X5cA5CcTeG3jBr7MBeAcI\n+XlX+O9MiT8eR/jp1+B+rQO8DFwYE+/HwMlNHdM4n30a8COgK+47XOv/jqX++xMCrsR9z9cD/wZ6\n+nWH+eN0rv+cb/vpTwKr/LF6G9ijoe9K7P+gP263AV/7121ATiJ/c+BY3Pe/xB/DXwZ9fknKOSvo\nADrjiziJwE9fBlzghx+g/oQd+fJnNrQN4BngH/6fqy/wAfATP+8cXDHURUAm7iRzArAI2M1Puwr4\nX9T2FHfi6A4MwZ3YJ/l5p/gv/b6AADvjrlBCwGzgd0A2sBPwJTCxgePwgP8HOtT/M96OPzn6z7Ec\nd4LPBPYG1gG7xx4fP74T7iQXwp1MluKKAyLzNvp5TW33VuA5oCeQh0syf/bzDvfH8Vogy58EyoAe\ncT6b4JLRC8CJQGHM/HNoPBHcgTthDgQygAP9MRrqj9lpPoZe1Ce4xmL/My4xZPnXIT7GUf54DIj6\nro3ww9dQnwjOAt6LinV3f7xzmjqmcY7NNOBHUce0OGb+xcAMYJDf/j+Ax2L+Fx7y+40kuh/6zxw5\nqc+J+Z5dF7OPJdQngmv9/voCfXBJ9g+J/M1xieEQP9wD/4Oos70CD6Azvmg4EcwAfuuH6768NJEI\ncEVHFZF/Cj/tNGCqHz4HWBazr5eBc6PGQ/4LPtSPK3Bw1Px/A1f64VeBi+PEv3+c/fwa+GcDx+EB\n4PGo8W64eyaDge8B78Qs/w/g6tjjEzV/OTAOOBWYjEuGu+JOUM/5ZRrcLu7EuBV/IvTzxgNf+eHD\ncb9go/8Oa3D3buJ9vkHA36m/gnsbGBn1N4mbCPzfYhswJs42fw08E2d6U7FfCzyLTzRRy+zsP8OR\nQFbMvGuoTwR5fvuR78cfgfubOqYNHJdpNJ4I5uOvSvx4f9yv+kzq/xd2auT/q7tfpqCR78oS6v9/\nFgPHRs2bCCxJ5G+O+/H2EyA/kf/9jvqyewSpNRB3ad1cQ3G/VlaKyCYR2YT7R+wbtczyOOvcHrX8\nBtzJZGDUMquihstwJ2pwJ+rFDcQxILJNv93f4BJVQ+riUtVSH8cAv639Y7Z1Oq68vSFv4f5xD/XD\n04DD/OutqBgb2m4foAswO2reK356xHpVrY4ajz4u21HVYlW9UFVH+P1uxf2SbUpv3D2heMe4oWPf\nVOx/wV0BviYiX4rIlT7GRcAluJP+GhF5XEQGxPksJcCLuCQL7ofGI364JX+rxgwFnona1nzcD4To\n71Hd90ZEMkTkehFZLCJbcCd5cMcxEZEryIilflpEY3/zb+OuEpaKyFsiMj7BfXYodiMmRURkX9xJ\n+N0WrL4cd0XQO+YLG03jrPNHVX0k3sIJ7G9EA9O/UtWRzdjW4MiAiHTDFWt87bf1lqoe1cB6sZ8H\n3Mn+W8Bw4E+4oovTcb+M/x4VY9ztikjkl/geqrqiGZ+hSaq6XETuAB7zk7biTtyRfUefNNcB5bhj\n/HHMppbj7ifFWkcjsfsT+WXAZSKyJ/CmiMxU1Smq+ijwqIjk435A3ACcGWcfjwFXi8jbuEQ1NSqm\nxv5WjYn3d1wO/FBV34udISLD4qz3fVxR55G4JFCAKwqURvYR7Wu2f/hiiJ/WJFWdCZzgn+y7EHfl\nPLjxtToeuyJIMhHJF5HjgMdxl+Fzm7sNVV0JvAbc7LcXEpERInJYI6vdDfzaP+KIiBSIyCkJ7vJe\n4Jciso84O4vIUFxRTImIXCEiuf6X2p4+yTXkWBE5WESygT8AM1R1Oa5sfRcROVNEsvxrXxHZza+3\nGlf2H+0t4Bu4IrJi3A3RSbhy9I/8Mg1uV92ju/cAt4pIX39cBorIxASPSx0R6SEiv/fHJiQivXHl\n2DP8Ih8De4jIWBEJ436RA+DjuB+4RUQG+OM4XkRycL/CjxSR74pIpq9TMrap2EXkOB+L4G6o1gC1\nIjJKRCb4bZdTf/M2npdwJ8xrgSe0/lHnpv5WjVkN9BKRgqhpdwN/9N8pRKSPiJzQyDbycD+E1uOS\n65/i7CP2uxLtMeAqv5/euHtc/2oqcBHJFpHTRaRAVatwT4Q1dOw6NEsEyfO8iJTgfv38FrgFV5bd\nUmfhbtB+hvs19BSubDUuVX0G98vvcX85PQ84JpEdqeqTuDLiR3E3Lv+Le6qjBjgOGIt7YmgdLmkU\nNLAp/DauxhUJ7QOc4fdRAhyNK4r4GldMdQPuZiC4J6R298UH//XrfIF7MuQdP74Fd7P6PR9bItu9\nAleEMsMflzdwN1SbqxJXnv0G7gQxD3eyOicq1mv9/IXseCX4S2AuMNMfmxtwT/kswxVFXOanzwHG\nJBD7SD9eCkwH7lTVqf5zX4/7W63CFSf+Ot4HUtUK4D+4X96PRk1v6pg2SFU/x52Iv/R/ywG4hwae\nwxVjleCS5/6NbOYhXHHOCtz3f0bM/B2+KzGuA2YBn+CO+Yd+WiLOBJb4430+7gq007EKZcYYk+bs\nisAYY9KcJQJjjElzSUsEIhIWkQ9E5GNxDUD93k9/QES+EpE5/jU2WTEYY4xpWjIfH60AJqhqqX/0\n6l0RednPu1xVn0rivo0xxiQoaYlA3V3oUj8aqfbeojvTvXv31mHDhrVRZMYYkx5mz569TlX7NLVc\nUp8aEpEMXNs0OwN3qOoV4pqMHY+7YpiCa9agIs665+EalGLIkCH7LF26NHYRY4wxjRCR2apa1NRy\nSb1ZrKo1qjoW1ybLfr7G469x7cPsi6tlekUD605W1SJVLerTp8mEZowxpoVS8tSQqm7CVVefpKor\n1akA/kn86vTGGGNSJJlPDfURke5+OBc4CvhcRPr7aYJrvndesmIwxhjTtGQ+NdQfeNDfJwgB/1bV\nF0TkTRHpg2swag6u2rYxxmynqqqK4uJiysvLgw6l3QuHwwwaNIisrKwWrZ/Mp4Y+wXVgETt9QrL2\naYzpPIqLi8nLy2PYsGG4AgQTj6qyfv16iouLGT58eIu2YTWLjTHtUnl5Ob169bIk0AQRoVevXq26\ncrJEYIxptywJJKa1x6lTJ4Ip81dz57RFQYdhjDHtWqdOBO8uWsedU+P1+meMMU3r1i1uL6WtMmzY\nMNatWxfIvhvSqRNBYX6Y0opqSisa6t3RGGNMJ08ErgOl1Vvs8TNjTNt4/vnn2X///dl777058sgj\nWb16NQDXXHMNZ599NocccghDhw7lP//5D7/61a8YPXo0kyZNoqqqqm4bN954I6NHj2a//fZj0SJX\nfP3VV18xfvx4Ro8ezVVXXVW3bGlpKUcccQTjxo1j9OjRPPvss23+mTp15/WF+WHAJYIRfVJ3mWWM\naVu/f/5TPvt6S5tuc/cB+Vz9rT2avd7BBx/MjBkzEBHuvfdebrzxRm6++WYAFi9ezNSpU/nss88Y\nP348Tz/9NDfeeCMnnXQSL774IieeeCIABQUFzJ07l4ceeohLLrmEF154gYsvvpgLLriAs846izvu\nuKNuf+FwmGeeeYb8/HzWrVvHAQccwPHHH9+mN9I7+RWBSwRrtuzQpp0xxrRIcXExEydOZPTo0fzl\nL3/h008/rZt3zDHHkJWVxejRo6mpqWHSpEkAjB49miVLltQtd9ppp9W9T58+HYD33nuvbvqZZ55Z\nt6yq8pvf/Ia99tqLI488khUrVtRdhbSVtLgiWGVFQ8Z0aC355Z4sF110EZdeeinHH38806ZN45pr\nrqmbl5PjiqNDoRBZWVl1v9pDoRDV1fX3KqN/zTc0HPHII4+wdu1aZs+eTVZWFsOGDWvz2tad+oqg\nW04m3XIy7R6BMabNbN68mYEDBwLw4IMPtmgbTzzxRN37+PHjATjooIN4/PHHAXfyj95f3759ycrK\nYurUqSSjSf5OfUUA0Dc/x4qGjDEtUlZWxqBBg+rGL730Uq655hpOOeUUevTowYQJE/jqq6+avd2N\nGzey1157kZOTw2OPPQbA7bffzve//31uuOEGTjjhhLplTz/9dL71rW8xevRoioqK2HXXXVv/wWIk\ntWOatlJUVKSzZs1q0bqnTZ5BZU0tT19wYBtHZYxJpvnz57PbbrsFHUaHEe94tYuOadqDfgVhKxoy\nxphGdPpEECka6ghXPsYYE4ROnwgK88JU1tSysayq6YWNMSYNdfpE0K+gvlKZMcaYHXX6RGDNTBhj\nTOM6fSLom2dXBMYY05jOnwjqrgisLoExpnkyMjIYO3Yse+65J6eccgplZWUNLvujH/2Izz77LIXR\ntZ1OnwhyMjPo2TXbrgiMMc2Wm5vLnDlzmDdvHtnZ2dx9990NLnvvvfey++67pzC6ttPpEwFA37wc\nSwTGmFY55JBDWLRoEVu3buWb3/wmY8aMYc8996xrLuLwww+npRVfg5a0JiZEJAy8DeT4/TylqleL\nyHDgcaAXMBs4U1UrkxUHRCqVWdGQMR3Wy1fCqrltu81+o+GY6xNatLq6mpdffplJkybxyiuvMGDA\nAF588UXAtQXU0SXziqACmKCqY4CxwCQROQC4AbhVVXcGNgLnJjEGwNUlsCsCY0xzbdu2jbFjx1JU\nVMSQIUM499xzGT16NK+//jpXXHEF77zzDgUFBUGH2WpJuyJQV5W31I9m+ZcCE4Dv++kPAtcAdyUr\nDoDCgjDrSiuorqklMyMtSsOM6VwS/OXe1iL3CKLtsssufPjhh7z00ktcddVVHHHEEfzud78LJL62\nktTWR0UkA1f8szNwB7AY2KSqkYa5i4GBDax7HnAewJAhQ1oVR2F+DrUK60or6yqYGWNMS3z99df0\n7NmTM844g+7du3PvvfcGHVKrJTURqGoNMFZEugPPAAm3n6qqk4HJ4FofbU0chXn1HdRYIjDGtMbc\nuXO5/PLL6zqfueuupBZopERK+iNQ1U0iMhUYD3QXkUx/VTAIWJHs/VszE8aYligtLd1h2sSJE5k4\nceIO06dNm5aCiJIjaQXmItLHXwkgIrnAUcB8YCrwHb/Y2cCzyYohIlKpbI0lAmOM2UEyrwj6Aw/6\n+wQh4N+q+oKIfAY8LiLXAR8B9yUxBgB6dc0hIyTWd7ExxsSRzKeGPgH2jjP9S2C/ZO03noyQ+Epl\nVpfAmI5EVeN26G6219r+VtLmWcq++VaXwJiOJBwOs379eutUqgmqyvr16wmHW/4gTKfvvD6iMC+H\nJeu3Bh2GMSZBgwYNori4mLVr1wYdSrsXDocZNGhQi9dPm0TQryDM+19tCDoMY0yCsrKyGD58eNBh\npIW0KRoqzA+zeVsV5VU1QYdijDHtStokgr551lOZMcbEkzaJoL5SmT05ZIwx0dImERTmW+1iY4yJ\nJ30SgfVdbIwxcaVNIsjPzSScFbJEYIwxMdImEYgIhfnWU5kxxsRKm0QArnjI2hsyxpjtpVciKAhb\nC6TGGBMjvRKBb3jO2i4xxph66ZUI8sNsq6phS3l10wsbY0yaSK9E4CuVWfGQMcbUS69EUNfMhD05\nZIwxEemVCPLrO7E3xhjjpGUisEplxhhTL60SQW52BvnhTLtHYIwxUdIqEYC7KrCiIWOMqddkIhCR\nLiLyfyJyjx8fKSLHJbDeYBGZKiKficinInKxn36NiKwQkTn+dWzrP0bi+hVYMxPGGBMtka4q/wnM\nBsb78RXAk8ALTaxXDVymqh+KSB4wW0Re9/NuVdWbWhJwa/XNC7N4zbogdm2MMe1SIkVDI1T1RqAK\nQFXLAGlqJVVdqaof+uESYD4wsBWxtonC/BzWlFRQW2u1i40xBhJLBJUikgsogIiMAJpVtiIiw4C9\ngff9pAtF5BMRuV9EejSwznkiMktEZq1du7Y5u2tUv4Iw1bXK+q2VbbZNY4zpyBJJBFcDrwCDReQR\nYArwq0R3ICLdgKeBS1R1C3AXMAIYC6wEbo63nqpOVtUiVS3q06dPortrUl/roMYYY7bT5D0CVX1d\nRD4EDsAVCV2sqgkVsotIFi4JPKKq//HbWx01/x6avtfQpgrz6zux33NgQSp3bYwx7VIiTw2dBFSr\n6ouq+gJQLSInJrCeAPcB81X1lqjp/aMWOwmY1/ywW846sTfGmO0lVDSkqpsjI6q6CVdc1JSDgDOB\nCTGPit4oInNF5BPgG8AvWhJ4S/XuloOIFQ0ZY0xEIo+PxksWiRQpvUv8p4teSmCfSZOVEaJX1xxL\nBMYY4yVyRTBLRG4RkRH+dQuuXkGH1a/AEoExxkQkkgguAiqBJ/yrAvhZMoNKtsI8q11sjDERiRTx\nbAWuTEEsKdM3P8yc5ZuCDsMYY9qFJhOBiOwC/BIYFr28qk5IXljJ1S8/zPqtlVRW15KdmXbt7hlj\nzHYSuVn8JHA3cC9Qk9xwUiNSl2BtaQUDu+cGHI0xxgQrkURQrap3JT2SFKrrqWxzuSUCY0zaS6Rc\n5HkR+amI9BeRnpFX0iNLokgisA5qjDEmsSuCs/375VHTFNip7cNJjehmJowxJt0l8tTQ8FQEkko9\numSTlSGsskdIjTEm4R7KrhKRyX48oR7K2rNQSOibF7aiIWOMIbF7BP/EVSg70I+vAK5LWkQpUpif\nw+oSSwTGGJO0Hsrau34FYVZttkRgjDEp6aGsPXJFQx3+YxhjTKsl8tRQbA9lBwHnJDOoVCjMD1NS\nUc3Wimq65iRyGIwxpnNq9AzoO5f5HDiZFvRQ1p71K6h/hHSnPt0CjsYYY4LTaCJQVRWRl1R1NPBi\nimJKicK8+p7KLBEYY9JZIvcIPhSRfZMeSYr1jdQutieHjDFpLpHC8f2B00VkKbAVVzykqrpXUiNL\nskjfxfbkkDEm3SWSCCYmPYoAdMvJpGt2hnVQY4xJe00WDanqUmAwMMEPlyWyXkdQmB+2SmXGmLSX\nSBMTVwNXAL/2k7KAfyUzqFQpzA+z2oqGjDFpLpFf9icBx+PuD6CqXwN5Ta0kIoNFZKqIfCYin4rI\nxX56TxF5XUQW+vcerfkArWHNTBhjTII1i1VVqa9Z3DXBbVcDl6nq7rg6CD8Tkd1x/R9PUdWRwBQC\n7A+5MN91Yu8+njHGpKdEEsG/ReQfQHcR+THwBnBPUyup6kpV/dAPlwDzgYHACcCDfrEHgRNbEnhb\nKMwPU1ldy6ayqqBCMMaYwDX41JCI5KhqhareJCJHAVuAUcDvVPX15uxERIYBewPvA4WqutLPWgUU\nNrDOecB5AEOGDGnO7hIW6alsdUk5PbpmJ2UfxhjT3jX2+Oh0YJyIPKyqZwLNOvlHiEg34GngElXd\n4lqtcHzN5bjlMqo6GZgMUFRUlJSym/qeyirYtV8y9mCMMe1fY4kgW0S+DxwoIifHzlTV/zS1cRHJ\nwiWBR6KWXy0i/VV1pYj0B9a0JPC2UHdFYE8OGWPSWGOJ4HzgdKA78K2YeQo0mgh8g3X3AfNV9Zao\nWc/h+kG+3r8/28yY20xf67vYGGMaTQT9VfUCEfnIF9M010HAmcBcEZnjp/0GlwD+LSLnAkuB77Zg\n220iJzODHl2yWGWJwBiTxhpLBL8GnsRdGTQ7EajquzTck9kRzd1eskQeITXGmHTVWCJYLyKvAcNF\n5LnYmap6fPLCSp3C/LC1QGqMSWuNJYJvAuOAh4GbUxNO6hXm5zB/5ZagwzDGmMA0mAhUtRKYISIH\nquraFMaUUv3yw6wrraC6ppbMjE7Rlp4xxjRLYxXKblPVS4D74z3r31mKhvrmh6lVWL+1su5xUmOM\nSSeNFQ097N9vSkUgQYmc/FdtLrdEYIxJS40VDc3272+JSB8/3OmKiPpFKpXZI6TGmDTVaKG4iFwj\nIuuABcAXIrJWRH6XmtBSo66ZiRJ7hNQYk54aTAQicimuUti+qtpTVXvg+i8+SER+kaoAk61Xtxwy\nQmLNTBhj0lZjVwRnAqep6leRCar6JXAGcFayA0uVjJDQp1uOFQ0ZY9JWY4kgS1XXxU709wmykhdS\n6rmeyqxoyBiTnhpLBJUtnNfh9LW+i40xaayxx0fHiEi8KrcCdKrnLPvlh5m5ZEPQYRhjTCAae3w0\nI5WBBKkwP4dNZVWUV9UQzkqbj22MMUBifRZ3en19XYI11gqpMSYNWSIgqlKZtUJqjElDjdUjyEll\nIEEqtNrFxpg01tgVwXQAEXm4kWU6hUjt4lX25JAxJg0ltfP6jqIgN4uczBBrrC6BMSYNJa3z+o5E\nRHyXlXZFYIxJP409Pvou8K6IzFLV+1IYUyAK83OsaMgYk5YSeWroYRH5uYg85V8XiUiTTUyIyP0i\nskZE5kVNu0ZEVojIHP86tlXRtyHXd7EVDRlj0k8iieBOYB//fieuH+O7EljvAWBSnOm3qupY/3op\n0UCTLVI0pLpDZ2zGGNOpNXaPIGJfVR0TNf6miHzc1Eqq+raIDGtpYKlWmJ9DWWUNJRXV5Ic7VZt6\nxhjTqESuCGpEZERkRER2Ampasc8LReQTX3TUo6GFROQ8EZklIrPWrk1+x2iFdbWL7T6BMSa9JJII\nLgemisg0EXkLeBO4rIX7uwsYAYwFVgI3N7Sgqk5W1SJVLerTp08Ld5e4+kpldp/AGJNemiwaUtUp\nIjISGOUnLVDVFp0tVXV1ZFhE7gFeaMl2kiG6E/s6Jatgw5cw9MCAojLGmORLqK0hVa1Q1U/8q8U/\nmUWkf9ToScC8hpZNtfq+i6MSwZt/gH8eC3OfCigqY4xJvkRuFreIiDwGHA70FpFi4GrgcBEZi6uQ\ntgT4SbL231xdsjPJC2du3wLp0v8BCs+cD116wogJgcVnjDHJkrREoKqnxZncriumFeaH64uGSla7\nYqFDL4cFr8DjZ8A5z8PAfYIN0hhj2lhCRUMispeIHC8iJ0deyQ4sCP3yw/VFQ8tnuPeRE+GMp6Br\nb3jkFFi3MLgAjTEmCZpMBCJyP3A/8G1cm0PfAo5LclyB6JufU180tGwGZIah/xjI6wdnPgMSgodP\nhi1fBxuoMca0oUSKhg5Q1d2THkk70M/XLq6tVULLpsPAIsjMdjN7jYDTn4IHjoN/fRt+8BLkNlgN\nwhhjOoxEioami0haJILC/DDVtcqGTRtg5Scw5IDtFxgwFk59BNYvgkdPhcqyYAI1xpg2lEgieAiX\nDBb4GsFzReSTZAcWhMgjpKWLZ4DWwJDxOy6002Fw8j2w/H146gdQU53iKI0xpm0lUjR0H3AmMBeo\nTW44wYpUKqtdOh0QGLxv/AX3OBHKboIXL4PnL4YT/g4iqQvUGGPaUCKJYK2qPpf0SNqBSCLosnIm\nFO4J4YKGF973R1C6Ft663j1RdNTvUxSlMca0rUQSwUci8ijwPFBX26ozdVUZ0Scvh0ypoefGj2Gf\nM5pe4fArYesaeO826NYXxv8s+UEaY0wbSyQR5OISwNFR0zpVV5URWRkhDuiykuyabTveKI5HBI69\nCcrWw6u/gS69Ycz3kh+oMca0oUYTgYhkAJ+o6q0piidwh+QshDLi3yiOJ5Thbh6XbYBnf+qaohh5\nVFJjNMaYttToU0OqWgPEayqi0xrHAlaH+kLBwMRXysyBUx+FvrvDv8+C4lnJC9AYY9pYIo+Pvici\nfxeRQ0RkXOSV9MiCoMqoyk+ZraOaXjZWOB/OeBq6FcIj34G1C1oeR1W5q8fwyZMwv9201G2M6aQS\nuUcw1r9fGzVNgc7XFOfGr8ivXs97VSM5qqaWrIyEmmKq162va4ri/onw8Elw7mtQMKjh5au2wbov\nXNJYM9+9r50PG5eARj2pe85LMOygFn0kY4xpSiId03wjFYG0C8tcQ3Mza0extqSCAd1zm7+NnsPd\nlcE/j3XtEv3wFddm0bovYO3n7rXGv29cgsupQCgTeu0M/UbD6O9Cn1GuWYsnzoAXfgHnv1vf3IUx\nxrShJhOBiBTg+hI41E96C7hWVTcnM7BALJtOVVY+C8sHsmpLecsSAbiT+WmPuURw+xioKKH+hJ8F\nvUe65irGnAp9dnWvnjvFP9EfexM8+l3431/h0F+2+KMZY0xDEikauh/Xk9h3/fiZwD+BztcU9bL3\nKe9fhJaEWt+J/bCD4fuPwyf/dif5PqOgz27uiiEjK/Ht7DIRdjse3v4L7Hmy25YxxrShRBLBCFX9\ndtT470VkTrICCszW9bBuAaFDToEv2qgT+xET2qZXs2NugMVvwou/dMVO1pyFMaYNJXI3dJuIHBwZ\nEZGDgG3JCykgy98HIHfEQWRlCKtae0XQlvIHwISrYPEU+LTT1eMzxgQskURwPnCHiCwRkaXA3/20\nzmXZdMjIJjRwHH3zXL8E7coN4SkYAAAaHklEQVR+57lOcl75NZR3vtszxpjgNJkIVPVjVR0D7AWM\nVtW9VfXj5IeWYstmwIBxkBXevqey9iKUAcfdBlvXwpQ/BB2NMaYTSeSpoRxcN5XDgEzx5dOqem0j\nq3UsVdvg64/qGo3rlx9m4ZrSgIOKY+A42PfH8MFkGHMaDNon6IiMMZ1AIkVDzwInANXA1qhXo0Tk\nfhFZIyLzoqb1FJHXRWShf28ffT2u+BBqq+raFyrMb4dFQxETfutqL79wsXWKY4xpE4kkgkGq+j1V\nvVFVb468EljvAWBSzLQrgSmqOhKY4seDt2y6ex+8H+A6sS8pr6assh2eaMMFcMz1sGquuzIwxphW\nSiQR/E9ERjd3w6r6NrAhZvIJwIN++EHgxOZuNymWzXDP+HfpCbiiIWijR0iTYfcTYeejYOofYfOK\noKMxxnRwiSSCg4HZbdRncaGqrvTDq4DChhYUkfNEZJaIzFq7dm0Ld5eA2hpY/sF2/Q8U1iWCdlo8\nJALfvAlqq+HlXwUdjTGmg0ukQtkxydixqqqIaCPzJwOTAYqKihpcrtXWzIeKzdv1PxDpxL7dJgKA\nHsPgsF/BlGthwcswKil/JmNMGkjk8dGl8V4t3N9qEekP4N/XtHA7bSdyf6AjXRFEjL/ItVP00uVQ\n2eT9e2OMiauZ7Sy32nPA2X74bNwTScFaNgPyBkD3IXWTuuVk0iU7o/3eI4jIzHZ1CzYvh2nXBx2N\nMaaDSloiEJHHgOnAKBEpFpFzgeuBo0RkIXCkHw/WshnuaiCq/R4Rad+PkEYbOh72PhOm3wGr5jW9\nvDHGxEjkHkGLqGpDXVwekax9Ntum5bClGIZcvMOswvycjpEIAI66Fha85Pot+OGrEEr1hZ4xpiNr\n8IwhIoNF5HEReUdEfiMiWVHz/pua8JLMd0QTfX8gwl0RtPOioYguPeHo66D4A/jwwaaXN8aYKI39\ndLwfmAZcBPQH3hKRXn7e0CTHlRrLpkN2HhTuscOsSNGQavIeWGpTY06DYYfAG1dDafD34I0xHUdj\niaCPqt6tqnNU9SLgTuBtERlBXXdbHdyyGa42cShjh1mF+WEqqmvZvK0qgMBaQAS+eQtUlsGrvw06\nGmNMB9JYIsgSkXBkRFX/BVwMvIq7QujYtm2ENZ9tV38gWn1dgg5SPATQZxc4+BKY+2/4clrQ0Rhj\nOojGEsG9wP7RE1T1DeAUXNeVHdvymYDGvT8A9XUJ2lUHNYk45DLoMRxeuBSqOljsxphANJgIVPVW\nVX0rzvSPgBeTGlUqLJsOoUwYGL8p534dpVJZrKxcOO4W2LAY3r016GiMMR1AS58zvLRNowjCshnQ\nfyxkd4k7u0+eKxpqdSf2QRgxAfb8Drx7C6xbGHQ0xph2rqWJoGP3nl5dAStmN1gsBBDOyqB7l6yO\nVzQUMfFPkJkLz14IpUlstM8Y0+G1NBF07KeGvp4DNRUN3iiO6NfMugS1tcq60grmrdjMlPmrefT9\nZbz26arWRtsyeYVw7I2wYhb8bRy8e5tLgMYYE6PBmsUiUkL8E74AuUmLKBXiNDQXT9/8cF3RUGlF\nNau3lLN6czmrS8pZtbnCjW8pZ9WWctZsqWBNSTlVNdsfMhF45eJDGdUvLykfpVFjTnX3QF67ytUv\nmHW/q4W8+wnbNalhjElvDSYCVQ3gzJUiy2ZAr5HQtXejixXm5fDeonXsefWrlFbs2FtZXk4mffNz\n6FcQZv+delKYH6ZffpjC/DCF+TnkhTM56Y7/cfNrC5h8VlGyPk3jeo+E7z8Bi6e6+gVPng1DDoRJ\nf4IBewcTkzGmXUlaW0PtVm0tLJ8Bux7X5KLf2WcQ26pq6N3NnewL83O2O9l3zWn68P340J245fUv\nmLN8E2MHd2+LT9AyI74B578DHz4Eb14Hkw93tZGP+B3kDwguLmNM4KQjNKFQVFSks2bNapuNrfkc\n7twfTrgT9j69bbbZiNKKag69cSp7DMjn4XP3b3qFVCjfDO/cAjPudI/QHnQxHPjzBp+gMsZ0TCIy\nW1WbLI5Iv2YqE7w/0Fa65WTy08NH8M7CdUxfvD4l+2xSuACO+j1cOBNGHg3T/gx/2wc+fsJdMRlj\n0koaJoIZ0LUv9NwpZbs844Ch9MsPc9NrC9pXI3Y9hsF3H4QfvOKeMnrmPLj3iPpWWY0xaSENE8H0\nHTqiSbZwVgYXHbEzs5duZNqCdvhM/9Dx8KM34cS7oWQl3D8RnjwHNra0R1JjTEeSXolgy9ewaWmT\n9QeS4btFgxnSswt/eXUBtbXt6KogIhSCsafBRbPhsCthwSvw933hzT9acZExnVx6JYJGOqJJtqyM\nEL84aiSfrdzCy/MCqmSWiOyu8I1fu4Sw+wnw9o3w/M8tGRjTiaVfIsjqCv32CmT3x48ZyMi+3bj5\n9QVU17TzE2vBQDh5Mhx6OXz0MLz4C0sGxnRSaZYIpsOgIsgIpvpERki47OhRfLl2K898tCKQGJpF\nBL7xWzj4Upj9ALz0S2hPN7uNMW0ikEQgIktEZK6IzBGRNqog0ITyLbB6XiD3B6JN3KOQvQYVcNsb\nC6morgk0loSIuEpnB10Ms+6Dly63ZGBMJxPkFcE3VHVsIpUd2kTxTNDaQO4PRBNxVwUrNm3jiZnL\nA40lYSJw5O9h/IUw8x545UpLBsZ0IulTNLRsBkiGKxoK2KEje7Pf8J787c1FbKvsAFcF4JLB0dfB\nAT+F9+927RZZMjCmUwgqESjwmojMFpHzUrLHZdOh32jICb4tPRHh8omjWFtSwYPTlwQdTuJEXD8H\n+58PM+6A1//PkoExnUBQieBgVR0HHAP8TEQOjV1ARM4TkVkiMmvt2lZWwqqpguJZgd8fiLbvsJ4c\nPqoPd7+1mC3lVUGHkzgRmHQ97Ptj+N/f4I1rLBkY08EFkghUdYV/XwM8A+wXZ5nJqlqkqkV9+vRp\n3Q5XfgLV2wK/PxDrl0ePYlNZFfe+81XQoTSPCBz7Fyg6F967Dd78gyUDYzqwlCcCEekqInmRYeBo\nYF5Sd5rihuYStefAAo7Zsx/3vfMlG7ZWBh1O84jAsTfBPufAOzfD1D8FHZExpoWCuCIoBN4VkY+B\nD4AXVfWVpO5x2XToMRzy+iV1Ny1x6VG7sK2qhrumLQo6lOYLheCbt8K4s1wN5GnXBx2RMaYFUl6z\nSlW/BMakcIfuiaGRR6dsl80xsjCPE/ceyEPTl3LuwTvRryAcdEjNEwrBcbe7WsfT/uyezDrs8qCj\nMsY0Q+d/fHT9Yihb1+6KhaL94shdqFXlb28uDDqUlgmF4Pi/uh7Ppl4Hb98UdETGmGbo/Img7v5A\n+3liKNbgnl343r6DeWLmcpatLws6nJYJZcAJd8Be33M3j9+9NeiIjDEJSoNEMANye7pO3NuxiyaM\nJCMk3Dbli6BDablQBpx4F+z5HfdY6Xt/DToiY0wC0iARpL4jmpYozA9z9oHDeOajFSxcXRJ0OC0X\nyoCT/gF7nOwqnL19E1SUBh2VMaYRnTsRlK6BDYvb9f2BaOcfNoKu2Znc8noHvioA17rryffA7ie6\nYqIbd4J/fRs+uAc2LQs6OmNMjM6dCOo6omm/9wei9eyazbkHD+fleauYW7w56HBaJyMTvvNPOPt5\n2PdHsOFL14z1baPhzgNhyrWw/AOo7SBtLRnTiXX+RJAZhv6pe1q1tX50yHC6d8niptcWBB1K64VC\nMPxQmPQn+PlHcOEs13Bdbg949za47yi4aRd45gL47Fmo6MBFYsZ0YMH00JIqu0yE7oMhMyfoSBKW\nF87i/MNGcP3Ln/PBVxvYb3jPoENqO71HuteBF8G2jbBoCnzxCix4CT5+FEJZMOxgGHWM+9v1GBZ0\nxMakBdEO0EZMUVGRzpqVmv5r2oNtlTUc+pepDO/VlSd+cgDSzm90t1pNNSx/H754Gb54Fdb5eyR9\ndoPdjoP9L4CuvYKN0ZgOSERmJ9LnS+cuGuqgcrMzuGjCznywZANvL1wXdDjJl5EJww5yxUYXzoSL\nPoSJf4ZufV07Rn8d654+qtwadKTGdEp2RdBOVVbXMuHmadTWKvsN70n3Ltn06JJN9y5ZdO+SRY+Y\n8W45mZ3zymHN5+7G8oIXoVs/OPxK2PvMwPqdNqYjSfSKwBJBO/beonXc+OoCNmytYNPWKkoqqhtc\nNitDKMjNpodPEgVdsujRJYshPbvw7X0G0b8gN4WRJ8GyGfD61bB8BvTa2fWjvNvx7b5+iDFBskTQ\nCVXV1LJ5WxWbyirZWFbFprIqNpZVbjfuhiv9cBWrS8oJiXD07oWcNX4YB+zUs+NeOajCgpdhyu9h\n7ecwcB846lp3g9kYswNLBAaA5RvK+Nf7S3li5nI2lVWxS2E3zhw/jJP3HkjXnA5avFJTDR8/5lo7\n3bLCtSx7xNXQb8+gIzOmXbFEYLZTXlXDcx9/zUPTlzBvxRbycjL59j6DOHP8UEb06RZ0eC1TtQ0+\nmOxuKJdvcQ3efeM30GNo0JEZ0y5YIjBxqSofLd/EQ/9bwotzV1JVoxwysjdnjR/GhF37khHqgMVG\n2za61k7f/wdoretP+ZDL7JFTk/YsEZgmrS2p4ImZy/jXjGWs2lLOwO65nHHAUL6372B6ds0OOrzm\n27zCFRfNeQSyu8FBF8MBF0B216AjMyYQlghMwqprann9s9U8NH0p079cT3ZmiG/tNYCzDxzKXoO6\nBx1e80U/cprVBcLdIbsLZOW68bpXrntld42Zl1v/3qUXDCpyw8Z0MJYITIt8sbqEh6Yv4T8frqCs\nsoad+3ajf0HY12OI1F/IokfX7O3qMvTsmk2X7Iz29UTSshnw6TNQWeruJ1SWQVWZG66KGa4sg9qq\n+NvJDMPQg2DnI92r90h7bNV0CJYITKtsKa/iP7OLeWfhOjb4x1E3bK1k87YGTpZAdkaIHl2ztksO\nvbrmMLBHLoN65DKwey6DenShd7fs9pUwImqqdkwSm4th8VRY9Aas912JFgyGnY+AEUfATodBuCDY\nuI1pgCUCkxQ1tcrmbS4pbCqr9O+uPsOGsko2bXXDG33dhrUlFTskj3BWqC4pRJLEoB5d/Hsufbrl\ntM9EsXEpLJ7iGsv78i2oLAHJgMH71SeG/mNdq6vGtAPtOhGIyCTgdiADuFdVr29seUsEHVtJeRUr\nNm2jeMM2976xjOKN2/yrjI1l2yeK7MwQg7rnMrBHLv3yw+RkhcjKCJGdGSLHv0fGI8M5mSGyM3ac\nHs4KkZeTRbdwJnnhTLIy2ugkXVMFxTPdlcKiN2Dlx256l14wYoIrQhoxwbWXFFFVDuWbYNummPeN\nDU+rLof8AVAwKOo1uP49nN82n8d0Su02EYhIBvAFcBRQDMwETlPVzxpaxxJB51ZaUc2KjdtYsWn7\nBFG8cRtrtlRQWVNLZXVt3Xtr5GSGyAtnkhd27TN1y3EJols4k7wcP90njW45mXTJziTeE7WxFyxZ\n5evpuepdeq18l56r3iG7YgMAZd2GkFFTQWblZjJqyhuNTbPz0NweEO6O5LoXmTlQsgo2L3dPRcXe\nx8gp2C5JaOQ93yeKrn2QmnKkqsw12ldZ6roOjQxXbo16lUQN+3lV5ZCR5e6TZOZs/54VjpkeZ5mM\nLP/hal3N8O3e/Yvocd1+PgoS2v4VyvDDGVHjEjMeNb/ujxX1R5OYgbjLRE0T8e+hqOE4743N225b\nMdut21/MsIh7Ai6zZU/xJZoIgqhauh+wSFW/BBCRx4ETgAYTgencuuVkMqpfHqP65TW5rKpSXasu\nMVTXUlVTS4V/jySKyLTK6lrKq2rZWlFNSXkVpRXVlFRUU1JeTWl5/bRlG8rcNL9cbYt/Gw0Evodw\nCrvLUg4Lfczum5eyVXPZTFc2a9cd3jfRjc3alRK6UFOeAVvqtyYCmT4LqYJQS2820Z91DJD1DJB1\nDKxZx4Cy9QxY9TkD5D16iOsfuiUFaxWaRRlhthJmm4QpI0wl2WRSTQ5VZFNFNpXkUEVO1HuI9l+8\n3JEtOPIBRh18UlL3EUQiGAgsjxovBvaPXUhEzgPOAxgyZEhqIjPtnoiQlSFkZYTomoT+hlSVbVU1\nlJZXs6W8mrLKHRv6a+gievvJhwBnUKtKdY1SXVNLVa1/r1Gqa2uprlGqamqpjpleVePXqXXz6n4/\nCogfi/yoXO9fc/2E7Oqt5FetJq9iFfkVq+lSuYGqUA4VGV2oCuVSmdGFilAuFaHIey6VoS5USJga\nyUSBWlX/w1ybPsWrEqKGzNoKMmsrydTK7d4zaisBoVYEJfIKUYugEqJWQQmhkfkqqAi1GkIlhKJ+\nH7UItYT8lUJI3bhQi/gYYpcTapHayBVk9Cdxw4L6wch4/VISGdLaqGUVoj4Ffv26cY1sp3a79SVq\n31L35dHtp2+3fuRv7sbH9tmlqb9Cq7XbxmZUdTIwGVzRUMDhmDQhInTJdkVCfa343aSJIB5vWAEM\njhof5KcZY4wJQBCJYCYwUkSGi0g2cCrwXABxGGOMIYCiIVWtFpELgVdxj4/er6qfpjoOY4wxTiD3\nCFT1JeClIPZtjDFme1YF0hhj0pwlAmOMSXOWCIwxJs1ZIjDGmDTXIVofFZG1wNIWrt4bWNeG4bQ1\ni691LL7Wsfharz3HOFRV+zS1UIdIBK0hIrMSaXQpKBZf61h8rWPxtV5HiLEpVjRkjDFpzhKBMcak\nuXRIBJODDqAJFl/rWHytY/G1XkeIsVGd/h6BMcaYxqXDFYExxphGWCIwxpg012kSgYhMEpEFIrJI\nRK6MMz9HRJ7w898XkWEpjG2wiEwVkc9E5FMRuTjOMoeLyGYRmeNfv0tVfH7/S0Rkrt/3Dh1Ei/NX\nf/w+EZFxKYxtVNRxmSMiW0TkkphlUnr8ROR+EVkjIvOipvUUkddFZKF/79HAumf7ZRaKyNkpjO8v\nIvK5//s9IyLdG1i30e9CEuO7RkRWRP0Nj21g3Ub/15MY3xNRsS0RkTkNrJv049fmVLXDv3DNWS8G\ndgKygY+B3WOW+Slwtx8+FXgihfH1B8b54TzgizjxHQ68EOAxXAL0bmT+scDLuB79DgDeD/BvvQpX\nUSaw4wccCowD5kVNuxG40g9fCdwQZ72ewJf+vYcf7pGi+I4GMv3wDfHiS+S7kMT4rgF+mcDfv9H/\n9WTFFzP/ZuB3QR2/tn51liuC/YBFqvqlqlYCjwMnxCxzAvCgH34KOEJEWtLHd7Op6kpV/dAPlwDz\ncX03dyQnAA+pMwPoLiL9A4jjCGCxqra0pnmbUNW3gQ0xk6O/Yw8CJ8ZZdSLwuqpuUNWNwOvApFTE\np6qvqWqkE+YZuN4BA9HA8UtEIv/rrdZYfP688V3gsbbeb1A6SyIYCCyPGi9mxxNt3TL+n2Ez0Csl\n0UXxRVJ7A+/HmT1eRD4WkZdFZI+UBuZ6yn5NRGaLyHlx5idyjFPhVBr+Bwzy+AEUqupKP7wKKIyz\nTHs5jj/EXeHF09R3IZku9EVX9zdQtNYejt8hwGpVXdjA/CCPX4t0lkTQIYhIN+Bp4BJV3RIz+0Nc\ncccY4G/Af1Mc3sGqOg44BviZiBya4v03yXdtejzwZJzZQR+/7agrI2iXz2aLyG+BauCRBhYJ6rtw\nFzACGAusxBW/tEen0fjVQLv/X4rVWRLBCmBw1PggPy3uMiKSCRQA61MSndtnFi4JPKKq/4mdr6pb\nVLXUD78EZIlI71TFp6or/Psa4BncJXi0RI5xsh0DfKiqq2NnBH38vNWR4jL/vibOMoEeRxE5BzgO\nON0nqx0k8F1IClVdrao1qloL3NPAfoM+fpnAycATDS0T1PFrjc6SCGYCI0VkuP/VeCrwXMwyzwGR\nJzS+A7zZ0D9CW/NlivcB81X1lgaW6Re5ZyEi++H+NilJVCLSVUTyIsO4m4rzYhZ7DjjLPz10ALA5\nqhgkVRr8JRbk8YsS/R07G3g2zjKvAkeLSA9f9HG0n5Z0IjIJ+BVwvKqWNbBMIt+FZMUXfc/ppAb2\nm8j/ejIdCXyuqsXxZgZ5/Fol6LvVbfXCPdXyBe6Jgt/6adfivvQAYVyRwiLgA2CnFMZ2MK6Y4BNg\njn8dC5wPnO+XuRD4FPcUxAzgwBTGt5Pf78c+hsjxi45PgDv88Z0LFKX479sVd2IviJoW2PHDJaSV\nQBWunPpc3D2nKcBC4A2gp1+2CLg3at0f+u/hIuAHKYxvEa58PfIdjDxFNwB4qbHvQorie9h/tz7B\nndz7x8bnx3f4X09FfH76A5HvXNSyKT9+bf2yJiaMMSbNdZaiIWOMMS1kicAYY9KcJQJjjElzlgiM\nMSbNWSIwxpg0Z4nApBURKfXvw0Tk+2287d/EjP+vLbdvTLJYIjDpahjQrETga5U2ZrtEoKoHNjMm\nYwJhicCkq+uBQ3yb8b8QkQzfXv9M3+jZT6Cun4N3ROQ54DM/7b++QbFPI42Kicj1QK7f3iN+WuTq\nQ/y25/l26r8Xte1pIvKUuH4CHklVi7jGRGvqF44xndWVuLbvjwPwJ/TNqrqviOQA74nIa37ZccCe\nqvqVH/+hqm4QkVxgpog8rapXisiFqjo2zr5OxjWkNgbo7dd528/bG9gD+Bp4DzgIeLftP64xDbMr\nAmOco3FtKc3BNRHeCxjp530QlQQAfi4ikaYsBkct15CDgcfUNai2GngL2Ddq28XqGlqbgyuyMial\n7IrAGEeAi1R1uwbgRORwYGvM+JHAeFUtE5FpuHasWqoiargG+580AbArApOuSnDdhka8ClzgmwtH\nRHbxrUfGKgA2+iSwK67bzoiqyPox3gG+5+9D9MF1g/hBm3wKY9qA/fow6eoToMYX8TwA3I4rlvnQ\n37BdS/yuJl8BzheR+cACXPFQxGTgExH5UFVPj5r+DDAe1yKlAr9S1VU+kRgTOGt91Bhj0pwVDRlj\nTJqzRGCMMWnOEoExxqQ5SwTGGJPmLBEYY0yas0RgjDFpzhKBMcakuf8HY6Dvq+MNxRAAAAAASUVO\nRK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x113d70c50>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(np.arange(ITERS), lam_diff, label='Lambda')\n",
    "plt.plot(np.arange(ITERS), psi_diff, label=\"Psi\")\n",
    "plt.legend()\n",
    "plt.xlabel('Iteration')\n",
    "plt.ylabel('L2 norm of Difference')\n",
    "plt.title('Difference between Successive Iterations')\n",
    "pass"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python [conda env:py36]",
   "language": "python",
   "name": "conda-env-py36-py"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}