/* Copyright (C) 2001-2021 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.,  1305 Grant Avenue - Suite 200, Novato,

   CA 94945, U.S.A., +1(415)492-9861, for further information.

*/

/* CIE color rendering operators */

#include "math_.h"

#include "ghost.h"

#include "oper.h"

#include "gsstruct.h"

#include "gscspace.h"

#include "gscolor2.h"

#include "gscrd.h"

#include "gscrdp.h"

#include "estack.h"

#include "ialloc.h"

#include "idict.h"

#include "idparam.h"

#include "igstate.h"

#include "icie.h"

#include "iparam.h"

#include "ivmspace.h"

#include "store.h"    /* for make_null */

/* Forward references */

static int zcrd1_proc_params(const gs_memory_t *mem, os_ptr op, ref_cie_render_procs * pcprocs);

static int zcrd1_params(os_ptr op, gs_cie_render * pcrd,

                         ref_cie_render_procs * pcprocs, gs_memory_t * mem);

/* - currentcolorrendering <dict> */

static int

zcurrentcolorrendering(i_ctx_t *i_ctx_p)

{

    os_ptr op = osp;

    push(1);

    *op = istate->colorrendering.dict;

    return 0;

}

/* <dict> .buildcolorrendering1 <crd> */

static int

zbuildcolorrendering1(i_ctx_t *i_ctx_p)

{

    os_ptr op = osp;

    gs_memory_t *mem = gs_gstate_memory(igs);

    int code;

    es_ptr ep = esp;

    gs_cie_render *pcrd;

    ref_cie_render_procs procs;

    check_read_type(*op, t_dictionary);

    check_dict_read(*op);

    code = gs_cie_render1_build(&pcrd, mem, ".buildcolorrendering1");

    if (code < 0)

        return code;

    code = zcrd1_params(op, pcrd, &procs, mem);

    if (code < 0 ) {

        rc_free_struct(pcrd, ".buildcolorrendering1");

        esp = ep;

        return code;

    }

    /****** FIX refct ******/

    /*rc_decrement(pcrd, ".buildcolorrendering1"); *//* build sets rc = 1 */

    istate->colorrendering.dict = *op;

    make_istruct_new(op, a_readonly, pcrd);

    return (esp == ep ? 0 : o_push_estack);

}

/* <dict> .builddevicecolorrendering1 <crd> */

static int

zbuilddevicecolorrendering1(i_ctx_t *i_ctx_p)

{

    os_ptr op = osp;

    gs_memory_t *mem = gs_gstate_memory(igs);

    dict_param_list list;

    gs_cie_render *pcrd = 0;

    int code;

    check_type(*op, t_dictionary);

    code = dict_param_list_read(&list, op, NULL, false, iimemory);

    if (code < 0)

        return code;

    code = gs_cie_render1_build(&pcrd, mem, ".builddevicecolorrendering1");

    if (code >= 0) {

        code = param_get_cie_render1(pcrd, (gs_param_list *) & list,

                                     gs_currentdevice(igs));

        if (code >= 0) {

            /****** FIX refct ******/

            /*rc_decrement(pcrd, ".builddevicecolorrendering1"); *//* build sets rc = 1 */

        }

    }

    iparam_list_release(&list);

    if (code < 0) {

        rc_free_struct(pcrd, ".builddevicecolorrendering1");

        return code;

    }

    istate->colorrendering.dict = *op;

    make_istruct_new(op, a_readonly, pcrd);

    return 0;

}

/* <dict> <crd> .setcolorrendering1 - */

static int

zsetcolorrendering1(i_ctx_t *i_ctx_p)

{

    os_ptr op = osp;

    es_ptr ep = esp;

    ref_cie_render_procs procs;

    int code;

    check_type(op[-1], t_dictionary);

    check_stype(*op, st_cie_render1);

    code = zcrd1_proc_params(imemory, op - 1, &procs);

    if (code < 0)

        return code;

    code = gs_setcolorrendering(igs, r_ptr(op, gs_cie_render));

    if (code < 0)

        return code;

    if (gs_cie_cs_common(igs) != 0 &&

        (code = cie_cache_joint(i_ctx_p, &procs, gs_cie_cs_common(igs), igs)) < 0

        )

        return code;

    istate->colorrendering.dict = op[-1];

    istate->colorrendering.procs = procs;

    pop(2);

    return (esp == ep ? 0 : o_push_estack);

}

