cadCAD-0.4.20-example.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Collecting cadCAD\n",
      "  Using cached cadCAD-0.4.20-py3-none-any.whl (66 kB)\n",
      "Processing /Users/bscholtz/Library/Caches/pip/wheels/5d/3a/c0/2de7dc84ea4b39a4705dfbd284e4c14dc1fe0011b7369161b0/fn-0.4.3-py3-none-any.whl\n",
      "Processing /Users/bscholtz/Library/Caches/pip/wheels/85/46/45/fe0554d6310d4198baa18d48f080b5da2a53afe7872a6dc98d/funcy-1.14-py2.py3-none-any.whl\n",
      "Processing /Users/bscholtz/Library/Caches/pip/wheels/83/12/43/a7b335eac30d213e8cbe768bff7456ef4bc76351769864dd79/pathos-0.2.6-py3-none-any.whl\n",
      "Collecting pandas\n",
      "  Using cached pandas-1.1.0-cp38-cp38-macosx_10_9_x86_64.whl (10.6 MB)\n",
      "Processing /Users/bscholtz/Library/Caches/pip/wheels/f4/c7/82/0bcdc3506e6a8365963db24f477b81693438d4b9165c0757b7/multiprocess-0.70.10-py3-none-any.whl\n",
      "Processing /Users/bscholtz/Library/Caches/pip/wheels/f2/af/b4/36b9ce8f3b271919abb82400d1548f534578420ab5ef36fb20/pox-0.2.8-py3-none-any.whl\n",
      "Processing /Users/bscholtz/Library/Caches/pip/wheels/32/6c/df/4335813452c1ff652a854b78093302f0caf2359f49f3417fb4/ppft-1.6.6.2-py3-none-any.whl\n",
      "Processing /Users/bscholtz/Library/Caches/pip/wheels/93/7f/7d/78ec535a4340ef2696aad8b17fe8bb063d56301bd62881b069/dill-0.3.2-py3-none-any.whl\n",
      "Collecting pytz>=2017.2\n",
      "  Using cached pytz-2020.1-py2.py3-none-any.whl (510 kB)\n",
      "Collecting numpy>=1.15.4\n",
      "  Using cached numpy-1.19.1-cp38-cp38-macosx_10_9_x86_64.whl (15.3 MB)\n",
      "Collecting python-dateutil>=2.7.3\n",
      "  Using cached python_dateutil-2.8.1-py2.py3-none-any.whl (227 kB)\n",
      "Collecting six>=1.7.3\n",
      "  Using cached six-1.15.0-py2.py3-none-any.whl (10 kB)\n",
      "Installing collected packages: fn, funcy, dill, multiprocess, pox, six, ppft, pathos, pytz, numpy, python-dateutil, pandas, cadCAD\n",
      "  Attempting uninstall: fn\n",
      "    Found existing installation: fn 0.4.3\n",
      "    Uninstalling fn-0.4.3:\n",
      "      Successfully uninstalled fn-0.4.3\n",
      "  Attempting uninstall: funcy\n",
      "    Found existing installation: funcy 1.14\n",
      "    Uninstalling funcy-1.14:\n",
      "      Successfully uninstalled funcy-1.14\n",
      "  Attempting uninstall: dill\n",
      "    Found existing installation: dill 0.3.2\n",
      "    Uninstalling dill-0.3.2:\n",
      "      Successfully uninstalled dill-0.3.2\n",
      "  Attempting uninstall: multiprocess\n",
      "    Found existing installation: multiprocess 0.70.10\n",
      "    Uninstalling multiprocess-0.70.10:\n",
      "      Successfully uninstalled multiprocess-0.70.10\n",
      "  Attempting uninstall: pox\n",
      "    Found existing installation: pox 0.2.8\n",
      "    Uninstalling pox-0.2.8:\n",
      "      Successfully uninstalled pox-0.2.8\n",
      "  Attempting uninstall: six\n",
      "    Found existing installation: six 1.15.0\n",
      "    Uninstalling six-1.15.0:\n",
      "      Successfully uninstalled six-1.15.0\n",
      "  Attempting uninstall: ppft\n",
      "    Found existing installation: ppft 1.6.6.2\n",
      "    Uninstalling ppft-1.6.6.2:\n",
      "      Successfully uninstalled ppft-1.6.6.2\n",
      "  Attempting uninstall: pathos\n",
      "    Found existing installation: pathos 0.2.6\n",
      "    Uninstalling pathos-0.2.6:\n",
      "      Successfully uninstalled pathos-0.2.6\n",
      "  Attempting uninstall: pytz\n",
      "    Found existing installation: pytz 2020.1\n",
      "    Uninstalling pytz-2020.1:\n",
      "      Successfully uninstalled pytz-2020.1\n",
      "  Attempting uninstall: numpy\n",
      "    Found existing installation: numpy 1.19.1\n",
      "    Uninstalling numpy-1.19.1:\n",
      "      Successfully uninstalled numpy-1.19.1\n",
      "  Attempting uninstall: python-dateutil\n",
      "    Found existing installation: python-dateutil 2.8.1\n",
      "    Uninstalling python-dateutil-2.8.1:\n",
      "      Successfully uninstalled python-dateutil-2.8.1\n",
      "  Attempting uninstall: pandas\n",
      "    Found existing installation: pandas 1.1.0\n",
      "    Uninstalling pandas-1.1.0:\n",
      "      Successfully uninstalled pandas-1.1.0\n",
      "  Attempting uninstall: cadCAD\n",
      "    Found existing installation: cadCAD 0.4.20\n",
      "    Uninstalling cadCAD-0.4.20:\n",
      "      Successfully uninstalled cadCAD-0.4.20\n",
      "Successfully installed cadCAD-0.4.20 dill-0.3.2 fn-0.4.3 funcy-1.14 multiprocess-0.70.10 numpy-1.19.1 pandas-1.1.0 pathos-0.2.6 pox-0.2.8 ppft-1.6.6.2 python-dateutil-2.8.1 pytz-2020.1 six-1.15.0\n",
      "Name: cadCAD\n",
      "Version: 0.4.20\n",
      "Summary: cadCAD: a differential games based simulation software package for research, validation, and         Computer Aided Design of economic systems\n",
      "Home-page: https://github.com/cadCAD-org/cadCAD\n",
      "Author: Joshua E. Jodesty\n",
      "Author-email: joshua@block.science\n",
      "License: LICENSE.txt\n",
      "Location: /usr/local/lib/python3.8/site-packages\n",
      "Requires: funcy, pandas, pathos, fn\n",
      "Required-by: \n"
     ]
    }
   ],
   "source": [
    "!python3 -m pip install cadCAD --force-reinstall\n",
    "!python3 -m pip show cadCAD"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [],
   "source": [
    "from cadCAD.configuration.utils import config_sim\n",
    "from cadCAD.configuration import Experiment\n",
    "from cadCAD.engine import ExecutionMode, ExecutionContext\n",
    "from cadCAD.engine import Executor\n",
    "\n",
    "import pandas as pd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [],
   "source": [
    "system_params = {\n",
    "    'add': [1]\n",
    "}\n",
    "\n",
    "initial_state = {\n",
    "    'a': 0\n",
    "}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [],
   "source": [
    "def state_add_to_a(params, substep, state_history, previous_state, policy_input, **kwargs):\n",
    "    return 'a', previous_state['a'] + params[0]['add']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [],
   "source": [
    "psubs = [\n",
    "    {\n",
    "        'policies': {},\n",
    "        'variables': {\n",
    "            'a': state_add_to_a\n",
    "        }\n",
    "    }\n",
    "]\n",
    "\n",
    "config = {\n",
    "    \"N\": 1,\n",
    "    \"T\": range(100),\n",
    "    \"M\": system_params\n",
    "}\n",
    "\n",
    "c = config_sim(config)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {},
   "outputs": [],
   "source": [
    "from cadCAD import configs\n",
    "del configs[:] # Clear any prior configs"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {},
   "outputs": [],
   "source": [
    "exp = Experiment()\n",
    "exp.append_configs(\n",
    "    initial_state = initial_state,\n",
    "    partial_state_update_blocks = psubs,\n",
    "    sim_configs = c\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {},
   "outputs": [],
   "source": [
    "exec_mode = ExecutionMode()\n",
    "local_mode_ctx = ExecutionContext(context=exec_mode.local_mode)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {},
   "outputs": [],
   "source": [
    "simulation = Executor(exec_context=local_mode_ctx, configs=configs)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "                  ___________    ____\n",
      "  ________ __ ___/ / ____/   |  / __ \\\n",
      " / ___/ __` / __  / /   / /| | / / / /\n",
      "/ /__/ /_/ / /_/ / /___/ ___ |/ /_/ /\n",
      "\\___/\\__,_/\\__,_/\\____/_/  |_/_____/\n",
      "by cadCAD\n",
      "\n",
      "Execution Mode: local_proc\n",
      "Configuration Count: 1\n",
      "Dimensions of the first simulation: (Timesteps, Params, Runs, Vars) = (100, 1, 1, 1)\n",
      "Execution Method: local_simulations\n",
      "SimIDs   : [0]\n",
      "SubsetIDs: [0]\n",
      "Ns       : [0]\n",
      "ExpIDs   : [0]\n",
      "Execution Mode: single_threaded\n",
      "Total execution time: 0.01s\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "<AxesSubplot:>"
      ]
     },
     "execution_count": 41,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy86wFpkAAAACXBIWXMAAAsTAAALEwEAmpwYAAAhUElEQVR4nO3deVxVdf7H8dfXBRFEEVcUERfcxVTcWi2tMXM0s5qaxmy1frPULL8US0vbtGWmmqUap2VspmkDXFLLzLQ9cykBURRxA1FwAZEd7vf3B7d5+DMt5V643HPfz8eDB/eeey/nczrw7nA8vK+x1iIiIs7SyNcDiIiI9yncRUQcSOEuIuJACncREQdSuIuIOFATXw8A0LZtWxsTE+PrMURE/MqmTZsOW2vbne6xBhHuMTExbNy40ddjiIj4FWPM3jM9ptMyIiIOpHAXEXEghbuIiAMp3EVEHEjhLiLiQD8a7saYV4wxecaYtJOWRRhjVhtjdro/t3YvN8aYPxtjMo0xKcaYIXU5vIiInN7ZHLn/Exh3yrIEYI21NhZY474PcCUQ6/6YDrzgnTFFRORc/Gi4W2s/AY6esngSsMh9exFw9UnLX7M1vgLCjTGRXppVRMQxSiqqmP/eNrKPldTJ16/tOfcO1tpc9+2DQAf37c7A/pOel+1e9j3GmOnGmI3GmI35+fm1HENExP98kXmYcc9+yt8/zmJtRt3kn8d/oWqttcaYc37HD2vtQmAhQHx8vN4xREQcr7C0kvkrt/Hmhv3EtAnhzekjGdm9TZ2sq7bhfsgYE2mtzXWfdslzL88Bupz0vCj3MhGRgLY6/RCzl6SSX1TOXZd053djexHctHGdra+2p2WWAdPct6cBS09afrP7qpmRQOFJp29ERALO4RPl/OaNb7jztY20Dgliya8uYNaVfes02OEsjtyNMW8Ao4G2xphs4CFgAfC2MeZ2YC9wvfvpK4HxQCZQAtxaBzOLiDR41lqWfnuAee9upbi8mj9c3ou7R/egaeP6+fOiHw13a+2NZ3hozGmea4FfeTqUiIg/O1BQyuwlaXy0PY/B0eE8OSWO2A5h9TpDg6j8FRFxApfL8p+v97Hgve1UuywPTujHtPNjaNzI1PssCncRES/YfbiYhKQU1u8+yoU92zL/moF0iQjx2TwKdxERD1RVu3jps908s3oHQU0a8eSUOK6Lj8KY+j9aP5nCXUSkltIPHGdmUgqpOYVc0a8Dj1w9gA4tg309FqBwFxE5Z+VV1fz1o0xeWLeL8JCmPH/TEK4c0NHnR+snU7iLiJyDTXuPMTMphcy8E1wzpDNzrupH69AgX4/1PQp3EZGzUFxexdMfZPDPL/bQqVVzFt02nEt6tfP1WGekcBcR+RGf7sxnVnIq2cdKuXlUV2aM60OLZg07Phv2dCIiPlRYUsljK9N5e2M23duG8vZdoxjeLcLXY50VhbuIyGms2nqQOUvSOFJcwS9H9+CeMbF13gfjTQp3EZGT5BeVM3fZVlak5tIvsiWv3DKMAZ1b+Xqsc6ZwFxGhpugreXMODy9Pp7Symvt+0pvpF3evt6Ivb1O4i0jAyz5Wwv2L0/hkRz5Du7bmiSlx9GzfwtdjeUThLiIBy+Wy/Hv9Xp54bzsWmDexP1NHdqWRD4q+vE3hLiIBaVf+CRKSUtiw5xgXxbbl8cm+LfryNoW7iASUymoX//g0i2c/3Enzpo15+rpBTBnSuUFVB3iDwl1EAkZaTiEzk1LYeuA44/p35OGr+9M+rGEUfXmbwl1EHK+sspq/fLSTFz/OonVIEC/cNIQrB0b6eqw6pXAXEUfbuOcoM5JSyMov5tqhUcy+qi/hIQ2v6MvbFO4i4kjF5VU8+f52XvtqL51aNee124ZzcQMu+vI2hbuIOM7HO/K5PzmVA4WlTBsVw30/6U1oAy/68rbA2loRcbSCkgoeWb6NpM3Z9GgXSuLdoxja1T+KvrxN4S4ijrAyNZcHl6ZRUFLJry/tya8v6+lXRV/epnAXEb+Wd7yMB5du5f2tBxnQuSWLbhtO/07+V/TlbQp3EfFL1lre2ZTNo8vTKatyMXNcH+68qBtN/LToy9sU7iLid/YfLeH+xal8uvMww2MiWDBlIN3b+XfRl7cp3EXEb1S7LK99uYenVmVggEcm9eemEc4o+vI2hbuI+IXMvCJmJKaweV8Bl/Rqx+PXDKRzeHNfj9VgKdxFpEGrrHbx94938ec1mYQ0a8wzPxvE1ec5r+jL2xTuItJgpeUUcl9iCttyj3NVXCTzJvanbYtmvh7LLyjcRaTBKaus5tkPd/KPT7NoExrEwqlDuaJ/R1+P5Vc8CndjzO+AOwALpAK3ApHAm0AbYBMw1Vpb4eGcIhIg1mcdISE5ld2Hi7lhWBdmje9Lq+ZNfT2W36n1BaHGmM7APUC8tXYA0Bi4AXgCeMZa2xM4BtzujUFFxNmKyiqZvSSVny38iiqXi9fvGMGCKXEK9lry9LRME6C5MaYSCAFygcuAn7sfXwTMBV7wcD0i4mBrt+fxwOJUco+XcfuF3fjDFb0ICdJZY0/U+r+etTbHGPM0sA8oBT6g5jRMgbW2yv20bKDz6V5vjJkOTAeIjo6u7Rgi4seOFlfwyPJ0Fn+TQ2z7FiT9z/kMiW7t67EcodbhboxpDUwCugEFwDvAuLN9vbV2IbAQID4+3tZ2DhHxP9ZaVqTm8tDSrRSWVnLvmFh+eWkPmjUJ3KIvb/Pk956xwG5rbT6AMSYZuAAIN8Y0cR+9RwE5no8pIk5x6HgZs5eksTr9EHFRrfj3HSPoG9nS12M5jifhvg8YaYwJoea0zBhgI7AWuJaaK2amAUs9HVJE/J+1lrc27OexlduoqHLxwPi+3HpBjIq+6ogn59zXG2MSgc1AFfANNadZVgBvGmMedS972RuDioj/2nekhITkFL7YdYQR3SJ4YkocMW1DfT2Wo3n0z9HW2oeAh05ZnAUM9+TriogzVLssr36+m6c/yKBJo0Y8NnkANw6LVtFXPdC1RiJSJ3Ycqin6+nZ/AZf1ac9jkwcQ2UpFX/VF4S4iXlVR5eKFdbv469qdhAU35bkbzmPioE4q+qpnCncR8Zot+wuYmZTC9oNFTDqvEw9O6EcbFX35hMJdRDxWWlHNn1Zn8PJnu2kfFsxLN8cztl8HX48V0BTuIuKRL3YdZlZyKnuPlHDj8Ghmje9Dy2D1wfiawl1EauV4WSXzV27nja/3ER0Rwn/uGMH5Pdv6eixxU7iLyDlbs+0QDyxOI6+ojDsv6sbvL+9N8yBVBzQkCncROWtHTpQz7910lm05QO8OYbw4dSjndQn39VhyGgp3EflR1lqWbTnAvHfTKSqr5LdjY/nl6J4ENVF1QEOlcBeRH3SwsIzZS1L5cFseg7qE8+SUOHp3DPP1WPIjFO4iclrWWt7csJ/HV2yj0uVi9lV9ufWCbjRWdYBfULiLyPfsOVxMQnIKX2UdZVT3NiyYMpCubVT05U8U7iLyX9Uuyyuf7eaPqzNo2qgR868ZyA3Duqg6wA8p3EUEgIyDRcxI3MKW7ELG9u3Ao1cPoGOrYF+PJbWkcBcJcOVV1Ty/dhfPr8ukZXBT/nLjYCbERepo3c8p3EUC2Df7jjEzKYUdh04weXBn5kzoR0RokK/HEi9QuIsEoJKKKv74wQ5e+Xw3HVsG88ot8VzWR0VfTqJwFwkwX2QeJiE5lX1HS/jFyGhmjutDmIq+HEfhLhIgCksrmb9yG29u2E+3tqG8OX0kI7u38fVYUkcU7iIB4IOtB5m9JI3DJ8q565Lu/G5sL4KbqujLyRTuIg52+EQ5c5dtZXlKLn06hvHStHjiosJ9PZbUA4W7iANZa1nybQ7z3k2npLyaP1zei7tH96BpYxV9BQqFu4jDHCgo5YHFqazNyGdwdE3RV2wHFX0FGoW7iEO4XJbXv97HE+9tp9pleXBCP6adH6OirwClcBdxgKz8EyQkp/L17qNc2LMt868ZSJeIEF+PJT6kcBfxY1XVLl76bDfPrN5BsyaNePLaOK4bGqXqAFG4i/ir9APHmZG0hbSc41zRrwOPXD2ADi1V9CU1FO4ifqa8qpq/fpTJC+t2ER7SlOdvGsKVAzrqaF3+H4W7iB/ZtPcYMxK3sCu/mClDopgzoS/hISr6ku9TuIv4geLyKp5alcGiL/fQqVVz/nnrMEb3bu/rsaQB8yjcjTHhwEvAAMACtwEZwFtADLAHuN5ae8yT9YgEsk935jMrOZXsY6VMG9WV+8b1oUUzHZfJD/P0z9WeA9631vYBBgHbgARgjbU2Fljjvi8i56iwpJL73tnC1Je/JqhxI965exTzJg1QsMtZqfV3iTGmFXAxcAuAtbYCqDDGTAJGu5+2CFgHzPRkSJFA837aQeYsTeNocQW/HN2De8bEquhLzoknhwDdgHzgVWPMIGATcC/QwVqb637OQeC07wBgjJkOTAeIjo72YAwR58grKmPusq2sTD1Iv8iWvHrLMAZ0buXrscQPeRLuTYAhwG+steuNMc9xyikYa601xtjTvdhauxBYCBAfH3/a54gECmstyZtzeHh5OqWV1dz3k95Mv7i7ir6k1jwJ92wg21q73n0/kZpwP2SMibTW5hpjIoE8T4cUcbLsYyXcvziNT3bkE9+1NQumxNGzfQtfjyV+rtbhbq09aIzZb4zpba3NAMYA6e6PacAC9+elXplUxGFcLsu/vtrLE+9vB2DexP5MHdmVRir6Ei/w9J/dfwO8bowJArKAW6m5AudtY8ztwF7geg/XIeI4mXknSEhKYePeY1wUW1P0FdVaRV/iPR6Fu7X2WyD+NA+N8eTrijhVZbWLhZ9k8dyanTRv2pg/XjeIa4Z0VnWAeJ0umBWpJ2k5hcxITCE99zjjB3Zk3sQBtAtr5uuxxKEU7iJ1rKyymj+v2cnfP8kiIjSIF38xhHEDIn09ljicwl2kDm3Yc5SZiSlkHS7muqFRzL6qH61Cmvp6LAkACneROnCivIqn3t/Oa1/tpXN4c/59+wgujG3r67EkgCjcRbxsXUYeDyxO40BhKdNGxXDfT3oTqj4YqWf6jhPxkoKSCh5enk7y5hx6tm9B4t3nM7Rra1+PJQFK4S7iIWst76Ud5MGlaRSUVPLrS3vymzE9adZERV/iOwp3EQ/kHS9jztI0Vm09xMDOrXjtthH069TS12OJKNxFasNayzsbs3l0RTrlVS5mjuvDnRd1o4mKvqSBULiLnKP9R0uYlZzKZ5mHGR4TwYIpA+neTkVf0rAo3EXOUrXLsuiLPTy1KoNGBh65egA3DY9W0Zc0SAp3kbOw81ARM5NS2LyvgNG92/H45IF0Cm/u67FEzkjhLvIDKqtdvLhuF3/5KJPQZo159mfnMem8Tir6kgZP4S5yBqnZhdyXuIXtB4uYEBfJ3In9adtCRV/iHxTuIqcoq6zmmQ938I9PsmjbohkLpw7liv4dfT2WyDlRuIucZH3WERKSU9l9uJgbhnVh1vi+tGquoi/xPwp3EaCorJIn38/gX1/tpUtEc16/YwQX9FTRl/gvhbsEvLUZeTyQnEru8TJuv7Abf7iiFyFB+tEQ/6bvYAlYR4sreGR5Oou/ySG2fQuS/ud8hkSr6EucQeEuAcday4rUXB5aupXC0kruGRPLry7toaIvcRSFuwSUQ8fLmL0kjdXph4iLasW/7xhB30gVfYnzKNwlIFhreWvDfh5buY2KKhf3j+/DbReo6EucS+EujrfvSAkJySl8sesII7pF8MSUOGLahvp6LJE6pXAXx6p2WV79fDdPf5BBk0aNeHzyQG4Y1kVFXxIQFO7iSDsOFTEjMYVv9xcwpk97Hp08gMhWKvqSwKFwF0epqHLxwrpd/HXtTsKCm/LcDecxcZCKviTwKNzFMbbsL2BGYgoZh4qYOKgTD/20H21U9CUBSuEufq+0opo/rc7g5c920z4smJdujmdsvw6+HkvEpxTu4te+3HWEhOQU9h4p4ecjokm4sg8tg1X0JaJwF790vKySBe9t5z/r99G1TQj/uXME5/dQ0ZfIdxTu4nc+TD/E7CVp5BWVcedF3fj95b1pHqTqAJGTeRzuxpjGwEYgx1o7wRjTDXgTaANsAqZaays8XY/IkRPlzHs3nWVbDtCnYxh/nzqUQV3CfT2WSIPkjb+9vhfYdtL9J4BnrLU9gWPA7V5YhwQway1Lv81h7J8+5r20XH43thfLfn2hgl3kB3gU7saYKOAq4CX3fQNcBiS6n7IIuNqTdUhgyy0s5Y5FG7n3zW+JbhPKinsu4t6xsQQ1USeMyA/x9LTMs8AMIMx9vw1QYK2tct/PBjqf7oXGmOnAdIDo6GgPxxCncbksb27Yz/yV26hyWWZf1ZdbL+hGY1UHiJyVWoe7MWYCkGet3WSMGX2ur7fWLgQWAsTHx9vaziHOs+dwMQnJKXyVdZQLerZh/uQ4otuE+HosEb/iyZH7BcBEY8x4IBhoCTwHhBtjmriP3qOAHM/HlEBQVe3ilc9388cPdhDUpBFPTBnI9fFdVB0gUgu1Dndr7SxgFoD7yP1/rbU3GWPeAa6l5oqZacBSz8cUp9t+8DgzE1PYkl3I2L4deGzyADq0DPb1WCJ+qy6uc58JvGmMeRT4Bni5DtYhDlFeVc3f1u7i+bWZtGrelL/cOJgJcZE6WhfxkFfC3Vq7Dljnvp0FDPfG1xVn27zvGDMTU9iZd4LJgzvz4IR+tA4N8vVYIo6gv1CVeldSUcXTq3bw6he7iWwZzKu3DOPSPu19PZaIoyjcpV59nnmYhOQU9h8tZerIrswY15swFX2JeJ3CXepFYWklj6/Yxlsb99OtbShvTR/JiO5tfD2WiGMp3KXOrdp6kDlL0jhSXMHdl/Tgt2NjCW6qoi+RuqRwlzqTX1TO3GVbWZGaS5+OYbw8bRgDo1r5eiyRgKBwF6+z1rL4mxweXp5OSXk1/3tFL+66pAdNG6sPRqS+KNzFq3IKSnlgcSrrMvIZEh3Ok9fG0bN92I+/UES8SuEuXuFyWV5fv5cF723HAnN/2o+po2JU9CXiIwp38VhW/gkSklL5es9RLopty+OTB9IlQkVfIr6kcJdaq6p28Y9Pd/PMhzsIbtKIp66N49qhUaoOEGkAFO5SK+kHjjMjaQtpOccZ178jD0/qT3sVfYk0GAp3OSdlldX89aNMXvx4F+EhQbxw0xCuHBjp67FE5BQKdzlrm/YeZUZiCrvyi5kyJIo5E/oSHqKiL5GGSOEuP6q4vIqnVmWw6Ms9dGrVnEW3DeeSXu18PZaI/ACFu/ygT3bkMys5lQOFpdw8siv3jetDi2b6thFp6PRTKqdVWFLJoyvSeWdTNt3bhfL2XaMYFhPh67FE5Cwp3OV73k/LZc7SrRwtruCXo3twzxgVfYn4G4W7/FdeURkPLd3Ke2kH6d+pJf+8dRj9O6noS8QfKdwFay1Jm3N4ZHk6pZXVzBjXmzsv6q6iLxE/pnAPcNnHSrh/cRqf7MhnWExrFkyJo0e7Fr4eS0Q8pHAPUC6X5bUv9/DkqgwM8PCk/vxiRFcaqehLxBEU7gEoM+8ECUkpbNx7jIt7tePxyQOIaq2iLxEnUbgHkMpqFws/yeK5D3cS0qwxf7xuENcM6ayiLxEHUrgHiLScQmYkppCee5zxAzsyb+IA2oU18/VYIlJHFO4OV1ZZzZ/X7OTvn2QRERrEi78YyrgBHX09lojUMYW7g23Yc5SZiSlkHS7muqFRzL6qH61Cmvp6LBGpBwp3BzpRXsWT72/ntS/3EtW6Of+6fTgXxaroSySQKNwdZl1GHg8sTuNAYSm3nB/DfT/pTaiKvkQCjn7qHaKgpIKHl6eTvDmHHu1CSbx7FEO7quhLJFAp3P2ctZb30g7y4NI0Ckoq+c1lPfn1ZT1p1kRFXyKBrNbhbozpArwGdAAssNBa+5wxJgJ4C4gB9gDXW2uPeT6qnCrveBlzlqaxaushBnZuxWu3jaBfp5a+HktEGgBPjtyrgD9YazcbY8KATcaY1cAtwBpr7QJjTAKQAMz0fFT5jrWWdzZl8+jydMqrXCRc2Yc7LuxGExV9iYhbrcPdWpsL5LpvFxljtgGdgUnAaPfTFgHrULh7zf6jJcxKTuWzzMMMj4lgwZSBdFfRl4icwivn3I0xMcBgYD3QwR38AAepOW1zutdMB6YDREdHe2MMR6t2WRZ9sYenVmXQyMAjk/pzk4q+ROQMPA53Y0wLIAn4rbX2+Mk9JdZaa4yxp3udtXYhsBAgPj7+tM+RGjsPFTEzKYXN+woY3bsdj00eSOfw5r4eS0QaMI/C3RjTlJpgf91am+xefMgYE2mtzTXGRAJ5ng4ZqCqrXby4bhd/+SiTkGaNefZn5zHpvE4q+hKRH+XJ1TIGeBnYZq3900kPLQOmAQvcn5d6NGGASs0u5L7ELWw/WMSEuEjmTuxP2xYq+hKRs+PJkfsFwFQg1RjzrXvZ/dSE+tvGmNuBvcD1Hk0YYMoqq3lm9Q7+8WkW7cKasXDqUK7or6IvETk3nlwt8xlwpvMDY2r7dQPZV1lHSEhKYc+REm4c3oWEK/vSqrmKvkTk3OkvVBuAorJKFry3ndfX76NLRHP+c8cIzu/Z1tdjiYgfU7j72EfbD/HA4jQOHS/jjgu78fsrehESpN0iIp5RivjI0eIKHn53K0u+PUCvDi14/qbzGRzd2tdjiYhDKNzrmbWWd1NymbtsK8dLK7l3TCy/urQnQU1UHSAi3qNwr0cHC8uYvSSND7cdYlBUK564cwR9OqroS0S8T+FeD6y1vLVhP4+t2Ealy8UD4/ty24XdaKzqABGpIwr3Orb3SDEJSal8mXWEkd0jWHBNHDFtQ309log4nMK9jlS7LK9+vpunP8igaaNGzL9mID+L76KiLxGpFwr3OpBxsIgZSSls2V/AmD7teXTyACJbqehLROqPwt2LKqpcPL8uk7+tzSQsuCl/vnEwP42LVNGXiNQ7hbuXfLu/gJmJKWQcKmLSeZ146Kf9iQgN8vVYIhKgFO4eKq2o5k+rM3j5s920Dwvm5WnxjOl72vcnERGpNwp3D3yx6zAJSansO1rCz0dEk3BlH1oGq+hLRHxP4V4Lx8sqmb9yG298vZ+YNiG8cedIRvVo4+uxRET+S+F+jj5MP8QDS1LJLypn+sXd+d3YXjQPauzrsURE/h+F+1k6cqKcue+m8+6WA/TpGMY/bo4nLirc12OJiJyWwv1HWGtZtuUAc5dt5UR5Fb+/vBd3X9JDRV8i0qAp3H/AgYJSZi9J46PteQyODufJKXHEdgjz9VgiIj9K4X4aLpfljQ37mL9yO9Uuy5wJ/bjl/BgVfYmI31C4n2LP4WISklP4KusoF/Rsw/zJcUS3CfH1WCIi50Th7lZV7eKVz3fzxw92ENSkEU9OieO6+ChVB4iIX1K4A9sPHmdGYgop2YVc3q8Dj149gA4tg309lohIrQV0uJdXVfO3jzJ5ft0uwkOa8refD2H8wI46WhcRvxew4b553zFmJqawM+8E1wzuzJwJ/Witoi8RcYiAC/eSiiqeXrWDV7/YTWTLYF69dRiX9m7v67FERLwqoML9s52HSUhOIftYKVNHdmXGuN6EqehLRBwoIMK9sLSSx1ak8/bGbLq1DeWt6SMZ0V1FXyLiXI4P91VbDzJnSRpHiiu4+5Ie/HZsLMFNVfQlIs7m2HDPLypn7rKtrEjNpU/HMF6eNoyBUa18PZaISL1wXLhba1n8TQ4PL0+npLya/72iF3dd0oOmjVX0JSKBw1HhnlNQyv3JqXy8I58h0eE8eW0cPdur6EtEAk+dhLsxZhzwHNAYeMlau6Au1vMdl8vy7/V7eeK97Vhg7k/7MXWUir5EJHB5PdyNMY2BvwGXA9nABmPMMmtturfXBbAr/wQJSSls2HOMi2Lb8vjkgXSJUNGXiAS2ujhyHw5kWmuzAIwxbwKTAK+H+9sb9jN7aRrBTRrx1LVxXDtURV8iIlA34d4Z2H/S/WxgxKlPMsZMB6YDREdH12pF3dqFMqZPe+ZN7E97FX2JiPyXz/5B1Vq7EFgIEB8fb2vzNYbFRDAsJsKrc4mIOEFdXB+YA3Q56X6Ue5mIiNSTugj3DUCsMaabMSYIuAFYVgfrERGRM/D6aRlrbZUx5tfAKmouhXzFWrvV2+sREZEzq5Nz7tbalcDKuvjaIiLy4/Q3+SIiDqRwFxFxIIW7iIgDKdxFRBzIWFurvx/y7hDG5AN7a/nytsBhL47jD7TNgUHbHBg82eau1tp2p3ugQYS7J4wxG6218b6eoz5pmwODtjkw1NU267SMiIgDKdxFRBzICeG+0NcD+IC2OTBomwNDnWyz359zFxGR73PCkbuIiJxC4S4i4kB+He7GmHHGmAxjTKYxJsHX89QFY0wXY8xaY0y6MWarMeZe9/IIY8xqY8xO9+fWvp7Vm4wxjY0x3xhjlrvvdzPGrHfv67fcddKOYYwJN8YkGmO2G2O2GWNGBcA+/p37ezrNGPOGMSbYafvZGPOKMSbPGJN20rLT7ldT48/ubU8xxgzxZN1+G+4nvRH3lUA/4EZjTD/fTlUnqoA/WGv7ASOBX7m3MwFYY62NBda47zvJvcC2k+4/ATxjre0JHANu98lUdec54H1rbR9gEDXb7th9bIzpDNwDxFtrB1BTD34DztvP/wTGnbLsTPv1SiDW/TEdeMGTFfttuHPSG3FbayuA796I21GstbnW2s3u20XU/NB3pmZbF7mftgi42icD1gFjTBRwFfCS+74BLgMS3U9x2va2Ai4GXgaw1lZYawtw8D52awI0N8Y0AUKAXBy2n621nwBHT1l8pv06CXjN1vgKCDfGRNZ23f4c7qd7I+7OPpqlXhhjYoDBwHqgg7U21/3QQaCDr+aqA88CMwCX+34boMBaW+W+77R93Q3IB151n4p6yRgTioP3sbU2B3ga2EdNqBcCm3D2fv7OmfarVzPNn8M9oBhjWgBJwG+ttcdPfszWXM/qiGtajTETgDxr7SZfz1KPmgBDgBestYOBYk45BeOkfQzgPs88iZr/sXUCQvn+6QvHq8v96s/hHjBvxG2MaUpNsL9urU12Lz703a9s7s95vprPyy4AJhpj9lBzqu0yas5Hh7t/fQfn7etsINtau959P5GasHfqPgYYC+y21uZbayuBZGr2vZP383fOtF+9mmn+HO4B8Ubc7vPNLwPbrLV/OumhZcA09+1pwNL6nq0uWGtnWWujrLUx1OzTj6y1NwFrgWvdT3PM9gJYaw8C+40xvd2LxgDpOHQfu+0DRhpjQtzf499ts2P380nOtF+XATe7r5oZCRSedPrm3Flr/fYDGA/sAHYBD/h6njraxgup+bUtBfjW/TGemvPQa4CdwIdAhK9nrYNtHw0sd9/uDnwNZALvAM18PZ+Xt/U8YKN7Py8BWjt9HwPzgO1AGvAvoJnT9jPwBjX/plBJzW9ot59pvwKGmisAdwGp1FxJVOt1q35ARMSB/Pm0jIiInIHCXUTEgRTuIiIOpHAXEXEghbuIiAMp3EVEHEjhLiLiQP8HmwhparpVdVEAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "raw_system_events, tensor_field, sessions = simulation.execute()\n",
    "simulation_result = pd.DataFrame(raw_system_events)\n",
    "\n",
    "simulation_result['a'].plot()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "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.8.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}