1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
|
/* -*- c++ -*- */
/*
* Copyright 2015 Free Software Foundation, Inc.
*
* This 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.
*
* This software 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 this software; see the file COPYING. If not, write to
* the Free Software Foundation, Inc., 51 Franklin Street,
* Boston, MA 02110-1301, USA.
*/
#ifndef INCLUDED_fec_mtrx_impl_H
#define INCLUDED_fec_mtrx_impl_H
#include <string>
#include <gnuradio/fec/fec_mtrx.h>
#include <gsl/gsl_blas.h>
#include <gsl/gsl_linalg.h>
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_permutation.h>
#include <gsl/gsl_randist.h>
namespace gr {
namespace fec {
namespace code {
class fec_mtrx_impl : public fec_mtrx
{
protected:
//! Constructor
fec_mtrx_impl();
//! Codeword length n
unsigned int d_n;
//! Information word length k
unsigned int d_k;
//! Number of rows in the matrix read in from alist file
unsigned int d_num_rows;
//! Number of columns in the matrix read in from alist file
unsigned int d_num_cols;
//! GSL matrix structure for the parity check matrix
matrix_sptr d_H_sptr;
//! Flag for whether or not the parity bits come first or last
bool d_par_bits_last;
public:
//! Returns the parity check matrix H (needed by decoder)
const gsl_matrix* H() const;
//! Get the codeword length n
unsigned int n() const;
//! Get the information word length k
unsigned int k() const;
//! Subtract matrices using mod2 operations
void
add_matrices_mod2(gsl_matrix* result, const gsl_matrix*, const gsl_matrix*) const;
//! Multiply matrices using mod2 operations
void
mult_matrices_mod2(gsl_matrix* result, const gsl_matrix*, const gsl_matrix*) const;
//! Invert a square matrix using mod2 operations
gsl_matrix* calc_inverse_mod2(const gsl_matrix*) const;
/*!
* \brief Get Boolean for whether or not parity bits come first or last
* \details
* The decoder will need to know if the parity bits are
* coming first or last
*/
bool parity_bits_come_last() const;
virtual ~fec_mtrx_impl();
};
} // namespace code
} // namespace fec
} // namespace gr
#endif /* INCLUDED_fec_mtrx_impl_H */
|
