Make plots show per-round-LER; Add windowed soft info decoding vs whole decoding comparison

src/param-exploration/main.tex

@@ -11,6 +11,7 @@
 \usepackage{subcaption}
 \usepackage{xcolor}
 \usepackage{pgfplots}
+\usepackage{pgfplotstable}
 \pgfplotsset{compat=newest}
 \usepackage{acro}
 \usepackage{braket}
@@ -67,7 +68,6 @@
 \section{Whole vs Windowed decoding}
 \label{sec:Whole vs Windowed decoding}
 
-
 \begin{figure}[H]
     \centering
 
@@ -82,7 +82,7 @@
                         cells={align=left},
                     },
                     enlargelimits=false,
-                    ymin=1e-2, ymax=1,
+                    ymin=1e-3, ymax=2.5e-1,
                     grid=both,
                     legend pos = south east,
                     xtick={0.001,0.0015,...,0.004},
@@ -90,12 +90,17 @@
                     xticklabel style={/pgf/number format/precision=4},
                     scaled x ticks=false,
                     xlabel={Physical error rate},
-                    ylabel={LER},
+                    ylabel={Per-round-LER},
                 ]
                 \foreach \W/\col in {3/scol2,4/scol1,5/scol0} {
                     \edef\temp{\noexpand
-                        \addplot+[mark=o, line width=1pt, densely dashed, \col]
+                        \addplot+[
-                        table[col sep=comma, x=physical_p, y=LER_whole]
+                            mark=o, line width=1pt, densely dashed, \col,
+                        ]
+                        table[
+                            col sep=comma, x=physical_p,
+                            y expr={1 - (1-\noexpand\thisrow{LER_whole})^(1/12)}
+                        ]
                         {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=minimum_sum,+experiment=windowed_vs_whole_param_exploration,+system.F=1,+system.W=\W/2026-03-29_23-39-18/LERs.csv};
                     }
                     \temp
@@ -106,7 +111,10 @@
                 \foreach \W/\col in {3/scol2,4/scol1,5/scol0} {
                     \edef\temp{\noexpand
                         \addplot+[mark=o, line width=1pt, \col]
-                        table[col sep=comma, x=physical_p, y=LER_sw]
+                        table[
+                            col sep=comma, x=physical_p,
+                            y expr={1 - (1-\noexpand\thisrow{LER_sw})^(1/12)}
+                        ]
                         {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=minimum_sum,+experiment=windowed_vs_whole_param_exploration,+system.F=1,+system.W=\W/2026-03-29_23-39-18/LERs.csv};
                     }
                     \temp
@@ -128,7 +136,7 @@
                         cells={align=left},
                     },
                     enlargelimits=false,
-                    ymin=1e-2, ymax=1,
+                    ymin=1e-3, ymax=2.5e-1,
                     grid=both,
                     legend pos = south east,
                     xtick={0.001,0.0015,...,0.004},
@@ -136,12 +144,15 @@
                     xticklabel style={/pgf/number format/precision=4},
                     scaled x ticks=false,
                     xlabel={Physical error rate},
-                    ylabel={LER},
+                    ylabel={Per-round-LER},
                 ]
                 \foreach \W/\col in {3/scol2,4/scol1,5/scol0} {
                     \edef\temp{\noexpand
                         \addplot+[mark=o, line width=1pt, densely dashed, \col]
-                        table[col sep=comma, x=physical_p, y=LER_whole]
+                        table[
+                            col sep=comma, x=physical_p,
+                            y expr={1 - (1-\noexpand\thisrow{LER_whole})^(1/12)}
+                        ]
                         {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=product_sum,+experiment=windowed_vs_whole_param_exploration,+system.F=1,+system.W=\W/2026-03-29_23-39-18/LERs.csv};
                     }
                     \temp
@@ -152,7 +163,10 @@
                 \foreach \W/\col in {3/scol2,4/scol1,5/scol0} {
                     \edef\temp{\noexpand
                         \addplot+[mark=o, line width=1pt, \col]
-                        table[col sep=comma, x=physical_p, y=LER_sw]
+                        table[
+                            col sep=comma, x=physical_p,
+                            y expr={1 - (1-\noexpand\thisrow{LER_sw})^(1/12)}
+                        ]
                         {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=product_sum,+experiment=windowed_vs_whole_param_exploration,+system.F=1,+system.W=\W/2026-03-29_23-39-18/LERs.csv};
                     }
                     \temp
@@ -165,7 +179,7 @@
     \end{subfigure}
 
     \caption{BP simulations for the [[144,12,12]]-BB code with $F=1,
-    n_\text{iter,BP}=32$}
+    n_\text{iter,BP}=32,n_\text{rounds}=12$}
 \end{figure}
 
 \begin{figure}[H]
@@ -182,7 +196,7 @@
                         cells={align=left},
                     },
                     enlargelimits=false,
