Skip to content

Commit

Permalink
experiments
Browse files Browse the repository at this point in the history
  • Loading branch information
@Yuriyzabegaev
Yuriyzabegaev committed Dec 26, 2023
1 parent 12108cb commit b48234c
Show file tree
Hide file tree
Showing 6 changed files with 607 additions and 222 deletions.
395 changes: 384 additions & 11 deletions examples/1/results_1_poro.ipynb
View file

Large diffs are not rendered by default.

31 changes: 15 additions & 16 deletions examples/1/results_1_thermal.ipynb
View file

Large diffs are not rendered by default.

39 changes: 21 additions & 18 deletions examples/2_source_location/results_2_thermal.ipynb
View file

Large diffs are not rendered by default.

143 changes: 73 additions & 70 deletions examples/3_exploration_sensitivity/results_3_2.ipynb
View file
Original file line number Diff line number Diff line change
Expand Up @@ -148,19 +148,19 @@
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>Standard exploration</th>\n",
" <th>Standard exploration $\\varepsilon = 0.5$</th>\n",
" <td>71.675894</td>\n",
" <td>69.395625</td>\n",
" <td>67.614741</td>\n",
" </tr>\n",
" <tr>\n",
" <th>No exploration</th>\n",
" <th>No exploration $\\varepsilon = 0$</th>\n",
" <td>102.419440</td>\n",
" <td>70.188748</td>\n",
" <td>65.645999</td>\n",
" </tr>\n",
" <tr>\n",
" <th>More exploration</th>\n",
" <th>More exploration $\\varepsilon = 0.7$</th>\n",
" <td>71.297287</td>\n",
" <td>69.732346</td>\n",
" <td>68.168296</td>\n",
Expand All @@ -176,11 +176,11 @@
"</div>"
],
"text/plain": [
" Worst Mean Best\n",
"Standard exploration 71.675894 69.395625 67.614741\n",
"No exploration 102.419440 70.188748 65.645999\n",
"More exploration 71.297287 69.732346 68.168296\n",
"Random choice 88.472429 85.484517 83.025788"
" Worst Mean Best\n",
"Standard exploration $\\varepsilon = 0.5$ 71.675894 69.395625 67.614741\n",
"No exploration $\\varepsilon = 0$ 102.419440 70.188748 65.645999\n",
"More exploration $\\varepsilon = 0.7$ 71.297287 69.732346 68.168296\n",
"Random choice 88.472429 85.484517 83.025788"
]
},
"execution_count": 2,
Expand All @@ -191,9 +191,9 @@
"source": [
"data = {\n",
" # \"Standard exploration\": load_data('../1/poro_dynamic', n_newest=NUM_SIMS),\n",
" \"Standard exploration\": load_data(\"poro_eps_default\", n_newest=NUM_SIMS),\n",
" \"No exploration\": load_data(\"poro_eps0\", n_newest=NUM_SIMS),\n",
" \"More exploration\": load_data(\"poro_eps_big\", n_newest=NUM_SIMS),\n",
" r\"Standard exploration $\\varepsilon = 0.5$\": load_data(\"poro_eps_default\", n_newest=NUM_SIMS),\n",
" r\"No exploration $\\varepsilon = 0$\": load_data(\"poro_eps0\", n_newest=NUM_SIMS),\n",
" r\"More exploration $\\varepsilon = 0.7$\": load_data(\"poro_eps_big\", n_newest=NUM_SIMS),\n",
" \"Random choice\": load_data('../1/poro_lrandom', n_newest=NUM_SIMS)\n",
" # \"GP\": load_data(\"poro_gp\", n_newest=NUM_SIMS),\n",
"}\n",
Expand Down Expand Up @@ -372,19 +372,19 @@
" <td>109.033182</td>\n",
" </tr>\n",
" <tr>\n",
" <th>Standard exploration</th>\n",
" <th>Standard exploration $\\varepsilon = 0.5$</th>\n",
" <td>122.746576</td>\n",
" <td>114.867089</td>\n",
" <td>107.933468</td>\n",
" </tr>\n",
" <tr>\n",
" <th>No exploration</th>\n",
" <th>No exploration $\\varepsilon = 0$</th>\n",
" <td>134.638330</td>\n",
" <td>114.571609</td>\n",
" <td>107.830565</td>\n",
" </tr>\n",
" <tr>\n",
" <th>More exploration</th>\n",
" <th>More exploration $\\varepsilon = 0.7$</th>\n",
" <td>135.631345</td>\n",
" <td>117.405578</td>\n",
" <td>108.233092</td>\n",
Expand All @@ -400,12 +400,12 @@
"</div>"
],
"text/plain": [
" Worst Mean Best\n",
"Standard from 1 124.957382 116.552458 109.033182\n",
"Standard exploration 122.746576 114.867089 107.933468\n",
"No exploration 134.638330 114.571609 107.830565\n",
"More exploration 135.631345 117.405578 108.233092\n",
"Random choice 202.020621 182.638595 143.524424"
" Worst Mean Best\n",
"Standard from 1 124.957382 116.552458 109.033182\n",
"Standard exploration $\\varepsilon = 0.5$ 122.746576 114.867089 107.933468\n",
"No exploration $\\varepsilon = 0$ 134.638330 114.571609 107.830565\n",
"More exploration $\\varepsilon = 0.7$ 135.631345 117.405578 108.233092\n",
"Random choice 202.020621 182.638595 143.524424"
]
},
"execution_count": 3,
Expand All @@ -416,9 +416,9 @@
"source": [
"data = {\n",
" \"Standard from 1\": load_data('../1/thermal_dynamic', n_newest=NUM_SIMS),\n",
" \"Standard exploration\": load_data(\"thermal_eps_default\", n_newest=NUM_SIMS),\n",
" \"No exploration\": load_data(\"thermal_eps0\", n_newest=NUM_SIMS),\n",
" \"More exploration\": load_data(\"thermal_eps_big\", n_newest=NUM_SIMS),\n",
" r\"Standard exploration $\\varepsilon = 0.5$\": load_data(\"thermal_eps_default\", n_newest=NUM_SIMS),\n",
" r\"No exploration $\\varepsilon = 0$\": load_data(\"thermal_eps0\", n_newest=NUM_SIMS),\n",
" r\"More exploration $\\varepsilon = 0.7$\": load_data(\"thermal_eps_big\", n_newest=NUM_SIMS),\n",
" \"Random choice\": load_data('../1/thermal_random', n_newest=NUM_SIMS)\n",
" # \"GP\": load_data(\"thermal_gp\", n_newest=NUM_SIMS),\n",
"}\n",
Expand Down Expand Up @@ -456,69 +456,69 @@
"data": {
"text/html": [
"<style type=\"text/css\">\n",
"#T_cac7b_row0_col1, #T_cac7b_row0_col3, #T_cac7b_row1_col2, #T_cac7b_row1_col4, #T_cac7b_row1_col5, #T_cac7b_row2_col0 {\n",
"#T_fc087_row0_col1, #T_fc087_row0_col3, #T_fc087_row1_col2, #T_fc087_row1_col4, #T_fc087_row1_col5, #T_fc087_row2_col0 {\n",
" textbf: --rwrap;\n",
"}\n",
"</style>\n",
"<table id=\"T_cac7b\">\n",
"<table id=\"T_fc087\">\n",
" <thead>\n",
" <tr>\n",
" <th class=\"blank level0\" >&nbsp;</th>\n",
" <th id=\"T_cac7b_level0_col0\" class=\"col_heading level0 col0\" colspan=\"3\">Mandel problem</th>\n",
" <th id=\"T_cac7b_level0_col3\" class=\"col_heading level0 col3\" colspan=\"3\">Thermal problem</th>\n",
" <th id=\"T_fc087_level0_col0\" class=\"col_heading level0 col0\" colspan=\"3\">Mandel's problem</th>\n",
" <th id=\"T_fc087_level0_col3\" class=\"col_heading level0 col3\" colspan=\"3\">Non-isothermal flow problem</th>\n",
" </tr>\n",
" <tr>\n",
" <th class=\"blank level1\" >&nbsp;</th>\n",
" <th id=\"T_cac7b_level1_col0\" class=\"col_heading level1 col0\" >Worst</th>\n",
" <th id=\"T_cac7b_level1_col1\" class=\"col_heading level1 col1\" >Mean</th>\n",
" <th id=\"T_cac7b_level1_col2\" class=\"col_heading level1 col2\" >Best</th>\n",
" <th id=\"T_cac7b_level1_col3\" class=\"col_heading level1 col3\" >Worst</th>\n",
" <th id=\"T_cac7b_level1_col4\" class=\"col_heading level1 col4\" >Mean</th>\n",
" <th id=\"T_cac7b_level1_col5\" class=\"col_heading level1 col5\" >Best</th>\n",
" <th id=\"T_fc087_level1_col0\" class=\"col_heading level1 col0\" >Worst</th>\n",
" <th id=\"T_fc087_level1_col1\" class=\"col_heading level1 col1\" >Mean</th>\n",
" <th id=\"T_fc087_level1_col2\" class=\"col_heading level1 col2\" >Best</th>\n",
" <th id=\"T_fc087_level1_col3\" class=\"col_heading level1 col3\" >Worst</th>\n",
" <th id=\"T_fc087_level1_col4\" class=\"col_heading level1 col4\" >Mean</th>\n",
" <th id=\"T_fc087_level1_col5\" class=\"col_heading level1 col5\" >Best</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th id=\"T_cac7b_level0_row0\" class=\"row_heading level0 row0\" >Standard exploration</th>\n",
" <td id=\"T_cac7b_row0_col0\" class=\"data row0 col0\" >71.68</td>\n",
" <td id=\"T_cac7b_row0_col1\" class=\"data row0 col1\" >69.40</td>\n",
" <td id=\"T_cac7b_row0_col2\" class=\"data row0 col2\" >67.61</td>\n",
" <td id=\"T_cac7b_row0_col3\" class=\"data row0 col3\" >122.75</td>\n",
" <td id=\"T_cac7b_row0_col4\" class=\"data row0 col4\" >114.87</td>\n",
" <td id=\"T_cac7b_row0_col5\" class=\"data row0 col5\" >107.93</td>\n",
" <th id=\"T_fc087_level0_row0\" class=\"row_heading level0 row0\" >Standard exploration $\\varepsilon = 0.5$</th>\n",
" <td id=\"T_fc087_row0_col0\" class=\"data row0 col0\" >71.68</td>\n",
" <td id=\"T_fc087_row0_col1\" class=\"data row0 col1\" >69.40</td>\n",
" <td id=\"T_fc087_row0_col2\" class=\"data row0 col2\" >67.61</td>\n",
" <td id=\"T_fc087_row0_col3\" class=\"data row0 col3\" >122.75</td>\n",
" <td id=\"T_fc087_row0_col4\" class=\"data row0 col4\" >114.87</td>\n",
" <td id=\"T_fc087_row0_col5\" class=\"data row0 col5\" >107.93</td>\n",
" </tr>\n",
" <tr>\n",
" <th id=\"T_cac7b_level0_row1\" class=\"row_heading level0 row1\" >No exploration</th>\n",
" <td id=\"T_cac7b_row1_col0\" class=\"data row1 col0\" >102.42</td>\n",
" <td id=\"T_cac7b_row1_col1\" class=\"data row1 col1\" >70.19</td>\n",
" <td id=\"T_cac7b_row1_col2\" class=\"data row1 col2\" >65.65</td>\n",
" <td id=\"T_cac7b_row1_col3\" class=\"data row1 col3\" >134.64</td>\n",
" <td id=\"T_cac7b_row1_col4\" class=\"data row1 col4\" >114.57</td>\n",
" <td id=\"T_cac7b_row1_col5\" class=\"data row1 col5\" >107.83</td>\n",
" <th id=\"T_fc087_level0_row1\" class=\"row_heading level0 row1\" >No exploration $\\varepsilon = 0$</th>\n",
" <td id=\"T_fc087_row1_col0\" class=\"data row1 col0\" >102.42</td>\n",
" <td id=\"T_fc087_row1_col1\" class=\"data row1 col1\" >70.19</td>\n",
" <td id=\"T_fc087_row1_col2\" class=\"data row1 col2\" >65.65</td>\n",
" <td id=\"T_fc087_row1_col3\" class=\"data row1 col3\" >134.64</td>\n",
" <td id=\"T_fc087_row1_col4\" class=\"data row1 col4\" >114.57</td>\n",
" <td id=\"T_fc087_row1_col5\" class=\"data row1 col5\" >107.83</td>\n",
" </tr>\n",
" <tr>\n",
" <th id=\"T_cac7b_level0_row2\" class=\"row_heading level0 row2\" >More exploration</th>\n",
" <td id=\"T_cac7b_row2_col0\" class=\"data row2 col0\" >71.30</td>\n",
" <td id=\"T_cac7b_row2_col1\" class=\"data row2 col1\" >69.73</td>\n",
" <td id=\"T_cac7b_row2_col2\" class=\"data row2 col2\" >68.17</td>\n",
" <td id=\"T_cac7b_row2_col3\" class=\"data row2 col3\" >135.63</td>\n",
" <td id=\"T_cac7b_row2_col4\" class=\"data row2 col4\" >117.41</td>\n",
" <td id=\"T_cac7b_row2_col5\" class=\"data row2 col5\" >108.23</td>\n",
" <th id=\"T_fc087_level0_row2\" class=\"row_heading level0 row2\" >More exploration $\\varepsilon = 0.7$</th>\n",
" <td id=\"T_fc087_row2_col0\" class=\"data row2 col0\" >71.30</td>\n",
" <td id=\"T_fc087_row2_col1\" class=\"data row2 col1\" >69.73</td>\n",
" <td id=\"T_fc087_row2_col2\" class=\"data row2 col2\" >68.17</td>\n",
" <td id=\"T_fc087_row2_col3\" class=\"data row2 col3\" >135.63</td>\n",
" <td id=\"T_fc087_row2_col4\" class=\"data row2 col4\" >117.41</td>\n",
" <td id=\"T_fc087_row2_col5\" class=\"data row2 col5\" >108.23</td>\n",
" </tr>\n",
" <tr>\n",
" <th id=\"T_cac7b_level0_row3\" class=\"row_heading level0 row3\" >Random choice</th>\n",
" <td id=\"T_cac7b_row3_col0\" class=\"data row3 col0\" >88.47</td>\n",
" <td id=\"T_cac7b_row3_col1\" class=\"data row3 col1\" >85.48</td>\n",
" <td id=\"T_cac7b_row3_col2\" class=\"data row3 col2\" >83.03</td>\n",
" <td id=\"T_cac7b_row3_col3\" class=\"data row3 col3\" >202.02</td>\n",
" <td id=\"T_cac7b_row3_col4\" class=\"data row3 col4\" >182.64</td>\n",
" <td id=\"T_cac7b_row3_col5\" class=\"data row3 col5\" >143.52</td>\n",
" <th id=\"T_fc087_level0_row3\" class=\"row_heading level0 row3\" >Random choice</th>\n",
" <td id=\"T_fc087_row3_col0\" class=\"data row3 col0\" >88.47</td>\n",
" <td id=\"T_fc087_row3_col1\" class=\"data row3 col1\" >85.48</td>\n",
" <td id=\"T_fc087_row3_col2\" class=\"data row3 col2\" >83.03</td>\n",
" <td id=\"T_fc087_row3_col3\" class=\"data row3 col3\" >202.02</td>\n",
" <td id=\"T_fc087_row3_col4\" class=\"data row3 col4\" >182.64</td>\n",
" <td id=\"T_fc087_row3_col5\" class=\"data row3 col5\" >143.52</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n"
],
"text/plain": [
"<pandas.io.formats.style.Styler at 0x7f4cb3517510>"
"<pandas.io.formats.style.Styler at 0x7f3f211dad50>"
]
},
"execution_count": 4,
Expand All @@ -527,7 +527,7 @@
}
],
"source": [
"df = pd.concat((df_mandel, df_thermal), join='inner', axis=1, keys=(\"Mandel problem\", \"Thermal problem\"))\n",
"df = pd.concat((df_mandel, df_thermal), join='inner', axis=1, keys=(\"Mandel's problem\", \"Non-isothermal flow problem\"))\n",
"df = df.style.highlight_min(axis=0, props=\"textbf:--rwrap;\").format(precision=2)\n",
"df"
]
Expand All @@ -547,20 +547,23 @@
"name": "stdout",
"output_type": "stream",
"text": [
"\\begin{tabular}{lrrrrrr}\n",
" & \\multicolumn{3}{r}{Mandel problem} & \\multicolumn{3}{r}{Thermal problem} \\\\\n",
"\\begin{tabular}{l | rrr | rrr}\n",
"\\toprule\n",
" & \\multicolumn{3}{r}{Mandel's problem} & \\multicolumn{3}{r}{Non-isothermal flow problem} \\\\\n",
" & Worst & Mean & Best & Worst & Mean & Best \\\\\n",
"Standard exploration & 71.68 & \\textbf{69.40} & 67.61 & \\textbf{122.75} & 114.87 & 107.93 \\\\\n",
"No exploration & 102.42 & 70.19 & \\textbf{65.65} & 134.64 & \\textbf{114.57} & \\textbf{107.83} \\\\\n",
"More exploration & \\textbf{71.30} & 69.73 & 68.17 & 135.63 & 117.41 & 108.23 \\\\\n",
"\\midrule\n",
"Standard exploration $\\varepsilon = 0.5$ & 71.68 & \\textbf{69.40} & 67.61 & \\textbf{122.75} & 114.87 & 107.93 \\\\\n",
"No exploration $\\varepsilon = 0$ & 102.42 & 70.19 & \\textbf{65.65} & 134.64 & \\textbf{114.57} & \\textbf{107.83} \\\\\n",
"More exploration $\\varepsilon = 0.7$ & \\textbf{71.30} & 69.73 & 68.17 & 135.63 & 117.41 & 108.23 \\\\\n",
"Random choice & 88.47 & 85.48 & 83.03 & 202.02 & 182.64 & 143.52 \\\\\n",
"\\bottomrule\n",
"\\end{tabular}\n",
"\n"
]
}
],
"source": [
"print(df.to_latex())"
"print(df.to_latex(column_format=' l | rrr | rrr ', hrules=True))"
]
}
],
Expand All @@ -580,7 +583,7 @@
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.11.7"
"version": "3.11.2"
}
},
"nbformat": 4,
Expand Down
157 changes: 82 additions & 75 deletions examples/4_extended_solver_space/results_4_thermal.ipynb
View file

Large diffs are not rendered by default.

64 changes: 32 additions & 32 deletions examples/spe10_data/spe10.ipynb
View file

Large diffs are not rendered by default.

0 comments on commit b48234c

Please sign in to comment.