#!/usr/bin/env python
#
# Copyright 2014 Free Software Foundation, Inc.
#
# This file is part of GNU Radio
#
# SPDX-License-Identifier: GPL-3.0-or-later
#
#
from __future__ import print_function
from __future__ import absolute_import
from __future__ import division
from __future__ import unicode_literals
from gnuradio import gr, blocks, digital
from . import fec_python as fec
from .bitflip import *
from .threaded_decoder import threaded_decoder
from .capillary_threaded_decoder import capillary_threaded_decoder
class extended_decoder(gr.hier_block2):
#solution to log_(1-2*t)(1-2*.0335) = 1/taps where t is thresh (syndrome density)
#for i in numpy.arange(.1, .499, .01):
#print(str(log((1-(2 * .035)), (1-(2 * i)))) + ':' + str(i);)
garbletable = {
0.310786835319:0.1,
0.279118162802:0.11,
0.252699589071:0.12,
0.230318516016:0.13,
0.211108735347:0.14,
0.194434959095:0.15,
0.179820650401:0.16,
0.166901324951:0.17,
0.15539341766:0.18,
0.145072979886:0.19,
0.135760766313:0.2,
0.127311581396:0.21,
0.119606529806:0.22,
0.112547286766:0.23,
0.106051798775:0.24,
0.10005101381:0.25,
0.0944863633098:0.26,
0.0893078003966:0.27,
0.084472254501:0.28,
0.0799424008658:0.29,
0.0756856701944:0.3,
0.0716734425668:0.31,
0.0678803831565:0.32,
0.0642838867856:0.33,
0.0608636049994:0.34,
0.0576010337489:0.35,
0.0544791422522:0.36,
0.0514820241933:0.37,
0.0485945507251:0.38,
0.0458019998183:0.39,
0.0430896262596:0.4,
0.0404421166935:0.41,
0.0378428350972:0.42,
0.0352726843274:0.43,
0.0327082350617:0.44,
0.0301183562535:0.45,
0.0274574540266:0.46,
0.0246498236897:0.47,
0.0215448131298:0.48,
0.0177274208353:0.49,
}
def __init__(self, decoder_obj_list, threading, ann=None, puncpat='11',
integration_period=10000, flush=None, rotator=None):
gr.hier_block2.__init__(self, "extended_decoder",
gr.io_signature(1, 1, gr.sizeof_float),
gr.io_signature(1, 1, gr.sizeof_char))
self.blocks=[]
self.ann=ann
self.puncpat=puncpat
self.flush=flush
if(type(decoder_obj_list) == list):
if(type(decoder_obj_list[0]) == list):
gr.log.info("fec.extended_decoder: Parallelism must be 1.")
raise AttributeError
else:
# If it has parallelism of 0, force it into a list of 1
decoder_obj_list = [decoder_obj_list,]
message_collector_connected=False
##anything going through the annihilator needs shifted, uchar vals
if fec.get_decoder_input_conversion(decoder_obj_list[0]) == "uchar" or \
fec.get_decoder_input_conversion(decoder_obj_list[0]) == "packed_bits":
self.blocks.append(blocks.multiply_const_ff(48.0))
if fec.get_shift(decoder_obj_list[0]) != 0.0:
self.blocks.append(blocks.add_const_ff(fec.get_shift(decoder_obj_list[0])))
elif fec.get_decoder_input_conversion(decoder_obj_list[0]) == "packed_bits":
self.blocks.append(blocks.add_const_ff(128.0))
if fec.get_decoder_input_conversion(decoder_obj_list[0]) == "uchar" or \
fec.get_decoder_input_conversion(decoder_obj_list[0]) == "packed_bits":
self.blocks.append(blocks.float_to_uchar());
const_index = 0; #index that corresponds to mod order for specinvert purposes
if not self.flush:
flush = 10000;
else:
flush = self.flush;
if self.ann: #ann and puncpat are strings of 0s and 1s
cat = fec.ULLVector();
for i in fec.read_big_bitlist(ann):
cat.append(i);
synd_garble = .49
idx_list = list(self.garbletable.keys())
idx_list.sort()
for i in idx_list:
if 1.0 / self.ann.count('1') >= i:
synd_garble = self.garbletable[i]
print('using syndrom garble threshold ' + str(synd_garble) + 'for conv_bit_corr_bb')
print('ceiling: .0335 data garble rate')
self.blocks.append(fec.conv_bit_corr_bb(cat, len(puncpat) - puncpat.count('0'),
len(ann), integration_period, flush, synd_garble))
if self.puncpat != '11':
self.blocks.append(fec.depuncture_bb(len(puncpat), read_bitlist(puncpat), 0))
if fec.get_decoder_input_conversion(decoder_obj_list[0]) == "packed_bits":
self.blocks.append(blocks.uchar_to_float())
self.blocks.append(blocks.add_const_ff(-128.0))
self.blocks.append(digital.binary_slicer_fb())
self.blocks.append(blocks.unpacked_to_packed_bb(1,0))
if(len(decoder_obj_list) > 1):
if(fec.get_history(decoder_obj_list[0]) != 0):
gr.log.info("fec.extended_decoder: Cannot use multi-threaded parallelism on a decoder with history.")
raise AttributeError
if threading == 'capillary':
self.blocks.append(capillary_threaded_decoder(decoder_obj_list,
fec.get_decoder_input_item_size(decoder_obj_list[0]),
fec.get_decoder_output_item_size(decoder_obj_list[0])))
elif threading == 'ordinary':
self.blocks.append(threaded_decoder(decoder_obj_list,
fec.get_decoder_input_item_size(decoder_obj_list[0]),
fec.get_decoder_output_item_size(decoder_obj_list[0])))
else:
self.blocks.append(fec.decoder(decoder_obj_list[0],
fec.get_decoder_input_item_size(decoder_obj_list[0]),
fec.get_decoder_output_item_size(decoder_obj_list[0])))
if fec.get_decoder_output_conversion(decoder_obj_list[0]) == "unpack":
self.blocks.append(blocks.packed_to_unpacked_bb(1, gr.GR_MSB_FIRST));
self.connect((self, 0), (self.blocks[0], 0));
self.connect((self.blocks[-1], 0), (self, 0));
for i in range(len(self.blocks) - 1):
self.connect((self.blocks[i], 0), (self.blocks[i+1], 0));