CRZ_PERF_Tables_MD-11/V4/MD11-ECON-CRZ v4.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import scipy as sp\n",
    "from IPython.display import display, Markdown\n",
    "\n",
    "skip_plots = True"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Define the bounds of our data:\n",
    "- Weights\n",
    "- Indexes corresponding to weights for easier lookup\n",
    "- Altitudes\n",
    "- Indexes corresponding to altitudes for easier lookup\n",
    "- Number of points in the linspaces used for fitting and plotting"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "weights         = [140,  150, 160, 170, 180, 190, 200, 210, 220, 230, 240, 250, 260, 270, 280, 290]\n",
    "weights_index   = [  0,    1,   2,   3,   4,   5,   6,   7,   8,   9,  10,  11,  12,  13,  14,  15]\n",
    "alts            = [250,  270, 290, 310, 330, 350, 370, 390, 410, 430]\n",
    "alts_index      = [  0,    1,   2,   3,   4,   5,   6,   7,   8,   9]\n",
    "ext_alts        = [110, 130, 150, 170, 190, 210, 230]\n",
    "ext_alts_index  = [  0,   1,   2,   3,   4,   5,   6]\n",
    "tip_tanks       = [ 60,   90]\n",
    "tip_tanks_index = [  0,    1]\n",
    "points_for_fit  = 1000"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Transform functions for altitudes, weights, and mach to their respective indices and back."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "def weight2index(weight: int) -> int:\n",
    "  \"\"\"Return the weight as a weight index.\"\"\"\n",
    "  return int((weight - 140) / 10)\n",
    "\n",
    "def index2weight(index: int) -> int:\n",
    "  \"\"\"Return the weight index as a weight.\"\"\"\n",
    "  return weights[index]\n",
    "\n",
    "def alt2index(alt: int) -> int:\n",
    "  \"\"\"Return the altitude as an altitude index.\"\"\"\n",
    "  return int((alt - 250) / 20)\n",
    "\n",
    "def index2alt(index: int) -> int:\n",
    "  \"\"\"Return the altitude index as an altitude.\"\"\"\n",
    "  return alts[index]\n",
    "\n",
    "def mach2index(mach: float) -> int:\n",
    "  \"\"\"Return the mach as a mach index\"\"\"\n",
    "  if mach == 60: return 0\n",
    "  if mach == 90: return 1\n",
    "  return -1\n",
    "\n",
    "def index2mach(index: int) -> float:\n",
    "  \"\"\"Return the mach index as a mach\"\"\"\n",
    "  return tip_tanks[index]\n",
    "\n",
    "def tip2mach(tip: int) -> float:\n",
    "  \"\"\"Return old name for tip\"\"\"\n",
    "  if tip == 60: return 0.85\n",
    "  if tip == 90: return 0.87\n",
    "  return -1"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Very basic IAS to TAS conversion based on a speed gain of 2% per 1000ft."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "def ias2tas(ias: int, alt: int) -> float:\n",
    "  \"\"\"Return TAS from IAS and altitude.\"\"\"\n",
    "  return ias + (alt/10.0 * (ias * 0.02))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "MMO mapping for altitudes based on the 320/365 limit"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "def alt2mmo60(alt: int) -> float:\n",
    "  \"\"\"Return MMO based on altitude <60% tip tank\"\"\"\n",
    "  dict = {\n",
    "    110: 0.587,\n",
    "    130: 0.608,\n",
    "    150: 0.631,\n",
    "    170: 0.655,\n",
    "    190: 0.680,\n",
    "    210: 0.706,\n",
    "    230: 0.733,\n",
    "    250: 0.761,\n",
    "    270: 0.791,\n",
    "    290: 0.822,\n",
    "    310: 0.850,\n",
    "    330: 0.850,\n",
    "    350: 0.850,\n",
    "    370: 0.850,\n",
    "    390: 0.850,\n",
    "    410: 0.850,\n",
    "    430: 0.850,\n",
    "  }\n",
    "  return dict[alt]\n",
    "\n",
    "def alt2mmo90(alt: int) -> float:\n",
    "  \"\"\"Return MMO based on altitude >90% tip tank\"\"\"\n",
    "  dict = {\n",
    "    110: 0.667,\n",
    "    130: 0.691,\n",
    "    150: 0.716,\n",
    "    170: 0.742,\n",
    "    190: 0.770,\n",
    "    210: 0.798,\n",
    "    230: 0.828,\n",
    "    250: 0.859,\n",
    "    270: 0.870,\n",
    "    290: 0.870,\n",
    "    310: 0.870,\n",
    "    330: 0.870,\n",
    "    350: 0.870,\n",
    "    370: 0.870,\n",
    "    390: 0.870,\n",
    "    410: 0.870,\n",
    "    430: 0.870,\n",
    "  }\n",
    "  return dict[alt]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Parsing function to read in the CRZ PERF data files and construct a key value list for all known tables.\n",
    "\n",
    "See `Format description.md` for a description of the file format."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "matrix = {}\n",
    "\n",
    "def parseFile(file_path: str) -> dict:\n",
    "  \"\"\"\n",
    "  Parse file and return as matrix.\n",
    "\n",
    "  Matrix has altitude keys and (Fuel Flow, IAS, mach) as values\n",
    "  \"\"\"\n",
    "  matrix = {}\n",
    "  file = open(file_path, 'r')\n",
    "  lines = file.readlines()\n",
    "  \n",
    "  count = 0\n",
    "  for line in lines:\n",
    "    count += 1\n",
    "    if count == 1:\n",
    "      continue\n",
    "    fields = line.split(\";\")\n",
    "    alt = int(fields[0]) / 100\n",
    "    ff = [int(item.strip()) for item in fields[1].split(\",\")]\n",
    "    ias = [int(item.strip()) for item in fields[2].split(\",\")]\n",
    "    mach = [float(item.strip()) for item in fields[3].split(\",\")]\n",
    "\n",
    "    inner = {}\n",
    "    for i in range(len(ff)):\n",
    "      _ias = ias[i] if len(ias) > 1 else ias[0]\n",
    "      _mach = mach[i] if len(mach) > 1 else mach[0]\n",
    "      inner[weights[i]] = (ff[i], _ias, _mach)\n",
    "    matrix[alt] = inner\n",
    "  return matrix\n",
    "\n",
    "def parseCSV(file_path: str) -> tuple[np.ndarray, np.ndarray]:\n",
    "  x = []\n",
    "  y = []\n",
    "\n",
    "  file = open(file_path, 'r')\n",
    "  lines = file.readlines()\n",
    "  for line in lines:\n",
    "    fields = line.split(\",\")\n",
    "    average = (float(fields[0]) + float(fields[1])) / 2\n",
    "    # We don't really care for anything other than 90 to 60\n",
    "    # so we filter them out to have less data to slog through\n",
    "    if average > 91 or average < 59: continue\n",
    "    x.append(float(fields[2]))\n",
    "    y.append(average)\n",
    "  return np.array(x), np.array(y)\n",
    "\n",
    "matrix[\"LRC\"] = parseFile(\"../data/LRC.data\")\n",
    "matrix[82] = parseFile(\"../data/.82.data\")\n",
    "matrix[83] = parseFile(\"../data/.83.data\")\n",
    "matrix[84] = parseFile(\"../data/.84.data\")\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Basic parametrized functions used for the fitting."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "def func_quad(x: np.ndarray, a: float, b: float, c: float) -> np.ndarray:\n",
    "  \"\"\"Return f(x) = ax^2 + bx + c.\"\"\"\n",
    "  return a * x**2 + b * x + c\n",
    "\n",
    "def func_root(x: np.ndarray, a: float, b: float, c: float) -> np.ndarray:\n",
    "  \"\"\"Return f(x) = a*sqrt(x) + bx + c.\"\"\"\n",
    "  return a * np.sqrt(x) + b * x + c\n",
    "\n",
    "def func_lin(x: np.ndarray, a: float, b: float) -> np.ndarray:\n",
    "  \"\"\"Return f(x) = a*x + b.\"\"\"\n",
    "  return a * x + b"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Quick lookup of $LRC mach$ from our matrix."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "def get_lrc_mach(weight: int, alt: int) -> float:\n",
    "  \"\"\"Return LRC mach at given altitude and weight.\"\"\"\n",
    "  return matrix[\"LRC\"][alt][weight][2]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Generate \"base\" lists.\n",
    "\n",
    "Base lists are extracted from the matrix for a given weight, altitude, and expected $MRC mach$.\n",
    "\n",
    "Returned values are base lists for mach, mileage, colour (for the plot). Additional returned value is the $MRC mileage$ and $mMRC mach$.\n",
    "\n",
    "$MRC mileage$ is calculated as the value required to get $LRC mileage$ if $LRC mileage$ is a $1\\%$ drop from $MRC mileage$."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "def generate_base_lists(weight: int, alt: int, mrc_mach: float) -> tuple[np.ndarray, np.ndarray, list[str], float, float]:\n",
    "  \"\"\"\n",
    "  For a given weight, altitude, and MRC mach, return tuple containing\n",
    "  - mach_base: [MRC mach, LRC mach, .82, .83, .84]\n",
    "  - mileage_base: [MRC mileage, LRC mileage, .82 mileage, .83 mileage, .84 mileage]\n",
    "  - col: Colour list for plotting, green is MRC; orange is LRC; .82, .83, .84 are blue; unknown are red\n",
    "  - mrc_mileage: Calculated MRC mileage\n",
    "  - mrc_mach: Pass through of MRC mach\n",
    "  \"\"\"\n",
    "  raw = []\n",
    "  col = []\n",
    "  try:\n",
    "    for spd in range(82, 85):\n",
    "      pt = matrix[spd][alt][weight]\n",
    "      raw.append((pt[0], ias2tas(pt[1], alt), pt[2]))\n",
    "  except KeyError:\n",
    "    pass\n",
    "\n",
    "  lrc = matrix[\"LRC\"][alt][weight]\n",
    "  raw.append((lrc[0], ias2tas(lrc[1], alt), lrc[2], \"LRC\"))\n",
    "  raw.sort(key=lambda i: i[2])\n",
    "\n",
    "  spd_base = np.array([x[1] for x in raw])\n",
    "  ff_base = np.array([x[0] * 3 for x in raw])\n",
    "  mileage_base = spd_base / ff_base * 1000\n",
    "  index = next((i for i, item in enumerate(raw) if len(item) == 4 and item[3] == 'LRC'), -10)\n",
    "  lrc_mileage = mileage_base[index]\n",
    "  mrc_mileage = lrc_mileage * 100.0/99.0\n",
    "\n",
    "  raw.append((-1, -1, mrc_mach, \"MRC\"))\n",
    "  raw.sort(key=lambda i: i[2])\n",
    "  for i in range(len(raw)):\n",
    "     if len(raw[i]) == 4 and raw[i][3] == \"MRC\":\n",
    "        mileage_base = np.insert(mileage_base, i, mrc_mileage)\n",
    "        break\n",
    "\n",
    "  for x in raw:\n",
    "    if len(x) < 4:\n",
    "      col.append(\"blue\")\n",
    "    elif x[3] == \"MRC\":\n",
    "      col.append(\"green\")\n",
    "    elif x[3] == \"LRC\":\n",
    "      col.append(\"orange\")\n",
    "    else:\n",
    "      col.append(\"red\")\n",
    "\n",
    "  mach_base = np.array([x[2] for x in raw])\n",
    "  return mach_base, mileage_base, col, mrc_mileage, mrc_mach"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Iteratively generate potential $MRC mach$ numbers for a given weight and altitude.\n",
    "\n",
    "Use $MRC mach$, $LRC mach$, $.82$, $.83$, and $.84$ as the $x$ and mileage as the $y$ to fit a quadratic function.\n",
    "\n",
    "Returned values are the $MRC mach$ value that has the lowest absolute error to the `max()` of the fitted function, the x axis linspace, and the corresponding fitted function.\n",
    "\n",
    "Optionally plot all generate curves.\n",
    "\n",
    "Optionally print altitude, weight, $MRC mach$, calculated $MRC mileage$, `max()` $MRC mileage$, and the absolute difference of the latter two."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "def find_best_mrc_mach(weight: int, alt: int, show_intermediate_graphs: bool = False, debug_output: bool = False) -> tuple[float, np.ndarray, np.ndarray]:\n",
    "  \"\"\"\n",
    "  For a given weight and altitude, return tuple containing\n",
    "  - best_mrc_mach: MRC mach closest to the max() of the mach/mileage curve\n",
    "  - x: x axis values for the mach/mileage curve\n",
    "  - best_func: mach/mileage curve\n",
    "\n",
    "  Optional parameters:\n",
    "  - show_intermediate_graphs: Default False. Used to plot all curves.\n",
    "  - debug_output: Default False. Used for debug\n",
    "  \"\"\"\n",
    "  lrc_mach = get_lrc_mach(weight, alt)\n",
    "  rng = np.linspace(lrc_mach - 0.1, lrc_mach, 100, False)\n",
    "  best_mrc_mach = 0\n",
    "  best_func = 0\n",
    "  last_mrc_mileage_error = 10000\n",
    "\n",
    "  for mrc_mach in rng:\n",
    "    mach, mileage, colour, mrc_mileage, _ = generate_base_lists(weight, alt, mrc_mach)\n",
    "    x = np.linspace(mrc_mach - 0.1, 0.95, points_for_fit)\n",
    "\n",
    "    popt, _ = sp.optimize.curve_fit(func_quad, mach, mileage, maxfev=100000)\n",
    "    mileage_mach_func = func_quad(x, *popt)\n",
    "\n",
    "    if debug_output:\n",
    "      print(f\"{alt=}, {weight=}, {mrc_mach=}\", mrc_mileage, max(mileage_mach_func), abs(mrc_mileage - max(mileage_mach_func)))\n",
    "\n",
    "    if abs(mrc_mileage - max(mileage_mach_func)) < last_mrc_mileage_error:\n",
    "      last_mrc_mileage_error = abs(mrc_mileage - max(mileage_mach_func))\n",
    "      best_mrc_mach = round(mrc_mach, 3)\n",
    "      best_func = mileage_mach_func\n",
    "    else:\n",
    "      return best_mrc_mach, x, best_func\n",
    "\n",
    "    if show_intermediate_graphs and not skip_plots:\n",
    "      plt.figure(figsize=(20,10))\n",
    "      plt.title(f\"mach/Mileage for {weight=} and {mrc_mach=}\")\n",
    "      plt.plot(x, mileage_mach_func)\n",
    "      plt.scatter(mach, mileage, c=colour)\n",
    "      plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Generate all curves for all weight and altitude combinations.\n",
    "\n",
    "Use the mach/mileage curves to generate CI/mach curves for all weight and altitude combinations by fitting a root function using $0$, $200$, and $999$ as the $x$ and $MRC mach$, $LRC mach$, and $0.85$ as the corresponding $y$."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/markdown": [
       "**Found: 0.534 MRC and 0.615 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "C:\\Users\\llego\\AppData\\Local\\Packages\\PythonSoftwareFoundation.Python.3.11_qbz5n2kfra8p0\\LocalCache\\local-packages\\Python311\\site-packages\\scipy\\optimize\\_minpack_py.py:1010: OptimizeWarning: Covariance of the parameters could not be estimated\n",
      "  warnings.warn('Covariance of the parameters could not be estimated',\n"
     ]
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.534 MRC and 0.615 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.55 MRC and 0.633 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.55 MRC and 0.633 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.57 MRC and 0.653 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.57 MRC and 0.653 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.591 MRC and 0.674 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.591 MRC and 0.674 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.612 MRC and 0.694 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.612 MRC and 0.694 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.634 MRC and 0.714 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.634 MRC and 0.714 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.659 MRC and 0.734 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.659 MRC and 0.734 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.695 MRC and 0.758 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.695 MRC and 0.758 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.732 MRC and 0.781 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.732 MRC and 0.781 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.762 MRC and 0.799 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.762 MRC and 0.799 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.783 MRC and 0.812 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.783 MRC and 0.812 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.798 MRC and 0.821 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.798 MRC and 0.821 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.801 MRC and 0.824 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.801 MRC and 0.824 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.806 MRC and 0.827 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.806 MRC and 0.827 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.81 MRC and 0.829 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.81 MRC and 0.829 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.812 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.812 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.555 MRC and 0.639 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.555 MRC and 0.639 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.577 MRC and 0.661 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.577 MRC and 0.661 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.598 MRC and 0.683 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.598 MRC and 0.683 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.623 MRC and 0.705 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.623 MRC and 0.705 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.647 MRC and 0.727 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.647 MRC and 0.727 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.683 MRC and 0.752 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.683 MRC and 0.752 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.725 MRC and 0.778 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.725 MRC and 0.778 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.759 MRC and 0.798 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.759 MRC and 0.798 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.782 MRC and 0.812 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.782 MRC and 0.812 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.798 MRC and 0.821 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.798 MRC and 0.821 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.802 MRC and 0.825 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.802 MRC and 0.825 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.807 MRC and 0.828 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.807 MRC and 0.828 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.81 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.81 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.81 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.81 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.81 MRC and 0.831 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.81 MRC and 0.831 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.814 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.814 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.583 MRC and 0.668 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.583 MRC and 0.668 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.606 MRC and 0.692 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.606 MRC and 0.692 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.632 MRC and 0.716 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.632 MRC and 0.716 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.667 MRC and 0.742 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.667 MRC and 0.742 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.711 MRC and 0.77 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.711 MRC and 0.77 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.75 MRC and 0.793 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.75 MRC and 0.793 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.78 MRC and 0.811 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.78 MRC and 0.811 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.795 MRC and 0.821 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.795 MRC and 0.821 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.802 MRC and 0.825 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.802 MRC and 0.825 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.803 MRC and 0.828 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.803 MRC and 0.828 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.809 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.809 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.809 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.809 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.807 MRC and 0.831 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.807 MRC and 0.831 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.811 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.811 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.811 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.811 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.816 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.816 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.614 MRC and 0.7 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.614 MRC and 0.7 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.647 MRC and 0.727 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.647 MRC and 0.727 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.691 MRC and 0.757 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.691 MRC and 0.757 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.737 MRC and 0.785 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.737 MRC and 0.785 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.776 MRC and 0.808 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.776 MRC and 0.808 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.795 MRC and 0.819 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.795 MRC and 0.819 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.802 MRC and 0.825 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.802 MRC and 0.825 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.809 MRC and 0.828 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.809 MRC and 0.828 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.811 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.811 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.811 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.811 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.813 MRC and 0.831 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.813 MRC and 0.831 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.812 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.812 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.818 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.818 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.813 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.813 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.821 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.821 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.821 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.821 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.66 MRC and 0.737 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.66 MRC and 0.737 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.712 MRC and 0.771 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.712 MRC and 0.771 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.758 MRC and 0.798 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.758 MRC and 0.798 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.786 MRC and 0.815 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.786 MRC and 0.815 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.798 MRC and 0.824 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.798 MRC and 0.824 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.805 MRC and 0.828 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.805 MRC and 0.828 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.808 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.808 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.807 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.807 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.805 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.805 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.81 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.81 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.806 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.806 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.808 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.808 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.811 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.811 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.807 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.807 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.816 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.816 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.729 MRC and 0.782 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.729 MRC and 0.782 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.776 MRC and 0.809 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.776 MRC and 0.809 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.796 MRC and 0.822 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.796 MRC and 0.822 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.801 MRC and 0.827 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.801 MRC and 0.827 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.805 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.805 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.806 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.806 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.802 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.802 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.807 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.807 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.802 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.802 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.805 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.805 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.811 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.811 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.8 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.8 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.8 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.8 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.785 MRC and 0.814 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.785 MRC and 0.814 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.798 MRC and 0.824 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.798 MRC and 0.824 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.808 MRC and 0.829 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.808 MRC and 0.829 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.806 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.806 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.81 MRC and 0.831 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.81 MRC and 0.831 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.809 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.809 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.813 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.813 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.82 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.82 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.81 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.81 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.815 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.815 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.811 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.811 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.804 MRC and 0.827 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.804 MRC and 0.827 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.809 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.809 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.804 MRC and 0.831 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.804 MRC and 0.831 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.81 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.81 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.807 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.807 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.809 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.809 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.818 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.818 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.804 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.804 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.801 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.801 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.807 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.807 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.804 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.804 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.81 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.81 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.8 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.8 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.805 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.805 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.8 MRC and 0.838 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.8 MRC and 0.838 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.804 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.804 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.798 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.798 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.802 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.802 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.791 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.791 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.788 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.788 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.534 MRC and 0.615 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.534 MRC and 0.615 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.55 MRC and 0.633 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.55 MRC and 0.633 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.57 MRC and 0.653 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.57 MRC and 0.653 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.591 MRC and 0.674 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.591 MRC and 0.674 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.612 MRC and 0.694 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.612 MRC and 0.694 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.634 MRC and 0.714 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.634 MRC and 0.714 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.659 MRC and 0.734 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.659 MRC and 0.734 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.695 MRC and 0.758 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.695 MRC and 0.758 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.732 MRC and 0.781 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.732 MRC and 0.781 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.762 MRC and 0.799 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.762 MRC and 0.799 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.783 MRC and 0.812 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.783 MRC and 0.812 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.798 MRC and 0.821 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.798 MRC and 0.821 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.801 MRC and 0.824 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.801 MRC and 0.824 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.806 MRC and 0.827 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.806 MRC and 0.827 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.81 MRC and 0.829 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.81 MRC and 0.829 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.812 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.812 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.555 MRC and 0.639 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.555 MRC and 0.639 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.577 MRC and 0.661 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.577 MRC and 0.661 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.598 MRC and 0.683 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.598 MRC and 0.683 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.623 MRC and 0.705 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.623 MRC and 0.705 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.647 MRC and 0.727 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.647 MRC and 0.727 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.683 MRC and 0.752 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.683 MRC and 0.752 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.725 MRC and 0.778 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.725 MRC and 0.778 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.759 MRC and 0.798 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.759 MRC and 0.798 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.782 MRC and 0.812 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.782 MRC and 0.812 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.798 MRC and 0.821 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.798 MRC and 0.821 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.802 MRC and 0.825 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.802 MRC and 0.825 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.807 MRC and 0.828 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.807 MRC and 0.828 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.81 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.81 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.81 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.81 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.81 MRC and 0.831 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.81 MRC and 0.831 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.814 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.814 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.583 MRC and 0.668 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.583 MRC and 0.668 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.606 MRC and 0.692 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.606 MRC and 0.692 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.632 MRC and 0.716 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.632 MRC and 0.716 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.667 MRC and 0.742 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.667 MRC and 0.742 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.711 MRC and 0.77 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.711 MRC and 0.77 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.75 MRC and 0.793 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.75 MRC and 0.793 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.78 MRC and 0.811 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.78 MRC and 0.811 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.795 MRC and 0.821 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.795 MRC and 0.821 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.802 MRC and 0.825 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.802 MRC and 0.825 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.803 MRC and 0.828 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.803 MRC and 0.828 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.809 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.809 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.809 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.809 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.807 MRC and 0.831 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.807 MRC and 0.831 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.811 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.811 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.811 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.811 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.816 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.816 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.614 MRC and 0.7 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.614 MRC and 0.7 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.647 MRC and 0.727 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.647 MRC and 0.727 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.691 MRC and 0.757 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.691 MRC and 0.757 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.737 MRC and 0.785 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.737 MRC and 0.785 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.776 MRC and 0.808 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.776 MRC and 0.808 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.795 MRC and 0.819 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.795 MRC and 0.819 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.802 MRC and 0.825 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.802 MRC and 0.825 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.809 MRC and 0.828 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.809 MRC and 0.828 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.811 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.811 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.811 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.811 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.813 MRC and 0.831 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.813 MRC and 0.831 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.812 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.812 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.818 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.818 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.813 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.813 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.821 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.821 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.821 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.821 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.66 MRC and 0.737 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.66 MRC and 0.737 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.712 MRC and 0.771 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.712 MRC and 0.771 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.758 MRC and 0.798 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.758 MRC and 0.798 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.786 MRC and 0.815 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.786 MRC and 0.815 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.798 MRC and 0.824 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.798 MRC and 0.824 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.805 MRC and 0.828 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.805 MRC and 0.828 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.808 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.808 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.807 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.807 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.805 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.805 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.81 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.81 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.806 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.806 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.808 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.808 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.811 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.811 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.807 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.807 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.816 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.816 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.729 MRC and 0.782 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.729 MRC and 0.782 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.776 MRC and 0.809 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.776 MRC and 0.809 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.796 MRC and 0.822 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.796 MRC and 0.822 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.801 MRC and 0.827 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.801 MRC and 0.827 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.805 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.805 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.806 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.806 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.802 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.802 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.807 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.807 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.802 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.802 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.805 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.805 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.811 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.811 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.8 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.8 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.8 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.8 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.785 MRC and 0.814 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.785 MRC and 0.814 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.798 MRC and 0.824 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.798 MRC and 0.824 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.808 MRC and 0.829 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.808 MRC and 0.829 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.806 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.806 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.81 MRC and 0.831 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.81 MRC and 0.831 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.809 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.809 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.813 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.813 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.82 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.82 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.81 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.81 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.815 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.815 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.811 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.811 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.804 MRC and 0.827 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.804 MRC and 0.827 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.809 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.809 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.804 MRC and 0.831 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.804 MRC and 0.831 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.81 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.81 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.807 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.807 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.809 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.809 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.818 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.818 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.804 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.804 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.801 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.801 MRC and 0.83 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.807 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.807 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.804 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.804 MRC and 0.833 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.81 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.81 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.8 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.8 MRC and 0.836 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.805 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.805 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.8 MRC and 0.838 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.8 MRC and 0.838 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.804 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.804 MRC and 0.832 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.798 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.798 MRC and 0.834 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.802 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.802 MRC and 0.835 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.791 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.791 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Found: 0.788 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/markdown": [
       "**Sanity check: 0.788 MRC and 0.837 LRC**"
      ],
      "text/plain": [
       "<IPython.core.display.Markdown object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "curves = []\n",
    "for tips in tip_tanks:\n",
    "  _curves = []\n",
    "  for alt in alts:\n",
    "    __curves = []\n",
    "    for weight in weights:\n",
    "      try:\n",
    "        lrc_mach = get_lrc_mach(weight, alt)\n",
    "        best_mrc_mach, x, mileage_mach_func = find_best_mrc_mach(weight, alt)\n",
    "        display(Markdown(f\"**Found: {best_mrc_mach} MRC and {lrc_mach} LRC**\"))\n",
    "\n",
    "        _mach, mileage, colour, mrc_mileage, mrc_mach = generate_base_lists(weight, alt, best_mrc_mach)\n",
    "\n",
    "        if not skip_plots:\n",
    "          plt.figure(figsize=(20,10))\n",
    "\n",
    "          plt.subplot(1, 2, 1)\n",
    "          plt.title(f\"mach/Mileage for {alt=}, {weight=}, {mrc_mach=}\")\n",
    "          plt.xlabel(\"mach\")\n",
    "          plt.ylabel(\"Mileage in nmi/1000kg\")\n",
    "          plt.plot(x, mileage_mach_func)\n",
    "          plt.scatter(_mach, mileage, c=colour)\n",
    "\n",
    "        mmo = alt2mmo60(alt) if tips == 60 else alt2mmo90(alt)\n",
    "        if lrc_mach >= mmo:\n",
    "          mmo = lrc_mach\n",
    "        ci_mach_base = [best_mrc_mach, lrc_mach, mmo]\n",
    "        ci_base = [0, 200, 999]\n",
    "        x = np.linspace(0, 999, points_for_fit)\n",
    "\n",
    "        popt, _ = sp.optimize.curve_fit(func_root, ci_base, ci_mach_base, maxfev=10000)\n",
    "        ci_mach_func = func_root(x, *popt)\n",
    "\n",
    "        if not skip_plots:\n",
    "          plt.subplot(1, 2, 2)\n",
    "          plt.title(f\"CI/mach for {alt=}, {weight=}, {mrc_mach=}, {lrc_mach=} {tips=}%\\nwill be truncated during output to MMO if above\")\n",
    "          plt.xlabel(\"CI\")\n",
    "          plt.ylabel(\"mach\")\n",
    "          plt.scatter(ci_base, ci_mach_base)\n",
    "          plt.plot(x, ci_mach_func)\n",
    "          plt.show()\n",
    "\n",
    "        display(Markdown(f\"**Sanity check: {ci_mach_func[0]} MRC and {round(ci_mach_func[int(points_for_fit / 1000) * 200], 3)} LRC**\"))\n",
    "\n",
    "        __curves.append((ci_mach_func, popt))\n",
    "      except KeyError:\n",
    "        pass\n",
    "    _curves.append(__curves)\n",
    "  curves.append(_curves)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Diagnostic plots to see how the CI/mach curves change with altitude for a given weight."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "if not skip_plots:\n",
    "  x = np.linspace(0, 999, points_for_fit)\n",
    "\n",
    "  for tips in tip_tanks_index:\n",
    "    for weight in weights_index:\n",
    "      plt.figure(figsize=(20,10))\n",
    "      plt.title(f\"CI/mach for weight={index2weight(weight)}, tip={index2mach(tips)}%\")\n",
    "      plt.xlabel(\"CI\")\n",
    "      plt.ylabel(\"mach\")\n",
    "      for alt in alts_index:\n",
    "        try:\n",
    "          plt.plot(x, curves[tips][alt][weight][0], label=f\"alt={index2alt(alt)}\")\n",
    "        except IndexError:\n",
    "          pass\n",
    "      plt.legend()\n",
    "      plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Diagnostic plot to see how linear $MRC mach$ is for all altitudes at a given weight."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "if not skip_plots:\n",
    "  x = np.linspace(250, 430, 10)\n",
    "\n",
    "  for tips in tip_tanks_index:\n",
    "    for weight in weights_index:\n",
    "      plt.title(f\"MRC at weight={index2weight(weight)}, tips={index2mach(tips)}%\")\n",
    "      plt.xlabel(\"Altitude in FL\")\n",
    "      plt.ylabel(\"mach\")\n",
    "      for alt in alts_index:\n",
    "        try:\n",
    "          plt.scatter(x[alt], curves[tips][alt][weight][0][0], label=f\"alt={index2alt(alt)}\")\n",
    "        except IndexError:\n",
    "          pass\n",
    "      plt.legend()\n",
    "      plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Total fuel to tip tank function"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjIAAAHHCAYAAACle7JuAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAB0S0lEQVR4nO3ddXRU19rH8e/EQ0KCJZCgQYpTnCItEtwrUEoFKxUkxaUttEBxtwotDTUqVHAp7lZIgBQpLsE1wUIyc94/+pLbITOQQGyS32etrHtnnzPP7M1wyNOtJsMwDEREREQckFNaV0BERETkcSmREREREYelREZEREQclhIZERERcVhKZERERMRhKZERERERh6VERkRERByWEhkRERFxWEpkRERExGEpkRFJA+vWrcNkMrFu3bpH3rtz505q1KiBl5cXJpOJ8PDwFKlTnTp1qFOnTorEfhIff/wxJpMpravxWE6cOIHJZGLChAlpXRWRDEuJjGQaJpMpUT+JSS5GjRrF/PnzU7zOsbGxtGnThqtXrzJ58mS+++47ChYsmOKf+zAp0fbbt2/z8ccfJ+rP/kmdPXuWjz/+OMUSQkma1HqWJONySesKiKSW7777zur1t99+y8qVKxOUlyxZ8pGxRo0axUsvvUTr1q2Ts4oJHD16lJMnT/Lll1/y5ptvpuhnJVZKtP327dsMGzYMIEGv0IcffsigQYOS7bPOnj3LsGHDKFSoEOXLl0+2uPJ4UutZkoxLiYxkGq+99prV623btrFy5coE5enJxYsXAciWLVvaViQNubi44OKif6pSyq1bt/Dy8krraog8Ng0tifzHrVu36Nu3L/nz58fd3Z3ixYszYcIE/ntIvMlk4tatW3zzzTfxw1EdO3YE4OTJk3Tr1o3ixYvj6elJzpw5adOmDSdOnEhyXTp27Ejt2rUBaNOmDSaTKb63wt58lo4dO1KoUCGrMovFwpQpUyhdujQeHh7kzp2bt99+m2vXriW5Tg9rO0BYWBhNmjTBx8cHb29vgoOD2bZt20NjnjhxAj8/PwCGDRsWH/fjjz8GbM+RMZlM9OjRgx9++IHixYvj4eFBpUqV2LBhw0M/a926dVSpUgWATp06xX/WnDlzANi4cSNt2rShQIECuLu7kz9/fnr37s2dO3es4nTs2BFvb28iIyNp3bo13t7e+Pn50a9fP8xm80PrYBgGb731Fm5ubvz+++8P/XO5P79m8uTJFCxYEE9PT2rXrk1ERESC+w8ePMhLL71Ejhw58PDwoHLlyixcuNDqnjlz5mAymVi/fj3dunXD39+ffPnyxV9ftmwZtWvXJmvWrPj4+FClShXmzp1rFWP79u00btwYX19fsmTJQu3atdm8ebPVPfe/syNHjtCxY0eyZcuGr68vnTp14vbt2/H3Pervk0hi6D9zRP6fYRi0bNmStWvX0qVLF8qXL8+KFSvo378/kZGRTJ48Gfh3iOrNN9+katWqvPXWWwAUKVIE+Hdi7pYtW2jXrh358uXjxIkTfPbZZ9SpU4f9+/eTJUuWRNfn7bffJm/evIwaNYqQkBCqVKlC7ty5k9yut99+mzlz5tCpUydCQkI4fvw4M2bMICwsjM2bN+Pq6proWA9r+99//82zzz6Lj48PAwYMwNXVlS+++II6deqwfv16qlWrZjOmn58fn332Ge+++y7PP/88L7zwAgDlypV7aF3Wr1/Pzz//TEhICO7u7nz66ac0btyYHTt2UKZMGZvvKVmyJMOHD2fo0KG89dZbPPvsswDUqFEDgHnz5nH79m3effddcubMyY4dO5g+fTpnzpxh3rx5VrHMZjONGjWiWrVqTJgwgVWrVjFx4kSKFCnCu+++a/PzzWYznTt35ueff+aPP/6gWbNmD20j/DsEGh0dTffu3bl79y5Tp06lXr167Nu3L/7vw99//03NmjXJmzcvgwYNwsvLi19++YXWrVvz22+/8fzzz1vF7NatG35+fgwdOpRbt24B/yY5nTt3pnTp0gwePJhs2bIRFhbG8uXLad++PQBr1qyhSZMmVKpUiY8++ggnJydCQ0OpV68eGzdupGrVqlaf07ZtW4KCghg9ejS7d+/mq6++wt/fn7FjxwIP//skkmiGSCbVvXt347+PwPz58w3A+OSTT6zue+mllwyTyWQcOXIkvszLy8vo0KFDgpi3b99OULZ161YDML799tv4srVr1xqAsXbt2ofW8f598+bNsyqvXbu2Ubt27QT3d+jQwShYsGD8640bNxqA8cMPP1jdt3z58gTl9mI+yF7bW7dubbi5uRlHjx6NLzt79qyRNWtW47nnnntozEuXLhmA8dFHHyW49tFHHxkP/lMFGIDx119/xZedPHnS8PDwMJ5//vmHftbOnTsNwAgNDU1wzdb3N3r0aMNkMhknT56ML+vQoYMBGMOHD7e6t0KFCkalSpXiXx8/ftwAjPHjxxuxsbHGyy+/bHh6ehorVqx4aB3/+15PT0/jzJkz8eXbt283AKN3797xZcHBwUbZsmWNu3fvxpdZLBajRo0aRrFixeLLQkNDDcCoVauWERcXF19+/fp1I2vWrEa1atWMO3fuWNXDYrHE/2+xYsWMRo0axZfd/zMLCgoyGjRoEF92/zvr3LmzVaznn3/eyJkzp1WZvb9PIomloSWR/7d06VKcnZ0JCQmxKu/bty+GYbBs2bJHxvD09Iz//7GxsVy5coWiRYuSLVs2du/enex1fpR58+bh6+tLgwYNuHz5cvxPpUqV8Pb2Zu3atcnyOWazmT///JPWrVtTuHDh+PKAgADat2/Ppk2biIqKSpbPuq969epUqlQp/nWBAgVo1aoVK1aseOTwjj3//f5u3brF5cuXqVGjBoZhEBYWluD+d955x+r1s88+y7FjxxLcd+/ePdq0acPixYtZunQpDRs2THSdWrduTd68eeNfV61alWrVqrF06VIArl69ypo1a2jbti3R0dHx3/GVK1do1KgRhw8fJjIy0ipm165dcXZ2jn+9cuVKoqOjGTRoEB4eHlb33h/WCw8P5/Dhw7Rv354rV67Ef86tW7cIDg5mw4YNWCyWR/75XLlyJdn/LkjmpqElkf938uRJAgMDyZo1q1X5/VVMJ0+efGSMO3fuMHr0aEJDQ4mMjLSaW3Pjxo3krXAiHD58mBs3buDv72/z+v3JxE/q0qVL3L59m+LFiye4VrJkSSwWC6dPn6Z06dLJ8nkAxYoVS1D21FNPcfv2bS5dukSePHmSHPPUqVMMHTqUhQsXJphD9OD35+HhET+3577s2bPbnHs0evRobt68ybJly5K8V4+9dv7yyy8AHDlyBMMwGDJkCEOGDLEZ4+LFi1bJUFBQkNX1o0ePAtgdkoN//y4BdOjQwe49N27cIHv27PGvCxQoYHX9/rVr167h4+NjN45IUiiREUlGPXv2JDQ0lF69elG9enV8fX0xmUy0a9cuwX+tPgmTyWSVJN33YE+ExWLB39+fH374wWacB38RZ2Zms5kGDRpw9epVBg4cSIkSJfDy8iIyMpKOHTsm+P7+26PxKI0aNWL58uWMGzeOOnXqJOj1eBL369WvXz8aNWpk856iRYtavf5vz1NSP2f8+PF2l617e3tbvbb3Z2Tr767I41IiI/L/ChYsyKpVq4iOjrbqlTl48GD89fvs7TT766+/0qFDByZOnBhfdvfuXa5fv56sdc2ePbvNIYwHe42KFCnCqlWrqFmz5mP98rLFVtv9/PzIkiULhw4dSnDt4MGDODk5kT9//iTFfJT7PQT/9c8//5AlS5aHJmj2Pmvfvn38888/fPPNN7zxxhvx5StXrkxy3R70zDPP8M4779C8eXPatGnDH3/8kegl5fbaeX912v2hPFdXV+rXr/9Y9bs/wTYiIiJB0vPgPT4+Po/9ObY46q7Nkn5ojozI/2vatClms5kZM2ZYlU+ePBmTyUSTJk3iy7y8vGwmJ87Ozgn+a3P69OmPPWfDniJFinDw4EEuXboUX7Znz54Ey2Dbtm2L2WxmxIgRCWLExcU9VoJlq+3Ozs40bNiQBQsWWC01v3DhAnPnzqVWrVoPHUq4v5orKfXZunWr1byj06dPs2DBAho2bPjQ3pL7e6bYagNY9xYYhsHUqVMTXaeHqV+/Pj/99BPLly/n9ddfT3QP3fz5863muOzYsYPt27fH/3309/enTp06fPHFF5w7dy7B+//7d8Sehg0bkjVrVkaPHs3du3etrt3/86hUqRJFihRhwoQJ3Lx587E+xxZ7z5JIYqlHRuT/tWjRgrp16/LBBx9w4sQJnn76af78808WLFhAr169rJaFVqpUiVWrVjFp0iQCAwMJCgqiWrVqNG/enO+++w5fX19KlSrF1q1bWbVqFTlz5kzWunbu3JlJkybRqFEjunTpwsWLF/n8888pXbq01UTK2rVr8/bbbzN69GjCw8Np2LAhrq6uHD58mHnz5jF16lReeumlJH22vbZ/8sknrFy5klq1atGtWzdcXFz44osviImJYdy4cQ+N6enpSalSpfj555956qmnyJEjB2XKlHnonI0yZcrQqFEjq+XXQPwOwfYUKVKEbNmy8fnnn5M1a1a8vLyoVq0aJUqUoEiRIvTr14/IyEh8fHz47bffHmu/HXtat25NaGgob7zxBj4+PnzxxRePfE/RokWpVasW7777LjExMUyZMoWcOXMyYMCA+HtmzpxJrVq1KFu2LF27dqVw4cJcuHCBrVu3cubMGfbs2fPQz/Dx8WHy5Mm8+eabVKlShfbt25M9e3b27NnD7du3+eabb3BycuKrr76iSZMmlC5dmk6dOpE3b14iIyNZu3YtPj4+LFq0KMl/Jvb+PokkWlotlxJJaw8uvzYMw4iOjjZ69+5tBAYGGq6urkaxYsWM8ePHWy03NQzDOHjwoPHcc88Znp6eBhC/fPTatWtGp06djFy5chne3t5Go0aNjIMHDxoFCxa0WmL6pMuvDcMwvv/+e6Nw4cKGm5ubUb58eWPFihUJll/fN2vWLKNSpUqGp6enkTVrVqNs2bLGgAEDjLNnz8bfk9jl1/babhiGsXv3bqNRo0aGt7e3kSVLFqNu3brGli1bHhnTMAxjy5YtRqVKlQw3Nzerpdj2ll93797d+P77741ixYoZ7u7uRoUKFR7553nfggULjFKlShkuLi5WS7H3799v1K9f3/D29jZy5cpldO3a1dizZ0+C5dodOnQwvLy8EsR9sK7/XX79X59++qkBGP369bNbx/++d+LEiUb+/PkNd3d349lnnzX27NmT4P6jR48ab7zxhpEnTx7D1dXVyJs3r9G8eXPj119/jb/n/vLrnTt32vzMhQsXGjVq1DA8PT0NHx8fo2rVqsaPP/5odU9YWJjxwgsvGDlz5jTc3d2NggULGm3btjVWr16d4M/h0qVLVu+9//nHjx+PL3vY3yeRxDAZhmZdiYhjMZlMdO/ePcEwYEZy4sQJgoKCGD9+PP369Uvr6oikW5ojIyIiIg5LiYyIiIg4LCUyIiIi4rA0R0ZEREQclnpkRERExGEpkRERERGHleE3xLNYLJw9e5asWbNqK2wREREHYRgG0dHRBAYG4uRkv98lwycyZ8+efegZLyIiIpJ+nT59mnz58tm9nuETmfuH/50+fVrHxouIiDiIqKgo8ufPb3WIry0ZPpG5P5zk4+OjREZERMTBPGpaiCb7ioiIiMNSIiMiIiIOS4mMiIiIOCwlMiIiIuKwlMiIiIiIw1IiIyIiIg5LiYyIiIg4LCUyIiIi4rCUyIiIiIjDyvA7+4qIiEjyuxdn4butJzh59TYFc2Th9eqFcHNJ/f4RJTIiIiKSJCMWRzB700mrspFLD9D12SAGNy2VqnVRIiMiIiKJ1mzqOv4+dytBucWALzYcB0jVZEZzZEREROSRzBaDyiP+tEpinCxm6h7daXXflxuPcy/Okmr1UiIjIiIiD7U4PJIi7y/l8q3Y+DK/m9f49pehhP46jBb718eXWwz4buuJVKubhpZERETEJrPFoM1nm9h9OsqqvMaJcKYunoDfrevcdnXHhGF1/eTV26lWRyUyIiIiksDi8Eh6/hRulaI4Wcy8t/lHem75GScMDuYqSI9WAzmSq4DVewvmyJJq9VQiIyIiIla6zNnO6oOXrcr8o68wbdF4njkdAcCP5RoyrP5b3HX1SPD+16sXSo1qAkpkRERE5P+ZLQbBE9Zw4updq/Lax3YxafFEct6J4qabJ+836s7CUnVsxuhSq2Cq7iejREZERERsDiU5W8z03fgd3bb9CsB+/yC6txrE8Rx5bcYoHeDFkOZlUqG2/6NERkREJJOzNZQUEHWJaQvHUyVyPwDfVWjKJ/XeJMbFzWaMfNncWfJenZSuagJpuvw6OjqaXr16UbBgQTw9PalRowY7d/5vPbphGAwdOpSAgAA8PT2pX78+hw8fTsMai4iIZBxmi0GdcasTJDH1juxgaWgIVSL3E+WWhe4tBzKkYTe7SUz+7B5sGlQ/NaqcQJomMm+++SYrV67ku+++Y9++fTRs2JD69esTGRkJwLhx45g2bRqff/4527dvx8vLi0aNGnH37t1HRBYREZGHWbr3LMXeX2o1H8bVHMsHa77i69+Gk/1uNHvzFKV5x6ksKfms3Tj1SvixcWBwalTZJpNhGMajb0t+d+7cIWvWrCxYsIBmzZrFl1eqVIkmTZowYsQIAgMD6du3L/369QPgxo0b5M6dmzlz5tCuXbtEfU5UVBS+vr7cuHEDHx+fFGmLiIiII7F1VlK+GxeYvmAcFc4dAiC0UgtG1+nMPRdXmzFMwLRXKtDi6cAUqWNif3+n2RyZuLg4zGYzHh7Wy7Y8PT3ZtGkTx48f5/z589Sv/7+uKl9fX6pVq8bWrVvtJjIxMTHExMTEv46KirJ5n4iISGZjthi89Nkmwh7Y4K7RP1sYt3QqvjG3uOHuRf+mvfjzqep241TMn41579bA2cmU0lV+pDRLZLJmzUr16tUZMWIEJUuWJHfu3Pz4449s3bqVokWLcv78eQBy585t9b7cuXPHX7Nl9OjRDBs2LEXrLiIi4mhsrUpyi4tl8Lqv6bRrEQBhAcXp2WoAZ3xz2w4CTE/BXpjHkaZzZL777jsMwyBv3ry4u7szbdo0XnnlFZycHr9agwcP5saNG/E/p0+fTsYai4iIOBazxeCFmRvp8UASU+DaOX79oX98EjOryvO0fXWM3SSmUA5Pjo5qmq6SGEjj5ddFihRh/fr13Lp1i6ioKAICAnj55ZcpXLgwefLkAeDChQsEBATEv+fChQuUL1/ebkx3d3fc3d1TuuoiIiLp3tK9Z+kxN4wHz6JudmAjY5ZPI+u9O1zzyErfZr1ZU7Sq3Tj1ivvxdSf719NSujj92svLi4CAAK5du8aKFSto1aoVQUFB5MmTh9WrV8ffFxUVxfbt26le3f64nYiIiMDIJX/T7YEkxj3uHp+smMnMhWPJeu8OO/OWommnaQ9NYrrUCkq3SQykcY/MihUrMAyD4sWLc+TIEfr370+JEiXo1KkTJpOJXr168cknn1CsWDGCgoIYMmQIgYGBtG7dOi2rLSIikq4NWxRB6GbrVUlBVyOZuWAMpS4eB2DmM22Y9OxrmJ2cbcZwAma0r0jTcgE2r6cXaZrI3Lhxg8GDB3PmzBly5MjBiy++yMiRI3F1/Xep14ABA7h16xZvvfUW169fp1atWixfvjzBSicRERGxvyqp1d9rGbViJl6xd7mcxZc+zfqwoXAlu3Eq5PPl124108WqpEdJs31kUov2kRERkczA1qokj9i7fLxqFu32/gnA1gJlea95Py5mzWk3TpdaQQxpXiqFa/to6X4fGREREUkets5KKnr5FDMXjKH45VNYMDG9Rjum1myHxcGHkh6kREZERMRBmS0GwRPWWB0zAPDSvlUMX/kZWWJjuOSVjfea92NLofJ246SnDe6SSomMiIiIA7I1lJTl3h1GrPyMFyPWALCxYHl6t+jLZa/sNmOk9DEDqUGJjIiIiAMxWwzafLaJ3Q9M6C1+6QQz54+h6NUzmE1OTK7Vnk+faWN3KKl8Xh9+617LIXth/kuJjIiIiIOwucGdYdBuzwo+Xj0Lj7h7nPfOwXst+rO9QFm7cdLzBndJpURGRETEAYxc8jdfbjxhVeYVc5tRK2bS6sB6ANYFVaJP8z5czeJrN079kv581aFKSlY1VSmRERERSedsbXBX+sJRZiwYQ9C1c8SZnJjw3Bt8Ue0FDJPtTfszwnwYW5TIiIiIpFM2N7gzDF4LW8qQNV/ibo4jMqsfPVsOYHe+knbjOPKqpEdRIiMiIpIO2VqVlDXmFmOWTaPZoc0ArCxalf5Ne3Hd0/6GcdMzYC/MfymRERERSUfsrUoqd+4fZiwYS4EbF4h1cmZMnU7MrtwKTLZ7WQrl8GR1v7oZshfmv5TIiIiIpBP2ViV12rWQwWtDcbPEcdo3Nz1aDmBPYHG7cYJL+DO7Y8aZ0PswSmRERETSgRGLI5i9yXpCr++daMYvm0rDw9sAWPZUDQY2CSHKw9tunIw+lPQgJTIiIiJprNPX21j7zxWrsgqRB5m+cCz5oi4R4+zCyLpd+LZi80w/lPQgJTIiIiJpxGwxeGbkSi7dio0vMxkWuu74g/4bvsXVYuZEtgB6tBpIRJ6iduNkpA3ukkqJjIiISBpYHB5Jj5/Crcqy377BhKVTCD66E4BFJZ5lcOOe3HTPYjdOl1pBDGleKiWrmq4pkREREUllXeZsZ/XBy1ZlVU5HMG3heAJuXiHG2ZVh9d9i7tON7Q4lOQEz2lekabmAVKhx+qVERkREJJWYLQbBE9Zw4urd+DKTYeHdbb/SZ+P3uBgWjubIS49WAzngX9hunIy8wV1SKZERERFJBbY2uMt56zqTF0/kuRNhAPxeui4fNuzGbTdPmzEy6jEDT0KJjIiISAqyt8Fd9ZN7mbpoPP63rnHHxZ2hDd5mXtkGdoeSKuTz5dduNdUL8wAlMiIiIinEVi+Mk8VMzy0/E7LlJ5wNC//kLED3VgM57FfQbpzMPqH3YZTIiIiIpABbE3r9bl5lyuIJ1Dy5F4Bfytbno/rvcMfNw26cTzWh96GUyIiIiCSz5lPXE3HuplVZreNhTF48Eb/b17nl6sGHDbvxR5l6dmP4ebmy7YMGGkp6BCUyIiIiycTWqiRni5lem+bSfesvOGFwwK8QPVoN5GjO/Hbj1Cvhx9cdM+cGd0mlREZERCQZ2JoPkyfqMlMXjafamb8B+KF8Y4bX60qMq7vdOJntrKQnpURGRETkCdhblVTn6E4mLZlMjjtRRLt58n6jHiwqVdtunMx6VtKTUiIjIiLymGz1wriY4+i34Vve2fE7ABG5i9Cj5QBO5MhrN05wCX9md6ySwrXNmJTIiIiIPAZbq5ICoy4yfcE4Kp09CMCcis0ZXbczMS5uNmNog7snp0RGREQkCWxN6AWof3g7E5ZOJtvdm0S5ezGgSQjLi9e0G0fHDCQPJTIiIiKJZGsoydUcy6B1c+jy1wIAwgOK0bPlQE5ny2M3ztSXy9Oqgv2hJkk8JTIiIiKPYG9Cb77r55mxcCzlzx0G4KvKrRhbpyOxzq52Y2mDu+SlREZEROQhlu49S4+5YVgeKG90aAvjl03FJ+YW1z286de0N6uKVbMbx9kEM1+tSOMySmKSkxIZERERO0Yu+ZsvN56wKnOPu8fgtV/TcfdiAHYFliCk5QAiff3txtF8mJSjREZERMSGYYsiCN180qqs4LWzzFwwljIXjgLwebUXmfDs68Q52/51qlVJKU+JjIiIyH+YLQYvfbaJsAfmw7TYv55RK2aQ9d4drnr60KdZb9YVsb/3i3phUocSGRERkf9na1WSe2wMH63+kvZ7lgOwPV9p3mvRn/M+uezG0TEDqUeJjIiIZHr2ViUVuXKaGQvGUvLSCSyYmFG9LVNrtcfs5Gwzjo4ZSH1KZEREJFOztyrp+Yg1fPLnp3jF3uVSlmz0bt6XTUEV7MbRMQNpQ4mMiIhkWiMWRzB7k/WEXs97dxm26nPa7lsFwJYC5XivRT8ueeewG0dDSWlHiYyIiGRKnUO3s+aQ9VlJxS6dZOaCsTx15RRmkxNTa77CjOptsWgoKd1SIiMiIpmKzbOSDIM2+1YyfOUXeMbFcNErOyEt+7OtQDm7ceoV9+PrTlVTocbyMEpkREQk07C1KinLvTt88uenvPD3WgA2FKpA7+Z9ueKVzW6cLrWCGNK8VMpWVhJFiYyIiGQKXeZsZ/VB66GkkhePMWPBWIpcjSTO5MSkZ1/js2dewjA52YzhBMzQWUnpiu1vKpWYzWaGDBlCUFAQnp6eFClShBEjRmAY/8uVO3bsiMlksvpp3LhxGtZaREQcyb04C8+MWmmdxBgG7cOXMf/bvhS5Gsk575y0az+aT6u3tZvEVMyfjcOjmiqJSWfStEdm7NixfPbZZ3zzzTeULl2av/76i06dOuHr60tISEj8fY0bNyY0NDT+tbu7e1pUV0REHIytVUneMbcZvXw6LQ5uBGBN4cr0bdaba1l8bcbQMQPpW5omMlu2bKFVq1Y0a9YMgEKFCvHjjz+yY8cOq/vc3d3JkydPWlRRREQckNliEDxxDSeu3LUqL33+CDMXjKXQ9XPEOjkz7rkOfFW1td1emAr5fPm1W02tSkrH0nRoqUaNGqxevZp//vkHgD179rBp0yaaNGlidd+6devw9/enePHivPvuu1y5csVuzJiYGKKioqx+REQk81gcHkmx95daJzGGwRu7FvH79/0odP0cZ3z8aNt+LF9We8FuEtOpZiH+6FFLSUw6l6Y9MoMGDSIqKooSJUrg7OyM2Wxm5MiRvPrqq/H3NG7cmBdeeIGgoCCOHj3K+++/T5MmTdi6dSvOzgnX9Y8ePZphw4alZjNERCQdsHfMgM/dm4xdNo0m/2wB4M9iz9C/yXvc8MxqN1bXZ4P4oJlWJTkCk/HfmbWp7KeffqJ///6MHz+e0qVLEx4eTq9evZg0aRIdOnSw+Z5jx45RpEgRVq1aRXBwcILrMTExxMTExL+Oiooif/783LhxAx8fnxRri4iIpJ3lEefo/v1uzA+UP332EDMWjiP/jQvcc3JhdN1OhFZqCSbbvSxalZR+REVF4evr+8jf32naI9O/f38GDRpEu3btAChbtiwnT55k9OjRdhOZwoULkytXLo4cOWIzkXF3d9dkYBGRTGTp3rN0mxtmXWgYdNk5n0Hr5+BqMXPKNzc9Wg1kb8BTduNUzJ+Nee/W0FCSg0nTROb27ds4OVmPTTo7O2OxPHh01/+cOXOGK1euEBCgbFlEJLOb/9cZev26x6rM9040E5ZOpsGRfxeOLClek0FNQoh297IZQ6uSHFuaJjItWrRg5MiRFChQgNKlSxMWFsakSZPo3LkzADdv3mTYsGG8+OKL5MmTh6NHjzJgwACKFi1Ko0aN0rLqIiKShuytSqp45gDTF44jb/QlYpxdGRHcle/LN7E7lKRVSY4vTefIREdHM2TIEP744w8uXrxIYGAgr7zyCkOHDsXNzY07d+7QunVrwsLCuH79OoGBgTRs2JARI0aQO3fuRH1GYsfYRETEMSzde5buc8OsjhkwGRbe3v47/TZ8i4th4Vj2QHq0GsT+3IXtxtExA+lbYn9/p2kikxqUyIiIZBy2NrjLcfsGE5dMou6xXQAsKFmb9xt155Z7FpsxnE0w/RVN6E3vHGKyr4iISGKYLQYvfbaJsAeWVlc9HcG0hePIc/Mqd13c+Kj+2/xcrqHdoaSmZfIwvX1FDSVlIEpkREQkXbN1YrWTxUy3bfPovWkuzoaFIzny0b31IA75FbIZI7e3GxsHBePmkqb7wEoKUCIjIiLpkr0N7nLdusbkRRN59mQ4AL+VqceQBu9y283TZpzSAd4sea92SldX0ogSGRERSXeW7j1Lj7lhPLgZR40T4UxdPAG/W9e57erO0Abv8mvZ+nbjlAnMyuKQ51K2spKmlMiIiEi6MnLJ33y58YRVmZPFzHubf6Lnlp9wwuBQrgJ0bzWII7kK2I0TXMKf2R2rpHBtJa0pkRERkXRj2KIIQjdbr0ryj77C1MUTqH5qHwA/lmvIsPpvcdfVw26c6drgLtNQIiMiImnO3qqk547tYtKSSeS6fYObbp6836g7C0vVsRunfD4ffuumE6szEyUyIiKSpmytSnK2mOmz8Xu6b5sHwH7/ILq3GsTxHHntxtEGd5mTEhkREUkT9lYlBURdYtrC8VSJ3A/AdxWa8km9N4lxcbMZRydWZ25KZEREJNXZW5VU9+hOJi2eRPa70US7eTKocQhLSj5rN45OrBYlMiIikqpsHTPgYo6j/4ZveXvH7wDszVOUHi0Hciq77V4WnVgt9ymRERGRVGFvQm++GxeYvmAcFc4dAiC0UgtG1+nMPRdXm3HK5/Xht+6a0Cv/UiIjIiIpztaEXoCG/2xl/NIp+Mbc4oa7FwOavseKp2rYjVOvuB9fd6qaspUVh6JERkREUlSXOdtZffCyVZlbXCyD1oXSeddCAMIDnqJHq4Gc8c1tN079kv581UEb3Ik1JTIiIpJimk9dT8S5m1ZlBa6dY8bCsZQ7fwSAWVWeZ3ztN4h1tj2UpPkw8jBKZEREJNmZLQbBE9Zw4updq/KmBzcxZtk0fO7d5ppHVvo2682aovaHirQqSR5FiYyIiCQrW/Nh3OPu8eGar3g9bCkAO/OWIqRlf875+NmMoV4YSSwlMiIikizsbXAXdDWSmQvGUOricQBmPtOGybVeJc7Z9q+gCvl8+bVbTfXCSKIokRERkSdmb1VSy/3rGLViJt737nDF04fezfuyoXAlu3F0zIAklRIZERF5IrZWJXnE3uWjVbN4Ze+fAGzLX4aQFv25mDWnzRg6ZkAelxIZERF5LPYm9Ba5fJqZC8ZQ4vJJLJiYXqMd02q2w+zkbDOOhpLkSSiRERGRJLM3lPTivtWMWPkpWWJjuOSVjfea92NLofJ242goSZ6UEhkREUmSzqHbWHPoilWZ5727jFj5GS9FrAZgU8Gn6d28H5e8s9uMYQJmaihJkoESGRERSRSzxeCZkSu5dCvWqrz4pRPMnD+GolfPYDY5MblWez59pg0WO0NJhXJ4srpfXQ0lSbJQIiMiIo+0dO9Zus8Nsx5KMgxe3vsnw1Z9gUfcPc575+C9Fv3ZXqCs3TjBJfyZ3VHHDEjyUSIjIiIPNWJxBLM3nbQq84q5zcg/Z9J6/3oA1gVVok/zPlzN4mszhja4k5SiREZERGwyWwxe+mwTYQ9scFfqwjFmLBhD4WtniTM5MeG5N/ii2gsYJiebcXTMgKQkJTIiIpKAzVVJhsFr4csYsvpL3M2xnM2ai54tB7Arn/1VR9PVCyMpTImMiIjEs3fMQNaYW4xeNp3mhzYBsKpIFfo16811Tx+bcTShV1KLEhkREQH+ndDbY24YlgfKy547zIyFYyl4/TyxTs6Mqd2R2VVag8l2kqIJvZKalMiIiIjNCb0YBp12LWTw2lDcLHGc8fGnR6uBhAcWtxtHQ0mS2pTIiIhkYvYm9Prcvcn4pVNodHgbAMufqs6AJu8R5eFtM46flyvbPmigoSRJdUpkREQyKXtDSRUiDzJ94TjyRV0kxtmFkXW78G3F5naHkuoWz0Vop2opX2ERG5TIiIhkQraGkkyGhTd3zGfAhm9wtZg5kS2AHq0GEpGnqN04OitJ0poSGRGRTKZz6HbWHLpsVZbtThQTl0wm+OhOABaXeJbBjXsQ7e5lM4YTMENnJUk6oERGRCSTMFsMgies4cTVu1bllc/8zbSF4wmMvkyMsyvD6r/F3Kcb2x1K0gZ3kp4okRERyQRsbXBnMiy8u+1X+mz8HhfDwtEceenRaiAH/AvbjKFjBiQ9UiIjIpKB2dvgLuet60xePJHnToQB8EepOnzYsBu33LPYjFMhny+/dqupXhhJd5TIiIhkUDaPGQCeObWXqYsmkPvmVe64uDO0wdvMK9vA7lCSJvRKeqZERkQkA+oyZzurD1pP6HWymOmx9Rfe2/wjzoaFwznz063VIA77FbQZwwTM1IReSeeUyIiIZCD2JvT63bzGlMXjqXlyLwC/lK3PR/Xf4Y6bh804OitJHIXtM9cfYsOGDcTFxSUoj4uLY8OGDUmKZTabGTJkCEFBQXh6elKkSBFGjBiBYfyvI9QwDIYOHUpAQACenp7Ur1+fw4cPJ7XaIiIZ3uLwSIq+vzRBElPzRDhLQ3tS8+Rebrl60LtZHwY07WU3iQku4c+6AfWUxIhDSHKPTN26dTl37hz+/v5W5Tdu3KBu3bqYzeZExxo7diyfffYZ33zzDaVLl+avv/6iU6dO+Pr6EhISAsC4ceOYNm0a33zzDUFBQQwZMoRGjRqxf/9+PDxsP4QiIpmNraEkZ4uZ9zbNpcfWX3DC4IBfIXq0GsjRnPntxtFZSeJokpzIGIaBycaEsCtXruDlZXvjJHu2bNlCq1ataNasGQCFChXixx9/ZMeOHfGfNWXKFD788ENatWoFwLfffkvu3LmZP38+7dq1S2r1RUQyFHtDSbmjLzNt0QSqnY4AYO7TjRkW3JUYV3ebcTSUJI4q0YnMCy+8AIDJZKJjx464u//vYTCbzezdu5caNWok6cNr1KjBrFmz+Oeff3jqqafYs2cPmzZtYtKkSQAcP36c8+fPU79+/fj3+Pr6Uq1aNbZu3apERkQyNXtnJdU+totJiyeS804UN908GdyoB4tK1bYbJ7iEP7M7VknZyoqkkEQnMr6+vsC/vSRZs2bF09Mz/pqbmxvPPPMMXbt2TdKHDxo0iKioKEqUKIGzszNms5mRI0fy6quvAnD+/HkAcufObfW+3Llzx197UExMDDExMfGvo6KibN4nIuLIbJ2V5GKOo+/G73l3+68AROQuQo+WAziRI6/NGNrgTjKCRCcyoaGhwL/DP/369UvyMJItv/zyCz/88ANz586ldOnShIeH06tXLwIDA+nQocNjxRw9ejTDhg174rqJiKRHZovBS59tIuyBDe4Coy4ybeF4KkceAOCbis0YVbcLMS5uNuPomAHJKEzGf5cIpbL8+fMzaNAgunfvHl/2ySef8P3333Pw4EGOHTtGkSJFCAsLo3z58vH31K5dm/LlyzN16tQEMW31yOTPn58bN27g4+OTou0REUlJ9ja4Cz6ynYlLJpPt7k2i3L0Y0CSE5cVr2o2jCb3iCKKiovD19X3k7+8kT/a9cOEC/fr1Y/Xq1Vy8eJEH86CkrFq6ffs2Tk7WK8CdnZ2xWP4d8Q0KCiJPnjysXr06PpGJiopi+/btvPvuuzZjuru7W83fERFxdPaOGXA1xzJg/Td03TkfgPCAYvRsOZDT2fLYjKMJvZIRJTmR6dixI6dOnWLIkCEEBATYXMGUWC1atGDkyJEUKFCA0qVLExYWxqRJk+jcuTPw78TiXr168cknn1CsWLH45deBgYG0bt36sT9XRMRR2JvQm+/6eWYsHEf5c/8AMLtyK8bU6Uiss6vNOPWK+/F1p6opXFuR1JfkRGbTpk1s3LjRaqjncU2fPp0hQ4bQrVs3Ll68SGBgIG+//TZDhw6Nv2fAgAHcunWLt956i+vXr1OrVi2WL1+uPWREJMOzNaEXoNGhLYxfNhWfmFtc9/CmX9PerCpWzW4cnZUkGVmS58iUKlWKH374gQoVKqRUnZJVYsfYRETSk86h21lzyHqDO7e4WN5fO5uOuxcDsDuwOD1bDiTS199WCJyAGTorSRxUis2RmTJlCoMGDeKLL76gUKFCT1JHERF5gL0N7gpeO8uMBWMpe+EoAJ9Xe5EJz75OnLPtf8Yr5PPl1241NR9GMrwkJzIvv/wyt2/fpkiRImTJkgVXV+vx2KtXryZb5UREMhN7q5KaH9jA6OXTyXrvDlc9fejTrDfritjfwE5DSZKZPFaPjIiIJB97q5LcY2MYuuZLXg1fDsCOfKUIaTGA8z65bMbRUJJkRklOZB53ozoREUnI3qqkwlfOMHPBGEpeOoEFEzOrt2VKrfaYnZxtxtEGd5JZJTmRATh69CihoaEcPXqUqVOn4u/vz7Jly+KXUYuIyKPZW5XU+u+1jFwxE6/Yu1zKko3ezfuyKcj2AgsdMyCZndOjb7G2fv16ypYty/bt2/n999+5efMmAHv27OGjjz5K9gqKiGQ0ZovB8zM3JkhiPGLvMnbpVKYsnohX7F22FChH007T7CYxFfL5cmRUUyUxkqklOZEZNGgQn3zyCStXrsTN7X9neNSrV49t27Yla+VERDKaxeGRFH1/aYKzkopdOsnCb/rw8r6VWDAxuWZ7Xnt5BJe8c9iM06lmIf7oUUtDSZLpJXload++fcydOzdBub+/P5cvX7bxDhERsTehF8Ogzb5VDF/5OZ5xMVz0ys57LfqztWA5u7G6PhvEB820KkkEHiORyZYtG+fOnSMoKMiqPCwsjLx5bR8VLyKSmS2POEf373fz4El0We7dYcSfn/Li32sB2FCoAr2b9+WKVzabcbQqSSShJCcy7dq1Y+DAgcybNw+TyYTFYmHz5s3069ePN954IyXqKCLisJbuPUu3uWEJyktcPM7MBWMpcvUMZpMTE599jc+eeQnDZHvEX6uSRGxLciIzatQounfvTv78+TGbzZQqVQqz2Uz79u358MMPU6KOIiIOaeHuSEJ+CbcuNAxe2bOCj1bPwiPuHue8cxLSsj8785exGUOrkkQeLslnLd13+vRp9u3bx82bN6lQoQLFihVL7rolC521JCKpzd58GO+Y24xaMYOWBzYAsKZwZfo26821LL424+iYAcnMUuyspfvy589P/vz5H/ftIiIZkr1jBkpfOMqMBWMIunaOWCdnxj/3Bl9Wfd7uUJKOGRBJnCQnMi+++CJVq1Zl4MCBVuXjxo1j586dzJs3L9kqJyLiSLrM2c7qgw+s3jQMXg9bwodrvsLdHMcZHz9CWg5gd96SNmNoQq9I0iQ5kdmwYQMff/xxgvImTZowceLE5KiTiIhDsXditc/dm4xZNo2m/2wBYGXRavRr2osbnlltxtGEXpGkS3Iic/PmTauN8O5zdXUlKirKxjtERDIue0NJ5c79w4wFYylw4wL3nFwYXbcToZVagsl2kjJdE3pFHkuSE5myZcvy888/M3ToUKvyn376iVKlNJ4rIpnDwza46/zXQgatC8XNEscp39z0aDWQvQFP2YxTMLsHa/rXUy+MyGNKciIzZMgQXnjhBY4ePUq9evUAWL16NT/++KPmx4hIpmDvxGrfO9FMWDqFBke2/3vfUzUY1CSEKA9vm3HKBGZlcchzKVxbkYwtyYlMixYtmD9/PqNGjeLXX3/F09OTcuXKsWrVKmrXrp0SdRQRSTdGLvmbLzeeSFBeMfIA0xeMI2/0JWKcXRhRryvfV2hqdygpuIQ/sztWSeHaimR8SUpk4uLiGDVqFJ07d2bz5s0pVScRkXRp2KIIQjdbn1htMiy8teN3+q//FhfDwvHsAfRoNYi/cxexGUMb3IkkryRviOft7U1ERASFChVKoSolL22IJyJPymwxeOmzTQlOrM5x+wYTl0yi7rFdACws+RzvN+rBTfcsNuNoVZJI4qXYhnjBwcGsX7/eYRIZEZEnYW9VUpXTEUxfOI48N69y18WNj4Pf4qenG2lVkkgqS3Ii06RJEwYNGsS+ffuoVKkSXl5eVtdbtmyZbJUTEUkr9lYlmQwL3bbOo8+mH3A2LBzNkY/urQZy0D/IZpxCOTxZ3a+uemFEUkiSh5acnGxvpw1gMpkwmx88qD5taWhJRJLK3qqkXLeuMWnxJJ478e9p1r+VrsuQht247eZpM44m9Io8vhQbWrJYHny0RUQyjhGLI5i96WSC8uon9zB10QT8b13jtqs7Qxu8y69l69uNo6EkkdTx2IdGAty9excPD4/kqouISJqxN6HXyWImZMtPhGz+CScMDuUqQPdWgziSq4DNOBpKEkld9seJ7DCbzYwYMYK8efPi7e3NsWPHgH83yps9e3ayV1BEJKUtDo+k6PtLEyQx/tFX+OHnD+m1+UecMPipXENavTHJbhJTr7gf6wZol16R1JTkRGbkyJHMmTOHcePGWZ25VKZMGb766qtkrZyISErrMmc7PWysSnr2+G6Wzgmh+ql93HL14L3mfRnUJIS7rrZ7obvUCuLrTlVTvsIiYiXJQ0vffvsts2bNIjg4mHfeeSe+/Omnn+bgwYPJWjkRkZTUfOp6Is7dtCpztpjpvekHum2dhxMG+/2D6N5qEMdz5LUZwwmY0b4iTcsFpEKNReRBSU5kIiMjKVq0aIJyi8VCbGxsslRKRCQlmS0GwRPWcOLqXavyPFGXmbZoHFXP7Afg+/JNGFHvTWJc3W3G0QZ3ImkvyYlMqVKl2LhxIwULFrQq//XXX6lQoUKyVUxEJCXY2+CuztGdTFoymRx3ooh282RQ4xCWlHzWZgwdMyCSfiQ5kRk6dCgdOnQgMjISi8XC77//zqFDh/j2229ZvHhxStRRROSJ2dvgzsUcR78N3/LOjt8B2Je7CN1bDeJUdttDRRXy+fJrt5rqhRFJJ5K8IR7Axo0bGT58OHv27OHmzZtUrFiRoUOH0rBhw5So4xPRhngiYq8XJu+Ni0xfOJaKZw8BEFqpBaPrdOaei6vNOF1qBTGkeakUrq2IQOJ/fz9WIuNIlMiIZG5d5mxn9cHLCcobHN7GhCWT8Y25xQ13LwY0fY8VT9WwGUMTekVSX4rt7HvfX3/9xYEDB4B/581UqlTpcUOJiCQ7exN6Xc2xDF4bSuddCwEID3iKHq0GcsY3t804GkoSSd+SnMicOXOGV155hc2bN5MtWzYArl+/To0aNfjpp5/Ily9fctdRRCRJ7A0l5b9+nhkLxvL0+cMAzKryPONrv0Gss4aSRBxVkjfEe/PNN4mNjeXAgQNcvXqVq1evcuDAASwWC2+++WZK1FFEJNHsbXDX5OAmloSG8PT5w1zzyEqXF4cwql4Xm0mMCfi0fUUlMSIOIMlzZDw9PdmyZUuCpda7du3i2Wef5fbt28lawSelOTIimYO9oST3uHt8sGY2b4QtAeCvvCXp2XIA53z8bMbRWUki6UOKzZHJnz+/zY3vzGYzgYHaU0FEUt/SvWfpMTcMywPlha5GMnPBWEpf/PdMuE+feYlJtV4jztn2P33BJfyZ3bFKCtdWRJJTkhOZ8ePH07NnT2bOnEnlypWBfyf+vvfee0yYMCHZKygi8jAjFkcwe9PJBOUt969n1IoZeN+7wxVPH/o078v6wrYXJWiDOxHHleShpezZs3P79m3i4uJwcfk3D7r//728vKzuvXr1avLV9DFpaEkkYzJbDF76bFOCE6vdY2P4aPUs2u9ZAcD2/GUIadGPC1lz2YyjYwZE0qcUG1qaMmXKk9RLROSJ2VuVVOTKaWbOH0OJyyexYGJ6jZeZVvMVzE7ONuNMVy+MiMNLciLToUOHZPvwQoUKcfJkwi7hbt26MXPmTOrUqcP69eutrr399tt8/vnnyVYHEXEc9o4ZAHghYjWf/PkpWWJjuOSVjV7N+7G5UHmbcTShVyTjeOwN8ZLDzp07MZvN8a8jIiJo0KABbdq0iS/r2rUrw4cPj3+dJUuWVK2jiKQP9nphPO/dZfjKz2kTsQqATQWfpnfzflzyzm4zTr3ifnzdqWoK11ZEUkuaJjJ+ftbLH8eMGUORIkWoXbt2fFmWLFnIkydPaldNRNIRe8cMPHXpBDMXjKXYldOYTU5MqfkKM6u3xWJnKEkb3IlkPEneEC+l3Lt3j++//57OnTtjMv2vu/eHH34gV65clClThsGDBz9yn5qYmBiioqKsfkTEcTWbsi5hEmMYtN3zJwu+7UuxK6c5752D9u1GMr3mKzaTGCe0wZ1IRpWmPTL/NX/+fK5fv07Hjh3jy9q3b0/BggUJDAxk7969DBw4kEOHDvH777/bjTN69GiGDRuWCjUWkZRkthhUG7mSy7es963yirnNJ39+yvP71wGwPqgivZv35WoWX5txdFaSSMaW5OXXa9eupW7dujavzZw5k+7duz9WRRo1aoSbmxuLFi2ye8+aNWsIDg7myJEjFClSxOY9MTExxMTExL+Oiooif/78Wn4t4kAWh0fS46fwBOUlLx5jxoKxFLkaSZzJiYnPvc7n1V7EMNnuXNZQkojjSuzy6yQPLb3wwgvs2rUrQfnUqVMZPHhwUsMBcPLkSVatWvXIs5qqVasGwJEjR+ze4+7ujo+Pj9WPiDgGs8XghZkbEyYxhsGrYUuZ/21filyN5GzWXLzcfgyfPdPGZhKjoSSRzOOxdvZt0qQJGzZsoESJEgBMnDiR4cOHs2TJkseqRGhoKP7+/jRr1uyh94WHhwMQEBDwWJ8jIumXvWMGssbcYvTyGTQ/uBGAVUWq0K9Zb6572v6PFG1wJ5K5JDmRefPNN7l69Sr169dn06ZN/Pzzz4waNYqlS5dSs2bNJFfAYrEQGhpKhw4d4ncKBjh69Chz586ladOm5MyZk71799K7d2+ee+45ypUrl+TPEZH0y94xA2XOH2HmgjEUvH6eWCdnxtTuyOwqrcGUMEnRMQMimdNjTfYdMGAAV65coXLlypjNZlasWMEzzzzzWBVYtWoVp06donPnzlblbm5urFq1iilTpnDr1i3y58/Piy++yIcffvhYnyMi6Y+9YwYwDDrsXsz7a2fjbo7jjI8/PVoNJDywuM046oURybwSNdl32rRpNssnTJjAc889R9Wq/9tcKiQkJPlqlwx01pJI+mRvgzufuzcZt2wqjf/ZCsCKYs/Qv2kvojy8bcbRMQMiGVNif38nKpEJCgpK1IeaTCaOHTuW+FqmAiUyIunLw44ZKH/2EDMWjCVf1EVinF0YVbcL31RsbnMoKVcWF7Z/2FC9MCIZVLIeGnn8+PFkq5iIZF72JvRiGLy58w8Grv8GV4uZk9ny0L3VICLyFLUZJ392DzYODE7x+opI+pduNsQTkYzN3oTebHeimLBkMvWP7gRgcYlnGdy4B9HuXjbj1Cvhx9cddVaSiPwryYmM2Wxmzpw5rF69mosXL2KxWP+31Zo1a5KtciLi+P4dStrM7tM3ElyrdGY/0xeOIzD6MjHOrgwP7soP5ZtoVZKIJFqSE5n33nuPOXPm0KxZM8qUKWN1LpKIyH8tDo8k5OdwLA/MxDMZFt7Z/ht9N3yHi2HhaI689Gg1kAP+hW3G0TEDImJPkhOZn376iV9++YWmTZumRH1EJIOwd2J1zlvXmbRkErWP7wbgj1J1+LBhN265Z7EdR8cMiMhDJDmRcXNzo2hR2xPwRETMFoPgCWs4cfVugmvVTu1j2qLx5L55lTsu7nxU/21+KdfA5lCSEzCjfUWaltNO3iJiX5ITmb59+zJ16lRmzJihYSURsWJvVZKTxUz3rb/Qa/OPOBsWDufMT/dWA/nHr5DNONrgTkQSK8mJzKZNm1i7di3Lli2jdOnSuLq6Wl3//fffk61yIuI47K1K8rt5jcmLJ1Dr5B4A5pWpz9AG73DHzSPBvZrQKyJJleREJlu2bDz//PMpURcRcUB2jxkAapwIZ+riCfjdus5tV3c+bNiN38vY3v9FE3pF5HEkOZEJDQ1NiXqIiAOyd8yAs8VMyOYf6bnlZ5wwOOBXiB6tBnI0Z36bcTShV0QelzbEE5HHYm9VUu7oy0xbNIFqpyMAmPt0I4YFv0WMq3uCe11MMO0VTegVkcf3WInMr7/+yi+//MKpU6e4d++e1bXdu3cnS8VEJH162Kqk2sd2MWnxRHLeieKmmyfvN+rBwlK1bcZpUjo3M16tpKEkEXkiTkl9w7Rp0+jUqRO5c+cmLCyMqlWrkjNnTo4dO0aTJk1Soo4ikk4sDo+k6PtLEyQxLuY4BqyfwzfzPiLnnSj+9i9M8w5T7CYxXZ8N4rPXKyuJEZEnluQemU8//ZRZs2bxyiuvMGfOHAYMGEDhwoUZOnQoV69eTYk6ikg60Dl0G2sOXUlQHhB1iekLx1E58gAA31Rsxqi6XYhxcUtwr/aGEZHkluRE5tSpU9SoUQMAT09PoqOjAXj99dd55plnmDFjRvLWUETSlNli8MzIlVy6FZvgWr0jO5i4ZDLZ70YT5ZaFgU1CWFails045fP68Fv3WuqFEZFkleShpTx58sT3vBQoUIBt27YBcPz4cQzjwbULIuLIlu49S9H3lyZIYlzNsXyw5iu+/m042e9GsydPMZp1mmY3iQku4c/8ns8qiRGRZJfoRKZevXpcv36devXqsXDhQgA6depE7969adCgAS+//LL2lxHJQEYsjqDb3LAES6vz3bjAvB8G0nXnfABmV27FS6+N43S2PDbjdH02iNkdq6RsZUUk0zIZiexGcXJy4vz58+TKlQuLxYKLy7+jUj/99BNbtmyhWLFivP3227i5JRwXT0tRUVH4+vpy48YNfHx80ro6Iunewza4a/TPFsYtnYpvzC1uuHvRr1lvVhZ7xmacYv5eLAl5DjeXJHf8iogk+vd3khMZf3//ZKtkalAiI5J49ja4c4uLZfC6r+m0axEAuwOL07PlQCJ9bf97MF3HDIjIE0rs7+8kTfbdv38/58+ff+g95cqVS0pIEUkHzBaDNp9tYreNXpgC184xY+FYyp0/AsDnVV9gwnNvEOec8J+PQjk8Wd2vrubCiEiqSVIiExwc/NAJvSaTCbPZ/MSVEpHUY+/EaoBmBzYyZvk0st67w1VPH/o2683aIrbnuwSX8NdcGBFJdUlKZLZv346fn19K1UVEUpm9E6vdY2MYsuYrXgtfBsCOfKUIaTGA8z65bMbRUJKIpJUkJTIFChRwuDkyIpLQvTgLzaZu4PClWwmuFb5yhpkLxlDy0gksmPi0ehsm13oVs5Nzgnv9vFzZ9kEDDSWJSJrRoZEimYjZYtBz7i6WRlyweb3V32sZtWImXrF3uZzFl97N+7IxqKLNe+sWz0Vop2opWV0RkUdKdCJTu3btdLe0WkQSz96KJACP2Lt8vGoW7fb+CcCWAuV4r0U/LnnnsBmrS60ghjQvlYK1FRFJnEQnMmvXrk3JeohICrkXZ6HZtA0cvphwGAmg6OVTzFwwhuKXT2HBxLSa7ZhWox0WG0NJOitJRNIbDS2JZGD2JvPe99K+VQxf+RlZYmO46JWd91r0Y2vBp23eWzF/Nua9W0PzYUQkXVEiI5IBmS0GwRPXcOLKXZvXs9y7w4iVn/FixBoANhSqQJ/mfbjslT3BvSZgmlYliUg6pURGJIN52FwYgOKXTjBz/hiKXj2D2eTExGdf47NnXsIwJTxKoEI+X37tVlO9MCKSbj1RInN/czyTSf/IiaS1h+3OC4Bh0G7PCj5ePQuPuHuc885JSMv+7MxfxubtmtArIo7gsU5zmz17NmXKlMHDwwMPDw/KlCnDV199ldx1E5FEWhweSdH3l9pNYrxjbjNt0XjGrJiBR9w91hauRNNO02wmMa5OJj5tX1FJjIg4hCT3yAwdOpRJkybRs2dPqlevDsDWrVvp3bs3p06dYvjw4cleSRGx781vdrDqwCW710tfOMqMBWMIunaOOJMT42p34Muqz9scSgqpW5T3GjyloSQRcRiJPv36Pj8/P6ZNm8Yrr7xiVf7jjz/Ss2dPLl++nKwVfFI6/Voysi5zdrD6oJ0kxjB4LWwpQ9Z8ibs5jjM+foS0HMDuvCUT3KoVSSKS3qTI6dcAsbGxVK5cOUF5pUqViIuLS2o4EXlMHy2MsJvE+Ny9yejl02l2aDMAK4tWo1/TXtzwzGp1n4sTTH5ZK5JExHElOZF5/fXX+eyzz5g0aZJV+axZs3j11VeTrWIiYl/n0O2sOWS797PcuX+YsWAsBW5c4J6TC2PqdOLryi3hgUn5GkYSkYzgsVYtzZ49mz///JNnnnkG+PdU7FOnTvHGG2/Qp0+f+PseTHZE5MmYLQbBE9Zw4qqN/WEMg85/LWTQulDcLHGc8s1Nj1YD2RvwlNVthXJ4srpfXSUwIpIhJDmRiYiIoGLFfw+RO3r0KAC5cuUiV65cRERExN+nJdkiyWt5xDm6f78bs41rvneimbB0Cg2ObAdg6VM1GNQkhCgPb6v7gkv4M7tjlVSorYhI6khyIqMzl0RS3/KIc7zz/W6b1ypGHmDawnHki7pEjLMLn9R7k+8qNLMaStLuvCKSUWlnX5F0zmwxbCYxJsNC1x1/0H/Dt7hazBzPHkCPVoP4O3cRq/u0O6+IZGSJSmReeOEF5syZg4+PDy+88MJD7/3999+TpWIi8q9yHy1LUJb99g0mLplMvWN/AbCw5HO836gHN92zWN3XoUZBhrW0vXOviEhGkKidfX19fePnvPj6+j70JykKFSqEyWRK8NO9e3cA7t69S/fu3cmZMyfe3t68+OKLXLhwIYlNFHFMZotBmY+WcSvWequnKqcjWBoaQr1jf3HXxY3BjXoQ0qJ/giQmuISfkhgRyfASvSHe8OHD6devH1myZHn0zYl06dIlzOb/TV2MiIigQYMGrF27ljp16vDuu++yZMkS5syZg6+vLz169MDJyYnNmzcn+jO0IZ44ouUR53j3+91WBz+aDAvvbvuVPhu/x8WwcDRHPrq3GshB/6AE7w8u4cfsjlVTr8IiIskssb+/E53IODs7c+7cOfz9/ZOtkg/q1asXixcv5vDhw0RFReHn58fcuXN56aWXADh48CAlS5Zk69at8Uu/H0WJjDiapXvP0m1umFVZzlvXmbx4Is+d+Lf8t9J1GdKwG7fdPBO8X4c9ikhGkOw7+ybxJIMku3fvHt9//z19+vTBZDKxa9cuYmNjqV+/fvw9JUqUoECBAklKZEQcycLdkYT8Em5VVv3kXqYuGo//rWvccXFnaIN3mFe2foIN7kzAzPYVaVouIPUqLCKSxpK0aikl94aZP38+169fp2PHjgCcP38eNzc3smXLZnVf7ty5OX/+vN04MTExxMTExL+OirJ9GrBIemK2GLT5bJPV6dVOFjM9t/xMyJafcDYsHMpVgO6tBnEkV4EE7y+QzZ21A4K1MklEMp0kJTJPPfXUI5OZq1evPlZFZs+eTZMmTQgMfLJ9LkaPHs2wYcOeKIZIalocHknPn8Kt5sP43bzK1EUTqHFqLwA/l23ARw3e5q6rR4L35/V1Y8Og+gnKRUQygyQlMsOGDUvyyqTEOHnyJKtWrbJaup0nTx7u3bvH9evXrXplLly4QJ48eezGGjx4sNUxCVFRUeTPnz/Z6yySHN78ZgerDlgf/FjreBiTF0/E7/Z1brl68EGj7swvXdfm+7N7OrN5cIPUqKqISLqUpESmXbt2KTLZNzQ0FH9/f5o1axZfVqlSJVxdXVm9ejUvvvgiAIcOHeLUqVNUr17dbix3d3fc3d2TvY4iye3Bgx+dLWZ6bZpL962/4ITBAb9CdG81iGM589l8f/YsLoQNbZRa1RURSZcSncik1PwYi8VCaGgoHTp0wMXlf9Xx9fWlS5cu9OnThxw5cuDj40PPnj2pXr26JvqKQzNbDF78dBPhZ/43HyZP1GWmLhpPtTN/A/B9+SaMqPcmMa62k/K82dzZrOEkEZG0X7W0atUqTp06RefOnRNcmzx5Mk5OTrz44ovExMTQqFEjPv300xSph0hqsDUfps7RnUxaMpkcd6KIdvNkcOOeLC75nN0YpQK8Wfpe7ZSvrIiIA0j0PjKOSvvISHpga1WSizmOfhu+5Z0d/84N25e7CD1aDeRkdvsT3ssEZmVxiP0kR0Qko0j2fWRE5PEs3XuWHnPDsPynLDDqItMXjKPS2YMAhFZqweg6nbnn4mo3TnAJf2Z3rJLCtRURcSxKZERS0IjFEczedNKqrP7h7UxYOplsd28S5e5F/ybvsaJ4DbsxTMC0VyrQ4ukn25pARCQjUiIjkgLMFoOXPttE2H+GklzNsQxaN4cufy0AIDzgKXq0HMCZbPa3E6iQz5dfu9XURnciInYokRFJZrYm9Oa7fp4ZC8dS/txhAL6s0ppxtTsQ62x/KKlTzUJ81KJ0CtdWRMSxKZERSSa2JvQCND60mXHLpuETc4vrHt70bdab1UWrPTRW12eD+KCZDn4UEXkUJTIiycDWhF73uHu8v3Y2HXYvAeCvvCUJadmfsz72N5V0Ambo4EcRkURTIiPyhEYu+ZsvN56wKit0NZIZC8dR5sJRAD6r9hITn32NOGf7j1zF/NmY924NzYcREUkCJTIiT2DYoghCN1uvSmqxfz2jVswg6707XPH0oU/zvqwvXMluDK1KEhF5fEpkRB6DrVVJ7rExfLT6S9rvWQ7A9vxlCGnRjwtZc9mNo1VJIiJPRomMSBLZWpVU5MppZiwYS8lLJ7BgYkb1tkyt1R6zk7PdOF1qBTGkuSb0iog8CSUyIknQZc52Vh+8bFX2fMQaPvnzU7xi73LJKxu9mvdjc6HydmNoQq+ISPJRIiOSCGaLQfCENZy4eje+zPPeXYav/Jw2EasA2FywHL2a9+eSd3a7cTShV0QkeSmREXkEW0NJxS6dZOaCsTx15RRmkxNTar7CzOptsdgZStKEXhGRlKFERsQOmxvcGQZt9q1k+Mov8IyL4YJ3Dt5r0Y9tBcrZjVM+rw+/da+lXhgRkRSgREbEBlu9MF4xt/nkz095fv86ANYHVaRPsz5c8cpmN0694n583alqitZVRCQzUyIj8gBbE3pLXjzGjAVjKXI1kjiTExOfe53Pq72IYXKyG6d+SX++6lAlpasrIpKpKZER+Y9mU9bx9/lb/yswDNrvWc5Hq2bhbo7lbNZchLTsz1/57B/mqPkwIiKpR4mMCP/Oh6k2ciWXb8XGl3nH3GbM8uk0P7gRgFVFqtCvWW+ue/rYjaNVSSIiqUuJjGR6i8Mj6fFTuFVZ6fNHmLlgLIWunyPWyZmxtTvwVZXnwWQ/QZmuXhgRkVSnREYyLXurkjrsXsz7a2fjbo7jjI8/PVsOICxvCbtxCuXwZHW/uuqFERFJA0pkJFOytSrJ5+5Nxi6bRpN/tgCwotgz9G/aiygPb7txgkv4M7ujJvSKiKQVJTKS6dhalfT02UPMWDiO/DcucM/JhVF1OzOnUgu7Q0ma0Csikj4okZFMw9YxAxgGXXbOZ9D6ObhazJzMloceLQeyL6CY3Tia0Csikn4okZFMwdZQku+daCYsnUyDIzv+vad4LQY36Um0u5fdOJrQKyKSviiRkQzN5oReoOKZA0xfOI680ZeIcXZleHBXfijfxO5QUq4sLmz/sKF6YURE0hklMpJhLd17lh5zw7D8p8xkWHh7++/02/AtLoaFY9kD6dFqEPtzF7YbJ392DzYODE75CouISJIpkZEMaeSSv/ly4wmrshy3bzBp8STqHN8FwPxStfmgYXduuWexG6deCT++7qizkkRE0islMpLhDFsUQejmk1ZlVU9HMG3hOPLcvMpdFzeG1n+HX8o10KokEREHp0RGMgyzxeClzzYR9p/5ME4WM922zaP3prk4GxYO58xP91YD+cevkN04WpUkIuI4lMhIhmBrVVKuW9eYsmgCtU7uAWBemfoMbfAOd9w87MbRqiQREceiREYcmr1VSTVOhDN18QT8bl3ntqs7Hzbsxu9l7E/Y1TEDIiKOSYmMOCxbq5KcLGbe2/wTPbf8hBMGB3MVpHurQRzNld9uHB0zICLiuJTIiEMasTiC2ZusJ/T6R19h2qLxPHM6AoC5TzdiWPBbxLi6242joSQREcemREYciq0JvQDPHdvFpCWTyHX7BjfdPHm/UQ8WlqptN46GkkREMgYlMuIwbE3odbaY6bvxO7pt+xWAv/0L06PVQI7nyGs3Tr3ifnzdSXvDiIhkBEpkxCHYOrE6IOoS0xaOp0rkfgC+rdCMkfW6EOPiZj9OrSCGNC+VonUVEZHUo0RG0r3mU9cTce6mVVndozuZtHgS2e9GE+WWhYFNQlhWopbdGE7AjPYVaVouIIVrKyIiqUmJjKRbZotB8IQ1nLh6N77M1RxL//Xf8tbOPwDYk6cYPVoN5HS2PHbjaIM7EZGMS4mMpEu25sPku3GB6QvGUeHcIQBmV27F2NoduefiajOGjhkQEcn4lMhIumJvg7uG/2xl/NIp+Mbc4oa7F/2a9WZlsWfsxqmQz5dfu9VUL4yISAanREbSDVu9MG5xsQxe9zWddi0CYHdgcXq2HEikr7/dOJrQKyKSeSiRkXTB1qqkAtfOMWPhWMqdPwLA51VfYMJzbxDnbPuvrSb0iohkPk5pXYHIyEhee+01cubMiaenJ2XLluWvv/6Kv96xY0dMJpPVT+PGjdOwxpKczBaDOuNWJ0himh7cxOI571Hu/BGuevrQ6aWPGFO3s90kpkI+Xw6PaqokRkQkk0nTHplr165Rs2ZN6taty7Jly/Dz8+Pw4cNkz57d6r7GjRsTGhoa/9rd3f6W8+I4bJ2V5B53jw/XfMXrYUsB2JGvFCEtBnDeJ5fdOBpKEhHJvNI0kRk7diz58+e3SlKCgoIS3Ofu7k6ePPaX14rjsXVWUtDVSGYuGEOpi8exYOLT6m2YXOtVzE7ONmOYgJkaShIRydTSdGhp4cKFVK5cmTZt2uDv70+FChX48ssvE9y3bt06/P39KV68OO+++y5XrlyxGzMmJoaoqCirH0k/zBaD52duTJDEtNy/jkXf9KLUxeNczuJLh7bDmPDcG3aTmEI5PDmioSQRkUzPZBiG8ejbUoaHhwcAffr0oU2bNuzcuZP33nuPzz//nA4dOgDw008/kSVLFoKCgjh69Cjvv/8+3t7ebN26FWfnhL/kPv74Y4YNG5ag/MaNG/j4+KRsg+ShbK1K8oi9y0erZvHK3j8B2FqgLO8178fFrDntxgku4c/sjlVSuLYiIpKWoqKi8PX1feTv7zRNZNzc3KhcuTJbtmyJLwsJCWHnzp1s3brV5nuOHTtGkSJFWLVqFcHBwQmux8TEEBMTE/86KiqK/PnzK5FJY51Dt7HmkHVPWpHLp5m5YAwlLp/EgolpNdsxrUY7LA8ZStIGdyIimUNiE5k0nSMTEBBAqVLWkzRLlizJb7/9Zvc9hQsXJleuXBw5csRmIuPu7q7JwOmI2WLwzMiVXLoVa1X+4r7VjFj5KVliY7jolZ33WvRja8Gn7cbRMQMiImJLmiYyNWvW5NChQ1Zl//zzDwULFrT7njNnznDlyhUCAjQ3Ir1bHB5Jj5/Crco8793lk5Wf8mLEGgA2FixP7xZ9ueyV3UaEf01XL4yIiNiRpolM7969qVGjBqNGjaJt27bs2LGDWbNmMWvWLABu3rzJsGHDePHFF8mTJw9Hjx5lwIABFC1alEaNGqVl1eUh7B0zUPzSCWbOH0PRq2cwm5yYVOtVPq3eBsNke855oRyerO5XV70wIiJiV5rOkQFYvHgxgwcP5vDhwwQFBdGnTx+6du0KwJ07d2jdujVhYWFcv36dwMBAGjZsyIgRI8idO3ei4id2jE2Sx/KIc3T/fjfm/xYaBi/v/ZNhq77AI+4e571zENJyADvyl7Ebp15xP77uVDXF6ysiIumTQ0z2TQ1KZFLP0r1n6TY3zKrMK+Y2o1bMpNWB9QCsLVyJvs36cDWLr9042uBOREQcYrKvZBwLd0cS8ku4VVmpC8eYsWAMha+dJc7kxPjabzCr6gt2h5J0VpKIiCSVEhl5IjbnwxgGr4UtZciar3A3xxKZ1Y+eLQewO19Ju3G0KklERB6HEhl5bLY2uMsac4sxy6bR7NBmAFYWrUa/pr244ZnVbhytShIRkcelREYeS5c52xOcWF323GFmLBxLwevnuefkwtg6HZlduRWYbPey+Hm5su2DBuqFERGRx6ZERpLEbDEInrCGE1fv/q/QMOi0ayGD14biZonjtG9uerQcwJ7A4nbj1Cvhx9cdtSpJRESejBIZSTRbQ0k+d28yfukUGh3eBsCyp2owsEkIUR7eduNoKElERJKLEhlJFFtDSRUiDzJ94TjyRV0kxtmFT+q9yXcVmtkdSiqY3YM1/etpKElERJKNEhl5KFtDSSbDwps75jNgwze4WsycyBZA91YD+TtPUbtxygRmZXHIc6lRZRERyUSUyIhdtoaSst++wYSlUwg+uhOARSWeZXDjntx0z2I3TnAJf2Z3rJLCtRURkcxIiYzYZGsoqfKZv5m2cDyB0ZeJcXbl4/pv8+PTjewOJZmAaZoPIyIiKUiJjFixN5T07rZf6bPxe1wMC0dz5KN7q4Ec9A+yG0cb3ImISGpQIiPxbA0l5bx1ncmLJ/LciX/PUPqtdF2GNOzGbTdPu3G0KklERFKLEhmxfcwA8MypvUxdNIHcN69yx8WdoQ3eYV7Z+naHkgrl8GR1v7rqhRERkVSjRCaTs9UL42Qx03PLz4Rs+Qlnw8I/OQvQvdVADvsVtBunXnE/vu6kDe5ERCR1KZHJxGxN6PW7eY0pi8dT8+ReAH4u24CPGrzNXVcP+3FqBTGkeakUrauIiIgtSmQyqWZT1vH3+VtWZTVPhDNl0QT8bl/nlqsHHzTqzvzSde3GcAJmtK9I03IBKVxbERER25TIZDJmi0G1kSu5fCs2vszZYqbXprl03/oLThgc8CtEj1YDOZozv904WpUkIiLpgRKZTGRxeCQ9fgq3KssdfZlpC8dT7czfAPxQvjHD63UlxtXdbhytShIRkfRCiUwmYG9VUp2jfzFxySRy3oki2s2TwY17srik/WMEtCpJRETSGyUyGdzSvWfpMTcMy3/KXMxx9Nv4He9s/w2AfbmL0KPVQE5mt9/LomMGREQkPVIik4GNWBzB7E0nrcoCoy4yfcE4Kp09CEBopRaMrtOZey6uNmPomAEREUnPlMhkQGaLwUufbSLsgaGk+oe3M2HpZLLdvUmUuxf9m7zHiuI17MbRhF4REUnvlMhkMLY2uHM1xzJw3Rze/GsBAOEBxejRciBnsuWxG0cTekVExBEokckg7E3ozXf9PDMWjqP8uX8A+LJKa8bV7kCss+2hpFxZXNj+YUP1woiIiENQIpMB2JrQC9Do0BbGL5uKT8wtrnt407dZb1YXrWY3Tv7sHmwcGJyylRUREUlGSmQc3Mglf/PlxhNWZW5xsby/djYddy8GYFdgCXq2GsBZH3+7ceqV8OPrjjorSUREHIsSGQc2bFEEoZutVyUVvHaWGQvGUvbCUQA+q/YSE599jThn21+1ViWJiIgjUyLjgOytSmp+YAOjl08n6707XPH0oW+zPqwrUtlunAr5fPm1W03NhxEREYelRMbB2FqV5B4bw9A1X/Jq+HIAtucrTUjL/lzImstuHJ1YLSIiGYESGQfSZc52Vh+8bFVW+MoZZi4YQ8lLJ7BgYkb1tkyt1R6zk7PNGDqxWkREMhIlMg7AbDEInrCGE1fvWpW3/nstI1fMxCv2LpeyZKNXi35sLlTebhxtcCciIhmNEpl0ztZQkkfsXYav/Jy2+1YBsLlgOXo178cl7xw2Y2hCr4iIZFRKZNIpexvcFbt0kpkLxvLUlVOYTU5MrfkKM6q3xWJnKKl8Xh9+615LvTAiIpIhKZFJh2xucGcYtNm3iuErP8czLoYL3jl4r0U/thUoZzdOveJ+fN1Je8OIiEjGpUQmnbF1YnWWe3f45M9PeeHvtQBsKFSB3s37csUrm9049Uv681WHKilZVRERkTSnRCYd6Ry6nTWHrFcllbh4nJkLxlLk6hniTE5MfO51Pq/2IobJyWYMzYcREZHMRIlMOmBzVZJh0H7Pcj5aNQt3cyxns+YipGV//spX2m4crUoSEZHMRolMGrO1Ksk75jajl0+nxcGNAKwuUoV+TXtxLYuv3TjT1QsjIiKZkBKZNGJvVVLp80eYsXAsQdfOEevkzNjaHZhdpbXdoaRCOTxZ3a+uemFERCRTUiKTBuytSnpj92I+WDsbd3McZ3z86dlyAGF5S9iNE1zCn9kdNaFXREQyLyUyqczWqiSfuzcZu2waTf7ZAsCKYs/Qv2kvojy8bcbQhF4REZF/2R6vSEWRkZG89tpr5MyZE09PT8qWLctff/0Vf90wDIYOHUpAQACenp7Ur1+fw4cPp2GNH4/ZYvD8zI0Jkpinzx5iyZz3aPLPFu45ufBx8Fu8/fwHdpOYCvl8OTKqqZIYERER0rhH5tq1a9SsWZO6deuybNky/Pz8OHz4MNmzZ4+/Z9y4cUybNo1vvvmGoKAghgwZQqNGjdi/fz8eHh5pWPvEszWhF8Ogy18LGLhuDm6WOE5my0OPlgPZF1DMbpxONQvxUQv7q5ZEREQyG5NhGMajb0sZgwYNYvPmzWzcuNHmdcMwCAwMpG/fvvTr1w+AGzdukDt3bubMmUO7du0e+RlRUVH4+vpy48YNfHx8krX+j2JvQq/vnWgmLJ1MgyM7AFhSvCaDmoQQ7e5lN1bXZ4P4oFmpFK2viIhIepHY399pOrS0cOFCKleuTJs2bfD396dChQp8+eWX8dePHz/O+fPnqV+/fnyZr68v1apVY+vWrWlR5URbHnGOp95fmiCJqXjmAEtDQ2hwZAcxzq582LAb3VsNspvEOAGftq+oJEZERMSGNB1aOnbsGJ999hl9+vTh/fffZ+fOnYSEhODm5kaHDh04f/48ALlz57Z6X+7cueOvPSgmJoaYmJj411FRUTbvS0lL956l29wwqzKTYeGtHb/Tf/23uBgWjmUPpEerQezPXdhuHG1wJyIi8nBpmshYLBYqV67MqFGjAKhQoQIRERF8/vnndOjQ4bFijh49mmHDhiVnNZNk4e5IQn4JtyrLcfsGE5dMou6xXQDML1WbDxp255Z7FpsxtCpJREQkcdJ0aCkgIIBSpayHTEqWLMmpU6cAyJMnDwAXLlywuufChQvx1x40ePBgbty4Ef9z+vTpFKh5QmaLwQszNyZIYqqejmBpaE/qHtvFXRc3BjQOoVfzfnaTGK1KEhERSbw07ZGpWbMmhw4dsir7559/KFiwIABBQUHkyZOH1atXU758eeDfoaLt27fz7rvv2ozp7u6Ou7t7itb7QbZWJZkMC922zqPPph9wNiwcyZGPbq0H8Y9fIbtxutQKYkhzzYURERFJrDRNZHr37k2NGjUYNWoUbdu2ZceOHcyaNYtZs2YBYDKZ6NWrF5988gnFihWLX34dGBhI69at07LqgP1VSbluXWPS4kk8d+LfeTK/lglmSIN3ueNme7m4EzCjfUWalgtI6SqLiIhkKGmayFSpUoU//viDwYMHM3z4cIKCgpgyZQqvvvpq/D0DBgzg1q1bvPXWW1y/fp1atWqxfPnyNN9DxuYxA0D1k3uYumgC/reucdvVnSENuvFb2WC7cTShV0RE5PGl6T4yqSEl9pEZueRvvtx4wqrMyWLmvc0/0XPLTzhhcDBXQbq3GsTRXPntxtGJ1SIiIrYl9ve3zlpKopFL9idIYvyjrzB18QSqn9oHwI/lGjKs/lvcdbXda1Qwuwdr+tdTL4yIiMgTUiKTBEv3nuPLjcetyp49vpvJiyeS6/YNbrp58n6j7iwsVcdujDKBWVkc8lwK11RERCRzUCKTSGaLwYcLIuJfO1vM9Nn4Pd23zQNgv38Q3VsN4niOvHZjBJfwZ3bHKileVxERkcxCiUwi7Th+lau37v37wjD48rcR1Dv27ynd31Zoxsh6XYhxcbP5Xm1wJyIikjKUyCTSxei7/3thMvFH6bpUPrOfQU1CWFqilt33aVWSiIhIylEik0j+Wa0n7i4qVZtNhcpzLYuvzfvVCyMiIpLylMgkUtWgHOTx8eB81P96ZuwlMRXy+fJrt5rqhREREUlhaXrWkiNxdjLxcctHHx/QpVYQf/SopSRGREQkFSiRSYLGZQL4/LWKZMvimuCat7sLn7avqLOSREREUpGGlpKocZkAGpTKw7ZjV9h69ApgUL1wLp4pklO9MCIiIqlMicxjcHYyUbNoLmoWzZXWVREREcnUNLQkIiIiDkuJjIiIiDgsJTIiIiLisJTIiIiIiMNSIiMiIiIOS4mMiIiIOCwlMiIiIuKwlMiIiIiIw1IiIyIiIg4rw+/saxgGAFFRUWlcExEREUms+7+37/8etyfDJzLR0dEA5M+fP41rIiIiIkkVHR2Nr6+v3esm41GpjoOzWCycPXuWrFmzEh0dTf78+Tl9+jQ+Pj5pXbVkFRUVlWHbBhm7fRm5bZCx25eR2wYZu30ZuW2QMdpnGAbR0dEEBgbi5GR/JkyG75FxcnIiX758AJhM/55O7ePj47Bf7KNk5LZBxm5fRm4bZOz2ZeS2QcZuX0ZuGzh++x7WE3OfJvuKiIiIw1IiIyIiIg4rUyUy7u7ufPTRR7i7u6d1VZJdRm4bZOz2ZeS2QcZuX0ZuG2Ts9mXktkHGb99/ZfjJviIiIpJxZaoeGREREclYlMiIiIiIw1IiIyIiIg5LiYyIiIg4rHSdyGzYsIEWLVoQGBiIyWRi/vz5Vtc//vhjSpQogZeXF9mzZ6d+/fps3749/vqJEyfo0qULQUFBeHp6UqRIET766CPu3btndY/JZErws23bNqvPmjdvHiVKlMDDw4OyZcuydOnSNG8fQKFChRLUfcyYMVb37N27l2effRYPDw/y58/PuHHjEtQludv3pG1bt26dze/FZDKxc+dOIH1/d//1zjvvYDKZmDJlilX51atXefXVV/Hx8SFbtmx06dKFmzdvWt2THr+7R7XN0Z+7R7UP0u9zlxztS8/P3qPa1rFjxwR1aty4sdU96fW5S472pfdnL6Wk60Tm1q1bPP3008ycOdPm9aeeeooZM2awb98+Nm3aRKFChWjYsCGXLl0C4ODBg1gsFr744gv+/vtvJk+ezOeff87777+fINaqVas4d+5c/E+lSpXir23ZsoVXXnmFLl26EBYWRuvWrWndujURERFp2r77hg8fblX3nj17xl+LioqiYcOGFCxYkF27djF+/Hg+/vhjZs2alaLte9K21ahRw6pN586d48033yQoKIjKlStbxUqP3919f/zxB9u2bSMwMDDBtVdffZW///6blStXsnjxYjZs2MBbb70Vfz29fnePapujP3ePat996fG5S472pednLzFta9y4sVWdfvzxR6vr6fW5S472pfdnL8UYDgIw/vjjj4fec+PGDQMwVq1aZfeecePGGUFBQfGvjx8/bgBGWFiY3fe0bdvWaNasmVVZtWrVjLfffjtRdU+Mx21fwYIFjcmTJ9t9z6effmpkz57diImJiS8bOHCgUbx48fjXKd2+5Pju7t27Z/j5+RnDhw+PL0vv392ZM2eMvHnzGhEREQm+p/379xuAsXPnzviyZcuWGSaTyYiMjDQMI31/dw9rmy2O9tw9qn2O8NwZRvJ8f+n12bPVtg4dOhitWrWy+x5Hee4M4/HaZ0t6ffaSU7rukUmKe/fuMWvWLHx9fXn66aft3nfjxg1y5MiRoLxly5b4+/tTq1YtFi5caHVt69at1K9f36qsUaNGbN26NXkqnwgPa9+YMWPImTMnFSpUYPz48cTFxcVf27p1K8899xxubm5WdT906BDXrl2Lvyct25eY727hwoVcuXKFTp06JbiWHr87i8XC66+/Tv/+/SldunSC61u3biVbtmxW/4Vbv359nJyc4ofY0ut396i22eJIz11i2+eoz11Svz9He/bWrVuHv78/xYsX59133+XKlStWdXLU5+6+h7XPFkd69h6Xwx8auXjxYtq1a8ft27cJCAhg5cqV5MqVy+a9R44cYfr06UyYMCG+zNvbm4kTJ1KzZk2cnJz47bffaN26NfPnz6dly5YAnD9/nty5c1vFyp07N+fPn0+5hv2/R7UvJCSEihUrkiNHDrZs2cLgwYM5d+4ckyZNiq97UFBQgrrfv5Y9e/Y0a19SvrvZs2fTqFGj+ANAIX1/d2PHjsXFxYWQkBCb18+fP4+/v79VmYuLCzly5IivW3r97h7Vtgc52nOXmPY58nOX1O/PkZ69xo0b88ILLxAUFMTRo0d5//33adKkCVu3bsXZ2dmhnzt4dPse5GjP3uNy+ESmbt26hIeHc/nyZb788kvatm3L9u3bE/xljYyMpHHjxrRp04auXbvGl+fKlYs+ffrEv65SpQpnz55l/Pjx8V9qWnpU+/5b93LlyuHm5sbbb7/N6NGj0/3W1In97s6cOcOKFSv45ZdfrMrT63e3a9cupk6dyu7du+NPXM8okto2R3vuEts+R33ukvr9Odqz165du/j/X7ZsWcqVK0eRIkVYt24dwcHBaVav5JKU9jnas/ckHH5oycvLi6JFi/LMM88we/ZsXFxcmD17ttU9Z8+epW7dutSoUcNqwpY91apV48iRI/Gv8+TJw4ULF6zuuXDhAnny5EmeRjxEYtr3X9WqVSMuLo4TJ04A9ut+/9rD7knp9iW2baGhoeTMmTNRD1l6+O42btzIxYsXKVCgAC4uLri4uHDy5En69u1LoUKF4ut18eJFq/fFxcVx9erVR34v96897J6Ual9i2nafIz53SWnffznKc5fU9jnas/egwoULkytXrvh6OepzZ8+D7bvPEZ+9J+HwicyDLBYLMTEx8a8jIyOpU6cOlSpVIjQ0FCenRzc5PDycgICA+NfVq1dn9erVVvesXLmS6tWrJ1/FE+nB9j0oPDwcJyen+F6N6tWrs2HDBmJjY+PvWblyJcWLFyd79uzx96SH9tlqm2EYhIaG8sYbb+Dq6vrIGOnhu3v99dfZu3cv4eHh8T+BgYH079+fFStWxNfr+vXr7Nq1K/59a9aswWKxUK1atfh70tt3l5i2geM+d4ltn626O8Jzl5T2OeKz96AzZ85w5cqV+Ho56nNnz4PtA8d99p5IWs82fpjo6GgjLCzMCAsLMwBj0qRJRlhYmHHy5Enj5s2bxuDBg42tW7caJ06cMP766y+jU6dOhru7uxEREWEYxr8z84sWLWoEBwcbZ86cMc6dOxf/c9+cOXOMuXPnGgcOHDAOHDhgjBw50nBycjK+/vrr+Hs2b95suLi4GBMmTDAOHDhgfPTRR4arq6uxb9++NG3fli1bjMmTJxvh4eHG0aNHje+//97w8/Mz3njjjfjPuH79upE7d27j9ddfNyIiIoyffvrJyJIli/HFF1+kaPuetG33rVq1ygCMAwcOJPiM9Prd2WJrZUjjxo2NChUqGNu3bzc2bdpkFCtWzHjllVfir6fH7y4xbXPk5y4x7UvPz11ytO++9PjsPaxt0dHRRr9+/YytW7cax48fN1atWmVUrFjRKFasmHH37t34GOn1uUuO9qX3Zy+lpOtEZu3atQaQ4KdDhw7GnTt3jOeff94IDAw03NzcjICAAKNly5bGjh074t8fGhpq8/3/zd/mzJljlCxZ0siSJYvh4+NjVK1a1Zg3b16Cuvzyyy/GU089Zbi5uRmlS5c2lixZkubt27Vrl1GtWjXD19fX8PDwMEqWLGmMGjXK6qE1DMPYs2ePUatWLcPd3d3ImzevMWbMmBRv35O27b5XXnnFqFGjhs3PSK/fnS22fllcuXLFeOWVVwxvb2/Dx8fH6NSpkxEdHW11T3r77hLTNkd+7hLTvvT83CVH++5Lj8/ew9p2+/Zto2HDhoafn5/h6upqFCxY0Ojatatx/vx5qxjp9blLjval92cvpZgMwzAetzdHREREJC1luDkyIiIiknkokRERERGHpURGREREHJYSGREREXFYSmRERETEYSmREREREYelREZEREQclhIZEUk1JpOJ+fPn271uGAZvvfUWOXLkwGQyER4eniyf+/HHH1O+fHm719etW4fJZOL69euP/RnJEUNEkk6JjEgmZDKZHvrz8ccf233viRMnkjXJ+K/ly5czZ84cFi9ezLlz5yhTpkyyf4YtNWrU4Ny5c/j6+qbK54lI8nFJ6wqISOo7d+5c/P//+eefGTp0KIcOHYov8/b2TotqcfToUQICAqhRo0aqfq6bm1u6PdlXRB5OPTIimVCePHnif3x9fTGZTPGv/f39mTRpEvny5cPd3Z3y5cuzfPny+PcGBQUBUKFCBUwmE3Xq1AFg586dNGjQgFy5cuHr60vt2rXZvXt3ouvUsWNHevbsyalTpzCZTBQqVAiAQoUKMWXKFKt7y5cvb9VrdP36dd588038/Pzw8fGhXr167NmzJ9Gf/eCw0Jw5c8iWLRsrVqygZMmSeHt707hxY6sE8FFu375NkyZNqFmzZnzcLVu2UL58eTw8PKhcuTLz589Psd4tkcxCiYyIWJk6dSoTJ05kwoQJ7N27l0aNGtGyZUsOHz4MwI4dOwBYtWoV586d4/fffwcgOjqaDh06sGnTJrZt20axYsVo2rQp0dHRif7c4cOHky9fPs6dO8fOnTsTXec2bdpw8eJFli1bxq5du6hYsSLBwcFcvXo1ia3/n9u3bzNhwgS+++47NmzYwKlTp+jXr1+i3nv9+nUaNGiAxWJh5cqVZMuWjaioKFq0aEHZsmXZvXs3I0aMYODAgY9dPxH5l4aWRMTKhAkTGDhwIO3atQNg7NixrF27lilTpjBz5kz8/PwAyJkzp9VwTL169azizJo1i2zZsrF+/XqaN2/+yM/19fUla9asODs7J2mYZ9OmTezYsYOLFy/i7u4e34b58+fz66+/8tZbbyU61n/Fxsby+eefU6RIEQB69OjB8OHDH/m+8+fP8/LLL1OsWDHmzp2Lm5sbAHPnzsVkMvHll1/i4eFBqVKliIyMpGvXro9VPxH5l3pkRCReVFQUZ8+epWbNmlblNWvW5MCBAw9974ULF+jatSvFihXD19cXHx8fbt68yalTp1KyyuzZs4ebN2+SM2dOvL2943+OHz/O0aNHHztulixZ4pMYgICAAC5evPjI9zVo0ICiRYvy888/xycxAIcOHaJcuXJ4eHjEl1WtWvWx6yci/1KPjIgkiw4dOnDlyhWmTp1KwYIFcXd3p3r16ty7d++J4jo5OWEYhlVZbGxs/P+/efMmAQEBrFu3LsF7s2XL9tif6+rqavXaZDIlqIctzZo147fffmP//v2ULVv2sT9fRBJHiYyIxPPx8SEwMJDNmzdTu3bt+PLNmzfH9x7c72Uwm81W7928eTOffvopTZs2BeD06dNcvnz5ievk5+dnNck2KiqK48ePx7+uWLEi58+fx8XFJX6CcFoaM2YM3t7eBAcHs27dOkqVKgVA8eLF+f7774mJiYkfAkvKPCARsU1DSyJipX///owdO5aff/6ZQ4cOMWjQIMLDw3nvvfcA8Pf3x9PTk+XLl3PhwgVu3LgBQLFixfjuu+84cOAA27dv59VXX8XT0/OJ61OvXj2+++47Nm7cyL59++jQoQPOzs7x1+vXr0/16tVp3bo1f/75JydOnGDLli188MEH/PXXX0/8+Y9jwoQJvPrqq9SrV4+DBw8C0L59eywWC2+99RYHDhxgxYoVTJgwAfi3t0dEHo8SGRGxEhISQp8+fejbty9ly5Zl+fLlLFy4kGLFigHg4uLCtGnT+OKLLwgMDKRVq1YAzJ49m2vXrlGxYkVef/11QkJC8Pf3f+L6DB48mNq1a9O8eXOaNWtG69atreaumEwmli5dynPPPUenTp146qmnaNeuHSdPniR37txP/PmPa/LkybRt25Z69erxzz//4OPjw6JFiwgPD6d8+fJ88MEHDB06FMBq3oyIJI3JSMygr4iIJLsffviBTp06cePGjWTpvRLJjDRHRkQklXz77bcULlyYvHnzsmfPHgYOHEjbtm2VxIg8ASUyIiKp5Pz58wwdOpTz588TEBBAmzZtGDlyZFpXS8ShaWhJREREHJYm+4qIiIjDUiIjIiIiDkuJjIiIiDgsJTIiIiLisJTIiIiIiMNSIiMiIiIOS4mMiIiIOCwlMiIiIuKwlMiIiIiIw/o/O1UJyKGk0B4AAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Parameters are [ 1.57613580e-02 -1.51221174e+02]\n"
     ]
    }
   ],
   "source": [
    "total2tip_x, total2tip_y = parseCSV(\"../data/md-11 tip to total 1_30h.csv\")\n",
    "\n",
    "popt, _ = sp.optimize.curve_fit(func_lin, total2tip_x, total2tip_y, maxfev=100000)\n",
    "total_tip_func = func_lin(total2tip_x, *popt)\n",
    "\n",
    "plt.title(f\"Total fuel to tip tank percent\")\n",
    "plt.xlabel(\"Total fuel in kg\")\n",
    "plt.ylabel(\"Tip Tank percent\")\n",
    "plt.scatter(total2tip_x, total2tip_y)\n",
    "plt.plot(total2tip_x, total_tip_func, c=\"r\")\n",
    "plt.show()\n",
    "\n",
    "print(\"Parameters are\", popt)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Output to files.\n",
    "\n",
    "See `Format description.md` for a description of the file format."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "for tips in tip_tanks_index:\n",
    "  raw = open(f\"curve_params_{index2mach(tips)}.data\", \"w\")\n",
    "  for alt in alts_index:\n",
    "    _raw = []\n",
    "    lut = open(f\"ci2mach_{index2alt(alt)}_{index2mach(tips)}.data\", \"w\")\n",
    "    for weight in weights_index:\n",
    "      try:\n",
    "        current_list = curves[tips][alt][weight][0]\n",
    "        mmo = curves[tips][alt][weight][0][999]\n",
    "        lut.writelines(\",\".join(str(round(min(x, mmo), 3)) for x in current_list.tolist()) + \"\\n\")\n",
    "        _raw.append(','.join(str(x) for x in curves[tips][alt][weight][1].tolist()))\n",
    "      except IndexError:\n",
    "        pass\n",
    "    lut.close()\n",
    "    raw.writelines(\";\".join(_raw) + \"\\n\")\n",
    "  raw.close()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Header file"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [],
   "source": [
    "for tips in tip_tanks_index:\n",
    "  lut = open(f\"ci2mach_{tip2mach(index2mach(tips))}.h\", \"w\")\n",
    "  lut.writelines(f\"#ifndef CI2MACH_{int(tip2mach(index2mach(tips)) * 100)}\\n\")\n",
    "  lut.writelines(f\"#define CI2MACH_{int(tip2mach(index2mach(tips)) * 100)}\\n\\n\")\n",
    "\n",
    "  lut.writelines(\"#ifndef NULL\\n\")\n",
    "  lut.writelines(\"#define NULL 0\\n\")\n",
    "  lut.writelines(\"#endif\\n\\n\")\n",
    "\n",
    "  lut.writelines(f\"class CI2Mach_{int(tip2mach(index2mach(tips)) * 100)} {{\\n\")\n",
    "  lut.writelines(\"public:\\n\")\n",
    "  lut.writelines(f\"  float values[{len(alts_index)}][{len(weights_index)}][1000] = {{\\n\")\n",
    "  for alt in alts_index:\n",
    "    lut.writelines(\"  {\\n\")\n",
    "    for weight in weights_index:\n",
    "      try:\n",
    "        current_list = curves[tips][alt][weight][0]\n",
    "        mmo = curves[tips][alt][weight][0][999]\n",
    "        lut.writelines(\"    {\" + \",\".join(\"{:.3f}f\".format(round(min(x, mmo), 3)) for x in current_list.tolist()) + \"},\\n\")\n",
    "      except IndexError:\n",
    "        lut.writelines(\"    {\" + \",\".join([\"-1\".format(x) for x in range(0, 1000)]) + \"},\\n\")\n",
    "    lut.writelines(\"  },\\n\")\n",
    "  lut.writelines(\"};\\n};\\n\\n\")\n",
    "\n",
    "  lut.writelines(f\"CI2Mach_{int(tip2mach(index2mach(tips)) * 100)}* ci2Mach_{int(tip2mach(index2mach(tips)) * 100)} = NULL;\\n\\n\")\n",
    "  lut.writelines(f\"void init_ci2mach_{int(tip2mach(index2mach(tips)) * 100)}() {{\\n\")\n",
    "  lut.writelines(f\"  ci2Mach_{int(tip2mach(index2mach(tips)) * 100)} = new CI2Mach_{int(tip2mach(index2mach(tips)) * 100)}();\\n\")\n",
    "  lut.writelines(\"}\\n\\n\");\n",
    "\n",
    "  lut.writelines(f\"void cleanup_ci2mach_{int(tip2mach(index2mach(tips)) * 100)}() {{\\n\")\n",
    "  lut.writelines(f\"  if (ci2Mach_{int(tip2mach(index2mach(tips)) * 100)} != NULL) {{\\n\")\n",
    "  lut.writelines(f\"    delete ci2Mach_{int(tip2mach(index2mach(tips)) * 100)};\\n\")\n",
    "  lut.writelines(f\"    ci2Mach_{int(tip2mach(index2mach(tips)) * 100)} = NULL;\\n\")\n",
    "  lut.writelines(\"  }\\n}\\n\\n\")\n",
    "\n",
    "  lut.writelines(\"#endif\")\n",
    "  lut.close()\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Demo area."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "alt=290, weight=250, ci=80: lut=0.8210969019894392, params=array([1.12456109e-03, 2.54815227e-05, 8.09000000e-01]), func=0.8210969019894392\n"
     ]
    }
   ],
   "source": [
    "mach = mach2index(90) # 0.87\n",
    "alt = 290\n",
    "weight = 250\n",
    "ci = 80\n",
    "\n",
    "lut = curves[mach][alt2index(alt)][weight2index(weight)][0][ci]\n",
    "params = curves[mach][alt2index(alt)][weight2index(weight)][1]\n",
    "func = func_root(ci, *curves[mach][alt2index(alt)][weight2index(weight)][1])\n",
    "\n",
    "print(f\"{alt=}, {weight=}, {ci=}: {lut=}, {params=}, {func=}\")"
   ]
  }
 ],
 "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.11.9"
  },
  "widgets": {
   "application/vnd.jupyter.widget-state+json": {
    "state": {},
    "version_major": 2,
    "version_minor": 0
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}