diff options
Diffstat (limited to 'gr-fec/python/fec/qa_ber_bf.py')
| -rw-r--r-- | gr-fec/python/fec/qa_ber_bf.py | 141 |
1 files changed, 141 insertions, 0 deletions
diff --git a/gr-fec/python/fec/qa_ber_bf.py b/gr-fec/python/fec/qa_ber_bf.py new file mode 100644 index 0000000000..1a60c68581 --- /dev/null +++ b/gr-fec/python/fec/qa_ber_bf.py @@ -0,0 +1,141 @@ +#!/usr/bin/env python +# +# Copyright 2014 Free Software Foundation, Inc. +# +# This file is part of GNU Radio +# +# GNU Radio is free software; you can redistribute it and/or modify +# it under the terms of the GNU General Public License as published by +# the Free Software Foundation; either version 3, or (at your option) +# any later version. +# +# GNU Radio is distributed in the hope that it will be useful, +# but WITHOUT ANY WARRANTY; without even the implied warranty of +# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +# GNU General Public License for more details. +# +# You should have received a copy of the GNU General Public License +# along with GNU Radio; see the file COPYING. If not, write to +# the Free Software Foundation, Inc., 51 Franklin Street, +# Boston, MA 02110-1301, USA. +# + +from gnuradio import gr, gr_unittest +import fec_swig as fec +import blocks_swig as blocks +import numpy, copy + +class test_ber_bf(gr_unittest.TestCase): + + def setUp(self): + self.tb = gr.top_block() + + def tearDown(self): + self.tb = None + + def test_000(self): + # Cause a single bit error out of 8*N bits + # using streaming mode + + mode = 0 + N = 10000 + data0 = numpy.random.randint(0, 256, N).tolist() + data1 = copy.deepcopy(data0) + data1[0] ^= 0x01 + + src0 = blocks.vector_source_b(data0) + src1 = blocks.vector_source_b(data1) + op = fec.ber_bf(mode) + dst = blocks.vector_sink_f() + + self.tb.connect(src0, (op,0)) + self.tb.connect(src1, (op,1)) + self.tb.connect(op, dst) + self.tb.run() + + data = dst.data() + expected_result = [numpy.log10(1.0/(8.0*N)),] + + self.assertFloatTuplesAlmostEqual(expected_result, data, 5) + + def test_001(self): + # Cause a single bit error out of 8*N bits + # using test mode + + mode = 1 + N = 1000 + data0 = numpy.random.randint(0, 256, N).tolist() + data1 = copy.deepcopy(data0) + data1[0] ^= 0x01 + + src0 = blocks.vector_source_b(data0) + src1 = blocks.vector_source_b(data1) + op = fec.ber_bf(mode, 1) + dst = blocks.vector_sink_f() + + self.tb.connect(src0, (op,0)) + self.tb.connect(src1, (op,1)) + self.tb.connect(op, dst) + self.tb.run() + + data = dst.data() + expected_result = [numpy.log10(1.0/(8.0*N)),] + + self.assertFloatTuplesAlmostEqual(expected_result, data, 5) + + def test_002(self): + # Cause 8 bit errors out of 8*N bits + # using test mode + + mode = 1 + N = 1000 + data0 = numpy.random.randint(0, 256, N).tolist() + data1 = copy.deepcopy(data0) + data1[0] ^= 0xFF + + src0 = blocks.vector_source_b(data0) + src1 = blocks.vector_source_b(data1) + op = fec.ber_bf(mode, 1, -2.0) + dst = blocks.vector_sink_f() + + self.tb.connect(src0, (op,0)) + self.tb.connect(src1, (op,1)) + self.tb.connect(op, dst) + self.tb.run() + + data = dst.data() + expected_result = [numpy.log10(8.0/(8.0*N)),] + + self.assertFloatTuplesAlmostEqual(expected_result, data, 5) + + def test_003(self): + # Cause a 8 bit errors out of 8*N bits + # using test mode + # Exit if BER < -2.0 + + mode = 1 + N = 1000 + data0 = numpy.random.randint(0, 256, N).tolist() + data1 = copy.deepcopy(data0) + data1[0] ^= 0xFF + + src0 = blocks.vector_source_b(data0) + src1 = blocks.vector_source_b(data1) + op = fec.ber_bf(mode, 10, -2.0) + dst = blocks.vector_sink_f() + + self.tb.connect(src0, (op,0)) + self.tb.connect(src1, (op,1)) + self.tb.connect(op, dst) + self.tb.run() + + data = dst.data() + expected_result = [-2.0,] + + print data + print expected_result + + self.assertFloatTuplesAlmostEqual(expected_result, data, 5) + +if __name__ == '__main__': + gr_unittest.run(test_ber_bf, "test_ber_bf.xml") |