| 1 |
/* |
| 2 |
* include needed system headers |
| 3 |
*/ |
| 4 |
#include <stdio.h> /* include standard i/o library */ |
| 5 |
#include <stdlib.h> /* include standard library */ |
| 6 |
#include <string.h> /* include string library */ |
| 7 |
#include <unistd.h> /* include unix standard library */ |
| 8 |
#include <sys/types.h> /* */ |
| 9 |
#include <sys/stat.h> /* */ |
| 10 |
#include <fcntl.h> /* */ |
| 11 |
#include <errno.h> /* error simbol definitions */ |
| 12 |
#include <time.h> /* system time definitions */ |
| 13 |
#include <math.h> /* math library */ |
| 14 |
|
| 15 |
void unpacking(int ndata, unsigned short *output); |
| 16 |
|
| 17 |
void dunpacker_(int *l_buffer, unsigned short *buffer, int *error, int *fd) |
| 18 |
|
| 19 |
{ |
| 20 |
int i; |
| 21 |
int nhead,nread; |
| 22 |
unsigned short word; |
| 23 |
|
| 24 |
*error = 0; |
| 25 |
|
| 26 |
nhead = 13; /* packed */ |
| 27 |
|
| 28 |
// printf("%i \n",*l_buffer); |
| 29 |
|
| 30 |
|
| 31 |
/* /\*** read the data ***\/ */ |
| 32 |
// if( nhead == 13) *l_buffer = (*l_buffer*13)/16+1; /* packed */ |
| 33 |
nread=read(*fd, buffer, 2**l_buffer); |
| 34 |
|
| 35 |
if( nread==0 ) { |
| 36 |
// printf("dunpacker: Unexpected end of file %d at %d \n",*fd,nread); |
| 37 |
*error=1; |
| 38 |
return ; |
| 39 |
} |
| 40 |
if ( nread != 2**l_buffer) { |
| 41 |
// printf("dunpacker: Error on reading: nread = %d \n",nread); |
| 42 |
*error=-1; |
| 43 |
return ; |
| 44 |
} |
| 45 |
|
| 46 |
|
| 47 |
/*** correct byte ***/ |
| 48 |
for(i=0; i<*l_buffer; i++) { |
| 49 |
word = buffer[i]; |
| 50 |
buffer[i] = ( (word&0x00ff) << 8 ) | ( (word&0xff00) >> 8 ); |
| 51 |
} |
| 52 |
|
| 53 |
/*** unpacking buffer data ***/ |
| 54 |
if( nhead == 13) unpacking(*l_buffer, buffer); |
| 55 |
|
| 56 |
// int num = (*l_buffer-13-1)/13*16+3;//??? |
| 57 |
// for( i=0; i<num; i++) printf("%i %x \n",i,buffer[i]);//??? |
| 58 |
|
| 59 |
/* /\*** Checksum ***\/ */ |
| 60 |
/* oi = buffer; */ |
| 61 |
/* word = 0; */ |
| 62 |
/* for( i=0; i<l_buffer; ++i) */ |
| 63 |
/* word = (word + *oi++) & 0x03ff; */ |
| 64 |
/* if( word == ( header[6] & 0x03ff ) ) */ |
| 65 |
/* /\* printf("Check Sum OK \n")*\/ ; */ |
| 66 |
/* else */ |
| 67 |
/* printf("Wrong Check Sum: %x, %x \n",word,header[6]&0x03ff); */ |
| 68 |
|
| 69 |
} |
| 70 |
|
| 71 |
/*========================================================= */ |
| 72 |
/*========================================================= */ |
| 73 |
/* |
| 74 |
* Routine to unpacking buffer file |
| 75 |
* |
| 76 |
*/ |
| 77 |
/*========================================================= */ |
| 78 |
/*========================================================= */ |
| 79 |
void unpacking(int nwords, unsigned short *output) { |
| 80 |
unsigned short out[10000]; |
| 81 |
int k; |
| 82 |
unsigned short *oi, *oo, *o1; |
| 83 |
|
| 84 |
for( k=0; k<10000; ++k) |
| 85 |
out[k]=0; |
| 86 |
oo=output; |
| 87 |
oi=out; |
| 88 |
for( k=0; k<nwords; ++k) { |
| 89 |
*oi++ = *oo; |
| 90 |
*oo++ = 0; |
| 91 |
} |
| 92 |
|
| 93 |
oo=output; |
| 94 |
oi=out; |
| 95 |
|
| 96 |
for(k=0; k<nwords; k+=13) { |
| 97 |
o1 = oi++; |
| 98 |
*oo++ = (*o1 & 0xfff8) >> 3; |
| 99 |
*oo++ = ( (*o1 & 0x0007) << 10 ) | ( (*oi & 0xffc0) >> 6); |
| 100 |
o1 = oi++; |
| 101 |
*oo++ = ( (*o1 & 0x003f) << 7 ) | ( (*oi & 0xfe00) >> 9); |
| 102 |
o1 = oi++; |
| 103 |
*oo++ = ( (*o1 & 0x01ff) << 4 ) | ( (*oi & 0xf000) >> 12); |
| 104 |
o1 = oi++; |
| 105 |
*oo++ = ( (*o1 & 0x0fff) << 1 ) | ( (*oi & 0x8000) >> 15); |
| 106 |
*oo++ = (*oi & 0x7ffc) >> 2; |
| 107 |
o1 = oi++; |
| 108 |
*oo++ = ( (*o1 & 0x0003) << 11 ) | ( (*oi & 0xffe0) >> 5); |
| 109 |
o1 = oi++; |
| 110 |
*oo++ = ( (*o1 & 0x001f) << 8 ) | ( (*oi & 0xff00) >> 8); |
| 111 |
o1 = oi++; |
| 112 |
*oo++ = ( (*o1 & 0x00ff) << 5 ) | ( (*oi & 0xf800) >> 11); |
| 113 |
o1 = oi++; |
| 114 |
*oo++ = ( (*o1 & 0x07ff) << 2 ) | ( (*oi & 0xc000) >> 14); |
| 115 |
*oo++ = (*oi & 0x3ffe) >> 1; |
| 116 |
o1 = oi++; |
| 117 |
*oo++ = ( (*o1 & 0x0001) << 12 ) | ( (*oi & 0xfff0) >> 4); |
| 118 |
o1 = oi++; |
| 119 |
*oo++ = ( (*o1 & 0x000f) << 9 ) | ( (*oi & 0xff80) >> 7); |
| 120 |
o1 = oi++; |
| 121 |
*oo++ = ( (*o1 & 0x007f) << 6 ) | ( (*oi & 0xfc00) >> 10); |
| 122 |
o1 = oi++; |
| 123 |
*oo++ = ( (*o1 & 0x03ff) << 3 ) | ( (*oi & 0xe000) >> 13); |
| 124 |
*oo++ = *oi++ & 0x1fff; |
| 125 |
} |
| 126 |
|
| 127 |
|
| 128 |
return; |
| 129 |
} |