/* Copyright (C) 2019-2023 Artifex Software, Inc.
   All Rights Reserved.

   This software is provided AS-IS with no warranty, either express or
   implied.

   This software is distributed under license and may not be copied,
   modified or distributed except as expressly authorized under the terms
   of the license contained in the file LICENSE in this distribution.

   Refer to licensing information at http://www.artifex.com or contact
   Artifex Software, Inc.,  39 Mesa Street, Suite 108A, San Francisco,
   CA 94129, USA, for further information.
*/

/* tifftop.c */
/* Top-level API implementation of "TIF" Language Interface */
#ifdef TIFF_INCLUDED

#include "pltop.h"
#include "gserrors.h"
#include "gxdevice.h"
#include "gsstate.h"
#include "strimpl.h"
#include "gscoord.h"
#include "gsicc_manage.h"
#include "gspaint.h"
#include "plmain.h"
#include "tiffio.h"
#if defined(SHARE_JPEG) && SHARE_JPEG==0
#include "jmemcust.h"
#endif
#include "gsmchunk.h"

#include <limits.h>

/* Forward decls */

/************************************************************/
/******** Language wrapper implementation (see pltop.h) *****/
/************************************************************/

typedef enum
{
    ii_state_identifying = 0,
    ii_state_tiff,
    ii_state_tiff_header,
    ii_state_tiff_decode,
    ii_state_flush
} ii_state;

/*
 * Tiff interpreter instance
 */
typedef struct tiff_interp_instance_s {
    gs_memory_t       *memory;
    gs_memory_t       *cmemory;
    gx_device         *dev;
    gx_device         *nulldev;

    gs_color_space    *gray;
    gs_color_space    *rgb;
    gs_color_space    *cmyk;

    /* Tiff parser state machine */
    ii_state           state;

    uint32_t           bpp;
    uint32_t           bpc;
    uint32_t           pal_bpc;
    uint32_t           cs;
    uint32_t           width;
    uint32_t           height;
    uint32_t           xresolution;
    uint32_t           yresolution;
    uint32_t           tile_height;
    uint32_t           tile_width;
    uint32_t           tiled;
    uint32_t           compression;
    uint32_t           photometric;
    uint8_t           *palette;

    uint32_t           raw_num_comps; /* As specified in the file */
    uint32_t           num_comps;     /* After processing */
    uint32_t           raw_byte_width;
    uint32_t           byte_width;

    gs_image_t         image;
    gs_image_enum     *penum;
    gs_gstate         *pgs;

    size_t             buffer_full;
    size_t             buffer_max;
    byte              *tiff_buffer;
    size_t             file_pos;
    TIFF              *handle;
    int                is_rgba;

    byte              *samples;
    byte              *proc_samples;
#if defined(SHARE_JPEG) && SHARE_JPEG==0
    jpeg_cust_mem_data jmem;
#endif
} tiff_interp_instance_t;

static int
tiff_detect_language(const char *s, int len)
{
    const byte *hdr = (const byte *)s;
    if (len >= 4) {
        if (hdr[0] == 'I' && hdr[1] == 'I' && (hdr[2] == 42 || hdr[2] == 43) && hdr[3] == 0)
            return 100; /* Intel (LSB) order */
        if (hdr[0] == 'M' && hdr[1] == 'M' && hdr[2] == 0 && (hdr[3] == 42 || hdr[3] == 43))
            return 100; /* Motorola (MSB) order */
    }

    return 0;
}

static const pl_interp_characteristics_t tiff_characteristics = {
    "TIFF",
    tiff_detect_language,
};

/* Get implementation's characteristics */
static const pl_interp_characteristics_t * /* always returns a descriptor */
tiff_impl_characteristics(const pl_interp_implementation_t *impl)     /* implementation of interpreter to alloc */
{
  return &tiff_characteristics;
}

static void
tiff_deallocate(tiff_interp_instance_t *tiff)
{
    if (tiff == NULL)
        return;

    rc_decrement_cs(tiff->gray, "tiff_deallocate");
    rc_decrement_cs(tiff->rgb, "tiff_deallocate");
    rc_decrement_cs(tiff->cmyk, "tiff_deallocate");

    if (tiff->pgs != NULL)
        gs_gstate_free_chain(tiff->pgs);
    gs_free_object(tiff->memory, tiff, "tiff_impl_allocate_interp_instance");
}

/* Deallocate a interpreter instance */
static int
tiff_impl_deallocate_interp_instance(pl_interp_implementation_t *impl)
{
    tiff_interp_instance_t *tiff = (tiff_interp_instance_t *)impl->interp_client_data;

    tiff_deallocate(tiff);
    impl->interp_client_data = NULL;

    return 0;
}

/* Do per-instance interpreter allocation/init. */
static int
tiff_impl_allocate_interp_instance(pl_interp_implementation_t *impl, gs_memory_t *mem)
{
    int code;
    tiff_interp_instance_t *tiff
        = (tiff_interp_instance_t *)gs_alloc_bytes(mem,
                                                  sizeof(tiff_interp_instance_t),
                                                  "tiff_impl_allocate_interp_instance");
    if (!tiff)
        return_error(gs_error_VMerror);
    memset(tiff, 0, sizeof(*tiff));

    tiff->memory = mem;
    tiff->pgs = gs_gstate_alloc(mem);
    if (tiff->pgs == NULL)
        goto failVM;

    /* Push one save level onto the stack to assuage the memory handling */
    code = gs_gsave(tiff->pgs);
    if (code < 0)
        goto fail;

    code = gsicc_init_iccmanager(tiff->pgs);
    if (code < 0)
        goto fail;

    tiff->gray = gs_cspace_new_ICC(mem, tiff->pgs, 1);
    tiff->rgb  = gs_cspace_new_ICC(mem, tiff->pgs, 3);
    tiff->cmyk = gs_cspace_new_ICC(mem, tiff->pgs, 4);

    impl->interp_client_data = tiff;

    return 0;

failVM:
    code = gs_note_error(gs_error_VMerror);
fail:
    (void)tiff_deallocate(tiff);
    return code;
}

