Coverage Report

Created: 2026-07-30 06:41

next uncovered line (L), next uncovered region (R), next uncovered branch (B)
/src/libjpeg-turbo.main/src/spng/zlib/inffast.c
Line
Count
Source
1
/* inffast.c -- fast decoding
2
 * Copyright (C) 1995-2026 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
404k
void ZLIB_INTERNAL inflate_fast(z_streamp strm, unsigned start) {
51
404k
    struct inflate_state FAR *state;
52
404k
    z_const unsigned char FAR *in;      /* local strm->next_in */
53
404k
    z_const unsigned char FAR *last;    /* have enough input while in < last */
54
404k
    unsigned char FAR *out;     /* local strm->next_out */
55
404k
    unsigned char FAR *beg;     /* inflate()'s initial strm->next_out */
56
404k
    unsigned char FAR *end;     /* while out < end, enough space available */
57
#ifdef INFLATE_STRICT
58
    unsigned dmax;              /* maximum distance from zlib header */
59
#endif
60
404k
    unsigned wsize;             /* window size or zero if not using window */
61
404k
    unsigned whave;             /* valid bytes in the window */
62
404k
    unsigned wnext;             /* window write index */
63
404k
    unsigned char FAR *window;  /* allocated sliding window, if wsize != 0 */
64
404k
    unsigned long hold;         /* local strm->hold */
65
404k
    unsigned bits;              /* local strm->bits */
66
404k
    code const FAR *lcode;      /* local strm->lencode */
67
404k
    code const FAR *dcode;      /* local strm->distcode */
68
404k
    unsigned lmask;             /* mask for first level of length codes */
69
404k
    unsigned dmask;             /* mask for first level of distance codes */
70
404k
    code const *here;           /* retrieved table entry */
71
404k
    unsigned op;                /* code bits, operation, extra bits, or */
72
                                /*  window position, window bytes to copy */
73
404k
    unsigned len;               /* match length, unused bytes */
74
404k
    unsigned dist;              /* match distance */
75
404k
    unsigned char FAR *from;    /* where to copy match from */
76
77
    /* copy state to local variables */
78
404k
    state = (struct inflate_state FAR *)strm->state;
79
404k
    in = strm->next_in;
80
404k
    last = in + (strm->avail_in - 5);
81
404k
    out = strm->next_out;
82
404k
    beg = out - (start - strm->avail_out);
83
404k
    end = out + (strm->avail_out - 257);
84
#ifdef INFLATE_STRICT
85
    dmax = state->dmax;
86
#endif
87
404k
    wsize = state->wsize;
88
404k
    whave = state->whave;
89
404k
    wnext = state->wnext;
90
404k
    window = state->window;
91
404k
    hold = state->hold;
92
404k
    bits = state->bits;
93
404k
    lcode = state->lencode;
94
404k
    dcode = state->distcode;
95
404k
    lmask = (1U << state->lenbits) - 1;
96
404k
    dmask = (1U << state->distbits) - 1;
97
98
    /* decode literals and length/distances until end-of-block or not enough
99
       input data or output space */
