def __init__(self):
gr.top_block.__init__(self)
##################################################
# Variables
##################################################
self.sigma = sigma = 0.4
self.max_iterations = max_iterations = 50
self.alist_file = alist_file = "/home/manu/repos/ldpc/gr-ldpc/python/alist-files/96.3.963"
self.samp_rate = samp_rate = 32000
self.encoder = _encoder = ldpc_encoder_bf(alist_file)
self.decoder = _decoder = ldpc_decoder_fb(alist_file, sigma, max_iterations)
##################################################
# Blocks
##################################################
self.generic_encoder_0 = fec.encoder(_encoder, gr.sizeof_char, gr.sizeof_float)
self.generic_decoder_0 = fec.decoder(_decoder, gr.sizeof_float, gr.sizeof_char)
self.blocks_file_source_0 = blocks.file_source(gr.sizeof_char*1, "/home/manu/Downloads/in.flac", False)
self.blocks_file_sink_0 = blocks.file_sink(gr.sizeof_char*1, "/home/manu/Downloads/out.flac")
self.blocks_file_sink_0.set_unbuffered(False)
self.blocks_add_xx_0 = blocks.add_vff(1)
self.analog_noise_source_x_0 = analog.noise_source_f(analog.GR_GAUSSIAN, 0, 0)
##################################################
# Connections
##################################################
self.connect((self.analog_noise_source_x_0, 0), (self.blocks_add_xx_0, 1))
self.connect((self.blocks_file_source_0, 0), (self.generic_encoder_0, 0))
self.connect((self.generic_encoder_0, 0), (self.blocks_add_xx_0, 0))
self.connect((self.blocks_add_xx_0, 0), (self.generic_decoder_0, 0))
self.connect((self.generic_decoder_0, 0), (self.blocks_file_sink_0, 0))
def __init__(self):
grc_wxgui.top_block_gui.__init__(self, title="Top Block")
_icon_path = "/usr/share/icons/hicolor/32x32/apps/gnuradio-grc.png"
self.SetIcon(wx.Icon(_icon_path, wx.BITMAP_TYPE_ANY))
##################################################
# Variables
##################################################
self.sigma = sigma = 0.4
self.max_iterations = max_iterations = 50
self.alist_file = alist_file = "/home/manu/repos/ldpc/gr-ldpc/python/alist-files/96.3.963"
self.samp_rate = samp_rate = 32000
self.encoder = encoder = ldpc_encoder_bf(alist_file)
self.decoder = decoder = ldpc_decoder_fb(alist_file, sigma, max_iterations)
##################################################
# Blocks
##################################################
self.generic_encoder_0 = fec.encoder(encoder, gr.sizeof_char, gr.sizeof_float)
self.generic_decoder_0 = fec.decoder(decoder, gr.sizeof_float, gr.sizeof_char)
self.blocks_file_source_0 = blocks.file_source(gr.sizeof_char*1, "/home/manu/Downloads/in.flac", False)
self.blocks_file_sink_0 = blocks.file_sink(gr.sizeof_char*1, "/home/manu/Downloads/out.flac")
self.blocks_file_sink_0.set_unbuffered(False)
self.blocks_add_xx_0 = blocks.add_vff(1)
self.analog_noise_source_x_0 = analog.noise_source_f(analog.GR_GAUSSIAN, 0, 0)
##################################################
# Connections
##################################################
self.connect((self.analog_noise_source_x_0, 0), (self.blocks_add_xx_0, 1))
self.connect((self.blocks_file_source_0, 0), (self.generic_encoder_0, 0))
self.connect((self.generic_encoder_0, 0), (self.blocks_add_xx_0, 0))
self.connect((self.blocks_add_xx_0, 0), (self.generic_decoder_0, 0))
self.connect((self.generic_decoder_0, 0), (self.blocks_file_sink_0, 0))
def __init__(self, decoder_list_0, input_size, output_size):
gr.hier_block2.__init__(
self,
"Threaded Decoder",
gr.io_signature(1, 1, input_size * 1),
gr.io_signature(1, 1, output_size * 1),
)
##################################################
# Parameters
##################################################
self.decoder_list_0 = decoder_list_0
##################################################
# Variables
##################################################
##################################################
# Blocks
##################################################
self.fec_deinterleave_0 = fec.deinterleave(
input_size, fec.get_decoder_input_size(decoder_list_0[0]))
self.generic_decoders_0 = []
for i in range(len(decoder_list_0)):
self.generic_decoders_0.append(
fec.decoder(decoder_list_0[i], input_size, output_size))
self.fec_interleave_0 = fec.interleave(
output_size, fec.get_decoder_output_size(decoder_list_0[0]))
##################################################
# Connections
##################################################
for i in range(len(decoder_list_0)):
self.connect((self.fec_deinterleave_0, i),
(self.generic_decoders_0[i], 0))
for i in range(len(decoder_list_0)):
self.connect((self.generic_decoders_0[i], 0),
(self.fec_interleave_0, i))
self.connect((self, 0), (self.fec_deinterleave_0, 0))
self.connect((self.fec_interleave_0, 0), (self, 0))
def __init__(self, decoder_list_0, input_size, output_size):
gr.hier_block2.__init__(
self, "Threaded Decoder",
gr.io_signature(1, 1, input_size*1),
gr.io_signature(1, 1, output_size*1),
)
##################################################
# Parameters
##################################################
self.decoder_list_0 = decoder_list_0
##################################################
# Variables
##################################################
##################################################
# Blocks
##################################################
self.fec_deinterleave_0 = fec.deinterleave(input_size, fec.get_decoder_input_size(decoder_list_0[0]))
self.generic_decoders_0 = [];
for i in range(len(decoder_list_0)):
self.generic_decoders_0.append(fec.decoder(decoder_list_0[i], input_size, output_size))
self.fec_interleave_0 = fec.interleave(output_size, fec.get_decoder_output_size(decoder_list_0[0]))
##################################################
# Connections
##################################################
for i in range(len(decoder_list_0)):
self.connect((self.fec_deinterleave_0, i), (self.generic_decoders_0[i], 0))
for i in range(len(decoder_list_0)):
self.connect((self.generic_decoders_0[i], 0), (self.fec_interleave_0, i))
self.connect((self, 0), (self.fec_deinterleave_0, 0))
self.connect((self.fec_interleave_0, 0), (self, 0))
def __init__(self, decoder_list_0, input_size, output_size):
gr.hier_block2.__init__(
self, "Capillary Threaded Decoder",
gr.io_signature(1, 1, input_size*1),
gr.io_signature(1, 1, output_size*1),
)
##################################################
# Parameters
##################################################
self.decoder_list_0 = decoder_list_0
##################################################
# Variables
##################################################
##################################################
# Blocks
##################################################
self.fec_deinterleaves_0 = [];
for i in range(int(math.log(len(decoder_list_0), 2))):
for j in range(int(math.pow(2, i))):
self.fec_deinterleaves_0.append(fec.deinterleave(input_size, fec.get_decoder_input_size(decoder_list_0[0])))
self.generic_decoders_0 = [];
for i in range(len(decoder_list_0)):
self.generic_decoders_0.append(fec.decoder(decoder_list_0[i], input_size, output_size))
self.fec_interleaves_0 = [];
for i in range(int(math.log(len(decoder_list_0), 2))):
for j in range(int(math.pow(2, i))):
self.fec_interleaves_0.append(fec.interleave(output_size, fec.get_decoder_output_size(decoder_list_0[0])))
##################################################
# Connections
##################################################
rootcount = 0;
branchcount = 1;
for i in range(int(math.log(len(decoder_list_0), 2)) - 1):
for j in range(int(math.pow(2, i))):
self.connect((self.fec_deinterleaves_0[rootcount], 0), (self.fec_deinterleaves_0[branchcount], 0))
self.connect((self.fec_deinterleaves_0[rootcount], 1), (self.fec_deinterleaves_0[branchcount + 1], 0))
rootcount += 1;
branchcount += 2;
codercount = 0;
for i in range(len(decoder_list_0)/2):
self.connect((self.fec_deinterleaves_0[rootcount], 0), (self.generic_decoders_0[codercount], 0))
self.connect((self.fec_deinterleaves_0[rootcount], 1), (self.generic_decoders_0[codercount + 1], 0))
rootcount += 1;
codercount += 2;
rootcount = 0;
branchcount = 1;
for i in range(int(math.log(len(decoder_list_0), 2)) - 1):
for j in range(int(math.pow(2, i))):
self.connect((self.fec_interleaves_0[branchcount], 0), (self.fec_interleaves_0[rootcount], 0))
self.connect((self.fec_interleaves_0[branchcount + 1], 0), (self.fec_interleaves_0[rootcount], 1))
rootcount += 1;
branchcount += 2;
codercount = 0;
for i in range(len(decoder_list_0)/2):
self.connect((self.generic_decoders_0[codercount], 0), (self.fec_interleaves_0[rootcount], 0))
self.connect((self.generic_decoders_0[codercount + 1], 0), (self.fec_interleaves_0[rootcount], 1))
rootcount += 1;
codercount += 2;
if ((len(self.decoder_list_0)) > 1):
self.connect((self, 0), (self.fec_deinterleaves_0[0], 0))
self.connect((self.fec_interleaves_0[0], 0), (self, 0))
else:
self.connect((self, 0), (self.generic_decoders_0[0], 0))
self.connect((self.generic_decoders_0[0], 0), (self, 0))
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_interface",
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
message_collector_connected = False
##anything going through the annihilator needs shifted, uchar vals
if fec.get_conversion(
decoder_obj_list[0]) == "uchar" or fec.get_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_conversion(decoder_obj_list[0]) == "packed_bits":
self.blocks.append(blocks.add_const_ff(128.0))
if fec.get_conversion(
decoder_obj_list[0]) == "uchar" or fec.get_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 read_big_bitlist(ann):
cat.append(i)
synd_garble = .49
idx_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 corr_bb'
print 'ceiling: .0335 data garble rate'
self.blocks.append(
fec.corr_bb(cat,
len(puncpat) - puncpat.count('0'), len(ann),
integration_period, flush, synd_garble))
#print puncpat
if self.puncpat != '11':
self.blocks.append(
fec.reinflate_bb(0, read_bitlist(puncpat), puncpat.count('0'),
len(puncpat)))
if fec.get_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):
assert fec.get_history(decoder_obj_list[0]) == 0
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_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))
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_interface",
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;
message_collector_connected=False;
##anything going through the annihilator needs shifted, uchar vals
if fec.get_conversion(decoder_obj_list[0]) == "uchar" or fec.get_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_conversion(decoder_obj_list[0]) == "packed_bits":
self.blocks.append(blocks.add_const_ff(128.0));
if fec.get_conversion(decoder_obj_list[0]) == "uchar" or fec.get_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 read_big_bitlist(ann):
cat.append(i);
synd_garble = .49;
idx_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 corr_bb'
print 'ceiling: .0335 data garble rate'
self.blocks.append(fec.corr_bb(cat, len(puncpat) - puncpat.count('0'), len(ann), integration_period, flush, synd_garble));
#print puncpat
if self.puncpat != '11':
self.blocks.append(fec.reinflate_bb(0, read_bitlist(puncpat), puncpat.count('0'), len(puncpat)));
if fec.get_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):
assert fec.get_history(decoder_obj_list[0]) == 0;
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_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));
def __init__(self, decoder_list_0, input_size, output_size):
gr.hier_block2.__init__(
self,
"Capillary Threaded Decoder",
gr.io_signature(1, 1, input_size * 1),
gr.io_signature(1, 1, output_size * 1),
)
##################################################
# Parameters
##################################################
self.decoder_list_0 = decoder_list_0
##################################################
# Variables
##################################################
##################################################
# Blocks
##################################################
self.fec_deinterleaves_0 = []
for i in range(int(math.log(len(decoder_list_0), 2))):
for j in range(int(math.pow(2, i))):
self.fec_deinterleaves_0.append(
fec.deinterleave(input_size, fec.get_decoder_input_size(decoder_list_0[0]))
)
self.generic_decoders_0 = []
for i in range(len(decoder_list_0)):
self.generic_decoders_0.append(fec.decoder(decoder_list_0[i], input_size, output_size))
self.fec_interleaves_0 = []
for i in range(int(math.log(len(decoder_list_0), 2))):
for j in range(int(math.pow(2, i))):
self.fec_interleaves_0.append(
fec.interleave(output_size, fec.get_decoder_output_size(decoder_list_0[0]))
)
##################################################
# Connections
##################################################
rootcount = 0
branchcount = 1
for i in range(int(math.log(len(decoder_list_0), 2)) - 1):
for j in range(int(math.pow(2, i))):
self.connect((self.fec_deinterleaves_0[rootcount], 0), (self.fec_deinterleaves_0[branchcount], 0))
self.connect((self.fec_deinterleaves_0[rootcount], 1), (self.fec_deinterleaves_0[branchcount + 1], 0))
rootcount += 1
branchcount += 2
codercount = 0
for i in range(len(decoder_list_0) / 2):
self.connect((self.fec_deinterleaves_0[rootcount], 0), (self.generic_decoders_0[codercount], 0))
self.connect((self.fec_deinterleaves_0[rootcount], 1), (self.generic_decoders_0[codercount + 1], 0))
rootcount += 1
codercount += 2
rootcount = 0
branchcount = 1
for i in range(int(math.log(len(decoder_list_0), 2)) - 1):
for j in range(int(math.pow(2, i))):
self.connect((self.fec_interleaves_0[branchcount], 0), (self.fec_interleaves_0[rootcount], 0))
self.connect((self.fec_interleaves_0[branchcount + 1], 0), (self.fec_interleaves_0[rootcount], 1))
rootcount += 1
branchcount += 2
codercount = 0
for i in range(len(decoder_list_0) / 2):
self.connect((self.generic_decoders_0[codercount], 0), (self.fec_interleaves_0[rootcount], 0))
self.connect((self.generic_decoders_0[codercount + 1], 0), (self.fec_interleaves_0[rootcount], 1))
rootcount += 1
codercount += 2
if (len(self.decoder_list_0)) > 1:
self.connect((self, 0), (self.fec_deinterleaves_0[0], 0))
self.connect((self.fec_interleaves_0[0], 0), (self, 0))
else:
self.connect((self, 0), (self.generic_decoders_0[0], 0))
self.connect((self.generic_decoders_0[0], 0), (self, 0))