/*
 * Get the allocator with which to allocate a device
 */
static gs_memory_t *
tiff_impl_get_device_memory(pl_interp_implementation_t *impl)
{
    tiff_interp_instance_t *tiff = (tiff_interp_instance_t *)impl->interp_client_data;

    return tiff->dev ? tiff->dev->memory : NULL;
}

#if 0 /* UNUSED */
static int
tiff_impl_set_param(pl_interp_implementation_t *impl,
                   pl_set_param_type           type,
                   const char                 *param,
                   const void                 *val)
{
    tiff_interp_instance_t *tiff = (tiff_interp_instance_t *)impl->interp_client_data;

    /* No params set here */
    return 0;
}

static int
tiff_impl_add_path(pl_interp_implementation_t *impl,
                  const char                 *path)
{
    tiff_interp_instance_t *tiff = (tiff_interp_instance_t *)impl->interp_client_data;

    /* No paths to add */
    return 0;
}

static int
tiff_impl_post_args_init(pl_interp_implementation_t *impl)
{
    tiff_interp_instance_t *tiff = (jpg_interp_instance_t *)impl->interp_client_data;

    /* No post args processing */
    return 0;
}
#endif

/* Prepare interp instance for the next "job" */
static int
tiff_impl_init_job(pl_interp_implementation_t *impl,
                  gx_device                  *device)
{
    tiff_interp_instance_t *tiff = (tiff_interp_instance_t *)impl->interp_client_data;

    tiff->dev = device;
    tiff->state = ii_state_identifying;
    tiff->buffer_full = 0;

    return 0;
}

#if 0 /* UNUSED */
static int
tiff_impl_run_prefix_commands(pl_interp_implementation_t *impl,
                             const char                 *prefix)
{
    tiff_interp_instance_t *tiff = (tiff_interp_instance_t *)impl->interp_client_data;

    return 0;
}

static int
tiff_impl_process_file(pl_interp_implementation_t *impl, const char *filename)
{
    tiff_interp_instance_t *tiff = (tiff_interp_instance_t *)impl->interp_client_data;

    return 0;
}
#endif

/* Do any setup for parser per-cursor */
static int                      /* ret 0 or +ve if ok, else -ve error code */
tiff_impl_process_begin(pl_interp_implementation_t * impl)
{
    return 0;
}

/* Ensure we have 'required' bytes to read, and further ensure
 * that we have no UEL's within those bytes. */
static int
ensure_bytes(tiff_interp_instance_t *jpg, stream_cursor_read *pr, int required)
{
    int n;
    const uint8_t *p = pr->ptr+1;
    const uint8_t *q;
    int avail;

    /* Find out how many bytes we need to check */
    n = pr->limit - pr->ptr;
    if (n > required)
        n = required;

    /* Make sure there are no UELs in that block */
    q = p + n;
    while (p != q) {
        while (p != q && *p != '\033')
            p++;
        if (p == q)
            break;
        avail = pr->limit - pr->ptr;
        if (memcmp(p, "\033%-12345X", min(avail, 9)) == 0) {
            /* At least a partial match to a UEL */
            return avail < 9 ? gs_error_NeedInput : gs_error_InterpreterExit;
        }
        p++;
    }

    /* If we have enough bytes, great, if not, get some more */
    return (n < required) ? gs_error_NeedInput : 0;
}

static int
flush_to_uel(stream_cursor_read *pr)
{
    const uint8_t *p = pr->ptr+1;
    const uint8_t *q = pr->limit+1;
    int avail;

    while (p != q) {
        while (p != q && *p != '\033')
            p++;
        if (p == q)
            break;
        avail = pr->limit - pr->ptr;
        if (memcmp(p, "\033%-12345X", min(avail, 9)) == 0) {
            /* At least a partial match to a UEL. Bin everything to
             * the start of the match. */
            pr->ptr = p-1;
            if (avail == 9) /* Complete match. Exit! */
                return gs_error_InterpreterExit;
            /* Partial match. Get more data. */
            return gs_error_NeedInput;
        }
        p++;
    }

    pr->ptr = pr->limit;

    return 0;
}

static int
bytes_until_uel(const stream_cursor_read *pr)
{
    const uint8_t *p = pr->ptr+1;
    const uint8_t *q = pr->limit+1;
    int avail;

    while (p != q) {
        while (p != q && *p != '\033')
            p++;
        if (p == q)
            break;
        avail = q - p;
        if (memcmp(p, "\033%-12345X", min(avail, 9)) == 0) {
            /* At least a partial match to a UEL. Everything up to
             * the start of the match is up for grabs. */
            return p - (pr->ptr+1);
        }
        p++;
    }

    return pr->limit - pr->ptr;
}

static tmsize_t tifsReadProc(thandle_t  tiff_,
                             void      *buf,
                             tmsize_t   size)
{
    tiff_interp_instance_t *tiff = (tiff_interp_instance_t *)tiff_;
    tmsize_t available = tiff->buffer_full - tiff->file_pos;
    if (available > size)
        available = size;

    memcpy(buf, &tiff->tiff_buffer[tiff->file_pos], available);
    tiff->file_pos += available;

    return available;
}

static tmsize_t tifsWriteProc(thandle_t  tiff_,
                              void      *buf,
                              tmsize_t   size)
{
    return 0;
}

static toff_t tifsSeekProc(thandle_t tiff_, toff_t offset, int whence)
{
    tiff_interp_instance_t *tiff = (tiff_interp_instance_t *)tiff_;
    size_t pos = tiff->file_pos;

    /* toff_t is unsigned, which kind of implies they'll never use
     * SEEK_CUR, or SEEK_END, which makes you wonder why they include
     * whence at all, but... */
    if (whence == 1) { /* SEEK_CURR */
        offset += pos;
    } else if (whence == 2) { /* SEEK_END */
        offset += tiff->buffer_full;
    }
    /* else assume SEEK_SET */

    /* Clamp (Don't check against 0 as toff_t is unsigned) */
    if (offset > tiff->buffer_full)
        offset = tiff->buffer_full;

    tiff->file_pos = offset;

    return offset;
}