/* <dict> <crd> .setdevicecolorrendering1 - */

static int

zsetdevicecolorrendering1(i_ctx_t *i_ctx_p)

{

    os_ptr op = osp;

    int code;

    ref_cie_render_procs procs;

    check_type(op[-1], t_dictionary);

    check_stype(*op, st_cie_render1);

    code = gs_setcolorrendering(igs, r_ptr(op, gs_cie_render));

    if (code < 0)

        return code;

    refset_null((ref *)&procs, sizeof(procs) / sizeof(ref));

    if (gs_cie_cs_common(igs) != 0 &&

        (code = cie_cache_joint(i_ctx_p, &procs, gs_cie_cs_common(igs), igs)) < 0

        )

        return code;

    istate->colorrendering.dict = op[-1];

    refset_null((ref *)&istate->colorrendering.procs,

                sizeof(istate->colorrendering.procs) / sizeof(ref));

    pop(2);

    return 0;

}

/* Get ColorRenderingType 1 procedures from the PostScript dictionary. */

static int

zcrd1_proc_params(const gs_memory_t *mem,

                  os_ptr op, ref_cie_render_procs * pcprocs)

{

    int code;

    ref *pRT;

    code = dict_proc3_param(mem, op, "EncodeLMN", &pcprocs->EncodeLMN);

    if (code < 0)

        return code;

    code = dict_proc3_param(mem, op, "EncodeABC", &pcprocs->EncodeABC);

    if (code < 0)

        return code;

    code = dict_proc3_param(mem, op, "TransformPQR", &pcprocs->TransformPQR);

    if (code < 0)

        return code;

    if (code == 1)

        return gs_note_error(gs_error_undefined);

    if (dict_find_string(op, "RenderTable", &pRT) > 0) {

        const ref *prte;

        int size;

        int i;

        check_read_type(*pRT, t_array);

        size = r_size(pRT);

        if (size < 5)

            return_error(gs_error_rangecheck);

        prte = pRT->value.const_refs;

        for (i = 5; i < size; i++)

            check_proc_only(prte[i]);

        make_const_array(&pcprocs->RenderTableT, a_readonly | r_space(pRT),

                         size - 5, prte + 5);

    } else

        make_null(&pcprocs->RenderTableT);

    return 0;

}

/* Get ColorRenderingType 1 parameters from the PostScript dictionary. */

static int

zcrd1_params(os_ptr op, gs_cie_render * pcrd,

             ref_cie_render_procs * pcprocs, gs_memory_t * mem)

{

    int code;

    int ignore;

    gx_color_lookup_table *const prtl = &pcrd->RenderTable.lookup;

    ref *pRT;

    if ((code = dict_int_param(op, "ColorRenderingType", 1, 1, 0, &ignore)) < 0)

        return code;

    if ((code = zcrd1_proc_params(mem, op, pcprocs)) < 0)

        return code;

    if ((code = dict_matrix3_param(mem, op, "MatrixLMN", &pcrd->MatrixLMN)) < 0)

        return code;

    if ((code = dict_range3_param(mem, op, "RangeLMN", &pcrd->RangeLMN)) < 0)

        return code;

    if ((code = dict_matrix3_param(mem, op, "MatrixABC", &pcrd->MatrixABC)) < 0)

        return code;

    if ((code = dict_range3_param(mem, op, "RangeABC", &pcrd->RangeABC)) < 0)

        return code;

    if ((code = cie_points_param(mem, op, &pcrd->points)) < 0)

        return code;

    if ((code = dict_matrix3_param(mem, op, "MatrixPQR", &pcrd->MatrixPQR)) < 0)

        return code;

    if ((code = dict_range3_param(mem,op, "RangePQR", &pcrd->RangePQR)) < 0)

        return code;

    if (dict_find_string(op, "RenderTable", &pRT) > 0) {

        const ref *prte;

        check_read_type(*pRT, t_array);

        prte = pRT->value.const_refs;

        /* Finish unpacking and checking the RenderTable parameter. */

        check_type_only(prte[4], t_integer);

        if (!(prte[4].value.intval == 3 || prte[4].value.intval == 4))

            return_error(gs_error_rangecheck);

        prtl->n = 3;

        prtl->m = prte[4].value.intval;

        if (r_size(pRT) != prtl->m + 5)

            return_error(gs_error_rangecheck);

        code = cie_table_param(pRT, prtl, mem);

        if (code < 0)

            return code;

    } else {

        prtl->table = 0;

    }

    pcrd->EncodeLMN = Encode_default;

    pcrd->EncodeABC = Encode_default;

    pcrd->TransformPQR = TransformPQR_default;

    pcrd->RenderTable.T = RenderTableT_default;

    return 0;

}

