Coverage Report

Created: 2025-06-10 07:27

/src/ghostpdl/base/scfd.c
Line
Count
Source (jump to first uncovered line)
1
/* Copyright (C) 2001-2023 Artifex Software, Inc.
2
   All Rights Reserved.
3
4
   This software is provided AS-IS with no warranty, either express or
5
   implied.
6
7
   This software is distributed under license and may not be copied,
8
   modified or distributed except as expressly authorized under the terms
9
   of the license contained in the file LICENSE in this distribution.
10
11
   Refer to licensing information at http://www.artifex.com or contact
12
   Artifex Software, Inc.,  39 Mesa Street, Suite 108A, San Francisco,
13
   CA 94129, USA, for further information.
14
*/
15
16
17
/* CCITTFax decoding filter */
18
#include "stdio_.h"   /* includes std.h */
19
#include "memory_.h"
20
#include "gdebug.h"
21
#include "strimpl.h"
22
#include "scf.h"
23
#include "scfx.h"
24
25
/* ------ CCITTFaxDecode ------ */
26
27
private_st_CFD_state();
28
29
14.8M
#define CFD_BUFFER_SLOP 4
30
31
static inline int
32
get_run(stream_CFD_state *ss, stream_cursor_read *pr, const cfd_node decode[],
33
    int initial_bits, int min_bits, int *runlen, const char *str);
34
35
static inline int
36
invert_data(stream_CFD_state *ss, stream_cursor_read *pr, int *rlen, byte black_byte);
37
38
static inline int
39
skip_data(stream_CFD_state *ss, stream_cursor_read *pr, int rlen);
40
41
/* Set default parameter values. */
42
static void
43
s_CFD_set_defaults(register stream_state * st)
44
253k
{
45
253k
    stream_CFD_state *const ss = (stream_CFD_state *) st;
46
47
253k
    s_CFD_set_defaults_inline(ss);
48
253k
}
49
50
/* Initialize CCITTFaxDecode filter */
51
static int
52
s_CFD_init(stream_state * st)
53
253k
{
54
253k
    stream_CFD_state *const ss = (stream_CFD_state *) st;
55
253k
    int raster = ss->raster =
56
253k
        ROUND_UP((ss->Columns + 7) >> 3, ss->DecodedByteAlign);
57
253k
    byte white = (ss->BlackIs1 ? 0 : 0xff);
58
59
253k
    if (raster < 0)
60
1
        return ERRC;
61
62
253k
    if (ss->Columns <= 0)
63
0
        return ERRC;
64
65
253k
    s_hcd_init_inline(ss);
66
    /* Because skip_white_pixels can look as many as 4 bytes ahead, */
67
    /* we need to allow 4 extra bytes at the end of the row buffers. */
68
253k
    ss->lbufstart = gs_alloc_bytes(st->memory, raster + CFD_BUFFER_SLOP * 2, "CFD lbuf");
69
253k
    ss->lprev = 0;
70
253k
    if (ss->lbufstart == 0)
71
0
        return ERRC;   /****** WRONG ******/
72
253k
    ss->lbuf = ss->lbufstart + CFD_BUFFER_SLOP;
73
253k
    memset(ss->lbufstart, 0xaa, CFD_BUFFER_SLOP);
74
253k
    memset(ss->lbuf, white, raster);
75
253k
    memset(ss->lbuf + raster, 0xaa, CFD_BUFFER_SLOP);  /* for Valgrind */
76
253k
    if (ss->K != 0) {
77
253k
        ss->lprevstart = gs_alloc_bytes(st->memory, raster + CFD_BUFFER_SLOP * 2, "CFD lprev");
78
253k
        if (ss->lprevstart == 0)
79
0
            return ERRC; /****** WRONG ******/
80
253k
        ss->lprev = ss->lprevstart + CFD_BUFFER_SLOP;
81
        /* Clear the initial reference line for 2-D encoding. */
82
253k
        memset(ss->lprev, white, raster);
83
        /* Ensure that the scan of the reference line will stop. */
84
253k
        memset(ss->lprev + raster, 0xaa, CFD_BUFFER_SLOP);
85
253k
        memset(ss->lprevstart, 0xaa, CFD_BUFFER_SLOP);
86
253k
    }
87
253k
    ss->k_left = min(ss->K, 0);
88
253k
    ss->run_color = 0;
89
253k
    ss->damaged_rows = 0;
90
253k
    ss->skipping_damage = false;
91
253k
    ss->cbit = 0;
92
253k
    ss->uncomp_run = 0;
93
253k
    ss->rows_left = (ss->Rows <= 0 || ss->EndOfBlock ? -1 : ss->Rows);
94
253k
    ss->row = 0;
95
253k
    ss->rpos = ss->wpos = -1;
96
253k
    ss->eol_count = 0;
97
253k
    ss->invert = white;
98
253k
    ss->min_left = 1;
99
253k
    return 0;
100
253k
}
101
102
/* Release the filter. */
103
static void
104
s_CFD_release(stream_state * st)
105
253k
{
106
253k
    stream_CFD_state *const ss = (stream_CFD_state *) st;
107
108
253k
    gs_free_object(st->memory, ss->lprevstart, "CFD lprev(close)");
109
253k
    gs_free_object(st->memory, ss->lbufstart, "CFD lbuf(close)");
110
253k
}
111
112
/* Declare the variables that hold the state. */
113
#define cfd_declare_state\
114
25.9M
        hcd_declare_state;\
115
25.9M
        register byte *q;\
116
25.9M
        int qbit
117
/* Load the state from the stream. */
118
#define cfd_load_state()\
119
44.3M
        hcd_load_state(),\
120
0
        q = ss->lbuf + ss->wpos, qbit = ss->cbit