static int tifsCloseProc(thandle_t tiff_)
{
    return 0;
}

static toff_t tifsSizeProc(thandle_t tiff_)
{
    tiff_interp_instance_t *tiff = (tiff_interp_instance_t *)tiff_;

    return tiff->buffer_full;
}

#if defined(SHARE_JPEG) && SHARE_JPEG==0
static void *gs_j_mem_alloc(j_common_ptr cinfo, size_t size)
{
    gs_memory_t *mem = (gs_memory_t *)(GET_CUST_MEM_DATA(cinfo)->priv);

    return(gs_alloc_bytes(mem, size, "JPEG allocation"));
}

static void gs_j_mem_free(j_common_ptr cinfo, void *object, size_t size)
{
    gs_memory_t *mem = (gs_memory_t *)(GET_CUST_MEM_DATA(cinfo)->priv);

    gs_free_object(mem, object, "JPEG free");
}

static long gs_j_mem_init (j_common_ptr cinfo)
{
    gs_memory_t *mem = (gs_memory_t *)(GET_CUST_MEM_DATA(cinfo)->priv);
    gs_memory_t *cmem = NULL;

    if (gs_memory_chunk_wrap(&(cmem), mem) < 0) {
        return (-1);
    }

    (void)jpeg_cust_mem_set_private(GET_CUST_MEM_DATA(cinfo), cmem);

    return 0;
}

static void gs_j_mem_term (j_common_ptr cinfo)
{
    gs_memory_t *cmem = (gs_memory_t *)(GET_CUST_MEM_DATA(cinfo)->priv);
    gs_memory_t *mem = gs_memory_chunk_target(cmem);

    gs_memory_chunk_release(cmem);

    (void)jpeg_cust_mem_set_private(GET_CUST_MEM_DATA(cinfo), mem);
}

static void *
tiff_jpeg_mem_callback(thandle_t tiff_)
{
    tiff_interp_instance_t *tiff = (tiff_interp_instance_t *)tiff_;

    (void)jpeg_cust_mem_init(&tiff->jmem, (void *)tiff->memory,
                             gs_j_mem_init, gs_j_mem_term, NULL,
                            gs_j_mem_alloc, gs_j_mem_free,
                            gs_j_mem_alloc, gs_j_mem_free, NULL);

    return &tiff->jmem;
}
#endif /* SHARE_JPEG == 0 */

static int
guess_pal_depth(int n, uint16_t *rmap, uint16_t *gmap, uint16_t *bmap)
{
    int i;
    for (i = 0; i < n; i++) {
        if (rmap[i] >= 256 || gmap[i] >= 256 || bmap[i] >= 256)
            return 16;
    }
    return 8;
}

static void
blend_alpha(tiff_interp_instance_t *tiff, size_t n, int nc, int planar)
{
    int i = tiff->raw_byte_width * nc;
    byte *p = tiff->samples;
    const byte *q = tiff->samples + i;

    switch (tiff->bpc)
    {
    case 1:
        p += i*8;
        do
        {
            byte a = *--q;
            *--p = ( a     & 1)*255;
            *--p = ((a>>1) & 1)*255;
            *--p = ((a>>2) & 1)*255;
            *--p = ((a>>3) & 1)*255;
            *--p = ((a>>4) & 1)*255;
            *--p = ((a>>5) & 1)*255;
            *--p = ((a>>6) & 1)*255;
            *--p = ((a>>7) & 1)*255;
        }
        while (--i);
        break;
    case 2:
        p += i*4;
        do
        {
            byte a = *--q;
            *--p = ( a     & 3)*0x55;
            *--p = ((a>>1) & 3)*0x55;
            *--p = ((a>>2) & 3)*0x55;
            *--p = ((a>>3) & 3)*0x55;
        }
        while (--i);
        break;
    case 4:
        p += i*2;
        do
        {
            byte a = *--q;
            *--p = ( a     & 15)*0x11;
            *--p = ((a>>1) & 15)*0x11;
            *--p = ((a>>2) & 15)*0x11;
            *--p = ((a>>3) & 15)*0x11;
        }
        while (--i);
        break;
    default:
        break;
    }

    nc--;
    p = tiff->samples;
    if (planar == PLANARCONFIG_CONTIG)
    {
        q = (const byte *)tiff->samples;
        while (n--) {
            byte a = q[nc];
            for (i = nc; i > 0; i--) {
                int c = *q++ * a + 255*(255-a);
                c += (c>>7);
                *p++ = c>>8;
            }
            q++;
        }
    }
    else
    {
        int next_comp = tiff->raw_byte_width;
        int alpha_offset = nc * next_comp;
        while (n--) {
            byte a = p[alpha_offset];
            for (i = nc; i > 0; i--) {
                int c = *p * a + 255*(255-a);
                c += (c>>7);
                *p = c>>8;
                p += next_comp;
            }
            p -= alpha_offset;
            p++;
        }
    }
}

/* Calulate (a*b*c+d) safely */
static uint32_t
safe_mla(const gs_memory_t *mem, int *code, uint32_t a, uint32_t b, uint32_t c, uint32_t d)
{
    /* UINT_MAX < b*a means overflow, but we can't calculate that... */
    if (UINT_MAX/b < a)
        goto fail;
    a *= b;
    if (UINT_MAX/c < a)
        goto fail;
    a *= c;
    if (UINT_MAX-c < d)
        goto fail;

    return a+d;

fail:
    emprintf(mem, "Numeric overflow!\n");
    *code = gs_error_rangecheck;

    return 0;
}

static size_t
size_mla(const gs_memory_t *mem, int *code, size_t a, size_t b, size_t c, size_t d)
{
    /* SIZE_MAX < b*a means overflow, but we can't calculate that... */
    if (SIZE_MAX/b < a)
        goto fail;
    a *= b;
    if (SIZE_MAX/c < a)
        goto fail;
    a *= c;
    if (SIZE_MAX-c < d)
        goto fail;

    return a+d;

fail:
    emprintf(mem, "Numeric overflow!\n");
    *code = gs_error_rangecheck;
    return 0;
}

