/proc/self/cwd/external/zlib/inffast.c
Line | Count | Source |
1 | | /* inffast.c -- fast decoding |
2 | | * Copyright (C) 1995-2017 Mark Adler |
3 | | * For conditions of distribution and use, see copyright notice in zlib.h |
4 | | */ |
5 | | |
6 | | #include "zutil.h" |
7 | | #include "inftrees.h" |
8 | | #include "inflate.h" |
9 | | #include "inffast.h" |
10 | | |
11 | | #ifdef ASMINF |
12 | | # pragma message("Assembler code may have bugs -- use at your own risk") |
13 | | #else |
14 | | |
15 | | /* |
16 | | Decode literal, length, and distance codes and write out the resulting |
17 | | literal and match bytes until either not enough input or output is |
18 | | available, an end-of-block is encountered, or a data error is encountered. |
19 | | When large enough input and output buffers are supplied to inflate(), for |
20 | | example, a 16K input buffer and a 64K output buffer, more than 95% of the |
21 | | inflate execution time is spent in this routine. |
22 | | |
23 | | Entry assumptions: |
24 | | |
25 | | state->mode == LEN |
26 | | strm->avail_in >= 6 |
27 | | strm->avail_out >= 258 |
28 | | start >= strm->avail_out |
29 | | state->bits < 8 |
30 | | |
31 | | On return, state->mode is one of: |
32 | | |
33 | | LEN -- ran out of enough output space or enough available input |
34 | | TYPE -- reached end of block code, inflate() to interpret next block |
35 | | BAD -- error in block data |
36 | | |
37 | | Notes: |
38 | | |
39 | | - The maximum input bits used by a length/distance pair is 15 bits for the |
40 | | length code, 5 bits for the length extra, 15 bits for the distance code, |
41 | | and 13 bits for the distance extra. This totals 48 bits, or six bytes. |
42 | | Therefore if strm->avail_in >= 6, then there is enough input to avoid |
43 | | checking for available input while decoding. |
44 | | |
45 | | - The maximum bytes that a single length/distance pair can output is 258 |
46 | | bytes, which is the maximum length that can be coded. inflate_fast() |
47 | | requires strm->avail_out >= 258 for each loop to avoid checking for |
48 | | output space. |
49 | | */ |
50 | | void ZLIB_INTERNAL inflate_fast(strm, start) |
51 | | z_streamp strm; |
52 | | unsigned start; /* inflate()'s starting value for strm->avail_out */ |
53 | 129k | { |
54 | 129k | struct inflate_state FAR *state; |
55 | 129k | z_const unsigned char FAR *in; /* local strm->next_in */ |
56 | 129k | z_const unsigned char FAR *last; /* have enough input while in < last */ |
57 | 129k | unsigned char FAR *out; /* local strm->next_out */ |
58 | 129k | unsigned char FAR *beg; /* inflate()'s initial strm->next_out */ |
59 | 129k | unsigned char FAR *end; /* while out < end, enough space available */ |
60 | | #ifdef INFLATE_STRICT |
61 | | unsigned dmax; /* maximum distance from zlib header */ |
62 | | #endif |
63 | 129k | unsigned wsize; /* window size or zero if not using window */ |
64 | 129k | unsigned whave; /* valid bytes in the window */ |
65 | 129k | unsigned wnext; /* window write index */ |
66 | 129k | unsigned char FAR *window; /* allocated sliding window, if wsize != 0 */ |
67 | 129k | unsigned long hold; /* local strm->hold */ |
68 | 129k | unsigned bits; /* local strm->bits */ |
69 | 129k | code const FAR *lcode; /* local strm->lencode */ |
70 | 129k | code const FAR *dcode; /* local strm->distcode */ |
71 | 129k | unsigned lmask; /* mask for first level of length codes */ |
72 | 129k | unsigned dmask; /* mask for first level of distance codes */ |
73 | 129k | code const *here; /* retrieved table entry */ |
74 | 129k | unsigned op; /* code bits, operation, extra bits, or */ |
75 | | /* window position, window bytes to copy */ |
76 | 129k | unsigned len; /* match length, unused bytes */ |
77 | 129k | unsigned dist; /* match distance */ |
78 | 129k | unsigned char FAR *from; /* where to copy match from */ |
79 | | |
80 | | /* copy state to local variables */ |
81 | 129k | state = (struct inflate_state FAR *)strm->state; |
82 | 129k | in = strm->next_in; |
83 | 129k | last = in + (strm->avail_in - 5); |
84 | 129k | out = strm->next_out; |
85 | 129k | beg = out - (start - strm->avail_out); |
86 | 129k | end = out + (strm->avail_out - 257); |
87 | | #ifdef INFLATE_STRICT |
88 | | dmax = state->dmax; |
89 | | #endif |
90 | 129k | wsize = state->wsize; |
91 | 129k | whave = state->whave; |
92 | 129k | wnext = state->wnext; |
93 | 129k | window = state->window; |
94 | 129k | hold = state->hold; |
95 | 129k | bits = state->bits; |
96 | 129k | lcode = state->lencode; |
97 | 129k | dcode = state->distcode; |
98 | 129k | lmask = (1U << state->lenbits) - 1; |
99 | 129k | dmask = (1U << state->distbits) - 1; |
100 | | |
101 | | /* decode literals and length/distances until end-of-block or not enough |
102 | | input data or output space */ |
103 | 8.77M | do { |
104 | 8.77M | if (bits < 15) { |
105 | 3.31M | hold += (unsigned long)(*in++) << bits; |
106 | 3.31M | bits += 8; |
107 | 3.31M | hold += (unsigned long)(*in++) << bits; |
108 | 3.31M | bits += 8; |
109 | 3.31M | } |
110 | 8.77M | here = lcode + (hold & lmask); |
111 | 9.40M | dolen: |
112 | 9.40M | op = (unsigned)(here->bits); |
113 | 9.40M | hold >>= op; |
114 | 9.40M | bits -= op; |
115 | 9.40M | op = (unsigned)(here->op); |
116 | 9.40M | if (op == 0) { /* literal */ |
117 | 6.83M | Tracevv((stderr, here->val >= 0x20 && here->val < 0x7f ? |
118 | 6.83M | "inflate: literal '%c'\n" : |
119 | 6.83M | "inflate: literal 0x%02x\n", here->val)); |
120 | 6.83M | *out++ = (unsigned char)(here->val); |
121 | 6.83M | } |
122 | 2.56M | else if (op & 16) { /* length base */ |
123 | 1.93M | len = (unsigned)(here->val); |
124 | 1.93M | op &= 15; /* number of extra bits */ |
125 | 1.93M | if (op) { |
126 | 826k | if (bits < op) { |
127 | 1.91k | hold += (unsigned long)(*in++) << bits; |
128 | 1.91k | bits += 8; |
129 | 1.91k | } |
130 | 826k | len += (unsigned)hold & ((1U << op) - 1); |
131 | 826k | hold >>= op; |
132 | 826k | bits -= op; |
133 | 826k | } |
134 | 1.93M | Tracevv((stderr, "inflate: length %u\n", len)); |
135 | 1.93M | if (bits < 15) { |
136 | 1.17M | hold += (unsigned long)(*in++) << bits; |
137 | 1.17M | bits += 8; |
138 | 1.17M | hold += (unsigned long)(*in++) << bits; |
139 | 1.17M | bits += 8; |
140 | 1.17M | } |
141 | 1.93M | here = dcode + (hold & dmask); |
142 | 2.51M | dodist: |
143 | 2.51M | op = (unsigned)(here->bits); |
144 | 2.51M | hold >>= op; |
145 | 2.51M | bits -= op; |
146 | 2.51M | op = (unsigned)(here->op); |
147 | 2.51M | if (op & 16) { /* distance base */ |
148 | 1.93M | dist = (unsigned)(here->val); |
149 | 1.93M | op &= 15; /* number of extra bits */ |
150 | 1.93M | if (bits < op) { |
151 | 18.5k | hold += (unsigned long)(*in++) << bits; |
152 | 18.5k | bits += 8; |
153 | 18.5k | if (bits < op) { |
154 | 106 | hold += (unsigned long)(*in++) << bits; |
155 | 106 | bits += 8; |
156 | 106 | } |
157 | 18.5k | } |
158 | 1.93M | dist += (unsigned)hold & ((1U << op) - 1); |
159 | | #ifdef INFLATE_STRICT |
160 | | if (dist > dmax) { |
161 | | strm->msg = (char *)"invalid distance too far back"; |
162 | | state->mode = BAD; |
163 | | break; |
164 | | } |
165 | | #endif |
166 | 1.93M | hold >>= op; |
167 | 1.93M | bits -= op; |
168 | 1.93M | Tracevv((stderr, "inflate: distance %u\n", dist)); |
169 | 1.93M | op = (unsigned)(out - beg); /* max distance in output */ |
170 | 1.93M | if (dist > op) { /* see if copy from window */ |
171 | 203k | op = dist - op; /* distance back in window */ |
172 | 203k | if (op > whave) { |
173 | 106 | if (state->sane) { |
174 | 106 | strm->msg = |
175 | 106 | (char *)"invalid distance too far back"; |
176 | 106 | state->mode = BAD; |
177 | 106 | break; |
178 | 106 | } |
179 | | #ifdef INFLATE_ALLOW_INVALID_DISTANCE_TOOFAR_ARRR |
180 | | if (len <= op - whave) { |
181 | | do { |
182 | | *out++ = 0; |
183 | | } while (--len); |
184 | | continue; |
185 | | } |
186 | | len -= op - whave; |
187 | | do { |
188 | | *out++ = 0; |
189 | | } while (--op > whave); |
190 | | if (op == 0) { |
191 | | from = out - dist; |
192 | | do { |
193 | | *out++ = *from++; |
194 | | } while (--len); |
195 | | continue; |
196 | | } |
197 | | #endif |
198 | 106 | } |
199 | 203k | from = window; |
200 | 203k | if (wnext == 0) { /* very common case */ |
201 | 13.1k | from += wsize - op; |
202 | 13.1k | if (op < len) { /* some from window */ |
203 | 753 | len -= op; |
204 | 12.6k | do { |
205 | 12.6k | *out++ = *from++; |
206 | 12.6k | } while (--op); |
207 | 753 | from = out - dist; /* rest from output */ |
208 | 753 | } |
209 | 13.1k | } |
210 | 190k | else if (wnext < op) { /* wrap around window */ |
211 | 56.6k | from += wsize + wnext - op; |
212 | 56.6k | op -= wnext; |
213 | 56.6k | if (op < len) { /* some from end of window */ |
214 | 1.26k | len -= op; |
215 | 66.4k | do { |
216 | 66.4k | *out++ = *from++; |
217 | 66.4k | } while (--op); |
218 | 1.26k | from = window; |
219 | 1.26k | if (wnext < len) { /* some from start of window */ |
220 | 359 | op = wnext; |
221 | 359 | len -= op; |
222 | 14.6k | do { |
223 | 14.6k | *out++ = *from++; |
224 | 14.6k | } while (--op); |
225 | 359 | from = out - dist; /* rest from output */ |
226 | 359 | } |
227 | 1.26k | } |
228 | 56.6k | } |
229 | 133k | else { /* contiguous in window */ |
230 | 133k | from += wnext - op; |
231 | 133k | if (op < len) { /* some from window */ |
232 | 9.98k | len -= op; |
233 | 787k | do { |
234 | 787k | *out++ = *from++; |
235 | 787k | } while (--op); |
236 | 9.98k | from = out - dist; /* rest from output */ |
237 | 9.98k | } |
238 | 133k | } |
239 | 2.21M | while (len > 2) { |
240 | 2.00M | *out++ = *from++; |
241 | 2.00M | *out++ = *from++; |
242 | 2.00M | *out++ = *from++; |
243 | 2.00M | len -= 3; |
244 | 2.00M | } |
245 | 203k | if (len) { |
246 | 94.0k | *out++ = *from++; |
247 | 94.0k | if (len > 1) |
248 | 61.0k | *out++ = *from++; |
249 | 94.0k | } |
250 | 203k | } |
251 | 1.72M | else { |
252 | 1.72M | from = out - dist; /* copy direct from output */ |
253 | 47.3M | do { /* minimum length is three */ |
254 | 47.3M | *out++ = *from++; |
255 | 47.3M | *out++ = *from++; |
256 | 47.3M | *out++ = *from++; |
257 | 47.3M | len -= 3; |
258 | 47.3M | } while (len > 2); |
259 | 1.72M | if (len) { |
260 | 1.22M | *out++ = *from++; |
261 | 1.22M | if (len > 1) |
262 | 1.00M | *out++ = *from++; |
263 | 1.22M | } |
264 | 1.72M | } |
265 | 1.93M | } |
266 | 579k | else if ((op & 64) == 0) { /* 2nd level distance code */ |
267 | 579k | here = dcode + here->val + (hold & ((1U << op) - 1)); |
268 | 579k | goto dodist; |
269 | 579k | } |
270 | 4 | else { |
271 | 4 | strm->msg = (char *)"invalid distance code"; |
272 | 4 | state->mode = BAD; |
273 | 4 | break; |
274 | 4 | } |
275 | 2.51M | } |
276 | 637k | else if ((op & 64) == 0) { /* 2nd level length code */ |
277 | 628k | here = lcode + here->val + (hold & ((1U << op) - 1)); |
278 | 628k | goto dolen; |
279 | 628k | } |
280 | 9.03k | else if (op & 32) { /* end-of-block */ |
281 | 9.02k | Tracevv((stderr, "inflate: end of block\n")); |
282 | 9.02k | state->mode = TYPE; |
283 | 9.02k | break; |
284 | 9.02k | } |
285 | 6 | else { |
286 | 6 | strm->msg = (char *)"invalid literal/length code"; |
287 | 6 | state->mode = BAD; |
288 | 6 | break; |
289 | 6 | } |
290 | 9.40M | } while (in < last && out < end); |
291 | | |
292 | | /* return unused bytes (on entry, bits < 8, so in won't go too far back) */ |
293 | 129k | len = bits >> 3; |
294 | 129k | in -= len; |
295 | 129k | bits -= len << 3; |
296 | 129k | hold &= (1U << bits) - 1; |
297 | | |
298 | | /* update state and return */ |
299 | 129k | strm->next_in = in; |
300 | 129k | strm->next_out = out; |
301 | 129k | strm->avail_in = (unsigned)(in < last ? 5 + (last - in) : 5 - (in - last)); |
302 | 129k | strm->avail_out = (unsigned)(out < end ? |
303 | 101k | 257 + (end - out) : 257 - (out - end)); |
304 | 129k | state->hold = hold; |
305 | 129k | state->bits = bits; |
306 | 129k | return; |
307 | 129k | } |
308 | | |
309 | | /* |
310 | | inflate_fast() speedups that turned out slower (on a PowerPC G3 750CXe): |
311 | | - Using bit fields for code structure |
312 | | - Different op definition to avoid & for extra bits (do & for table bits) |
313 | | - Three separate decoding do-loops for direct, window, and wnext == 0 |
314 | | - Special case for distance > 1 copies to do overlapped load and store copy |
315 | | - Explicit branch predictions (based on measured branch probabilities) |
316 | | - Deferring match copy and interspersed it with decoding subsequent codes |
317 | | - Swapping literal/length else |
318 | | - Swapping window/direct else |
319 | | - Larger unrolled copy loops (three is about right) |
320 | | - Moving len -= 3 statement into middle of loop |
321 | | */ |
322 | | |
323 | | #endif /* !ASMINF */ |