121
/* Store the state back in the stream. */
122
#define cfd_store_state()\
123
44.3M
        hcd_store_state(),\
124
16.4M
        ss->wpos = q - ss->lbuf, ss->cbit = qbit
125
126
/* Macros to get blocks of bits from the input stream. */
127
/* Invariants: 0 <= bits_left <= bits_size; */
128
/* bits [bits_left-1..0] contain valid data. */
129
130
196k
#define avail_bits(n) hcd_bits_available(n)
131
49.6M
#define ensure_bits(n, outl) hcd_ensure_bits(n, outl)
132
53.6M
#define peek_bits(n) hcd_peek_bits(n)
133
196k
#define peek_var_bits(n) hcd_peek_var_bits(n)
134
44.5M
#define skip_bits(n) hcd_skip_bits(n)
135
136
/* Get a run from the stream. */
137
#ifdef DEBUG
138
#  define IF_DEBUG(expr) expr
139
#else
140
196k
#  define IF_DEBUG(expr) DO_NOTHING
141
#endif
142
143
static inline int get_run(stream_CFD_state *ss, stream_cursor_read *pr, const cfd_node decode[],
144
           int initial_bits, int min_bits, int *runlen, const char *str)
145
4.25M
{
146
4.25M
    cfd_declare_state;
147
4.25M
    const cfd_node *np;
148
4.25M
    int clen;
149
4.25M
    cfd_load_state();
150
151
4.25M
    HCD_ENSURE_BITS_ELSE(initial_bits) {
152
        /* We might still have enough bits for the specific code. */
153
261
        if (bits_left < min_bits)
154
82
            goto outl;
155
179
        np = &decode[hcd_peek_bits_left() << (initial_bits - bits_left)];
156
179
        if ((clen = np->code_length) > bits_left)
157
95
            goto outl;
158
84
        goto locl;
159
179
    }
160
4.25M
    np = &decode[peek_bits(initial_bits)];
161
4.25M
    if ((clen = np->code_length) > initial_bits) {
162
196k
            IF_DEBUG(uint init_bits = peek_bits(initial_bits));
163
196k
            if (!avail_bits(clen)) goto outl;
164
196k
            clen -= initial_bits;
165
196k
            skip_bits(initial_bits);
166
196k
            ensure_bits(clen, outl);                /* can't goto outl */
167
196k
            np = &decode[np->run_length + peek_var_bits(clen)];
168
196k
            if_debug4('W', "%s xcode=0x%x,%d rlen=%d\n", str,
169
196k
                      (init_bits << np->code_length) +
170
196k
                        peek_var_bits(np->code_length),
171
196k
                      initial_bits + np->code_length,
172
196k
                      np->run_length);
173
196k
            skip_bits(np->code_length);
174
4.05M
    } else {
175
4.05M
locl:
176
4.05M
       if_debug4('W', "%s code=0x%x,%d rlen=%d\n", str,
177
4.05M
                      peek_var_bits(clen), clen, np->run_length);
178
4.05M
       skip_bits(clen);
179
4.05M
    }
180
4.25M
    *runlen = np->run_length;
181
182
4.25M
    cfd_store_state();
183
4.25M
    return(0);
184
191
outl:
185
191
    cfd_store_state();
186
191
    return(-1);
187
4.25M
}
188
189
190
/* Skip data bits for a white run. */
191
/* rlen is either less than 64, or a multiple of 64. */
192
static inline int skip_data(stream_CFD_state *ss, stream_cursor_read *pr, int rlen)
193
1.29M
{
194
1.29M
    cfd_declare_state;
195
1.29M
    cfd_load_state();
196
1.29M
    (void)rlimit;
197
198
1.29M
    if ( (qbit -= rlen) < 0 )
199
956k
    {
200
956k
        q -= qbit >> 3, qbit &= 7;
201
956k
        if ( rlen >= 64 ) {
202
66.3k
            cfd_store_state();
203
66.3k
            return(-1);
204
66.3k
        }
205
956k
    }
206
1.29M
    cfd_store_state();
207
1.23M
    return(0);
208
1.29M
}
209
210
211
/* Invert data bits for a black run. */
212
/* If rlen >= 64, execute makeup_action: this is to handle */
213
/* makeup codes efficiently, since these are always a multiple of 64. */
214
215
static inline int invert_data(stream_CFD_state *ss, stream_cursor_read *pr, int *rlen, byte black_byte)
216
12.8M
{
217
12.8M
    byte *qlim = ss->lbuf + ss->raster + CFD_BUFFER_SLOP;
218
12.8M
    cfd_declare_state;
219
12.8M
    cfd_load_state();
220
12.8M
    (void)rlimit;
221
222
12.8M
    if (q >= qlim || q < ss->lbufstart) {
223
0
        return(-1);
224
0
    }
225
226
12.8M
    if ( (*rlen) > qbit )
227
7.85M
    {
228
7.85M
        if (q + ((*rlen - qbit) >> 3) > qlim) {
229
0
            return(-1);
230
0
        }
231
232
7.85M
        if (q >= ss->lbuf) {
233
7.63M
          *q++ ^= (1 << qbit) - 1;
234
7.63M
        }
235
218k
        else {
236
218k
            q++;
237
218k
        }
238
7.85M
        (*rlen) -= qbit;
239
240
7.85M
        if (q + ((*rlen) >> 3) >= qlim) {
241
1
            return(-1);
242
1
        }
243
244
7.85M
        switch ( (*rlen) >> 3 )
245
7.85M
        {
246
29.6k
          case 7:         /* original rlen possibly >= 64 */
247
29.6k
                  if ( (*rlen) + qbit >= 64 ) {
248
12.9k
                    goto d;
249
12.9k
                  }
250
16.7k
                  *q++ = black_byte;
251
47.9k
          case 6: *q++ = black_byte;
252
88.8k
          case 5: *q++ = black_byte;
253
150k
          case 4: *q++ = black_byte;
254
262k
          case 3: *q++ = black_byte;
255
666k
          case 2: *q++ = black_byte;
256
1.75M
          case 1: *q = black_byte;
257
1.75M
              (*rlen) &= 7;
258
1.75M
              if ( !(*rlen) ) {
259
336k
                  qbit = 0;
260
336k
                  break;
261
336k
              }
262
1.42M
              q++;
263
7.44M
          case 0:                 /* know rlen != 0 */
264
7.44M
              qbit = 8 - (*rlen);
265
7.44M
              *q ^= 0xff << qbit;
266
7.44M
              break;
267
59.8k
          default:        /* original rlen >= 64 */
268
72.7k
d:            memset(q, black_byte, (*rlen) >> 3);
269
72.7k
              q += (*rlen) >> 3;
270
72.7k
              (*rlen) &= 7;
271
72.7k
              if ( !(*rlen) ) {
272
19.2k
                qbit = 0;
273
19.2k
                q--;
274
19.2k
              }
275
53.4k
              else {
276
53.4k
                qbit = 8 - (*rlen);
277
53.4k
                *q ^= 0xff << qbit;
278
53.4k
              }
279
72.7k
              cfd_store_state();
280
72.7k
              return(-1);
281
7.85M
          }
282
7.85M
    }
283
4.94M
    else {
284
4.94M
            qbit -= (*rlen),
285
4.94M
            *q ^= ((1 << (*rlen)) - 1) << qbit;
286
287
4.94M
    }
288
12.8M
    cfd_store_state();
289
12.7M
    return(0);
290
12.8M
}
291
292
293
/* Buffer refill for CCITTFaxDecode filter */
294
static int cf_decode_eol(stream_CFD_state *, stream_cursor_read *);
295
static int cf_decode_1d(stream_CFD_state *, stream_cursor_read *);
296
static int cf_decode_2d(stream_CFD_state *, stream_cursor_read *);
297
static int cf_decode_uncompressed(stream_CFD_state *, stream_cursor_read *);
298
static int
299
s_CFD_process(stream_state * st, stream_cursor_read * pr,
300
              stream_cursor_write * pw, bool last)