static int
do_tiff_decode(tiff_interp_instance_t *tiff)
{
    int code = 0;
    int tx, ty;
    short planar;
    float f, scale;
    gs_color_space *cs;
    unsigned int used[GS_IMAGE_MAX_COMPONENTS];
    gs_string plane_data[GS_IMAGE_MAX_COMPONENTS];
    int invert = 0;
    int alpha = 0;

    TIFFGetField(tiff->handle, TIFFTAG_COMPRESSION, &tiff->compression);
    if (tiff->compression == COMPRESSION_JPEG) {
        TIFFSetField(tiff->handle, TIFFTAG_JPEGCOLORMODE, JPEGCOLORMODE_RGB);
    }
    TIFFGetField(tiff->handle, TIFFTAG_PHOTOMETRIC, &tiff->photometric);
    if (tiff->photometric == PHOTOMETRIC_LOGL ||
        tiff->photometric == PHOTOMETRIC_LOGLUV) {
        TIFFSetField(tiff->handle, TIFFTAG_SGILOGDATAFMT, SGILOGDATAFMT_8BIT);
    }

    TIFFGetField(tiff->handle, TIFFTAG_IMAGEWIDTH, &tiff->width);
    TIFFGetField(tiff->handle, TIFFTAG_IMAGELENGTH, &tiff->height);
    TIFFGetField(tiff->handle, TIFFTAG_TILEWIDTH, &tiff->tile_width);
    TIFFGetField(tiff->handle, TIFFTAG_TILELENGTH, &tiff->tile_height);
    TIFFGetField(tiff->handle, TIFFTAG_BITSPERSAMPLE, &tiff->bpc);
    TIFFGetField(tiff->handle, TIFFTAG_SAMPLESPERPIXEL, &tiff->num_comps);
    tiff->raw_num_comps = tiff->num_comps;
    TIFFGetField(tiff->handle, TIFFTAG_PLANARCONFIG, &planar);
    f = 0;
    TIFFGetField(tiff->handle, TIFFTAG_XRESOLUTION, &f);
    tiff->xresolution = (uint32_t)(f+0.5);
    f = 0;
    TIFFGetField(tiff->handle, TIFFTAG_YRESOLUTION, &f);
    tiff->yresolution = (uint32_t)(f+0.5);

    if (tiff->xresolution == 0)
        tiff->yresolution = tiff->xresolution;
    if (tiff->yresolution == 0)
        tiff->xresolution = tiff->yresolution;
    if (tiff->xresolution == 0)
        tiff->xresolution = tiff->yresolution = 72;
    if (tiff->width == 0 || tiff->height == 0 || tiff->bpc == 0 || tiff->num_comps == 0 ||
        !(planar == PLANARCONFIG_CONTIG || planar == PLANARCONFIG_SEPARATE)) {
        emprintf(tiff->memory, "Unsupported TIFF format\n");
        return gs_error_unknownerror;
    }

    tiff->tiled = TIFFIsTiled(tiff->handle);

    if (!tiff->tiled) {
        tiff->tile_width = tiff->width;
        tiff->tile_height = tiff->height;
    }

    if (tiff->tiled || planar == PLANARCONFIG_CONTIG) {
        tiff->byte_width = safe_mla(tiff->memory, &code, tiff->bpc, tiff->num_comps, tiff->tile_width, 7)>>3;
    } else {
        tiff->byte_width = safe_mla(tiff->memory, &code, tiff->bpc, 1, tiff->tile_width, 7)>>3;
    }
    if (code < 0) {
        emprintf(tiff->memory, "Unsupported: TIFF size overflow\n");
        goto fail_decode;
    }

    tiff->raw_byte_width = tiff->byte_width;
    if (tiff->photometric == PHOTOMETRIC_RGB && tiff->num_comps == 4)
    {
        /* RGBA, so alpha data */
        alpha = 1;
    }
    if (alpha && tiff->bpp < 8)
    {
        /* We need to expand the data to 8bpp to blend for alpha. */
        if (tiff->bpc != 1 && tiff->bpc != 2 && tiff->bpc != 4)
        {
            emprintf1(tiff->memory, "Unsupported: TIFF with alpha and bpc=%d\n", tiff->bpc);
            code = gs_error_unknownerror;
            goto fail_decode;
        }
        tiff->byte_width *= 8/tiff->bpc;
    }

    /* Allocate 'samples' to hold the raw samples values read from libtiff.
     * The exact size of this buffer depends on which of the multifarious
     * read routines we are using. (Tiled/RGBAImage/Scanlines) */
    if (tiff->compression == COMPRESSION_OJPEG ||
        tiff->photometric == PHOTOMETRIC_YCBCR) {
        size_t z = size_mla(tiff->memory, &code, sizeof(uint32_t), tiff->width, tiff->height, 0);
        if (code < 0) {
            emprintf(tiff->memory, "Unsupported: TIFF size overflow\n");
            goto fail_decode;
        }
        tiff->is_rgba = 1;
        tiff->samples = gs_alloc_bytes(tiff->memory, z, "tiff_image");
        tiff->tile_width = tiff->width;
        tiff->tile_height = tiff->height;
        tiff->byte_width = safe_mla(tiff->memory, &code, tiff->bpc, tiff->num_comps, tiff->tile_width, 7)>>3;
        if (code < 0) {
            emprintf(tiff->memory, "Unsupported: TIFF size overflow\n");
            goto fail_decode;
        }
    } else if (tiff->tiled) {
        tiff->samples = gs_alloc_bytes(tiff->memory, TIFFTileSize(tiff->handle), "tiff_tile");
    } else if (planar == PLANARCONFIG_SEPARATE) {
        tiff->samples = gs_alloc_bytes(tiff->memory, tiff->byte_width * tiff->num_comps, "tiff_scan");
    } else {
        tiff->samples = gs_alloc_bytes(tiff->memory, tiff->byte_width, "tiff_scan");
    }
    if (tiff->samples == NULL) {
        code = gs_error_VMerror;
        goto fail_decode;
    }
    tiff->proc_samples = tiff->samples;

    tiff->bpp = tiff->bpc * tiff->num_comps;
    switch(tiff->photometric) {
    case PHOTOMETRIC_MINISWHITE:
        invert = 1;
        /* Fall through */
    case PHOTOMETRIC_MINISBLACK:
        if (tiff->num_comps != 1) {
            emprintf1(tiff->memory, "Unsupported: TIFF with MINISBLACK with nc=%d\n", tiff->num_comps);
            code = gs_error_unknownerror;
            goto fail_decode;
        }
        break;
    case PHOTOMETRIC_RGB:
        if (tiff->num_comps == 4) {
            alpha = 1;
            tiff->num_comps = 3;
            tiff->bpp = tiff->bpp * 3/4;
            tiff->byte_width = tiff->byte_width * 3/4;
        } else if (tiff->num_comps != 3) {
            emprintf1(tiff->memory, "Unsupported: RGB TIFF nc=%d\n", tiff->num_comps);
            code = gs_error_unknownerror;
            goto fail_decode;
        }
        break;
    case PHOTOMETRIC_PALETTE:
    {
        uint16_t *rmap, *gmap, *bmap;
        int i, n = 1<<tiff->bpc;
        if (tiff->num_comps != 1) {
            emprintf1(tiff->memory, "Unsupported: Paletted TIFF with nc=%d\n", tiff->num_comps);
            code = gs_error_unknownerror;
            goto fail_decode;
        }
        if (tiff->bpc > 8) {
            emprintf1(tiff->memory, "Unsupported: Paletted TIFF with bpc=%d\n", tiff->bpc);
            code = gs_error_unknownerror;
            goto fail_decode;
        }
        if (!TIFFGetField(tiff->handle, TIFFTAG_COLORMAP, &rmap, &gmap, &bmap)) {
            emprintf(tiff->memory, "Unsupported: Paletted TIFF with bad palette\n");
            code = gs_error_unknownerror;
            goto fail_decode;
        }
        tiff->palette = gs_alloc_bytes(tiff->memory, 3*256, "palette");
        if (tiff->palette == NULL) {
            code = gs_error_VMerror;
            goto fail_decode;
        }
        memset(tiff->palette, 0, 3 * 256);
        if (guess_pal_depth(n, rmap, gmap, bmap) == 8) {
            for (i=0; i < n; i++) {
                tiff->palette[3*i+0] = rmap[i];
                tiff->palette[3*i+1] = gmap[i];
                tiff->palette[3*i+2] = bmap[i];
            }
        } else {
            for (i=0; i < n; i++) {
                tiff->palette[3*i+0] = rmap[i]*255/65535;
                tiff->palette[3*i+1] = gmap[i]*255/65535;
                tiff->palette[3*i+2] = bmap[i]*255/65535;
            }
        }
        tiff->pal_bpc = tiff->bpc;
        tiff->bpc = 8;
        tiff->num_comps = 3;
        tiff->raw_num_comps = 1;
        tiff->bpp = 24;
        tiff->byte_width = tiff->tile_width * 3;
        tiff->raw_byte_width = tiff->byte_width;
        /* Now we need to make a "proc_samples" area to store the
         * processed samples in. */
        if (tiff->is_rgba) {
            emprintf(tiff->memory, "Unsupported: Paletted TIFF with RGBA\n");
            code = gs_error_unknownerror;
            goto fail_decode;
        } else if (tiff->tiled) {
            size_t z = size_mla(tiff->memory, &code, tiff->tile_width, tiff->tile_height, 3, 0);
            if (code < 0) {
                emprintf(tiff->memory, "Unsupported: TIFF size overflow\n");
                goto fail_decode;
            }
            tiff->proc_samples = gs_alloc_bytes(tiff->memory, z, "tiff_tile");
        } else {
            size_t z = size_mla(tiff->memory, &code, tiff->tile_width, 1, 3, 0);
            tiff->proc_samples = gs_alloc_bytes(tiff->memory, z, "tiff_scan");
        }
        break;
    }
    case PHOTOMETRIC_MASK:
        if (tiff->num_comps != 1) {
            emprintf1(tiff->memory, "Unsupported: Mask TIFF with nc=%d\n", tiff->num_comps);
            code = gs_error_unknownerror;
            goto fail_decode;
        }
        break;
    case PHOTOMETRIC_SEPARATED:
        if (tiff->num_comps == 3 || tiff->num_comps == 4)
        {
            emprintf1(tiff->memory, "Unsupported: Separated TIFF with nc=%d\n", tiff->num_comps);
            break;
        }
    case PHOTOMETRIC_YCBCR:
    case PHOTOMETRIC_CIELAB:
    case PHOTOMETRIC_ICCLAB:
    case PHOTOMETRIC_ITULAB:
        if (tiff->num_comps != 3) {
            emprintf1(tiff->memory, "Unsupported: YUV/LAB TIFF with nc=%d\n", tiff->num_comps);
            code = gs_error_unknownerror;
            goto fail_decode;
        }
        break;
    case PHOTOMETRIC_CFA:
    default:
        emprintf(tiff->memory, "Unsupported TIFF\n");
        tiff->state = ii_state_flush;
        break;
    }
    switch(tiff->num_comps) {
    default:
    case 1:
        cs = tiff->gray;
        break;
    case 3:
        cs = tiff->rgb;
        break;
    case 4:
        cs = tiff->cmyk;
        break;
    }

    switch (tiff->bpc)
    {
    case 1:
    case 2:
    case 4:
    case 8:
    case 16:
        /* We can cope with all these. */
        break;
    default:
        emprintf1(tiff->memory, "Unsupported: TIFF with bpc=%d\n", tiff->bpc);
        code = gs_error_unknownerror;
        goto fail_decode;
    }

    /* Scale to fit, if too large. */
    scale = 1.0f;
    if (tiff->width * tiff->dev->HWResolution[0] > tiff->dev->width * tiff->xresolution)
        scale = ((float)tiff->dev->width * tiff->xresolution) / (tiff->width * tiff->dev->HWResolution[0]);
    if (scale * tiff->height * tiff->dev->HWResolution[1] > tiff->dev->height * tiff->yresolution)
        scale = ((float)tiff->dev->height * tiff->yresolution) / (tiff->height * tiff->dev->HWResolution[1]);

    code = gs_erasepage(tiff->pgs);
    if (code < 0)
        return code;

    for (ty = 0; ty < tiff->height; ty += tiff->tile_height) {
        for (tx = 0; tx < tiff->width; tx += tiff->tile_width) {
            int y, s;
            byte *row;
            float xext, xoffset, yext, yoffset;
            int tremx, tremy;

            tiff->penum = gs_image_enum_alloc(tiff->memory, "tiff_impl_process(penum)");
            if (tiff->penum == NULL)
                return_error(gs_error_VMerror);

            /* Centre - Extents and offsets are all calculated in points (1/72 of an inch) */
            xext = (((float)tiff->width - tx * 2) * 72 * scale / tiff->xresolution);
            xoffset = (tiff->dev->width * 72 / tiff->dev->HWResolution[0] - xext)/2;
            yext = (((float)tiff->height - ty * 2) * 72 * scale / tiff->yresolution);
            yoffset = (tiff->dev->height * 72 / tiff->dev->HWResolution[1] - yext)/2;

            gs_initmatrix(tiff->pgs);

            /* By default the ctm is set to:
             *   xres/72   0
             *   0         -yres/72
             *   0         dev->height * yres/72
             * i.e. it moves the origin from being top right to being bottom left.
             * We want to move it back, as without this, the image will be displayed
             * upside down.
             */
            code = gs_translate(tiff->pgs, 0.0, tiff->dev->height * 72 / tiff->dev->HWResolution[1]);
            if (code >= 0)
                code = gs_translate(tiff->pgs, xoffset, -yoffset);
            if (code >= 0)
                code = gs_scale(tiff->pgs, scale, -scale);
            if (code < 0)
               goto fail_decode;

            memset(&tiff->image, 0, sizeof(tiff->image));
            gs_image_t_init(&tiff->image, cs);
            tiff->image.BitsPerComponent = tiff->bpp/tiff->num_comps;
            if (alpha)
                tiff->image.BitsPerComponent = 8;
            tiff->image.Width = tiff->tile_width;
            tiff->image.Height = tiff->tile_height;

            tiff->image.ImageMatrix.xx = tiff->xresolution / 72.0f;
            tiff->image.ImageMatrix.yy = tiff->yresolution / 72.0f;
            if (invert) {
                tiff->image.Decode[0] = 1;
                tiff->image.Decode[1] = 0;
                tiff->image.Decode[2] = 1;
                tiff->image.Decode[3] = 0;
                tiff->image.Decode[4] = 1;
                tiff->image.Decode[5] = 0;
                tiff->image.Decode[6] = 1;
                tiff->image.Decode[7] = 0;
            }

            if (tiff->is_rgba) {
                size_t z = size_mla(tiff->memory, &code, tiff->tile_width, tiff->tile_height, 1, 0);
                if (code < 0)
                    goto fail_decode;
                if (TIFFReadRGBAImage(tiff->handle, tiff->width, tiff->height,
                                      (uint32_t *)tiff->samples, 0) == 0) {
                    code = gs_error_unknownerror;
                    goto fail_decode;
                }
                blend_alpha(tiff, z, 4, planar);
            } else if (tiff->tiled) {
                if (TIFFReadTile(tiff->handle, tiff->samples, tx, ty, 0, 0) == 0) {
                    code = gs_error_unknownerror;
                    goto fail_decode;
                }
            } else if (planar != PLANARCONFIG_CONTIG) {
                tiff->image.format = gs_image_format_component_planar;
            }

            if (!tiff->is_rgba && tiff->tiled) {
                if (tiff->palette) {
                    size_t n = size_mla(tiff->memory, &code, tiff->tile_width, tiff->tile_height, 1, 0);
                    byte *q = tiff->samples;
                    byte *p = tiff->proc_samples;
                    if (code < 0)
                        goto fail_decode;
                    while (n--) {
                        byte *v = &tiff->palette[3 * *q++];
                        p[0] = *v++;
                        p[1] = *v++;
                        p[2] = *v++;
                        p += 3;
                    }
                }
                if (alpha) {
                    blend_alpha(tiff, tiff->tile_width, tiff->num_comps, planar);
                }
            }

            code = gs_image_init(tiff->penum,
                                 &tiff->image,
                                 false,
                                 false,
                                 tiff->pgs);
            if (code < 0) {
                goto fail_decode;
            }

            tremx = tiff->width - tx;
            if (tremx > tiff->tile_width)
                tremx = tiff->tile_width;
            tremy = tiff->height - ty;
            if (tremy > tiff->tile_height)
                tremy = tiff->tile_height;
            {
                /* Make sure we won't overflow in the loop. */
                (void)size_mla(tiff->memory, &code, tiff->byte_width, tiff->tile_height, 1, 0);
                if (code < 0)
                    goto fail_decode;
            }
            for (y = 0; y < tremy; y++) {
                if (tiff->is_rgba) {
                    row = tiff->proc_samples + (size_t)tiff->byte_width * (tiff->tile_height-1-y);
                } else if (tiff->tiled) {
                    row = tiff->proc_samples + (size_t)tiff->byte_width * y;
                } else if (planar == PLANARCONFIG_CONTIG) {
                    row = tiff->proc_samples;
                    if (TIFFReadScanline(tiff->handle, tiff->samples, ty+y, 0) == 0) {
                        code = gs_error_unknownerror;
                        goto fail_decode;
                    }
                } else {
                    int span = tiff->raw_byte_width;
                    byte *in_row = tiff->samples;
                    row = tiff->proc_samples;
                    for (s = 0; s < tiff->raw_num_comps; s++) {
                         plane_data[s].data = row;
                         plane_data[s].size = span;
                         if (TIFFReadScanline(tiff->handle, in_row, ty+y, s) == 0) {
                             code = gs_error_unknownerror;
                             goto fail_decode;
                         }
                         row += span;
                         in_row += span;
                    }
                    row -= span * tiff->raw_num_comps;
                }

                if (tiff->bpc == 16)
                {
                    byte *p = row;
                    int n = tiff->tile_width * tiff->raw_num_comps;
                    while (n--)
                    {
                        byte b = p[0];
                        p[0] = p[1];
                        p[1] = b;
                        p += 2;
                    }
                }

                if (!tiff->tiled) {
                    if (tiff->palette) {
                        int n = tiff->tile_width;
                        const byte *q = tiff->samples;
                        byte *p = tiff->proc_samples;
                        switch (tiff->pal_bpc)
                        {
                        case 8:
                            while (n--) {
                                byte *v = &tiff->palette[3 * *q++];
                                p[0] = *v++;
                                p[1] = *v++;
                                p[2] = *v++;
                                p += 3;
                            }
                            break;
                        case 1:
                        {
                            int sh = 7;
                            while (n--) {
                                byte *v = &tiff->palette[3 * (((*q)>>sh) & 1)];
                                sh--;
                                if (sh < 0)
                                    sh = 7, q++;
                                p[0] = *v++;
                                p[1] = *v++;
                                p[2] = *v++;
                                p += 3;
                            }
                            break;
                        }
                        case 2:
                        {
                            int sh = 6;
                            while (n--) {
                                byte *v = &tiff->palette[3 * (((*q)>>sh) & 3)];
                                sh -= 2;
                                if (sh < 0)
                                    sh = 6, q++;
                                p[0] = *v++;
                                p[1] = *v++;
                                p[2] = *v++;
                                p += 3;
                            }
                            break;
                        }
                        case 4:
                        {
                            int sh = 4;
                            while (n--) {
                                byte *v = &tiff->palette[3 * (((*q)>>sh) & 15)];
                                sh ^= 4;
                                if (sh == 4)
                                    q++;
                                p[0] = *v++;
                                p[1] = *v++;
                                p[2] = *v++;
                                p += 3;
                            }
                            break;
                        }
                        }
                    }
                    if (alpha) {
                        blend_alpha(tiff, tiff->tile_width, tiff->raw_num_comps, planar);
                    }
                }

                if (tiff->image.format == gs_image_format_component_planar) {
                    code = gs_image_next_planes(tiff->penum, (gs_const_string *)&plane_data[0], used, false);
                } else {
                    code = gs_image_next(tiff->penum, row, tiff->byte_width, used);
                }
                if (code < 0) {
                    code = gs_error_unknownerror;
                    goto fail_decode;
                }
            }
            code = gs_image_cleanup_and_free_enum(tiff->penum, tiff->pgs);
            tiff->penum = NULL;
            if (code < 0)
                return code;
        }
    }
    tiff->state = ii_state_flush;
    (void)pl_finish_page(tiff->memory->gs_lib_ctx->top_of_system,
                         tiff->pgs, 1, true);
    return 0;

fail_decode:
    if (tiff->penum)
    {
        (void)gs_image_cleanup_and_free_enum(tiff->penum, tiff->pgs);
        tiff->penum = NULL;
    }
    tiff->state = ii_state_flush;

    return code;
}

