source: proto/parabix2/bitutil.py @ 486

Last change on this file since 486 was 371, checked in by cameron, 9 years ago

Fix EOF_mask (went 1 position too far)

File size: 4.3 KB
Line 
1#
2# bitutil.py
3#
4# Bitstream Utilities
5#
6# These are quick-and-dirty Python implementations of utilities
7# for demo purposes only.
8#
9# Robert D. Cameron
10# August 19, 2009
11#
12#----------------------------------------------------------------------------
13#
14# We use python's unlimited precision integers for unbounded bit streams.
15# This permits simple logical operations on the entire stream.
16# Assumption: bitstreams are little-endian (e.g., as on x86).
17#
18#----------------------------------------------------------------------------
19#
20import sys
21import codecs
22# Utility functions for demo purposes.
23#
24
25def writefile(filename, s, charcode='utf-8',):
26        f = codecs.open(filename, 'w', encoding=charcode)
27        f.write(s) # writes encoded bytes to file
28        f.close()
29
30def readfile(filename):
31        f = open(filename, 'r')
32        contents = f.read() 
33        f.close()
34        return contents
35
36def count_leading_zeroes(strm):
37        zeroes = 0
38        while (strm & 0xFFFFFFFF) == 0: 
39                zeroes += 32
40                strm >>= 32
41        while (strm & 1) == 0:
42                zeroes += 1
43                strm >>= 1
44        return zeroes
45
46def EOF_mask(strm_lgth):
47        mask = 1
48        mask <<= (strm_lgth)
49        mask -= 1
50        return mask
51
52def extract_bit(strm, pos):
53        bit = (strm >> pos) & 1
54        return bit
55
56def transpose_streams(s):
57        b = []
58        mask = 128
59        index = 0
60        while index < 8:
61                current = 0
62                cursor = 1
63                for byte in s:
64                        if  (ord(byte) & mask != 0):
65                                current += cursor
66                        cursor <<= 1
67                       
68                index+=1
69                mask>>=1
70                b.append(current)
71        return (b, cursor-1)  # basis streams and EOF mask
72
73def inverse_transpose(bitset, len):
74        bytestream=""
75        cursor = 1
76        for i in range(0, len):
77                byteval = 0
78                for j in range(0,8):
79                        if bitset[j] & cursor != 0:
80                                byteval += 128 >> j
81                bytestream += chr(byteval)
82                cursor += cursor
83        return bytestream
84                       
85def filter_bits(bitstream, delmask):
86        newstream = 0
87        cursor = 1
88        while delmask > 0:
89                if delmask & 1 == 0:
90                        if bitstream & 1 == 1:
91                                newstream += cursor
92                        cursor += cursor
93                delmask >>= 1
94                bitstream >>= 1
95        while bitstream > 0:
96                if bitstream & 1 == 1:
97                        newstream += cursor
98                cursor += cursor
99                bitstream >>= 1         
100        return newstream
101                       
102def filter_bytes(bytestream, delmask):
103        newstream=""
104        cursor = 1
105        for c in bytestream:
106                if delmask & cursor == 0:
107                        newstream += c
108                cursor += cursor
109        return newstream
110                       
111def merge_bytes(stream1, stream2):
112        s = ""
113        for i in range(len(stream1)):
114                s += stream1[i]
115                s += stream2[i]
116        return s
117
118def bitstream2string(stream, lgth):
119        str = ""
120        for i in range(lgth):
121                if stream & 1 == 1: str += '1'
122                else: str += '_'
123                stream >>= 1
124        return str
125
126def print_aligned_streams(stream_list):
127        """Print out a set of aligned streams."""
128        label_max = max([len(p[0]) for p in stream_list])
129        for p in stream_list:
130                print (p[0] + " "*(label_max - len(p[0]))) + ": " + p[1] 
131
132def print_aligned_u8_byte_streams(u8_byte_stream):
133        """Print out a set of 'encoding' aligned streams."""
134
135        # Set the system info to print utf-8
136        info = codecs.lookup('utf-8')
137        sys.stdout = info.streamwriter(sys.stdout) 
138
139        label_max = max([len(p[0]) for p in u8_byte_stream])
140        for p in u8_byte_stream:
141                sys.stdout.write(p[0] + " "*(label_max - len(p[0])) + ": ")
142               
143                for c in (p[1].decode('utf-8')):                        # for each unicode character
144                        u8_seq_len = len(c.encode('utf-8'))     # encode the unicode character as utf-8 and get the utf-8 sequence length
145                        if u8_seq_len == 1:
146                                sys.stdout.write(c)
147                        elif u8_seq_len == 2:
148                                sys.stdout.write('2_')                          # align 2 byte sequences with a trailing 2
149                        elif u8_seq_len == 3:
150                                sys.stdout.write('3__')                         # align 3 byte sequences with 2 trailing 3's
151                        elif u8_seq_len == 4:
152                                sys.stdout.write('4___')                        # align 2 byte sequences with 3 trailing 4's             
153                        else:
154                                sys.stdout.write('Error: Invalid UTF-8 character sequence of length: ' + str(u8_seq_len) + '.')                 
155                sys.stdout.write('\n')   
156
157def high_nybble_stream(bytes):
158        str=""
159        for b in bytes:
160                h = hex(ord(b))[-2]
161                if h == 'x': str += '0'
162                else: str += h
163        return str
164
165def low_nybble_stream(bytes):
166        str=""
167        for b in bytes:
168                str += hex(ord(b))[-1]
169        return str
170
171
172       
173#
174#  ScanThru is the core operation for parallel scanning.
175#  Given a bitstream Cursors marking current cursor positions,
176#  advance each cursor through occurrences of a character
177#  class identified by ScanStream.
178
179def ScanThru(Cursors, ScanStream):
180        return (Cursors + ScanStream) & ~ScanStream
181
182#
183# Advance all cursors by one position.
184def Advance(stream):
185        return stream + stream
186
Note: See TracBrowser for help on using the repository browser.