/* ------ Internal procedures ------ */

/* Load the joint caches. */

static int

    cie_exec_tpqr(i_ctx_t *),

    cie_post_exec_tpqr(i_ctx_t *),

    cie_tpqr_finish(i_ctx_t *);

int

cie_cache_joint(i_ctx_t *i_ctx_p, const ref_cie_render_procs * pcrprocs,

                const gs_cie_common *pcie, gs_gstate * pgs)

{

    const gs_cie_render *pcrd = gs_currentcolorrendering(pgs);

    gx_cie_joint_caches *pjc = gx_unshare_cie_caches(pgs);

    gs_ref_memory_t *imem = (gs_ref_memory_t *) gs_gstate_memory(pgs);

    ref pqr_procs;

    uint space;

    int code;

    int i;

    if (pcrd == 0)   /* cache is not set up yet */

        return 0;

    if (pjc == 0)   /* must already be allocated */

        return_error(gs_error_VMerror);

    if (r_has_type(&pcrprocs->TransformPQR, t_null)) {

        /*

         * This CRD came from a driver, not from a PostScript dictionary.

         * Resample TransformPQR in C code.

         */

        return gs_cie_cs_complete(pgs, true);

    }

    gs_cie_compute_points_sd(pjc, pcie, pcrd);

    code = ialloc_ref_array(&pqr_procs, a_readonly, 3 * (1 + 4 + 4 * 6),

                            "cie_cache_common");

    if (code < 0)

        return code;

    /* When we're done, deallocate the procs and complete the caches. */

    check_estack(3);

    code = cie_cache_push_finish(i_ctx_p, cie_tpqr_finish, imem, pgs);

    if (code < 0)

        return code;

    *++esp = pqr_procs;

    space = r_space(&pqr_procs);

    for (i = 0; i < 3; i++) {

        ref *p = pqr_procs.value.refs + 3 + (4 + 4 * 6) * i;

        const float *ppt = (float *)&pjc->points_sd;

        int j;

        make_array(pqr_procs.value.refs + i, a_readonly | a_executable | space,

                   4, p);

        make_array(p, a_readonly | space, 4 * 6, p + 4);

        p[1] = pcrprocs->TransformPQR.value.refs[i];

        make_oper(p + 2, 0, cie_exec_tpqr);

        make_oper(p + 3, 0, cie_post_exec_tpqr);

        for (j = 0, p += 4; j < 4 * 6; j++, p++, ppt++)

            make_real(p, *ppt);

    }

    return cie_prepare_cache3(i_ctx_p, &pcrd->RangePQR,

                              pqr_procs.value.const_refs,

                              pjc->TransformPQR.caches,

                              pjc, imem, "Transform.PQR");

}

/* Private operator to shuffle arguments for the TransformPQR procedure: */

/* v [ws wd bs bd] proc -> -mark- ws wd bs bd v proc + exec */

static int

cie_exec_tpqr(i_ctx_t *i_ctx_p)

{

    os_ptr op = osp;

    const ref *ppt = op[-1].value.const_refs;

    uint space = r_space(op - 1);

    int i;

    check_op(3);

    push(4);

    *op = op[-4];   /* proc */

    op[-1] = op[-6];    /* v */

    for (i = 0; i < 4; i++)

        make_const_array(op - 5 + i, a_readonly | space,

                         6, ppt + i * 6);

    make_mark(op - 6);

    return zexec(i_ctx_p);

}