static int
decode_all_tiffs(tiff_interp_instance_t *tiff)
{
    int code;

    tiff->nulldev = gs_currentdevice(tiff->pgs);
    rc_increment(tiff->nulldev);
    code = gs_setdevice_no_erase(tiff->pgs, tiff->dev);
    if (code < 0)
        return code;

    do
    {
        code = do_tiff_decode(tiff);
        if (code < 0)
            return code;
    }
    while (TIFFReadDirectory(tiff->handle));

    return 0;
}


static int
do_impl_process(pl_interp_implementation_t * impl, stream_cursor_read * pr, int eof)
{
    tiff_interp_instance_t *tiff = (tiff_interp_instance_t *)impl->interp_client_data;
    int code = 0;
    ii_state ostate = (ii_state)-1;
    int bytes_left = 0;

    while (tiff->state != ostate || pr->limit - pr->ptr != bytes_left)
    {
        ostate = tiff->state;
        bytes_left = pr->limit - pr->ptr;
        switch(tiff->state)
        {
        case ii_state_identifying:
        {
            const byte *hdr;
            /* Try and get us 4 bytes */
            code = ensure_bytes(tiff, pr, 4);
            if (code < 0)
                return code;
            hdr = pr->ptr+1;
            if (hdr[0] == 'I' && hdr[1] == 'I' && (hdr[2] == 42 || hdr[2] == 43) && hdr[3] == 0) {
                tiff->state = ii_state_tiff;
                break;
            }
            if (hdr[0] == 'M' && hdr[1] == 'M' && hdr[2] == 0 && (hdr[3] == 42 || hdr[3] == 43)) {
                tiff->state = ii_state_tiff;
                break;
            }
            tiff->state = ii_state_flush;
            break;
        }
        case ii_state_tiff:
        {
            /* Gather data into a buffer */
            int bytes = bytes_until_uel(pr);

            if (bytes == 0 && pr->limit - pr->ptr > 9) {
                /* No bytes until UEL, and there is space for a UEL in the buffer */
                tiff->state = ii_state_tiff_decode;
                tiff->file_pos = 0;
                break;
            }
            if (bytes == 0 && eof) {
                /* No bytes until UEL, and we are at eof */
                tiff->state = ii_state_tiff_decode;
                tiff->file_pos = 0;
                break;
            }

            if (tiff->buffer_full + bytes > tiff->buffer_max) {
                /* Need to expand our buffer */
                size_t proposed = tiff->buffer_full*2;
                if (proposed == 0)
                    proposed = 32768;
                while (proposed < tiff->buffer_full + bytes)
                    proposed *= 2;

                if (tiff->tiff_buffer == NULL) {
                    tiff->tiff_buffer = gs_alloc_bytes(tiff->memory, proposed, "tiff_buffer");
                    if (tiff->tiff_buffer == NULL) {
                        tiff->state = ii_state_flush;
                        break;
                    }
                } else {
                    void *new_buf = gs_resize_object(tiff->memory, tiff->tiff_buffer, proposed, "tiff_buffer");
                    if (new_buf == NULL) {
                        tiff->state = ii_state_flush;
                        break;
                    }
                    tiff->tiff_buffer = new_buf;
                }
                tiff->buffer_max = proposed;
            }

            memcpy(&tiff->tiff_buffer[tiff->buffer_full], pr->ptr+1, bytes);
            tiff->buffer_full += bytes;
            pr->ptr += bytes;
            break;
        }
        case ii_state_tiff_decode:
        {
            tiff->handle = TIFFClientOpen("dummy", "rm",
                                          (thandle_t)tiff,
                                          tifsReadProc,
                                          tifsWriteProc,
                                          tifsSeekProc,
                                          tifsCloseProc,
                                          tifsSizeProc,
                                          NULL,
                                          NULL);
            if (tiff->handle == NULL) {
                tiff->state = ii_state_flush;
                break;
            }

#if defined(SHARE_JPEG) && SHARE_JPEG==0
            TIFFSetJpegMemFunction(tiff->handle,
                                   &tiff_jpeg_mem_callback);
#endif

            code = decode_all_tiffs(tiff);
            if (code < 0)
            {
                tiff->state = ii_state_flush;
                break;
            }
            break;
        }
        default:
        case ii_state_flush:

            if (tiff->handle) {
                TIFFClose(tiff->handle);
                tiff->handle = NULL;
            }

            if (tiff->penum) {
                (void)gs_image_cleanup_and_free_enum(tiff->penum, tiff->pgs);
                tiff->penum = NULL;
            }

            if (tiff->proc_samples && tiff->proc_samples != tiff->samples) {
                gs_free_object(tiff->memory, tiff->proc_samples, "tiff_impl_process(samples)");
                tiff->proc_samples = NULL;
            }

            if (tiff->samples) {
                gs_free_object(tiff->memory, tiff->samples, "tiff_impl_process(samples)");
                tiff->samples = NULL;
                tiff->proc_samples = NULL;
            }

            if (tiff->palette) {
                gs_free_object(tiff->memory, tiff->palette, "tiff_impl_process(samples)");
                tiff->palette = NULL;
            }

            if (tiff->tiff_buffer) {
                gs_free_object(tiff->memory, tiff->tiff_buffer, "tiff_impl_process(tiff_buffer)");
                tiff->tiff_buffer = NULL;
                tiff->buffer_max = 0;
                tiff->buffer_full = 0;
            }
            /* We want to bin any data we get up to, but not including
             * a UEL. */
            return flush_to_uel(pr);
        }
    }

    return code;
}