-                    ymin=1e-2, ymax=1,
+                    ymin=1e-3, ymax=2.5e-1,
                     grid=both,
                     legend pos = south east,
                     xtick={0.001,0.0015,...,0.004},
@@ -190,12 +204,15 @@
                     xticklabel style={/pgf/number format/precision=4},
                     scaled x ticks=false,
                     xlabel={Physical error rate},
-                    ylabel={LER},
+                    ylabel={Per-round-LER},
                 ]
                 \foreach \F/\col in {1/scol2,2/scol1,3/scol0} {
                     \edef\temp{\noexpand
                         \addplot+[mark=o, line width=1pt, densely dashed, \col]
-                        table[col sep=comma, x=physical_p, y=LER_whole]
+                        table[
+                            col sep=comma, x=physical_p,
+                            y expr={1 - (1-\noexpand\thisrow{LER_whole})^(1/12)}
+                        ]
                         {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=minimum_sum,+experiment=windowed_vs_whole_param_exploration,+system.F=\F,+system.W=5/2026-03-29_23-39-18/LERs.csv};
                     }
                     \temp
@@ -206,8 +223,10 @@
                 \foreach \F/\col in {1/scol2,2/scol1,3/scol0} {
                     \edef\temp{\noexpand
                         \addplot+[mark=o, line width=1pt, \col]
-                        table[col sep=comma, x=physical_p, y=LER_sw]
+                        table[
-                        % {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=minimum_sum,+experiment=windowed_vs_whole_param_exploration,+system.F=\F,+system.W=5/2026-03-29_23-39-18/LERs.csv};
+                            col sep=comma, x=physical_p,
+                            y expr={1 - (1-\noexpand\thisrow{LER_sw})^(1/12)}
+                        ]
                         {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=minimum_sum,+experiment=windowed_vs_whole_param_exploration,+system.F=\F,+system.W=5/2026-03-29_23-39-18/LERs.csv};
                     }
                     \temp
@@ -229,7 +248,7 @@
                         cells={align=left},
                     },
                     enlargelimits=false,
-                    ymin=1e-2, ymax=1,
+                    ymin=1e-3, ymax=2.5e-1,
                     grid=both,
                     legend pos = south east,
                     xtick={0.001,0.0015,...,0.004},
@@ -237,12 +256,15 @@
                     xticklabel style={/pgf/number format/precision=4},
                     scaled x ticks=false,
                     xlabel={Physical error rate},
-                    ylabel={LER},
+                    ylabel={Per-round-LER},
                 ]
                 \foreach \F/\col in {1/scol2,2/scol1,3/scol0} {
                     \edef\temp{\noexpand
                         \addplot+[mark=o, line width=1pt, densely dashed, \col]
-                        table[col sep=comma, x=physical_p, y=LER_whole]
+                        table[
+                            col sep=comma, x=physical_p,
+                            y expr={1 - (1-\noexpand\thisrow{LER_whole})^(1/12)}
+                        ]
                         {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=product_sum,+experiment=windowed_vs_whole_param_exploration,+system.F=\F,+system.W=5/2026-03-29_23-39-18/LERs.csv};
                     }
                     \temp
@@ -253,7 +275,10 @@
                 \foreach \F/\col in {1/scol2,2/scol1,3/scol0} {
                     \edef\temp{\noexpand
                         \addplot+[mark=o, line width=1pt, \col]
-                        table[col sep=comma, x=physical_p, y=LER_sw]
+                        table[
+                            col sep=comma, x=physical_p,
+                            y expr={1 - (1-\noexpand\thisrow{LER_sw})^(1/12)}
+                        ]
                         {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=product_sum,+experiment=windowed_vs_whole_param_exploration,+system.F=\F,+system.W=5/2026-03-29_23-39-18/LERs.csv};
                     }
                     \temp
@@ -266,7 +291,7 @@
     \end{subfigure}
 
     \caption{BP simulations for the [[144,12,12]]-BB code with $W=5,
-    n_\text{iter,BP}=32$}
+    n_\text{iter,BP}=32,n_\text{rounds}=12$}
 \end{figure}
 
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
@@ -287,7 +312,7 @@
                         cells={align=left},
                     },
                     enlargelimits=false,