/* Remove extraneous values from the stack after executing */

/* the TransformPQR procedure.  -mark- ... v -> v */

static int

cie_post_exec_tpqr(i_ctx_t *i_ctx_p)

{

    os_ptr op = osp;

    uint count = ref_stack_counttomark(&o_stack);

    ref vref;

    if (count < 2)

        return_error(gs_error_unmatchedmark);

    vref = *op;

    ref_stack_pop(&o_stack, count - 1);

    *osp = vref;

    return 0;

}

/* Free the procs array and complete the joint caches. */

static int

cie_tpqr_finish(i_ctx_t *i_ctx_p)

{

    os_ptr op = osp;

    gs_gstate *pgs = r_ptr(op, gs_gstate);

    gs_cie_render *pcrd =

        (gs_cie_render *)gs_currentcolorrendering(pgs);  /* break const */

    int code;

    ifree_ref_array(op - 1, "cie_tpqr_finish");

    pcrd->TransformPQR = TransformPQR_from_cache;

    code = gs_cie_cs_complete(pgs, false);

    pop(2);

    return code;

}

/* Ws Bs Wd Bd Ps .transformPQR_scale_wb[012] Pd

   The default TransformPQR procedure is implemented in C, rather than

   PostScript, as a speed optimization.

   This TransformPQR implements a relative colorimetric intent by scaling

   the XYZ values relative to the white and black points.

*/

static int

ztpqr_scale_wb_common(i_ctx_t *i_ctx_p, int idx)

{

    os_ptr op = osp;

    double a[4], Ps; /* a[0] = ws, a[1] = bs, a[2] = wd, a[3] = bd */

    double result;

    int code;

    int i;

    code = real_param(op, &Ps);

    if (code < 0) return code;

    for (i = 0; i < 4; i++) {

        ref tmp;

        code = array_get(imemory, op - 4 + i, idx, &tmp);

        if (code >= 0)

            code = real_param(&tmp, &a[i]);

        if (code < 0) return code;

    }

    if (a[0] == a[1])

        return_error(gs_error_undefinedresult);

    result = a[3] + (a[2] - a[3]) * (Ps - a[1]) / (a[0] - a[1]);

    make_real(op - 4, result);

    pop(4);

    return 0;

}

/* Ws Bs Wd Bd Ps .TransformPQR_scale_wb0 Pd */

static int

ztpqr_scale_wb0(i_ctx_t *i_ctx_p)

{

    return ztpqr_scale_wb_common(i_ctx_p, 3);

}

/* Ws Bs Wd Bd Ps .TransformPQR_scale_wb2 Pd */

static int

ztpqr_scale_wb1(i_ctx_t *i_ctx_p)

{

    return ztpqr_scale_wb_common(i_ctx_p, 4);

}

/* Ws Bs Wd Bd Ps .TransformPQR_scale_wb2 Pd */

static int

ztpqr_scale_wb2(i_ctx_t *i_ctx_p)

{

    return ztpqr_scale_wb_common(i_ctx_p, 5);

}

/* ------ Initialization procedure ------ */

const op_def zcrd_l2_op_defs[] =

{

    op_def_begin_level2(),

    {"0currentcolorrendering", zcurrentcolorrendering},

    {"2.setcolorrendering1", zsetcolorrendering1},

    {"2.setdevicecolorrendering1", zsetdevicecolorrendering1},

    {"1.buildcolorrendering1", zbuildcolorrendering1},

    {"1.builddevicecolorrendering1", zbuilddevicecolorrendering1},

                /* Internal "operators" */

    {"3%cie_exec_tpqr", cie_exec_tpqr},

    {"2%cie_post_exec_tpqr", cie_post_exec_tpqr},

    {"1%cie_tpqr_finish", cie_tpqr_finish},

    {"5.TransformPQR_scale_WB0", ztpqr_scale_wb0},

    {"5.TransformPQR_scale_WB1", ztpqr_scale_wb1},

    {"5.TransformPQR_scale_WB2", ztpqr_scale_wb2},

    op_def_end(0)

};