From 1ef4eb0245942d15d15459e938e921d116a4cf96 Mon Sep 17 00:00:00 2001 From: Parvfect Date: Wed, 3 Apr 2024 11:48:35 +0100 Subject: [PATCH] progres? --- distracted_coupon_collector.py | 6 +++--- tanner_qspa.py | 6 +++--- 2 files changed, 6 insertions(+), 6 deletions(-) diff --git a/distracted_coupon_collector.py b/distracted_coupon_collector.py index ee2aee6..ba8b3ae 100644 --- a/distracted_coupon_collector.py +++ b/distracted_coupon_collector.py @@ -196,7 +196,7 @@ def simulate_reads(C, symbols, read_length, P, n_motifs, n_picks): return likelihood_arr -def decoding_errors_fer(k, n, dv, dc, P, H, G, GF, graph, C, symbols, n_motifs, n_picks, decoder=None, decoding_failures_parameter=10, max_iterations=50, iterations=50, uncoded=False, bec_decoder=False, label=None, code_class="", read_lengths=np.arange(1,20)): +def decoding_errors_fer(k, n, dv, dc, P, H, G, GF, graph, C, symbols, n_motifs, n_picks, decoder=None, decoding_failures_parameter=20, max_iterations=50, iterations=50, uncoded=False, bec_decoder=False, label=None, code_class="", read_lengths=np.arange(1,20)): frame_error_rate = [] max_iterations = max_iterations @@ -269,8 +269,8 @@ def run_fer(n_motifs, n_picks, dv, dc, k, n, L, M, ffdim, P, code_class="", iter with Profile() as prof: n_motifs, n_picks = 8, 4 dv, dc, ffdim, P = 3, 9, 67, 2 * 0.038860387943791645 - k, n = 180, 270 - L, M = 12, 36 + k, n = 22, 33 + L, M = 12, 51 read_length = 6 read_lengths = np.arange(7, 13) diff --git a/tanner_qspa.py b/tanner_qspa.py index d038f87..6239cb2 100644 --- a/tanner_qspa.py +++ b/tanner_qspa.py @@ -15,7 +15,7 @@ def initialize_vn_links(self, P): for vn in self.vns: vn_index = vn.identifier for cn_index in vn.links: - self.vn_links[(cn_index, vn_index)] = P[vn_index] + self.vn_links[(cn_index, vn_index)] = P[vn_index].copy() def initialize_cn_links(self): """ Initializes CN Links""" @@ -24,7 +24,7 @@ def initialize_cn_links(self): for vn_index in cn.links: self.cn_links[(cn_index, vn_index)] = np.zeros(67) - def qspa_decode(self, symbol_likelihood_arr, H, GF, max_iterations=5): + def qspa_decode(self, symbol_likelihood_arr, H, GF, max_iterations=20): """Decodes using QSPA """ self.GF = GF @@ -64,7 +64,7 @@ def qspa_decode(self, symbol_likelihood_arr, H, GF, max_iterations=5): prev_max_prob_codeword = max_prob_codeword iterations+=1 - print(f"Iteration {iterations}") + #print(f"Iteration {iterations}") print("Decoding does not converge") return max_prob_codeword