301
257k
{
302
257k
    stream_CFD_state *const ss = (stream_CFD_state *) st;
303
257k
    int wstop = ss->raster - 1;
304
257k
    int eol_count = ss->eol_count;
305
257k
    int k_left = ss->k_left;
306
257k
    int rows_left = ss->rows_left;
307
257k
    int status = 0;
308
309
#ifdef DEBUG
310
    const byte *rstart = pr->ptr;
311
    const byte *wstart = pw->ptr;
312
313
#endif
314
315
    /* Update the pointers we actually use, in case GC moved the buffer */
316
257k
    ss->lbuf = ss->lbufstart + CFD_BUFFER_SLOP;
317
257k
    ss->lprev = ss->lprevstart + CFD_BUFFER_SLOP;
318
319
7.86M
top:
320
#ifdef DEBUG
321
    {
322
        hcd_declare_state;
323
        hcd_load_state();
324
        (void)rlimit;
325
        if_debug8m('w', ss->memory,
326
                   "[w]CFD_process top: eol_count=%d, k_left=%d, rows_left=%d\n"
327
                   "    bits="PRI_INTPTR", bits_left=%d, read %u, wrote %u%s\n",
328
                   eol_count, k_left, rows_left,
329
                   (intptr_t) bits, bits_left,
330
                   (uint) (p - rstart), (uint) (pw->ptr - wstart),
331
                   (ss->skipping_damage ? ", skipping damage" : ""));
332
    }
333
#endif
334
7.86M
    if (ss->skipping_damage) { /* Skip until we reach an EOL. */
335
0
        hcd_declare_state;
336
0
        int skip;
337
0
        (void)rlimit;
338
339
0
        status = 0;
340
0
        do {
341
0
            switch ((skip = cf_decode_eol(ss, pr))) {
342
0
                default:  /* not EOL */
343
0
                    hcd_load_state();
344
0
                    skip_bits(-skip);
345
0
                    hcd_store_state();
346
0
                    continue;
347
0
                case 0: /* need more input */
348
0
                    goto out;
349
0
                case 1: /* EOL */
350
0
                    {   /* Back up over the EOL. */
351
0
                        hcd_load_state();
352
0
                        bits_left += run_eol_code_length;
353
0
                        hcd_store_state();
354
0
                    }
355
0
                    ss->skipping_damage = false;
356
0
            }
357
0
        }
358
0
        while (ss->skipping_damage);
359
0
        ss->damaged_rows++;
360
0
    }
361
    /*
362
     * Check for a completed input scan line.  This isn't quite as
363
     * simple as it seems, because we could have run out of input data
364
     * between a makeup code and a 0-length termination code, or in a
365
     * 2-D line before a final horizontal code with a 0-length second
366
     * run.  There's probably a way to think about this situation that
367
     * doesn't require a special check, but I haven't found it yet.
368
     */
369
7.86M
    if (ss->wpos == wstop && ss->cbit <= (-ss->Columns & 7) &&
370
7.86M
        (k_left == 0 ? !(ss->run_color & ~1) : ss->run_color == 0)
371
7.86M
        ) {     /* Check for completed data to be copied to the client. */
372
        /* (We could avoid the extra copy step for 1-D, but */
373
        /* it's simpler not to, and it doesn't cost much.) */
374
7.61M
        if (ss->rpos < ss->wpos) {
375
7.61M
            stream_cursor_read cr;
376
377
7.61M
            cr.ptr = ss->lbuf + ss->rpos;
378
7.61M
            cr.limit = ss->lbuf + ss->wpos;
379
7.61M
            status = stream_move(&cr, pw);
380
7.61M
            ss->rpos = cr.ptr - ss->lbuf;
381
7.61M
            if (status)
382
3.37k
                goto out;
383
7.61M
        }
384
7.60M
        ss->row++;
385
7.60M
        if (rows_left > 0 && --rows_left == 0)
386
21
            goto ck_eol;  /* handle EOD if it is present */
387
7.60M
        if (ss->K != 0) {
388
7.60M
            byte *prev_bits = ss->lprev;
389
7.60M
            byte *prev_start = ss->lprevstart;
390
391
7.60M
            ss->lprev = ss->lbuf;
392
7.60M
            ss->lprevstart = ss->lbufstart;
393
7.60M
            ss->lbuf = prev_bits;
394
7.60M
            ss->lbufstart = prev_start;
395
7.60M
            if (ss->K > 0)
396
162
                k_left = (k_left == 0 ? ss->K : k_left) - 1;
397
7.60M
        }
398
7.60M
        ss->rpos = ss->wpos = -1;
399
7.60M
        ss->eol_count = eol_count = 0;
400
7.60M
        ss->cbit = 0;
401
7.60M
        ss->invert = (ss->BlackIs1 ? 0 : 0xff);
402
7.60M
        memset(ss->lbuf, ss->invert, wstop + 1);
403
7.60M
        ss->run_color = 0;
404
        /*
405
         * If EndOfLine is true, we want to include the byte padding
406
         * in the string of initial zeros in the EOL.  If EndOfLine
407
         * is false, we aren't sure what we should do....
408
         */
409
7.60M
        if (ss->EncodedByteAlign && !ss->EndOfLine)
410
0
            ss->bits_left &= ~7;
411
7.60M
    }
412
    /* If we're between scan lines, scan for EOLs. */
413
7.86M
    if (ss->wpos < 0) {
414
            /*
415
             * According to Adobe, the decoder should always check for
416
             * the EOD sequence, regardless of EndOfBlock: the Red Book's
417
             * documentation of EndOfBlock is wrong.
418
             */
419
7.86M
ck_eol:
420
8.12M
        while ((status = cf_decode_eol(ss, pr)) > 0) {
421
512k
            if_debug0m('w', ss->memory, "[w]EOL\n");
422
            /* If we are in a Group 3 mixed regime, */
423
            /* check the next bit for 1- vs. 2-D. */
424
512k
            if (ss->K > 0) {
425
169
                hcd_declare_state;
426
169
                hcd_load_state();
427
169
                ensure_bits(1, out); /* can't fail */
428
169
                k_left = (peek_bits(1) ? 0 : 1);
429
169
                skip_bits(1);
430
169
                hcd_store_state();
431
169
            }
432
512k
            ++eol_count;
433
512k
            if (eol_count == (ss->K < 0 ? 2 : 6)) {
434
253k
                status = EOFC;
435
253k
                goto out;
436
253k
            }
437
512k
        }
438
7.60M
        if (rows_left == 0) {
439
20
            status = EOFC;
440
20
            goto out;
441
20
        }
442
7.60M
        if (status == 0)  /* input empty while scanning EOLs */
443
175
            goto out;
444
7.60M
        switch (eol_count) {
445
7.60M
            case 0:
446
7.60M
                if (ss->EndOfLine) { /* EOL is required, but none is present. */
447
0
                    status = ERRC;
448
0
                    goto check;
449
0
                }
450
7.60M
            case 1:
451
7.60M
                break;
452
0
            default:
453
0
                status = ERRC;
454
0
                goto check;
455
7.60M
        }
456
7.60M
    }
457
    /* Now decode actual data. */
458
7.60M
    if (k_left < 0) {
459
7.60M
        if_debug0m('w', ss->memory, "[w2]new row\n");
460
7.60M
        status = cf_decode_2d(ss, pr);
461
7.60M
    } else if (k_left == 0) {
462
6.43k
        if_debug0m('w', ss->memory, "[w1]new row\n");
463
6.43k
        status = cf_decode_1d(ss, pr);
464
6.43k
    } else {
465
158
        if_debug1m('w', ss->memory, "[w1]new 2-D row, k_left=%d\n", k_left);
466
158
        status = cf_decode_2d(ss, pr);
467
158
    }
468
7.60M
    if_debug3m('w', ss->memory, "[w]CFD status = %d, wpos = %d, cbit = %d\n",
469
7.60M
               status, ss->wpos, ss->cbit);
470
7.60M
  check:switch (status) {
471
7.60M
        case 1:   /* output full */
472
7.60M
            goto top;
473
93
        case ERRC:
474
            /* Check for special handling of damaged rows. */
475
93
            if (ss->damaged_rows >= ss->DamagedRowsBeforeError ||
476
93
                !(ss->EndOfLine && ss->K >= 0)
477
93
                )
478
93
                break;
479
            /* Substitute undamaged data if appropriate. */
480
            /****** NOT IMPLEMENTED YET ******/
481
0
            {
482
0
                ss->wpos = wstop;
483
0
                ss->cbit = -ss->Columns & 7;
484
0
                ss->run_color = 0;
485
0
            }
486
0
            ss->skipping_damage = true;
487
0
            goto top;
488
242
        default:
489
242
            ss->damaged_rows = 0; /* finished a good row */
490
7.60M
    }
491
257k
  out:ss->k_left = k_left;
492
257k
    ss->rows_left = rows_left;
493
257k
    ss->eol_count = eol_count;
494
257k
    if (ss->ErrsAsEOD && status < 0)
495
0
        return EOFC;
496
257k
    return status;
497
257k
}
498
499
/*
500
 * Decode a leading EOL, if any.
501
 * If an EOL is present, skip over it and return 1;
502
 * if no EOL is present, read no input and return -N, where N is the
503
 * number of initial bits that can be skipped in the search for an EOL;
504
 * if more input is needed, return 0.
505
 * Note that if we detected an EOL, we know that we can back up over it;
506
 * if we detected an N-bit non-EOL, we know that at least N bits of data
507
 * are available in the buffer.
508
 */