static int
tiff_impl_process(pl_interp_implementation_t * impl, stream_cursor_read * pr) {
    return do_impl_process(impl, pr, 0);
}

static int
tiff_impl_process_end(pl_interp_implementation_t * impl)
{
    return 0;
}

/* Not implemented */
static int
tiff_impl_flush_to_eoj(pl_interp_implementation_t *impl, stream_cursor_read *cursor)
{
    const byte *p = cursor->ptr;
    const byte *rlimit = cursor->limit;

    /* Skip to, but leave UEL in buffer for PJL to find later */
    for (; p < rlimit; ++p)
        if (p[1] == '\033') {
            uint avail = rlimit - p;

            if (memcmp(p + 1, "\033%-12345X", min(avail, 9)))
                continue;
            if (avail < 9)
                break;
            cursor->ptr = p;
            return 1;           /* found eoj */
        }
    cursor->ptr = p;
    return 0;                   /* need more */
}

/* Parser action for end-of-file */
static int
tiff_impl_process_eof(pl_interp_implementation_t *impl)
{
    stream_cursor_read r;

    r.ptr = NULL;
    r.limit = NULL;
    return do_impl_process(impl, &r, 1);
}

/* Report any errors after running a job */
static int
tiff_impl_report_errors(pl_interp_implementation_t *impl,          /* interp instance to wrap up job in */
                        int                         code,          /* prev termination status */
                        long                        file_position, /* file position of error, -1 if unknown */
                        bool                        force_to_cout  /* force errors to cout */
)
{
    return 0;
}