100
82.6M
    do {
101
82.6M
        if (bits < 15) {
102
23.5M
            hold += (unsigned long)(*in++) << bits;
103
23.5M
            bits += 8;
104
23.5M
            hold += (unsigned long)(*in++) << bits;
105
23.5M
            bits += 8;
106
23.5M
        }
107
82.6M
        here = lcode + (hold & lmask);
108
89.5M
      dolen:
109
89.5M
        op = (unsigned)(here->bits);
110
89.5M
        hold >>= op;
111
89.5M
        bits -= op;
112
89.5M
        op = (unsigned)(here->op);
113
89.5M
        if (op == 0) {                          /* literal */
114
23.8M
            Tracevv((stderr, here->val >= 0x20 && here->val < 0x7f ?
115
23.8M
                    "inflate:         literal '%c'\n" :
116
23.8M
                    "inflate:         literal 0x%02x\n", here->val));
117
23.8M
            *out++ = (unsigned char)(here->val);
118
23.8M
        }
119
65.6M
        else if (op & 16) {                     /* length base */
120
58.7M
            len = (unsigned)(here->val);
121
58.7M
            op &= 15;                           /* number of extra bits */
122
58.7M
            if (op) {
123
20.4M
                if (bits < op) {
124
2.77k
                    hold += (unsigned long)(*in++) << bits;
125
2.77k
                    bits += 8;
126
2.77k
                }
127
20.4M
                len += (unsigned)hold & ((1U << op) - 1);
128
20.4M
                hold >>= op;
129
20.4M
                bits -= op;
130
20.4M
            }
131
58.7M
            Tracevv((stderr, "inflate:         length %u\n", len));
132
58.7M
            if (bits < 15) {
133
19.8M
                hold += (unsigned long)(*in++) << bits;
134
19.8M
                bits += 8;
135
19.8M
                hold += (unsigned long)(*in++) << bits;
136
19.8M
                bits += 8;
137
19.8M
            }
138
58.7M
            here = dcode + (hold & dmask);
139
66.5M
          dodist:
140
66.5M
            op = (unsigned)(here->bits);
141
66.5M
            hold >>= op;
142
66.5M
            bits -= op;
143
66.5M
            op = (unsigned)(here->op);
144
66.5M
            if (op & 16) {                      /* distance base */
145
58.7M
                dist = (unsigned)(here->val);
146
58.7M
                op &= 15;                       /* number of extra bits */
147
58.7M
                if (bits < op) {
148
12.5k
                    hold += (unsigned long)(*in++) << bits;
149
12.5k
                    bits += 8;
150
12.5k
                    if (bits < op) {
151
308
                        hold += (unsigned long)(*in++) << bits;
152
308
                        bits += 8;
153
308
                    }
154
12.5k
                }
155
58.7M
                dist += (unsigned)hold & ((1U << op) - 1);
156
#ifdef INFLATE_STRICT
157
                if (dist > dmax) {
158
                    strm->msg = (z_const char *)
159
                        "invalid distance too far back";
160
                    state->mode = BAD;
161
                    break;
162
                }
163
#endif
164
58.7M
                hold >>= op;
165
58.7M
                bits -= op;
166
58.7M
                Tracevv((stderr, "inflate:         distance %u\n", dist));
167
58.7M
                op = (unsigned)(out - beg);     /* max distance in output */
168
58.7M
                if (dist > op) {                /* see if copy from window */
169
210k
                    op = dist - op;             /* distance back in window */
170
210k
                    if (op > whave) {
171
954
                        if (state->sane) {
172
954
                            strm->msg = (z_const char *)
173
954
                                "invalid distance too far back";
174
954
                            state->mode = BAD;
175
954
                            break;
176
954
                        }
177
#ifdef INFLATE_ALLOW_INVALID_DISTANCE_TOOFAR_ARRR
178
                        if (len <= op - whave) {
179
                            do {
180
                                *out++ = 0;
181
                            } while (--len);
182
                            continue;
183
                        }
184
                        len -= op - whave;
185
                        do {
186
                            *out++ = 0;
187
                        } while (--op > whave);
188
                        if (op == 0) {
189
                            from = out - dist;
190
                            do {
191
                                *out++ = *from++;
192
                            } while (--len);
193
                            continue;
194
                        }
195
#endif
196
954
                    }
197
209k
                    from = window;
198
209k
                    if (wnext == 0) {           /* very common case */
199
82.5k
                        from += wsize - op;
200
82.5k
                        if (op < len) {         /* some from window */
201
6.53k
                            len -= op;
202
393k
                            do {
203
393k
                                *out++ = *from++;
204
393k
                            } while (--op);
205
6.53k
                            from = out - dist;  /* rest from output */
206
6.53k
                        }
207
82.5k
                    }
208
127k
                    else if (wnext < op) {      /* wrap around window */
209
17.4k
                        from += wsize + wnext - op;
210
17.4k
                        op -= wnext;
211
17.4k
                        if (op < len) {         /* some from end of window */
212
395
                            len -= op;
213
41.5k
                            do {
214
41.5k
                                *out++ = *from++;
215
41.5k
                            } while (--op);
216
395
                            from = window;
217
395
                            if (wnext < len) {  /* some from start of window */
218
155
                                op = wnext;
219
155
                                len -= op;
220
9.02k
                                do {
221
9.02k
                                    *out++ = *from++;
222
9.02k
                                } while (--op);
223
155
                                from = out - dist;      /* rest from output */
224
155
                            }
225
395
                        }
226
17.4k
                    }
227
109k
                    else {                      /* contiguous in window */
228
109k
                        from += wnext - op;
229
109k
                        if (op < len) {         /* some from window */
230
8.45k
                            len -= op;
231
820k
                            do {
232
820k
                                *out++ = *from++;
233
820k
                            } while (--op);
234
8.45k
                            from = out - dist;  /* rest from output */
235
8.45k
                        }
236
109k
                    }
237
5.76M
                    while (len > 2) {
238
5.56M
                        *out++ = *from++;
239
5.56M
                        *out++ = *from++;
240
5.56M
                        *out++ = *from++;
241
5.56M
                        len -= 3;
242
5.56M
                    }
243
209k
                    if (len) {
244
131k
                        *out++ = *from++;
245
131k
                        if (len > 1)
246
103k
                            *out++ = *from++;
247
131k
                    }
248
209k
                }
249
58.5M
                else {
250
58.5M
                    from = out - dist;          /* copy direct from output */
251
3.51G
                    do {                        /* minimum length is three */
252
3.51G
                        *out++ = *from++;
253
3.51G
                        *out++ = *from++;
254
3.51G
                        *out++ = *from++;
255
3.51G
                        len -= 3;
256
3.51G
                    } while (len > 2);
257
58.5M
                    if (len) {
258
20.2M
                        *out++ = *from++;
259
20.2M
                        if (len > 1)
260
20.1M
                            *out++ = *from++;
261
20.2M
                    }
262
58.5M
                }
263
58.7M
            }
264
7.83M
            else if ((op & 64) == 0) {          /* 2nd level distance code */
265
7.83M
                here = dcode + here->val + (hold & ((1U << op) - 1));
266
7.83M
                goto dodist;
267
7.83M
            }
268
280
            else {
269
280
                strm->msg = (z_const char *)"invalid distance code";
270
280
                state->mode = BAD;
271
280
                break;
272
280
            }
273
66.5M
        }
274
6.93M
        else if ((op & 64) == 0) {              /* 2nd level length code */
275
6.84M
            here = lcode + here->val + (hold & ((1U << op) - 1));
276
6.84M
            goto dolen;
277
6.84M
        }
278
92.3k
        else if (op & 32) {                     /* end-of-block */
279
92.2k
            Tracevv((stderr, "inflate:         end of block\n"));
280
92.2k
            state->mode = TYPE;
281
92.2k
            break;
282
92.2k
        }
283
125
        else {
284
125
            strm->msg = (z_const char *)"invalid literal/length code";
285
125
            state->mode = BAD;
286
125
            break;
287
125
        }
288
89.5M
    } while (in < last && out < end);
