23 #ifndef INCLUDED_volk_8u_conv_k7_r2puppet_8u_H
24 #define INCLUDED_volk_8u_conv_k7_r2puppet_8u_H
36 static inline int parity(
int x,
unsigned char* Partab)
45 unsigned int endstate,
46 unsigned int tailsize,
47 unsigned char* decisions)
51 int d_numstates = (1 << 6);
52 int d_decision_t_size = d_numstates/8;
54 int d_framebits = nbits;
61 endstate = (endstate%d_numstates) << d_ADDSHIFT;
68 d += tailsize * d_decision_t_size ;
70 int dif = tailsize - (d_k - 1);
73 while(nbits-- > d_framebits - (d_k - 1)) {
75 dec.
t = &d[nbits * d_decision_t_size];
76 k = (dec.w[(endstate>>d_ADDSHIFT)/32] >> ((endstate>>d_ADDSHIFT)%32)) & 1;
78 endstate = (endstate >> 1) | (k << (d_k-2+d_ADDSHIFT));
81 data[((nbits+dif)%d_framebits)] = k;
90 dec.t = &d[nbits * d_decision_t_size];
92 k = (dec.w[(endstate>>d_ADDSHIFT)/32] >> ((endstate>>d_ADDSHIFT)%32)) & 1;
94 endstate = (endstate >> 1) | (k << (d_k-2+d_ADDSHIFT));
95 data[((nbits+dif)%d_framebits)] = k;
100 return retval >> d_ADDSHIFT;
106 #include <pmmintrin.h>
107 #include <emmintrin.h>
108 #include <xmmintrin.h>
109 #include <mmintrin.h>
116 static inline void volk_8u_conv_k7_r2puppet_8u_spiral(
unsigned char* syms,
unsigned char* dec,
unsigned int framebits) {
120 int d_numstates = (1 << 6);
122 static unsigned char* D;
123 static unsigned char* Y;
124 static unsigned char* X;
125 static unsigned int excess = 6;
126 static unsigned char* Branchtab;
127 static unsigned char Partab[256];
129 int d_polys[2] = {79, 109};
153 for(state=0;state < d_numstates/2;state++){
154 for(i=0; i<rate; i++){
155 Branchtab[i*d_numstates/2+state] = (d_polys[i] < 0) ^
parity((2*state) & abs(d_polys[i]), Partab) ? 255 : 0;
163 memset(X, 31, d_numstates);
165 volk_8u_x4_conv_k7_r2_8u_spiral(Y, X, syms, D, framebits/2 - excess, excess, Branchtab);
167 unsigned int min = X[0];
168 int i = 0, state = 0;
169 for(i = 0; i < (d_numstates); ++i) {
190 static inline void volk_8u_conv_k7_r2puppet_8u_generic(
unsigned char* syms,
unsigned char* dec,
unsigned int framebits) {
195 int d_numstates = (1 << 6);
197 static unsigned char* Y;
198 static unsigned char* X;
199 static unsigned char* D;
200 static unsigned int excess = 6;
201 static unsigned char* Branchtab;
202 static unsigned char Partab[256];
204 int d_polys[2] = {79, 109};
229 for(state=0;state < d_numstates/2;state++){
230 for(i=0; i<rate; i++){
231 Branchtab[i*d_numstates/2+state] = (d_polys[i] < 0) ^
parity((2*state) & abs(d_polys[i]), Partab) ? 255 : 0;
242 memset(X, 31, d_numstates);
244 volk_8u_x4_conv_k7_r2_8u_generic(Y, X, syms, D, framebits/2 - excess, excess, Branchtab);
246 unsigned int min = X[0];
247 int i = 0, state = 0;
248 for(i = 0; i < (d_numstates); ++i) {
unsigned char * t
Definition: volk_8u_conv_k7_r2puppet_8u.h:32
float min(float a, float b)
Definition: volk_8u_conv_k7_r2puppet_8u.h:30
VOLK_API size_t volk_get_alignment(void)
Get the machine alignment in bytes.
static int chainback_viterbi(unsigned char *data, unsigned int nbits, unsigned int endstate, unsigned int tailsize, unsigned char *decisions)
Definition: volk_8u_conv_k7_r2puppet_8u.h:43
__VOLK_DECL_BEGIN VOLK_API void * volk_malloc(size_t size, size_t alignment)
Allocate size bytes of data aligned to alignment.
static int parity(int x, unsigned char *Partab)
Definition: volk_8u_conv_k7_r2puppet_8u.h:36
unsigned int * w
Definition: volk_8u_conv_k7_r2puppet_8u.h:33