diff options
Diffstat (limited to 'gr-fec/lib/reed-solomon/exercise.c')
| -rw-r--r-- | gr-fec/lib/reed-solomon/exercise.c | 138 |
1 files changed, 70 insertions, 68 deletions
diff --git a/gr-fec/lib/reed-solomon/exercise.c b/gr-fec/lib/reed-solomon/exercise.c index 41ed4f2148..562a051f17 100644 --- a/gr-fec/lib/reed-solomon/exercise.c +++ b/gr-fec/lib/reed-solomon/exercise.c @@ -34,101 +34,103 @@ #elif defined(CCSDS) #define PRINTPARM printf("CCSDS (255,223):"); #else -#define PRINTPARM printf("(%d,%d):",rs->nn,rs->nn-rs->nroots); +#define PRINTPARM printf("(%d,%d):", rs->nn, rs->nn - rs->nroots); #endif /* Exercise the RS codec passed as an argument */ int EXERCISE( #if !defined(CCSDS) && !defined(FIXED) -void *p, + void* p, #endif -int trials){ + int trials) +{ #if !defined(CCSDS) && !defined(FIXED) - struct rs *rs = (struct rs *)p; + struct rs* rs = (struct rs*)p; #endif #if MAX_ARRAY - DTYPE block[MAX_ARRAY],tblock[MAX_ARRAY]; - unsigned int i; - int errors; - int errlocs[MAX_ARRAY]; - int derrlocs[MAX_ARRAY]; + DTYPE block[MAX_ARRAY], tblock[MAX_ARRAY]; + unsigned int i; + int errors; + int errlocs[MAX_ARRAY]; + int derrlocs[MAX_ARRAY]; #else - DTYPE block[NN],tblock[NN]; - unsigned int i; - int errors; - int errlocs[NN]; - int derrlocs[NROOTS]; + DTYPE block[NN], tblock[NN]; + unsigned int i; + int errors; + int errlocs[NN]; + int derrlocs[NROOTS]; #endif - int derrors; - int errval,errloc; - int erasures; - int decoder_errors = 0; + int derrors; + int errval, errloc; + int erasures; + int decoder_errors = 0; - while(trials-- != 0){ - /* Test up to the error correction capacity of the code */ - for(errors=0;(unsigned int)errors <= NROOTS/2;errors++){ + while (trials-- != 0) { + /* Test up to the error correction capacity of the code */ + for (errors = 0; (unsigned int)errors <= NROOTS / 2; errors++) { - /* Load block with random data and encode */ - for(i=0;i<NN-NROOTS;i++) - block[i] = rand() & NN; + /* Load block with random data and encode */ + for (i = 0; i < NN - NROOTS; i++) + block[i] = rand() & NN; #if defined(CCSDS) || defined(FIXED) - ENCODE_RS(&block[0],&block[NN-NROOTS]); + ENCODE_RS(&block[0], &block[NN - NROOTS]); #else - ENCODE_RS(rs,&block[0],&block[NN-NROOTS]); + ENCODE_RS(rs, &block[0], &block[NN - NROOTS]); #endif - /* Make temp copy, seed with errors */ - memcpy(tblock,block,sizeof(tblock)); - memset(errlocs,0,sizeof(errlocs)); - memset(derrlocs,0,sizeof(derrlocs)); - erasures=0; - for(i=0;i<(unsigned int)errors;i++){ - do { - errval = rand() & NN; - } while(errval == 0); /* Error value must be nonzero */ + /* Make temp copy, seed with errors */ + memcpy(tblock, block, sizeof(tblock)); + memset(errlocs, 0, sizeof(errlocs)); + memset(derrlocs, 0, sizeof(derrlocs)); + erasures = 0; + for (i = 0; i < (unsigned int)errors; i++) { + do { + errval = rand() & NN; + } while (errval == 0); /* Error value must be nonzero */ - do { - errloc = rand() % NN; - } while(errlocs[errloc] != 0); /* Must not choose the same location twice */ + do { + errloc = rand() % NN; + } while (errlocs[errloc] != + 0); /* Must not choose the same location twice */ - errlocs[errloc] = 1; + errlocs[errloc] = 1; #if FLAG_ERASURE - if(rand() & 1) /* 50-50 chance */ - derrlocs[erasures++] = errloc; + if (rand() & 1) /* 50-50 chance */ + derrlocs[erasures++] = errloc; #endif - tblock[errloc] ^= errval; - } + tblock[errloc] ^= errval; + } - /* Decode the errored block */ + /* Decode the errored block */ #if defined(CCSDS) || defined(FIXED) - derrors = DECODE_RS(tblock,derrlocs,erasures); + derrors = DECODE_RS(tblock, derrlocs, erasures); #else - derrors = DECODE_RS(rs,tblock,derrlocs,erasures); + derrors = DECODE_RS(rs, tblock, derrlocs, erasures); #endif - if(derrors != errors){ - PRINTPARM - printf(" decoder says %d errors, true number is %d\n",derrors,errors); - decoder_errors++; - } - for(i=0;i<(unsigned int)derrors;i++){ - if(errlocs[derrlocs[i]] == 0){ - PRINTPARM - printf(" decoder indicates error in location %d without error\n",i); - decoder_errors++; - } - } - if(memcmp(tblock,block,sizeof(tblock)) != 0){ - PRINTPARM - printf(" uncorrected errors! output ^ input:"); - decoder_errors++; - for(i=0;i<NN;i++) - printf(" %02x",tblock[i] ^ block[i]); - printf("\n"); - } + if (derrors != errors) { + PRINTPARM + printf(" decoder says %d errors, true number is %d\n", derrors, errors); + decoder_errors++; + } + for (i = 0; i < (unsigned int)derrors; i++) { + if (errlocs[derrlocs[i]] == 0) { + PRINTPARM + printf(" decoder indicates error in location %d without error\n", i); + decoder_errors++; + } + } + if (memcmp(tblock, block, sizeof(tblock)) != 0) { + PRINTPARM + printf(" uncorrected errors! output ^ input:"); + decoder_errors++; + for (i = 0; i < NN; i++) + printf(" %02x", tblock[i] ^ block[i]); + printf("\n"); + } + } } - } - return decoder_errors; + return decoder_errors; } |