From 3dd0863aae3560decf80b8cf9955ec5ceca02038 Mon Sep 17 00:00:00 2001 From: Andreas Tsouchlos Date: Mon, 30 Mar 2026 10:49:20 +0200 Subject: [PATCH] Make plots show per-round-LER; Add windowed soft info decoding vs whole decoding comparison --- src/param-exploration/main.tex | 241 +++++++++++++++++++++++++++------ 1 file changed, 202 insertions(+), 39 deletions(-) diff --git a/src/param-exploration/main.tex b/src/param-exploration/main.tex index f75197d..3cee21c 100644 --- a/src/param-exploration/main.tex +++ b/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] - table[col sep=comma, x=physical_p, y=LER_whole] + \addplot+[ + 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] - % {/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}; + 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=\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}