509
static int
510
cf_decode_eol(stream_CFD_state * ss, stream_cursor_read * pr)
511
8.12M
{
512
8.12M
    hcd_declare_state;
513
8.12M
    int zeros;
514
8.12M
    int look_ahead;
515
516
8.12M
    hcd_load_state();
517
17.8M
    for (zeros = 0; zeros < run_eol_code_length - 1; zeros++) {
518
17.2M
        ensure_bits(1, out);
519
17.2M
        if (peek_bits(1))
520
7.60M
            return -(zeros + 1);
521
9.69M
        skip_bits(1);
522
9.69M
    }
523
    /* We definitely have an EOL.  Skip further zero bits. */
524
512k
    look_ahead = (ss->K > 0 ? 2 : 1);
525
512k
    for (;;) {
526
512k
        ensure_bits(look_ahead, back);
527
512k
        if (peek_bits(1))
528
512k
            break;
529
0
        skip_bits(1);
530
0
    }
531
512k
    skip_bits(1);
532
512k
    hcd_store_state();
533
512k
    return 1;
534
0
  back:     /*
535
                                 * We ran out of data while skipping zeros.
536
                                 * We know we are at a byte boundary, and have just skipped
537
                                 * at least run_eol_code_length - 1 zeros.  However,
538
                                 * bits_left may be 1 if look_ahead == 2.
539
                                 */
540
0
    bits &= (1 << bits_left) - 1;
541
0
    bits_left += run_eol_code_length - 1;
542
0
    hcd_store_state();
543
175
  out:return 0;
544
0
}
545
546
/* Decode a 1-D scan line. */
547
static int
548
cf_decode_1d(stream_CFD_state * ss, stream_cursor_read * pr)
549
6.43k
{
550
6.43k
    cfd_declare_state;
551
6.43k
    byte black_byte = (ss->BlackIs1 ? 0xff : 0);
552
6.43k
    int end_bit = -ss->Columns & 7;
553
6.43k
    byte *stop = ss->lbuf - 1 + ss->raster;
554
6.43k
    int run_color = ss->run_color;
555
6.43k
    int status;
556
6.43k
    int bcnt;
557
6.43k
    (void)rlimit;
558
559
6.43k
    cfd_load_state();
560
6.43k
    if_debug1m('w', ss->memory, "[w1]entry run_color = %d\n", ss->run_color);
561
6.43k
    if (run_color > 0)
562
55
        goto db;
563
6.38k
    else
564
6.38k
        goto dw;
565
348k
#define q_at_stop() (q >= stop && (qbit <= end_bit || q > stop))
566
171k
    while (1)
567
171k
    {
568
171k
        run_color = 0;
569
171k
        if (q_at_stop())
570
289
            goto done;
571
572
177k
  dw:   /* Decode a white run. */
573
228k
        do {
574
228k
            cfd_store_state();
575
228k
            status = get_run(ss, pr, cf_white_decode, cfd_white_initial_bits,
576
228k
                             cfd_white_min_bits, &bcnt, "[w1]white");
577
228k
            cfd_load_state();
578
228k
            if (status < 0) {
579
                /* We already set run_color to 0 or -1. */
580
89
                status = 0;
581
89
                goto out;
582
89
            }
583
584
228k
            if (bcnt < 0) {   /* exceptional situation */
585
5
                switch (bcnt) {
586
0
                    case run_uncompressed: /* Uncompressed data. */
587
0
                        cfd_store_state();
588
0
                        bcnt = cf_decode_uncompressed(ss, pr);
589
0
                        if (bcnt < 0)
590
0
                            return bcnt;
591
0
                        cfd_load_state();
592
0
                        if (bcnt)
593
0
                            goto db;
594
0
                        else
595
0
                            continue; /* Restart decoding white */
596
                        /*case run_error: */
597
                        /*case run_zeros: *//* Premature end-of-line. */
598
5
                    default:
599
5
                        status = ERRC;
600
5
                        goto out;
601
5
                }
602
5
            }
603
604
228k
            cfd_store_state();
605
228k
            status = skip_data(ss, pr, bcnt);
606
228k
            cfd_load_state();
607
228k
            if (status < 0) {
608
                /* If we run out of data after a makeup code, */
609
                /* note that we are still processing a white run. */
610
51.2k
                run_color = -1;
611
51.2k
            }
612
228k
        } while (status < 0);
613
614
177k
        if (q_at_stop()) {
615
5.99k
            run_color = 0;    /* not inside a run */
616
6.28k
  done:
617
6.28k
            if (q > stop || qbit < end_bit)
618
9
                status = ERRC;
619
6.27k
            else
620
6.27k
                status = 1;
621
6.28k
            break;
622
5.99k
        }
623
171k
        run_color = 1;
624
625
171k
  db:   /* Decode a black run. */
626
171k
        do {
627
171k
            cfd_store_state();
628
171k
            status = get_run(ss, pr, cf_black_decode, cfd_black_initial_bits,
629
171k
                             cfd_black_min_bits, &bcnt, "[w1]black");
630
171k
            cfd_load_state();
631
171k
            if (status < 0) {
632
                /* We already set run_color to 1 or 2. */
633
55
                status = 0;
634
55
                goto out;
635
55
            }
636
637
171k
            if (bcnt < 0) {  /* All exceptional codes are invalid here. */
638
                /****** WRONG, uncompressed IS ALLOWED ******/
639
4
                status = ERRC;
640
4
                goto out;
641
4
            }
642
643
            /* Invert bits designated by black run. */
644
171k
            cfd_store_state();
645
171k
            status = invert_data(ss, pr, &bcnt, black_byte);
646
171k
            cfd_load_state();
647
171k
            if (status < 0) {
648
                /* If we run out of data after a makeup code, */
649
                /* note that we are still processing a black run. */
650
36
                run_color = 2;
651
36
            }
652
171k
        } while (status < 0);
653
171k
    }
654
655
6.43k
out:
656
6.43k
    cfd_store_state();
657
6.43k
    ss->run_color = run_color;
658
6.43k
    if_debug1m('w', ss->memory, "[w1]exit run_color = %d\n", run_color);
659
6.43k
    return status;
660
0
}
661
662
/* Decode a 2-D scan line. */
663
static int
664
cf_decode_2d(stream_CFD_state * ss, stream_cursor_read * pr)
665
7.60M
{
666
7.60M
    cfd_declare_state;
667
7.60M
    byte invert_white = (ss->BlackIs1 ? 0 : 0xff);
668
7.60M
    byte black_byte = ~invert_white;
669
7.60M
    byte invert = ss->invert;
670
7.60M
    int end_count = -ss->Columns & 7;
671
7.60M
    uint raster = ss->raster;
672
7.60M
    byte *q0 = ss->lbuf;
673
7.60M
    byte *prev_q01 = ss->lprev + 1;
674
7.60M
    byte *endptr = q0 - 1 + raster;
675
7.60M
    int init_count = raster << 3;
676
7.60M
    register int count;
677
7.60M
    int rlen;
678
7.60M
    int status;
679
680
7.60M
    cfd_load_state();
681
7.60M
    count = ((endptr - q) << 3) + qbit;
682
7.60M
    endptr[1] = 0xa0;   /* a byte with some 0s and some 1s, */
683
    /* to ensure run scan will stop */
684
7.60M
    if_debug1m('W', ss->memory, "[w2]raster=%d\n", raster);
685
7.60M
    switch (ss->run_color) {
686
6
        case -2:
687
6
            ss->run_color = 0;
688
6
            goto hww;
689
12
        case -1:
690
12
            ss->run_color = 0;
691
12
            goto hbw;
692
9
        case 1:
693
9
            ss->run_color = 0;
694
9
            goto hwb;
695
3
        case 2:
696
3
            ss->run_color = 0;
697
3
            goto hbb;
698
            /*case 0: */
699
7.60M
    }
700
701
    /* top of decode loop */
702
39.2M
    while (1)
703
39.2M
    {
704
39.2M
        if (count <= end_count) {
705
7.60M
            status = (count < end_count ? ERRC : 1);
706
7.60M
            goto out;
707
7.60M
        }
708
        /* If invert == invert_white, white and black have their */
709
        /* correct meanings; if invert == ~invert_white, */
710
        /* black and white are interchanged. */
711
39.2M
        if_debug1m('W', ss->memory, "[w2]%4d:\n", count);
712
#ifdef DEBUG
713
        /* Check the invariant between q, qbit, and count. */
714
        {
715
            int pcount = (endptr - q) * 8 + qbit;
716
717
            if (pcount != count)
718
                dmlprintf2(ss->memory, "[w2]Error: count=%d pcount=%d\n",
719
                           count, pcount);
720
        }
721
#endif
722
        /*
723
         * We could just use get_run here, but we can do better.  However,
724
         * we must be careful to handle the case where the very last codes
725
         * in the input stream are 1-bit "vertical 0" codes: we can't just
726
         * use ensure_bits(3, ...) and go to get more data if it fails.
727
         */
728
31.6M
        ensure_bits(3, out3);
729
58.9M
#define vertical_0 (countof(cf2_run_vertical) / 2)
730
31.6M
        switch (peek_bits(3)) {
731
23.4M
        default /*4..7*/ :  /* vertical(0) */
732
23.4M
            if (0) {
733
66
 out3:
734
                /* Unless it's a 1-bit "vertical 0" code, exit. */
735
66
                if (!(bits_left > 0 && peek_bits(1)))
736
50
                    goto out0;
737
66
            }
738
23.4M
            skip_bits(1);
739
23.4M
            rlen = vertical_0;
740
23.4M
            break;
741
2.97M
        case 2:   /* vertical(+1) */
742
2.97M
            skip_bits(3);
743
2.97M
            rlen = vertical_0 + 1;
744
2.97M
            break;
745
2.40M
        case 3:   /* vertical(-1) */
746
2.40M
            skip_bits(3);
747
2.40M
            rlen = vertical_0 - 1;
748
2.40M
            break;
749
1.05M
        case 1:   /* horizontal */
750
1.05M
            skip_bits(3);
751
1.05M
            if (invert == invert_white) {
752
                /* We handle horizontal decoding here, so that we can
753
                 * branch back into it if we run out of input data. */
754
                /* White, then black. */
755
671k
  hww:
756
679k
                do {
757
679k
                    cfd_store_state();
758
679k
                    status = get_run(ss, pr, cf_white_decode,
759
679k
                                     cfd_white_initial_bits,
760
679k
                                     cfd_white_min_bits, &rlen, " white");
761
679k
                    cfd_load_state();
762
679k
                    if (status < 0) {
763
6
                        ss->run_color = -2;
764
6
                        goto out0;
765
6
                    }
766
767
679k
                    if ((count -= rlen) < end_count) {
768
9
                        status = ERRC;
769
9
                        goto out;
770
9
                    }
771
679k
                    if (rlen < 0) goto rlen_lt_zero;
772
773
679k
                    cfd_store_state();
774
679k
                    status = skip_data(ss, pr, rlen);
775
679k
                    cfd_load_state();
776
679k
                } while (status < 0);
777
778
                /* Handle the second half of a white-black horizontal code. */
779
671k
  hwb:
780
679k
                do {
781
679k
                    cfd_store_state();
782
679k
                    status = get_run(ss, pr, cf_black_decode,
783
679k
                                     cfd_black_initial_bits,
784
679k
                                     cfd_black_min_bits, &rlen, " black");
785
679k
                    cfd_load_state();
786
679k
                    if (status < 0){
787
10
                        ss->run_color = 1;
788
10
                        goto out0;
789
10
                    }
790
791
679k
                    if ((count -= rlen) < end_count) {
792
6
                        status = ERRC;
793
6
                        goto out;
794
6
                    }
795
679k
                    if (rlen < 0) goto rlen_lt_zero;
796
797
679k
                    cfd_store_state();
798
679k
                    status = invert_data(ss, pr, &rlen, black_byte);
799
679k
                    cfd_load_state();
800
679k
                } while (status < 0);
801
671k
            } else {
802
                /* Black, then white. */
803
383k
  hbb:
804
386k
                do {
805
386k
                    cfd_store_state();
806
386k
                    status = get_run(ss, pr, cf_black_decode,
807
386k
                                     cfd_black_initial_bits,
808
386k
                                     cfd_black_min_bits, &rlen, " black");
809
386k
                    cfd_load_state();
810
386k
                    if (status < 0) {
811
3
                        ss->run_color = 2;
812
3
                        goto out0;
813
3
                    }
814
815
386k
                    if ((count -= rlen) < end_count) {
816
9
                        status = ERRC;
817
9
                        goto out;
818
9
                    }
819
386k
                    if (rlen < 0) goto rlen_lt_zero;
820
821
386k
                    cfd_store_state();
822
386k
                    status = invert_data(ss, pr, &rlen, black_byte);
823
386k
                    cfd_load_state();
824
386k
                }
825
386k
                while (status < 0);
826
827
                /* Handle the second half of a black-white horizontal code. */
828
383k
  hbw:
829
390k
                do {
830
390k
                    cfd_store_state();
831
390k
                    status = get_run(ss, pr, cf_white_decode,
832
390k
                                     cfd_white_initial_bits,
833
390k
                                     cfd_white_min_bits, &rlen, " white");
834
390k
                    cfd_load_state();
835
390k
                    if (status < 0) {
836
12
                        ss->run_color = -1;
837
12
                        goto out0;
838
12
                    }
839
840
390k
                    if ((count -= rlen) < end_count) {
841
16
                        status = ERRC;
842
16
                        goto out;
843
16
                    }
844
390k
                    if (rlen < 0) goto rlen_lt_zero;
845
846
390k
                    cfd_store_state();
847
390k
                    status = skip_data(ss, pr, rlen);
848
390k
                    cfd_load_state();
849
390k
                } while (status < 0);
850
383k
            }
851
1.05M
            continue; /* jump back to top of decode loop */
852
1.71M
        case 0:   /* everything else */
853
1.71M
            cfd_store_state();
854
1.71M
            status = get_run(ss, pr, cf_2d_decode, cfd_2d_initial_bits,
855
1.71M
                             cfd_2d_min_bits, &rlen, "[w2]");
856
1.71M
            cfd_load_state();
857
1.71M
            if (status < 0) {
858
16
                goto out0;
859
16
            }
860
861
            /* rlen may be run2_pass, run_uncompressed, or */
862
            /* 0..countof(cf2_run_vertical)-1. */
863
1.71M
            if (rlen < 0) {
864
452k
rlen_lt_zero:
865
452k
                switch (rlen) {
866
452k
                case run2_pass:
867
452k
                    break;
868
0
                case run_uncompressed:
869
0
                {
870
0
                    int which;
871
872
0
                    cfd_store_state();
873
0
                    which = cf_decode_uncompressed(ss, pr);
874
0
                    if (which < 0) {
875
0
                        status = which;
876
0
                        goto out;
877
0
                    }
878
0
                    cfd_load_state();
879
/****** ADJUST count ******/
880
0
                    invert = (which ? ~invert_white : invert_white);
881
0
                    continue; /* jump back to top of decode loop */
882
0
                }
883
15
                default:  /* run_error, run_zeros */
884
15
                    status = ERRC;
885
15
                    goto out;
886
452k
                }
887
452k
            }
888
31.6M
        }
889
        /* Interpreting the run requires scanning the */
890
        /* previous ('reference') line. */
891
30.5M
        {
892
30.5M
            int prev_count = count;
893
30.5M
            byte prev_data;
894
30.5M
            int dlen;
895
30.5M
            static const byte count_bit[8] =
896
30.5M
                                   {0x80, 1, 2, 4, 8, 0x10, 0x20, 0x40};
897
30.5M
            byte *prev_q = prev_q01 + (q - q0);
898
30.5M
            int plen;
899
900
30.5M
            if (!(count & 7))
901
10.0M
                prev_q++;   /* because of skip macros */
902
30.5M
            prev_data = prev_q[-1] ^ invert;
903
            /* Find the b1 transition. */
904
30.5M
            if ((prev_data & count_bit[prev_count & 7]) &&
905
30.5M
                (prev_count < init_count || invert != invert_white)
906
30.5M
                ) {     /* Look for changing white first. */
907
3.56M
                if_debug1m('W', ss->memory, " data=0x%x", prev_data);
908
3.56M
                skip_black_pixels(prev_data, prev_q,
909
3.56M
                                  prev_count, invert, plen);
910
3.56M
                if (prev_count < end_count)    /* overshot */
911
744
                    prev_count = end_count;
912
3.56M
                if_debug1m('W', ss->memory, " b1 other=%d", prev_count);
913
3.56M
            }
914
30.5M
            if (prev_count != end_count) {
915
30.5M
                if_debug1m('W', ss->memory, " data=0x%x", prev_data);
916
30.5M
                skip_white_pixels(prev_data, prev_q,
917
30.5M
                                  prev_count, invert, plen);
918
30.5M
                if (prev_count < end_count)    /* overshot */
919
17.6k
                    prev_count = end_count;
920
30.5M
                if_debug1m('W', ss->memory, " b1 same=%d", prev_count);
921
30.5M
            }
922
            /* b1 = prev_count; */
923
30.5M
            if (rlen == run2_pass) { /* Pass mode.  Find b2. */
924
452k
                if (prev_count != end_count) {
925
452k
                    if_debug1m('W', ss->memory, " data=0x%x", prev_data);
926
452k
                    skip_black_pixels(prev_data, prev_q,
927
452k
                                      prev_count, invert, plen);
928
452k
                    if (prev_count < end_count)  /* overshot */
929
8
                        prev_count = end_count;
930
452k
                }
931
                /* b2 = prev_count; */
932
452k
                if_debug2m('W', ss->memory, " b2=%d, pass %d\n",
933
452k
                           prev_count, count - prev_count);
934
30.1M
            } else {   /* Vertical coding. */
935
                /* Remember that count counts *down*. */
936
30.1M
                prev_count += rlen - vertical_0; /* a1 */
937
30.1M
                if_debug2m('W', ss->memory, " vertical %d -> %d\n",
938
30.1M
                           (int)(rlen - vertical_0), prev_count);
939
30.1M
            }
940
            /* Now either invert or skip from count */
941
            /* to prev_count, and reset count. */
942
30.5M
            if (invert == invert_white) { /* Skip data bits. */
943
18.9M
                q = endptr - (prev_count >> 3);
944
18.9M
                qbit = prev_count & 7;
945
18.9M
            } else {   /* Invert data bits. */
946
11.5M
                dlen = count - prev_count;
947
948
11.5M
                cfd_store_state();
949
11.5M
                (void)invert_data(ss, pr, &dlen, black_byte);
950
11.5M
                cfd_load_state();
951
11.5M
            }
952
30.5M
            count = prev_count;
953
30.5M
            if (rlen >= 0)    /* vertical mode */
954
30.1M
                invert = ~invert; /* polarity changes */
955
30.5M
        }
956
        /* jump back to top of decode loop */
957
30.5M
    }
958
959
97
  out0:status = 0;
960
    /* falls through */
961
7.60M
  out:cfd_store_state();
962
7.60M
    ss->invert = invert;
963
    /* Ignore an error (missing EOFB/RTC when EndOfBlock == true) */
964
    /* if we have finished all rows. */
965
7.60M
    if (status == ERRC && ss->Rows > 0 && ss->row > ss->Rows)
966
1
        status = EOFC;
967
7.60M
    return status;
968
97
}
969
970
#if 1       /*************** */
971
static int
972
cf_decode_uncompressed(stream_CFD_state * ss, stream_cursor_read * pr)
973
0
{
974
0
    return ERRC;
975
0
}
976
#else /*************** */
977
978
/* Decode uncompressed data. */
979
/* (Not tested: no sample data available!) */
980
/****** DOESN'T CHECK FOR OVERFLOWING SCAN LINE ******/
981
static int
982
cf_decode_uncompressed(stream * s)
983
{
984
    cfd_declare_state;
985
    const cfd_node *np;
986
    int clen, rlen;
987
988
    cfd_load_state();
989
    while (1) {
990
        ensure_bits(cfd_uncompressed_initial_bits, NOOUT);
991
        np = &cf_uncompressed_decode[peek_bits(cfd_uncompressed_initial_bits)];
992
        clen = np->code_length;
993
        rlen = np->run_length;
994
        if (clen > cfd_uncompressed_initial_bits) { /* Must be an exit code. */
995
            break;
996
        }
997
        if (rlen == cfd_uncompressed_initial_bits) {  /* Longest representable white run */
998
            if_debug1m('W', s->memory, "[wu]%d\n", rlen);
999
            if ((qbit -= cfd_uncompressed_initial_bits) < 0)
1000
                qbit += 8, q++;
1001
        } else {
1002
            if_debug1m('W', s->memory, "[wu]%d+1\n", rlen);
1003
            if (qbit -= rlen < 0)
1004
                qbit += 8, q++;
1005
            *q ^= 1 << qbit;
1006
        }
1007
        skip_bits(clen);
1008
    }
1009
    clen -= cfd_uncompressed_initial_bits;
1010
    skip_bits(cfd_uncompressed_initial_bits);
1011
    ensure_bits(clen, NOOUT);
1012
    np = &cf_uncompressed_decode[rlen + peek_var_bits(clen)];
1013
    rlen = np->run_length;
1014
    skip_bits(np->code_length);
1015
    if_debug1m('w', s->memory, "[wu]exit %d\n", rlen);
1016
    if (rlen >= 0) {    /* Valid exit code, rlen = 2 * run length + next polarity */
1017
        if ((qbit -= rlen >> 1) < 0)
1018
            qbit += 8, q++;
1019
        rlen &= 1;
1020
    }
1021
  out:        /******* WRONG ******/
1022
    cfd_store_state();
1023
    return rlen;
1024
}
1025
1026
#endif /*************** */
1027
1028
/* Stream template */
1029
const stream_template s_CFD_template =
1030
{&st_CFD_state, s_CFD_init, s_CFD_process, 1, 1, s_CFD_release,
1031
 s_CFD_set_defaults
1032
};