-                    ymin=1e-2, ymax=1,
+                    ymin=8e-4, ymax=2.5e-1,
                     grid=both,
                     legend pos = south east,
                     xtick={0.001,0.0015,...,0.004},
@@ -295,12 +320,15 @@
                     xticklabel style={/pgf/number format/precision=4},
                     scaled x ticks=false,
                     xlabel={Physical error rate},
-                    ylabel={LER},
+                    ylabel={Per-round-LER},
                 ]
                 \foreach \W/\col in {3/scol2,4/scol1,5/scol0} {
                     \edef\temp{\noexpand
                         \addplot+[mark=o, line width=1pt, densely dashed, \col]
-                        table[col sep=comma, x=physical_p, y=LER_hard]
+                        table[
+                            col sep=comma, x=physical_p,
+                            y expr={1 - (1-\noexpand\thisrow{LER_hard})^(1/12)}
+                        ]
                         {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=minimum_sum,+experiment=soft_vs_hard_param_exploration,+system.F=1,+system.W=\W/2026-03-30_00-06-26/LERs.csv};
                     }
                     \temp
@@ -311,7 +339,10 @@
                 \foreach \W/\col in {3/scol2,4/scol1,5/scol0} {
                     \edef\temp{\noexpand
                         \addplot+[mark=*, line width=1pt, \col]
-                        table[col sep=comma, x=physical_p, y=LER_soft]
+                        table[
+                            col sep=comma, x=physical_p,
+                            y expr={1 - (1-\noexpand\thisrow{LER_soft})^(1/12)}
+                        ]
                         {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=minimum_sum,+experiment=soft_vs_hard_param_exploration,+system.F=1,+system.W=\W/2026-03-30_00-06-26/LERs.csv};
                     }
                     \temp
@@ -333,7 +364,7 @@
                         cells={align=left},
                     },
                     enlargelimits=false,
-                    ymin=1e-2, ymax=1,
+                    ymin=8e-4, ymax=2.5e-1,
                     grid=both,
                     legend pos = south east,
                     xtick={0.001,0.0015,...,0.004},
@@ -341,12 +372,15 @@
                     xticklabel style={/pgf/number format/precision=4},
                     scaled x ticks=false,
                     xlabel={Physical error rate},
-                    ylabel={LER},
+                    ylabel={Per-round-LER},
                 ]
                 \foreach \W/\col in {3/scol2,4/scol1,5/scol0} {
                     \edef\temp{\noexpand
                         \addplot+[mark=o, line width=1pt, densely dashed, \col]
-                        table[col sep=comma, x=physical_p, y=LER_hard]
+                        table[
+                            col sep=comma, x=physical_p,
+                            y expr={1 - (1-\noexpand\thisrow{LER_hard})^(1/12)}
+                        ]
                         {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=product_sum,+experiment=soft_vs_hard_param_exploration,+system.F=1,+system.W=\W/2026-03-30_00-06-26/LERs.csv};
                     }
                     \temp
@@ -357,7 +391,10 @@
                 \foreach \W/\col in {3/scol2,4/scol1,5/scol0} {
                     \edef\temp{\noexpand
                         \addplot+[mark=*, line width=1pt, \col]
-                        table[col sep=comma, x=physical_p, y=LER_soft]
+                        table[
+                            col sep=comma, x=physical_p,
+                            y expr={1 - (1-\noexpand\thisrow{LER_soft})^(1/12)}
+                        ]
                         {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=product_sum,+experiment=soft_vs_hard_param_exploration,+system.F=1,+system.W=\W/2026-03-30_00-06-26/LERs.csv};
                     }
                     \temp
@@ -370,7 +407,7 @@
     \end{subfigure}
 
     \caption{BP simulations for the [[144,12,12]]-BB code with $F=1,
-    n_\text{iter,BP}=32$}
+    n_\text{iter,BP}=32,n_\text{rounds}=12$}
 \end{figure}
 
 \begin{figure}[H]
@@ -387,7 +424,7 @@
                         cells={align=left},
                     },
                     enlargelimits=false,