289
290
    /* return unused bytes (on entry, bits < 8, so in won't go too far back) */
291
404k
    len = bits >> 3;
292
404k
    in -= len;
293
404k
    bits -= len << 3;
294
404k
    hold &= (1U << bits) - 1;
295
296
    /* update state and return */
297
404k
    strm->next_in = in;
298
404k
    strm->next_out = out;
299
404k
    strm->avail_in = (unsigned)(in < last ? 5 + (last - in) : 5 - (in - last));
300
404k
    strm->avail_out = (unsigned)(out < end ?
301
211k
                                 257 + (end - out) : 257 - (out - end));
302
404k
    state->hold = hold;
303
404k
    state->bits = bits;
304
404k
    return;
305
404k
}
306
307
/*
308
   inflate_fast() speedups that turned out slower (on a PowerPC G3 750CXe):
309
   - Using bit fields for code structure
310
   - Different op definition to avoid & for extra bits (do & for table bits)
311
   - Three separate decoding do-loops for direct, window, and wnext == 0
312
   - Special case for distance > 1 copies to do overlapped load and store copy
313
   - Explicit branch predictions (based on measured branch probabilities)
314
   - Deferring match copy and interspersed it with decoding subsequent codes
315
   - Swapping literal/length else
316
   - Swapping window/direct else
317
   - Larger unrolled copy loops (three is about right)
318
   - Moving len -= 3 statement into middle of loop
319
 */
320
321
#endif /* !ASMINF */