src/notes/ltb_benchmark/andes_benchmark.ipynb

{
 "cells": [
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# ANDES Power Flow Benchmark"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Reference:\n",
    "\n",
    "1. H. Cui, F. Li and K. Tomsovic, \"Hybrid Symbolic-Numeric Framework for Power System Modeling and Analysis,\" in IEEE Transactions on Power Systems, vol. 36, no. 2, pp. 1373-1384, March 2021, doi: 10.1109/TPWRS.2020.3017019.\n",
    "1. R. D. Zimmerman, C. E. Murillo-Sánchez and R. J. Thomas, \"MATPOWER: Steady-State Operations, Planning, and Analysis Tools for Power Systems Research and Education,\" in IEEE Transactions on Power Systems, vol. 26, no. 1, pp. 12-19, Feb. 2011, doi: 10.1109/TPWRS.2010.2051168.\n",
    "keywords: {Power system planning;Steady-state;Power system analysis computing;Power system simulation;Load flow;Open source software;Computer languages;Packaging;Costs;Mathematical model;Load flow analysis;optimal power flow;optimization methods;power engineering;power engineering education;power system economics;power system simulation;power systems;simulation software;software tools},\n",
    "1. Power Grid Lib - Optimal Power Flow, https://github.com/power-grid-lib/pglib-opf"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import platform\n",
    "import os\n",
    "\n",
    "from tqdm import tqdm\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "\n",
    "import andes\n",
    "\n",
    "import pypower.api as pyp\n",
    "\n",
    "import datetime\n",
    "import time"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Last run time: 2024-03-26 17:35:23\n",
      "Python: 3.9.18; OS: Darwin\n",
      "andes:1.9.1.post24+g7a87ad5d\n"
     ]
    }
   ],
   "source": [
    "print(\"Last run time:\", datetime.datetime.now().strftime(\"%Y-%m-%d %H:%M:%S\"))\n",
    "\n",
    "# OS information\n",
    "print(f\"Python: {platform.python_version()}; OS: {platform.system()}\")\n",
    "\n",
    "print(f'andes:{andes.__version__}')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "andes.config_logger(stream_level=50)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "directory = '../../cases/matpower/'\n",
    "\n",
    "# List all files and directories within the specified directory\n",
    "entries = os.listdir(directory)\n",
    "\n",
    "# Get a list of files in the directory along with their sizes\n",
    "files_with_sizes = [(file, os.path.getsize(os.path.join(directory, file))) for file in os.listdir(directory) if os.path.isfile(os.path.join(directory, file))]\n",
    "\n",
    "# Sort the list of files by size (the second element of the tuple)\n",
    "sorted_files_with_sizes = sorted(files_with_sizes, key=lambda x: x[1])\n",
    "\n",
    "# If you just want a list of filenames sorted by size\n",
    "sorted_files = [file for file, size in sorted_files_with_sizes]\n",
    "\n",
    "sorted_files.remove('.DS_Store')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "100%|██████████| 42/42 [00:34<00:00,  1.22it/s]\n"
     ]
    }
   ],
   "source": [
    "n_cases = len(sorted_files)\n",
    "\n",
    "res_cols = ['case', 'v_diff_max', 'a_diff_max', 't_ad', 't_pp', 'c_ad', 'c_pp']\n",
    "res = pd.DataFrame(columns=res_cols)\n",
    "\n",
    "for case_i in tqdm(range(n_cases)):\n",
    "    pf_case = sorted_files[case_i]\n",
    "    sa = andes.load(directory + pf_case,\n",
    "                    setup=True, no_output=True, default_config=True)\n",
    "    # --- PYPOWER ---\n",
    "    ppc = andes.io.matpower.system2mpc(sa)\n",
    "    ppopt = pyp.ppoption(VERBOSE=0, OUT_ALL=0, PF_ALG=1, OPF_ALG_DC=200)\n",
    "\n",
    "    t_ad = time.time()\n",
    "    sa.PFlow.run()\n",
    "    s_ad = time.time() - t_ad\n",
    "\n",
    "    t_pp = time.time()\n",
    "    ppc_sol, _ = pyp.runpf(ppc, ppopt)\n",
    "    s_pp = time.time() - t_pp\n",
    "\n",
    "    if (sa.exit_code == 0) & (ppc_sol['success']):\n",
    "        v_ad = sa.Bus.v.v\n",
    "        a_ad = sa.Bus.a.v - sa.Bus.a.v[0]\n",
    "        v_pp = ppc_sol['bus'][:, 7]\n",
    "        a_pp = (ppc_sol['bus'][:, 8] - ppc_sol['bus'][0, 8]) * andes.shared.deg2rad\n",
    "\n",
    "        v_diff = np.abs(v_ad - v_pp)\n",
    "        a_diff = np.abs(a_ad - a_pp)\n",
    "    else:\n",
    "        v_diff = np.array([-1])\n",
    "        a_diff = np.array([-1])\n",
    "\n",
    "    row_res = {'case': pf_case.split('.')[0],\n",
    "               'v_diff_max': np.max(v_diff),\n",
    "               'a_diff_max': np.max(a_diff),\n",
    "               't_ad': s_ad,\n",
    "               't_pp': s_pp,\n",
    "               'c_ad': int(1 - sa.exit_code),\n",
    "               'c_pp': int(ppc_sol['success'])}\n",
    "    res.loc[case_i] = row_res"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>case</th>\n",
       "      <th>v_diff_max</th>\n",
       "      <th>a_diff_max</th>\n",
       "      <th>t_ad</th>\n",
       "      <th>t_pp</th>\n",
       "      <th>c_ad</th>\n",
       "      <th>c_pp</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>case18</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.005228</td>\n",
       "      <td>0.008196</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>case22</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>0.009135</td>\n",
       "      <td>0.016628</td>\n",
       "      <td>True</td>\n",
       "      <td>False</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>case14</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.003891</td>\n",
       "      <td>0.005336</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>case30</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.004930</td>\n",
       "      <td>0.007562</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>case_ieee30</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.004656</td>\n",
       "      <td>0.006019</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>case69</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>0.003248</td>\n",
       "      <td>0.020533</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>case24_ieee_rts</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.005172</td>\n",
       "      <td>0.007628</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>case39</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.003234</td>\n",
       "      <td>0.003471</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>case85</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>0.003196</td>\n",
       "      <td>0.016355</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>case118zh</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>0.003179</td>\n",
       "      <td>0.017203</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>case57</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.004299</td>\n",
       "      <td>0.006711</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>case141</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>0.003319</td>\n",
       "      <td>0.017596</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>case136ma</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>-1.000000</td>\n",
       "      <td>0.004004</td>\n",
       "      <td>0.018393</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>case89pegase</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.006109</td>\n",
       "      <td>0.010191</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>case118</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.005065</td>\n",
       "      <td>0.006915</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>case_ACTIVSg200</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.004998</td>\n",
       "      <td>0.006081</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>case145</td>\n",
       "      <td>0.000919</td>\n",
       "      <td>0.000290</td>\n",
       "      <td>0.006282</td>\n",
       "      <td>0.007770</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>case300</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000002</td>\n",
       "      <td>0.008165</td>\n",
       "      <td>0.013262</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>case_ACTIVSg500</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.007865</td>\n",
       "      <td>0.011466</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>case1354pegase</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000002</td>\n",
       "      <td>0.022147</td>\n",
       "      <td>0.025832</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>case1888rte</td>\n",
       "      <td>0.000001</td>\n",
       "      <td>0.000003</td>\n",
       "      <td>0.021216</td>\n",
       "      <td>0.019682</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>case1951rte</td>\n",
       "      <td>0.000001</td>\n",
       "      <td>0.000003</td>\n",
       "      <td>0.020124</td>\n",
       "      <td>0.025314</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>case2383wp</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000001</td>\n",
       "      <td>0.044913</td>\n",
       "      <td>0.059615</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>case3120sp</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.051568</td>\n",
       "      <td>0.075086</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>case3012wp</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.038621</td>\n",
       "      <td>0.042295</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>case2848rte</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000001</td>\n",
       "      <td>0.038467</td>\n",
       "      <td>0.049454</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>case2868rte</td>\n",
       "      <td>0.000001</td>\n",
       "      <td>0.000002</td>\n",
       "      <td>0.039894</td>\n",
       "      <td>0.054357</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>case2869pegase</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000002</td>\n",
       "      <td>0.064755</td>\n",
       "      <td>0.071706</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>case3375wp</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000003</td>\n",
       "      <td>0.031728</td>\n",
       "      <td>0.030560</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>case_ACTIVSg2000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000010</td>\n",
       "      <td>0.036846</td>\n",
       "      <td>0.035593</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30</th>\n",
       "      <td>case2737sop</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.043554</td>\n",
       "      <td>0.057491</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>31</th>\n",
       "      <td>case2736sp</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.051892</td>\n",
       "      <td>0.075590</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>32</th>\n",
       "      <td>case2746wop</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.041671</td>\n",
       "      <td>0.046079</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>33</th>\n",
       "      <td>case2746wp</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.042574</td>\n",
       "      <td>0.046418</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>34</th>\n",
       "      <td>case6468rte</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000002</td>\n",
       "      <td>0.114465</td>\n",
       "      <td>0.095043</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>35</th>\n",
       "      <td>case6470rte</td>\n",
       "      <td>0.000001</td>\n",
       "      <td>0.000002</td>\n",
       "      <td>0.101663</td>\n",
       "      <td>0.097198</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>36</th>\n",
       "      <td>case6495rte</td>\n",
       "      <td>0.000001</td>\n",
       "      <td>0.000002</td>\n",
       "      <td>0.092424</td>\n",
       "      <td>0.061160</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>37</th>\n",
       "      <td>case6515rte</td>\n",
       "      <td>0.000001</td>\n",
       "      <td>0.000003</td>\n",
       "      <td>0.271681</td>\n",
       "      <td>0.158288</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>38</th>\n",
       "      <td>case9241pegase</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000007</td>\n",
       "      <td>0.343820</td>\n",
       "      <td>0.260913</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>39</th>\n",
       "      <td>case8387pegase</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000005</td>\n",
       "      <td>0.202764</td>\n",
       "      <td>0.112890</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>40</th>\n",
       "      <td>case13659pegase</td>\n",
       "      <td>0.000001</td>\n",
       "      <td>0.000126</td>\n",
       "      <td>0.386642</td>\n",
       "      <td>0.344090</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>41</th>\n",
       "      <td>case_ACTIVSg10k</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000011</td>\n",
       "      <td>0.203224</td>\n",
       "      <td>0.154309</td>\n",
       "      <td>True</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                case  v_diff_max  a_diff_max      t_ad      t_pp   c_ad   c_pp\n",
       "0             case18    0.000000    0.000000  0.005228  0.008196   True   True\n",
       "1             case22   -1.000000   -1.000000  0.009135  0.016628   True  False\n",
       "2             case14    0.000000    0.000000  0.003891  0.005336   True   True\n",
       "3             case30    0.000000    0.000000  0.004930  0.007562   True   True\n",
       "4        case_ieee30    0.000000    0.000000  0.004656  0.006019   True   True\n",
       "5             case69   -1.000000   -1.000000  0.003248  0.020533  False  False\n",
       "6    case24_ieee_rts    0.000000    0.000000  0.005172  0.007628   True   True\n",
       "7             case39    0.000000    0.000000  0.003234  0.003471   True   True\n",
       "8             case85   -1.000000   -1.000000  0.003196  0.016355  False  False\n",
       "9          case118zh   -1.000000   -1.000000  0.003179  0.017203  False  False\n",
       "10            case57    0.000000    0.000000  0.004299  0.006711   True   True\n",
       "11           case141   -1.000000   -1.000000  0.003319  0.017596  False  False\n",
       "12         case136ma   -1.000000   -1.000000  0.004004  0.018393  False  False\n",
       "13      case89pegase    0.000000    0.000000  0.006109  0.010191   True   True\n",
       "14           case118    0.000000    0.000000  0.005065  0.006915   True   True\n",
       "15   case_ACTIVSg200    0.000000    0.000000  0.004998  0.006081   True   True\n",
       "16           case145    0.000919    0.000290  0.006282  0.007770   True   True\n",
       "17           case300    0.000000    0.000002  0.008165  0.013262   True   True\n",
       "18   case_ACTIVSg500    0.000000    0.000000  0.007865  0.011466   True   True\n",
       "19    case1354pegase    0.000000    0.000002  0.022147  0.025832   True   True\n",
       "20       case1888rte    0.000001    0.000003  0.021216  0.019682   True   True\n",
       "21       case1951rte    0.000001    0.000003  0.020124  0.025314   True   True\n",
       "22        case2383wp    0.000000    0.000001  0.044913  0.059615   True   True\n",
       "23        case3120sp    0.000000    0.000000  0.051568  0.075086   True   True\n",
       "24        case3012wp    0.000000    0.000000  0.038621  0.042295   True   True\n",
       "25       case2848rte    0.000000    0.000001  0.038467  0.049454   True   True\n",
       "26       case2868rte    0.000001    0.000002  0.039894  0.054357   True   True\n",
       "27    case2869pegase    0.000000    0.000002  0.064755  0.071706   True   True\n",
       "28        case3375wp    0.000000    0.000003  0.031728  0.030560   True   True\n",
       "29  case_ACTIVSg2000    0.000000    0.000010  0.036846  0.035593   True   True\n",
       "30       case2737sop    0.000000    0.000000  0.043554  0.057491   True   True\n",
       "31        case2736sp    0.000000    0.000000  0.051892  0.075590   True   True\n",
       "32       case2746wop    0.000000    0.000000  0.041671  0.046079   True   True\n",
       "33        case2746wp    0.000000    0.000000  0.042574  0.046418   True   True\n",
       "34       case6468rte    0.000000    0.000002  0.114465  0.095043   True   True\n",
       "35       case6470rte    0.000001    0.000002  0.101663  0.097198   True   True\n",
       "36       case6495rte    0.000001    0.000002  0.092424  0.061160   True   True\n",
       "37       case6515rte    0.000001    0.000003  0.271681  0.158288   True   True\n",
       "38    case9241pegase    0.000000    0.000007  0.343820  0.260913   True   True\n",
       "39    case8387pegase    0.000000    0.000005  0.202764  0.112890   True   True\n",
       "40   case13659pegase    0.000001    0.000126  0.386642  0.344090   True   True\n",
       "41   case_ACTIVSg10k    0.000000    0.000011  0.203224  0.154309   True   True"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "res[['c_ad', 'c_pp']] = res[['c_ad', 'c_pp']].astype(bool)\n",
    "res.round(6)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "ams",
   "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.9.18"
  },
  "orig_nbformat": 4,
  "vscode": {
   "interpreter": {
    "hash": "d2b3bf80176349caa68dc4a3c77bd06eaade8abc678330f7d1c813c53380e5d2"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}