-                    ymin=1e-2, ymax=1,
+                    ymin=8e-4, ymax=2.5e-1,
                     grid=both,
                     legend pos = south east,
                     xtick={0.001,0.0015,...,0.004},
@@ -395,12 +432,15 @@
                     xticklabel style={/pgf/number format/precision=4},
                     scaled x ticks=false,
                     xlabel={Physical error rate},
-                    ylabel={LER},
+                    ylabel={Per-round-LER},
                 ]
                 \foreach \F/\col in {1/scol2,2/scol1,3/scol0} {
                     \edef\temp{\noexpand
                         \addplot+[mark=o, line width=1pt, densely dashed, \col]
-                        table[col sep=comma, x=physical_p, y=LER_hard]
+                        table[
+                            col sep=comma, x=physical_p,
+                            y expr={1 - (1-\noexpand\thisrow{LER_hard})^(1/12)}
+                        ]
                         {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=minimum_sum,+experiment=soft_vs_hard_param_exploration,+system.F=\F,+system.W=5/2026-03-30_00-06-26/LERs.csv};
                     }
                     \temp
@@ -411,7 +451,10 @@
                 \foreach \F/\col in {1/scol2,2/scol1,3/scol0} {
                     \edef\temp{\noexpand
                         \addplot+[mark=*, line width=1pt, \col]
-                        table[col sep=comma, x=physical_p, y=LER_soft]
+                        table[
+                            col sep=comma, x=physical_p,
+                            y expr={1 - (1-\noexpand\thisrow{LER_soft})^(1/12)}
+                        ]
                         {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=minimum_sum,+experiment=soft_vs_hard_param_exploration,+system.F=\F,+system.W=5/2026-03-30_00-06-26/LERs.csv};
                     }
                     \temp
@@ -433,7 +476,7 @@
                         cells={align=left},
                     },
                     enlargelimits=false,
-                    ymin=1e-2, ymax=1,
+                    ymin=8e-4, ymax=2.5e-1,
                     grid=both,
                     legend pos = south east,
                     xtick={0.001,0.0015,...,0.004},
@@ -441,12 +484,15 @@
                     xticklabel style={/pgf/number format/precision=4},
                     scaled x ticks=false,
                     xlabel={Physical error rate},
-                    ylabel={LER},
+                    ylabel={Per-round-LER},
                 ]
                 \foreach \F/\col in {1/scol2,2/scol1,3/scol0} {
                     \edef\temp{\noexpand
                         \addplot+[mark=o, line width=1pt, densely dashed, \col]
-                        table[col sep=comma, x=physical_p, y=LER_hard]
+                        table[
+                            col sep=comma, x=physical_p,
+                            y expr={1 - (1-\noexpand\thisrow{LER_hard})^(1/12)}
+                        ]
                         {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=product_sum,+experiment=soft_vs_hard_param_exploration,+system.F=\F,+system.W=5/2026-03-30_00-06-26/LERs.csv};
                     }
                     \temp
@@ -457,7 +503,10 @@
                 \foreach \F/\col in {1/scol2,2/scol1,3/scol0} {
                     \edef\temp{\noexpand
                         \addplot+[mark=*, line width=1pt, \col]
-                        table[col sep=comma, x=physical_p, y=LER_soft]
+                        table[
+                            col sep=comma, x=physical_p,
+                            y expr={1 - (1-\noexpand\thisrow{LER_soft})^(1/12)}
+                        ]
                         {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=product_sum,+experiment=soft_vs_hard_param_exploration,+system.F=\F,+system.W=5/2026-03-30_00-06-26/LERs.csv};
                     }
                     \temp
@@ -470,7 +519,121 @@
     \end{subfigure}
 
     \caption{BP simulations for the [[144,12,12]]-BB code with $W=5,
-    n_\text{iter,BP}=32$}
+    n_\text{iter,BP}=32,n_\text{rounds}=12$}
+\end{figure}
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+\section{Soft-Information Decoding: Best Performance}
+\label{sec:Soft-Information Decoding: Best Performance}
+
+\begin{figure}[H]
+    \centering
+
+    \begin{subfigure}{0.5\textwidth}
+        \begin{tikzpicture}
+            \begin{axis}[
+                    width=\figwidth,
+                    height=\figheight,
+                    ymode=log,
+                    legend style={
+                        cells={anchor=west},
+                        cells={align=left},
+                    },
+                    enlargelimits=false,
+                    ymin=8e-4, ymax=2.5e-1,
+                    grid=both,
+                    legend pos = south east,
+                    xtick={0.001,0.0015,...,0.004},
+                    xticklabel style={/pgf/number format/fixed},
+                    xticklabel style={/pgf/number format/precision=4},
+                    scaled x ticks=false,
+                    xlabel={Physical error rate},
+                    ylabel={Per-round-LER},
+                ]
+                \addplot+[mark=o, line width=1pt, densely dashed, scol0]
+                table[
+                    col sep=comma, x=physical_p,
+                    y expr={1 - (1-\thisrow{LER_hard})^(1/12)}
+                ]
+                {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=minimum_sum,+experiment=soft_vs_hard_param_exploration,+system.F=1,+system.W=5/2026-03-30_00-06-26/LERs.csv};
+
+                \addlegendentryexpanded{Windowed, hard}
+
+                \addplot+[mark=*, line width=1pt, scol0]
+                table[
+                    col sep=comma, x=physical_p,
+                    y expr={1 - (1-\thisrow{LER_soft})^(1/12)}
+                ]
+                {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=minimum_sum,+experiment=soft_vs_hard_param_exploration,+system.F=1,+system.W=5/2026-03-30_00-06-26/LERs.csv};
+
+                \addlegendentryexpanded{Windowed, soft}
+
+                \addplot+[mark=o, line width=1pt, densely dashed, scol1]
+                table[
+                    col sep=comma, x=physical_p,
+                    y expr={1 - (1-\thisrow{LER_whole})^(1/12)}
+                ]
+                {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=minimum_sum,+experiment=windowed_vs_whole_param_exploration,+system.F=1,+system.W=5/2026-03-29_23-39-18/LERs.csv};
+
+                \addlegendentryexpanded{Whole window}
+            \end{axis}
+        \end{tikzpicture}
+        \caption{Min-Sum}
+    \end{subfigure}%
+    \begin{subfigure}{0.5\textwidth}
+        \begin{tikzpicture}
+            \begin{axis}[
+                    width=\figwidth,
+                    height=\figheight,
+                    ymode=log,
+                    legend style={
+                        cells={anchor=west},
+                        cells={align=left},
+                    },
+                    enlargelimits=false,
+                    ymin=8e-4, ymax=2.5e-1,
+                    grid=both,
+                    legend pos = south east,
+                    xtick={0.001,0.0015,...,0.004},
+                    xticklabel style={/pgf/number format/fixed},
+                    xticklabel style={/pgf/number format/precision=4},
+                    scaled x ticks=false,
+                    xlabel={Physical error rate},
+                    ylabel={Per-round-LER},
+                ]
+                \addplot+[mark=o, line width=1pt, densely dashed, scol0]
+                table[
+                    col sep=comma, x=physical_p,
+                    y expr={1 - (1-\thisrow{LER_hard})^(1/12)}
+                ]
+                {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=product_sum,+experiment=soft_vs_hard_param_exploration,+system.F=1,+system.W=5/2026-03-30_00-06-26/LERs.csv};
+
+                \addlegendentryexpanded{Windowed, hard}
+
+                \addplot+[mark=*, line width=1pt, scol0]
+                table[
+                    col sep=comma, x=physical_p,
+                    y expr={1 - (1-\thisrow{LER_soft})^(1/12)}
+                ]
+                {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=product_sum,+experiment=soft_vs_hard_param_exploration,+system.F=1,+system.W=5/2026-03-30_00-06-26/LERs.csv};
+
+                \addlegendentryexpanded{Windowed, soft}
+
+                \addplot+[mark=o, line width=1pt, densely dashed, scol1]
+                table[
+                    col sep=comma, x=physical_p,
+                    y expr={1 - (1-\thisrow{LER_whole})^(1/12)}
+                ]
+                {/home/andreas/workspace/private/ma-sw-results/outputs/+decoder.bp_method=product_sum,+experiment=windowed_vs_whole_param_exploration,+system.F=1,+system.W=5/2026-03-29_23-39-18/LERs.csv};
+
+                \addlegendentryexpanded{Whole window}
+            \end{axis}
+        \end{tikzpicture}
+        \caption{SPA}
+    \end{subfigure}
+
+    \caption{BP simulations for the [[144,12,12]]-BB code,
+    $W=5, F=1, n_\text{iter,BP} = 32$}
 \end{figure}
 
 \end{document}