/* Wrap up interp instance after a "job" */
static int
tiff_impl_dnit_job(pl_interp_implementation_t *impl)
{
    tiff_interp_instance_t *tiff = (tiff_interp_instance_t *)impl->interp_client_data;

    if (tiff->nulldev) {
        int code = gs_setdevice(tiff->pgs, tiff->nulldev);
        tiff->dev = NULL;
        rc_decrement(tiff->nulldev, "tiff_impl_dnit_job(nulldevice)");
        tiff->nulldev = NULL;
        return code;
    }
    return 0;
}

/* Parser implementation descriptor */
const pl_interp_implementation_t tiff_implementation = {
  tiff_impl_characteristics,
  tiff_impl_allocate_interp_instance,
  tiff_impl_get_device_memory,
  NULL, /* tiff_impl_set_param */
  NULL, /* tiff_impl_add_path */
  NULL, /* tiff_impl_post_args_init */
  tiff_impl_init_job,
  NULL, /* tiff_impl_run_prefix_commands */
  NULL, /* tiff_impl_process_file */
  tiff_impl_process_begin,
  tiff_impl_process,
  tiff_impl_process_end,
  tiff_impl_flush_to_eoj,
  tiff_impl_process_eof,
  tiff_impl_report_errors,
  tiff_impl_dnit_job,
  tiff_impl_deallocate_interp_instance,
  NULL, /* tiff_impl_reset */
  NULL  /* interp_client_data */
};
#endif /* TIFF_INCLUDED */