2017-08-01_mg_test_sprng.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import os\n",
    "import re\n",
    "from os.path import join\n",
    "from io import StringIO\n",
    "from glob import glob\n",
    "import shutil\n",
    "\n",
    "import sympy\n",
    "import pandas as pd\n",
    "from sympy import Symbol, sqrt, cos, pi, symbols\n",
    "import numpy as np\n",
    "import matplotlib.pylab as plt\n",
    "import pandas as pd\n",
    "from collections import OrderedDict\n",
    "\n",
    "import qnet\n",
    "from qnet.algebra import *\n",
    "\n",
    "import QDYN\n",
    "\n",
    "from src.notebook_plots_v1 import plot_bs_decay, display_hamiltonian, display_eq, show_summary_dicke\n",
    "from src.single_sided_network_v1 import network_slh\n",
    "from src.dicke_single_model_v2 import write_dicke_single_model, err_dicke_single\n",
    "from src.doit_actions_v1 import run_traj_prop\n",
    "\n",
    "from doit.tools import register_doit_as_IPython_magic"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "qnet.init_printing()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "attributes": {
     "classes": [],
     "id": "",
     "n": "4"
    },
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "DOIT_CONFIG = {\n",
    "    'backend': 'json',\n",
    "    'dep_file': '.doit_db/2017-08-01_mg_test_sprng.json',\n",
    "}\n",
    "register_doit_as_IPython_magic()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "$\n",
    "\\newcommand{ket}[1]{\\vert #1 \\rangle}\n",
    "\\newcommand{bra}[1]{\\langle #1 \\vert}\n",
    "\\newcommand{Op}[1]{\\hat{#1}}\n",
    "$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# How should the error be evaluated?"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Set up runfolders"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "SOURCE_RF = './data/dicke1_QSL_1traj_xxx/rf_2nodes_T1_50'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "DATA_ROOT = './data/dicke1_test_sprng'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "params_data_str = r'''\n",
    "# n_nodes     T  n_trajs      mcwf  seed\n",
    "        2  1.50   10000       True    12\n",
    "'''\n",
    "\n",
    "def read_bool(s):\n",
    "    \"\"\"convert string to boolean\"\"\"\n",
    "    true_vals = ['1', 'true']\n",
    "    false_vals = ['0', 'false']\n",
    "    s = s.lower()\n",
    "    if s in true_vals:\n",
    "        return True\n",
    "    elif s in false_vals:\n",
    "        return False\n",
    "    else:\n",
    "        raise ValueError(\"Cannot convert '%s'\" % s)\n",
    "\n",
    "params_df = pd.read_fwf(\n",
    "        StringIO(params_data_str), comment='#', header=1,\n",
    "        names=['n_nodes', 'T', 'n_trajs', 'mcwf', 'seed'],\n",
    "        converters={\n",
    "            'n_nodes': int, 'T': float, 'n_trajs': int,\n",
    "            'mcwf': read_bool, 'seed': int})\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>n_nodes</th>\n",
       "      <th>T</th>\n",
       "      <th>n_trajs</th>\n",
       "      <th>mcwf</th>\n",
       "      <th>seed</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2</td>\n",
       "      <td>1.5</td>\n",
       "      <td>10000</td>\n",
       "      <td>True</td>\n",
       "      <td>12</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   n_nodes    T  n_trajs  mcwf  seed\n",
       "0        2  1.5    10000  True    12"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "params_df"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def runfolder(row):\n",
    "    int_part, frac_part = (\"%.2f\" % row['T']).split('.')\n",
    "    if row['mcwf']:\n",
    "        label = 'mcwf_%s' % row['seed']\n",
    "    else:\n",
    "        label = 'non_herm'\n",
    "    rf = 'rf_%dnodes_T%s_%s_%s' % (\n",
    "        row['n_nodes'], int_part, frac_part, label)\n",
    "    return join(DATA_ROOT, rf)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def copy_optimized_pulses(from_rf, to_rf):\n",
    "    for file in glob(join(from_rf, 'pulse*.oct.dat')):\n",
    "        shutil.copy(file, to_rf)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def task_create_runfolder():\n",
    "    \"\"\"Create all necessary runfolders for the runs defined in params_df\"\"\"\n",
    "    slh_models = {}\n",
    "    runfolders = set()\n",
    "    for ind, row in params_df.iterrows():\n",
    "        rf = runfolder(row)\n",
    "        if rf in runfolders:\n",
    "            raise ValueError(\"Duplicate runfolder %s, row %d\" % (rf, ind))\n",
    "        runfolders.add(rf)\n",
    "        n_nodes = int(row['n_nodes'])\n",
    "        try:\n",
    "            slh = slh_models[n_nodes]\n",
    "        except KeyError:\n",
    "            slh = network_slh(\n",
    "                n_cavity=2, n_nodes=n_nodes, topology='driven_bs_fb')\n",
    "            slh_models[n_nodes] = slh\n",
    "        yield {\n",
    "            'name': str(rf),\n",
    "            'actions': [\n",
    "                (write_dicke_single_model, [slh, ], dict(\n",
    "                    rf=rf, T=row['T'], theta=0, nt=500, seed=row['seed'],\n",
    "                    kappa=1.0, E0_cycles=2, mcwf=row['mcwf'], non_herm=True,\n",
    "                    observables='last')),\n",
    "                (copy_optimized_pulses, [SOURCE_RF, rf])\n",
    "            ],\n",
    "            'targets': [join(rf, 'config'), join(rf, 'pulse1.oct.dat')],\n",
    "            'uptodate': [True, ] # up to date if targets exist\n",
    "        }\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "N_PROCS = 1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def task_propagate():\n",
    "    \"\"\"Run optimization for every runfolder from params_df\"\"\"\n",
    "    for ind, row in params_df.iterrows():\n",
    "        rf = runfolder(row)\n",
    "        n_procs = N_PROCS\n",
    "        if not row['mcwf']:\n",
    "            n_procs = 1\n",
    "        yield {\n",
    "            'name': str(rf),\n",
    "            'actions': [\n",
    "                (run_traj_prop, [rf, ], dict(n_trajs=row['n_trajs'], n_procs=n_procs))],\n",
    "            'file_dep': [join(rf, 'pulse1.oct.dat'),],\n",
    "            'targets': [join(rf, 'state_final.dat.1'),]}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      ".  create_runfolder:./data/dicke1_test_sprng/rf_2nodes_T1_50_mcwf_12\n"
     ]
    }
   ],
   "source": [
    "%doit create_runfolder"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "-- create_runfolder:./data/dicke1_test_sprng/rf_2nodes_T1_50_mcwf_12\n",
      ".  propagate:./data/dicke1_test_sprng/rf_2nodes_T1_50_mcwf_12\n"
     ]
    }
   ],
   "source": [
    "%doit propagate"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "with open(\"./data/dicke1_test_sprng/rf_2nodes_T1_50_mcwf_12/P_target_combined.dat\") as in_fh:\n",
    "    fid_table = pd.read_fwf(\n",
    "            in_fh, comment='#', header=1,\n",
    "            names=['T', 'Fid', 'Var'])\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>T</th>\n",
       "      <th>Fid</th>\n",
       "      <th>Var</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59969</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59970</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59971</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59972</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59973</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59974</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59975</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59976</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59977</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59978</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59979</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59980</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59981</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59982</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59983</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59984</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59985</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59986</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59987</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59988</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59989</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59990</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59991</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59992</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59993</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59994</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59995</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59996</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59997</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59998</th>\n",
       "      <td>1.5</td>\n",
       "      <td>0.999486</td>\n",
       "      <td>0.999486</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>59999 rows × 3 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "         T       Fid       Var\n",
       "0      1.5  0.999486  0.999486\n",
       "1      1.5  0.999486  0.999486\n",
       "2      1.5  0.999486  0.999486\n",
       "3      1.5  0.999486  0.999486\n",
       "4      1.5  0.999486  0.999486\n",
       "5      1.5  0.999486  0.999486\n",
       "6      1.5  0.999486  0.999486\n",
       "7      1.5  0.999486  0.999486\n",
       "8      1.5  0.999486  0.999486\n",
       "9      1.5  0.999486  0.999486\n",
       "10     1.5  0.999486  0.999486\n",
       "11     1.5  0.999486  0.999486\n",
       "12     1.5  0.999486  0.999486\n",
       "13     1.5  0.999486  0.999486\n",
       "14     1.5  0.999486  0.999486\n",
       "15     1.5  0.999486  0.999486\n",
       "16     1.5  0.999486  0.999486\n",
       "17     1.5  0.999486  0.999486\n",
       "18     1.5  0.999486  0.999486\n",
       "19     1.5  0.999486  0.999486\n",
       "20     1.5  0.999486  0.999486\n",
       "21     1.5  0.999486  0.999486\n",
       "22     1.5  0.999486  0.999486\n",
       "23     1.5  0.999486  0.999486\n",
       "24     1.5  0.999486  0.999486\n",
       "25     1.5  0.999486  0.999486\n",
       "26     1.5  0.999486  0.999486\n",
       "27     1.5  0.999486  0.999486\n",
       "28     1.5  0.999486  0.999486\n",
       "29     1.5  0.999486  0.999486\n",
       "...    ...       ...       ...\n",
       "59969  1.5  0.999486  0.999486\n",
       "59970  1.5  0.999486  0.999486\n",
       "59971  1.5  0.999486  0.999486\n",
       "59972  1.5  0.999486  0.999486\n",
       "59973  1.5  0.999486  0.999486\n",
       "59974  1.5  0.999486  0.999486\n",
       "59975  1.5  0.999486  0.999486\n",
       "59976  1.5  0.999486  0.999486\n",
       "59977  1.5  0.999486  0.999486\n",
       "59978  1.5  0.999486  0.999486\n",
       "59979  1.5  0.999486  0.999486\n",
       "59980  1.5  0.999486  0.999486\n",
       "59981  1.5  0.999486  0.999486\n",
       "59982  1.5  0.999486  0.999486\n",
       "59983  1.5  0.999486  0.999486\n",
       "59984  1.5  0.999486  0.999486\n",
       "59985  1.5  0.999486  0.999486\n",
       "59986  1.5  0.000000  0.000000\n",
       "59987  1.5  0.999486  0.999486\n",
       "59988  1.5  0.999486  0.999486\n",
       "59989  1.5  0.999486  0.999486\n",
       "59990  1.5  0.999486  0.999486\n",
       "59991  1.5  0.999486  0.999486\n",
       "59992  1.5  0.999486  0.999486\n",
       "59993  1.5  0.999486  0.999486\n",
       "59994  1.5  0.999486  0.999486\n",
       "59995  1.5  0.999486  0.999486\n",
       "59996  1.5  0.000000  0.000000\n",
       "59997  1.5  0.000000  0.000000\n",
       "59998  1.5  0.999486  0.999486\n",
       "\n",
       "[59999 rows x 3 columns]"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "fid_table"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "default",
   "language": "python",
   "name": "default"
  },
  "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.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}