home *** CD-ROM | disk | FTP | other *** search
/ PC World Komputer 1996 February / PCWK0296.iso / sharewar / dos / program / gs300sr1 / gs300sr1.exe / ZCIE.C < prev    next >
Encoding:
C/C++ Source or Header  |  1994-07-27  |  20.1 KB  |  585 lines
  1. /* Copyright (C) 1992, 1993, 1994 Aladdin Enterprises.  All rights reserved.
  2.   
  3.   This file is part of Aladdin Ghostscript.
  4.   
  5.   Aladdin Ghostscript is distributed with NO WARRANTY OF ANY KIND.  No author
  6.   or distributor accepts any responsibility for the consequences of using it,
  7.   or for whether it serves any particular purpose or works at all, unless he
  8.   or she says so in writing.  Refer to the Aladdin Ghostscript Free Public
  9.   License (the "License") for full details.
  10.   
  11.   Every copy of Aladdin Ghostscript must include a copy of the License,
  12.   normally in a plain ASCII text file named PUBLIC.  The License grants you
  13.   the right to copy, modify and redistribute Aladdin Ghostscript, but only
  14.   under certain conditions described in the License.  Among other things, the
  15.   License requires that the copyright notice and this notice be preserved on
  16.   all copies.
  17. */
  18.  
  19. /* zcie.c */
  20. /* CIE color operators */
  21. #include "ghost.h"
  22. #include "errors.h"
  23. #include "oper.h"
  24. #include "gsstruct.h"
  25. #include "gscspace.h"
  26. #include "gscolor2.h"
  27. #include "gscie.h"
  28. #include "estack.h"
  29. #include "ialloc.h"
  30. #include "idict.h"
  31. #include "idparam.h"
  32. #include "igstate.h"
  33. #include "isave.h"
  34. #include "store.h"        /* for make_null */
  35.  
  36. /* There are actually only two CIE-specific operators, */
  37. /* but CIE color dictionaries are so complex that */
  38. /* we handle the CIE case of setcolorspace here as well. */
  39.  
  40. /* Forward references */
  41. private int cache_colorrendering(P3(gs_cie_render *,
  42.   const ref_cie_render_procs *, gs_state *));
  43. private int cache_common(P4(gs_cie_common *, const ref_cie_procs *,
  44.   const ref_cie_render_procs *, gs_state *));
  45.  
  46. /* Allocator structure types for CIE structures */
  47. private_st_cie_abc();
  48. private_st_cie_a();
  49. private_st_cie_render();
  50.  
  51. /* Empty procedures */
  52. static ref empty_procs[3];
  53.  
  54. /* Original CIE color space types */
  55. extern const gs_color_space_type
  56.     gs_color_space_type_CIEBasedABC,
  57.     gs_color_space_type_CIEBasedA;
  58. /* Redefined CIE color space types (that load the cache when installed) */
  59. gs_color_space_type
  60.     cs_type_zCIEBasedABC,
  61.     cs_type_zCIEBasedA;
  62. private cs_proc_install_cspace(cs_install_zCIEBasedABC);
  63. private cs_proc_install_cspace(cs_install_zCIEBasedA);
  64.  
  65. /* Initialization */
  66. private void
  67. zcie_init(void)
  68. {
  69.     /* Make the null (default) transformation procedures. */
  70.     make_const_array(&empty_procs[0], a_readonly + a_executable, 0, NULL);
  71.     make_const_array(&empty_procs[1], a_readonly + a_executable, 0, NULL);
  72.     make_const_array(&empty_procs[2], a_readonly + a_executable, 0, NULL);
  73.  
  74.     /* Create the modified color space types. */
  75.     cs_type_zCIEBasedABC = gs_color_space_type_CIEBasedABC;
  76.     cs_type_zCIEBasedABC.install_cspace = cs_install_zCIEBasedABC;
  77.     cs_type_zCIEBasedA = gs_color_space_type_CIEBasedA;
  78.     cs_type_zCIEBasedA.install_cspace = cs_install_zCIEBasedA;
  79.  
  80. }
  81.  
  82. /* ------ CIE setcolorspace ------ */
  83.  
  84. /* Get a 3-element range parameter from a dictionary. */
  85. #define dict_range3_param(op, kstr, prange)\
  86.     dict_float_array_param(op, kstr, 6, (float *)prange, (float *)&Range3_default)
  87. #define range3_ok 6
  88.  
  89. /* Get a 3x3 matrix parameter from a dictionary. */
  90. #define dict_matrix3_param(op, kstr, pmat)\
  91.     dict_float_array_param(op, kstr, 9, (float *)pmat, (float *)&Matrix3_default)
  92. #define matrix3_ok 9
  93.  
  94. /* Get an array of procedures from a dictionary. */
  95. /* We know count <= 3. */
  96. private int
  97. dict_proc_array_param(const ref *pdict, const char _ds *kstr,
  98.   uint count, ref *pparray)
  99. {    ref *pvalue;
  100.     if ( dict_find_string(pdict, kstr, &pvalue) > 0 )
  101.     {    uint i;
  102.         check_array(*pvalue);
  103.         if ( r_size(pvalue) != count )
  104.             return_error(e_rangecheck);
  105.         for ( i = 0; i < count; i++ )
  106.         {    ref proc;
  107.             array_get(pvalue, (long)i, &proc);
  108.             check_proc(proc);
  109.         }
  110.         *pparray = *pvalue;
  111.     }
  112.     else
  113.         make_const_array(pparray, a_readonly | a_foreign,
  114.                  count, &empty_procs[0]);
  115.     return 0;
  116. }
  117.  
  118. /* Get 3 procedures from a dictionary. */
  119. #define dict_proc3_param(op, kstr, pparray)\
  120.     dict_proc_array_param(op, kstr, 3, pparray)
  121.         
  122. /* Shared code for getting WhitePoint and BlackPoint values. */
  123. private int
  124. cie_points_param(const ref *pdref, gs_cie_wb *pwb)
  125. {    int code;
  126.     if ( (code = dict_float_array_param(pdref, "WhitePoint", 3, (float *)&pwb->WhitePoint, NULL)) != 3 ||
  127.          (code = dict_float_array_param(pdref, "BlackPoint", 3, (float *)&pwb->BlackPoint, (float *)&BlackPoint_default)) != 3
  128.        )
  129.         return (code < 0 ? code : e_rangecheck);
  130.     if ( pwb->WhitePoint.u <= 0 ||
  131.          pwb->WhitePoint.v != 1 ||
  132.          pwb->WhitePoint.w <= 0 ||
  133.          pwb->BlackPoint.u < 0 ||
  134.          pwb->BlackPoint.v < 0 ||
  135.          pwb->BlackPoint.w < 0
  136.        )
  137.         return_error(e_rangecheck);
  138.     return 0;
  139. }
  140.  
  141. /* Common code for the CIEBasedA[BC] cases of setcolorspace. */
  142. private int
  143. cie_lmnp_param(const ref *pdref, gs_cie_common *pcie, ref_cie_procs *pcprocs)
  144. {    int code;
  145.     if ( (code = dict_range3_param(pdref, "RangeLMN", &pcie->RangeLMN)) != range3_ok ||
  146.          (code = dict_proc3_param(pdref, "DecodeLMN", &pcprocs->DecodeLMN)) < 0 ||
  147.          (code = dict_matrix3_param(pdref, "MatrixLMN", &pcie->MatrixLMN)) != matrix3_ok ||
  148.          (code = cie_points_param(pdref, &pcie->points)) < 0
  149.        )
  150.         return (code < 0 ? code : e_rangecheck);
  151.     pcie->DecodeLMN = DecodeLMN_default;
  152.     return 0;
  153. }
  154.  
  155. /* Get the parameters of a CIEBasedABC color space. */
  156. /* This doesn't actually implement setcolorspace, */
  157. /* since it has to be usable for the base color space */
  158. /* of Separation, Indexed, and Pattern spaces as well. */
  159. int
  160. zcolorspace_CIEBasedABC(const ref *pdref, gs_color_space *pcs,
  161.   ref_cie_procs *pcprocs)
  162. {    gs_memory_t *mem = gs_state_memory(igs);
  163.     gs_cie_abc *pcie;
  164.     int code;
  165.     check_type(*pdref, t_dictionary);
  166.     check_dict_read(*pdref);
  167.     rc_alloc_struct_0(pcie, gs_cie_abc, &st_cie_abc, mem,
  168.               return_error(e_VMerror),
  169.               "setcolorspace(CIEBasedABC)");
  170.     if ( (code = dict_range3_param(pdref, "RangeABC", &pcie->RangeABC)) != range3_ok ||
  171.          (code = dict_proc3_param(pdref, "DecodeABC", &pcprocs->Decode.ABC)) < 0 ||
  172.          (code = dict_matrix3_param(pdref, "MatrixABC", &pcie->MatrixABC)) != matrix3_ok ||
  173.          (code = cie_lmnp_param(pdref, &pcie->common, pcprocs)) < 0
  174.        )
  175.     {    rc_free_struct(pcie, mem, "setcolorspace(CIEBasedABC)");
  176.         return (code < 0 ? code : e_rangecheck);
  177.     }
  178.     pcie->DecodeABC = DecodeABC_default;
  179.     pcs->params.abc = pcie;
  180.     pcs->type = &cs_type_zCIEBasedABC;
  181.     return 0;    /* installation will load the caches */
  182. }
  183.  
  184. /* Get the parameters of a CIEBasedA color space. */
  185. /* See above. */
  186. int
  187. zcolorspace_CIEBasedA(const ref *pdref, gs_color_space *pcs,
  188.   ref_cie_procs *pcprocs)
  189. {    gs_memory_t *mem = gs_state_memory(igs);
  190.     gs_cie_a *pcie;
  191.     int code;
  192.     check_type(*pdref, t_dictionary);
  193.     check_dict_read(*pdref);
  194.     if ( (code = dict_proc_param(pdref, "DecodeA", &pcprocs->Decode.A)) < 0 )
  195.         return code;
  196.     rc_alloc_struct_0(pcie, gs_cie_a, &st_cie_a, mem,
  197.               return_error(e_VMerror),
  198.               "setcolorspace(CIEBasedA)");
  199.     if ( (code = dict_float_array_param(pdref, "RangeA", 2, (float *)&pcie->RangeA, (float *)&RangeA_default)) != 2 ||
  200.          (code = dict_float_array_param(pdref, "MatrixA", 3, (float *)&pcie->MatrixA, (float *)&MatrixA_default)) != 3 ||
  201.          (code = cie_lmnp_param(pdref, &pcie->common, pcprocs)) < 0
  202.        )
  203.     {    rc_free_struct(pcie, mem, "setcolorspace(CIEBasedA)");
  204.         return (code < 0 ? code : e_rangecheck);
  205.     }
  206.     pcie->DecodeA = DecodeA_default;
  207.     pcs->params.a = pcie;
  208.     pcs->type = &cs_type_zCIEBasedA;
  209.     return 0;    /* installation will load the caches */
  210. }
  211.  
  212. /* ------ CIE rendering dictionary ------ */
  213.  
  214. /* - currentcolorrendering <dict> */
  215. private int
  216. zcurrentcolorrendering(register os_ptr op)
  217. {    push(1);
  218.     *op = istate->colorrendering.dict;
  219.     return 0;
  220. }
  221.  
  222. /* <dict> setcolorrendering - */
  223. private int zsetcolorrendering_internal(P4(os_ptr, gs_cie_render *, ref_cie_render_procs *, gs_memory_t *));
  224. private int
  225. zsetcolorrendering(register os_ptr op)
  226. {    gs_memory_t *mem = gs_state_memory(igs);
  227.     int code;
  228.     es_ptr ep = esp;
  229.     gs_cie_render *pcie;
  230.     ref_cie_render_procs procs_old;
  231.     check_read_type(*op, t_dictionary);
  232.     check_dict_read(*op);
  233.     rc_alloc_struct_0(pcie, gs_cie_render, &st_cie_render, mem,
  234.               return_error(e_VMerror),
  235.               "setcolorrendering");
  236.     /* gs_setcolorrendering may refer to istate->colorrendering.procs. */
  237.     procs_old = istate->colorrendering.procs;
  238.     code = zsetcolorrendering_internal(op, pcie, &istate->colorrendering.procs, mem);
  239.     if ( code < 0 )
  240.     {    rc_free_struct(pcie, mem, "setcolorrendering");
  241.         istate->colorrendering.procs = procs_old;
  242.         esp = ep;
  243.         return code;
  244.     }
  245.     istate->colorrendering.dict = *op;
  246.     pop(1);
  247.     return (esp == ep ? 0 : o_push_estack);
  248. }
  249. private int
  250. zsetcolorrendering_internal(os_ptr op, gs_cie_render *pcie,
  251.   ref_cie_render_procs *pcprocs, gs_memory_t *mem)
  252. {    int code;
  253.     int discard;
  254.     ref *pRT;
  255.     if ( (code = dict_int_param(op, "ColorRenderingType", 1, 1, 0, &discard)) < 0 ||
  256.          (code = dict_matrix3_param(op, "MatrixLMN", &pcie->MatrixLMN)) != matrix3_ok ||
  257.          (code = dict_proc3_param(op, "EncodeLMN", &pcprocs->EncodeLMN)) < 0 ||
  258.          (code = dict_range3_param(op, "RangeLMN", &pcie->RangeLMN)) != range3_ok ||
  259.          (code = dict_matrix3_param(op, "MatrixABC", &pcie->MatrixABC)) != matrix3_ok ||
  260.          (code = dict_proc3_param(op, "EncodeABC", &pcprocs->EncodeABC)) < 0 ||
  261.          (code = dict_range3_param(op, "RangeABC", &pcie->RangeABC)) != range3_ok ||
  262.          (code = cie_points_param(op, &pcie->points)) < 0 ||
  263.          (code = dict_matrix3_param(op, "MatrixPQR", &pcie->MatrixPQR)) != matrix3_ok ||
  264.          (code = dict_range3_param(op, "RangePQR", &pcie->RangePQR)) != range3_ok ||
  265.          (code = dict_proc3_param(op, "TransformPQR", &pcprocs->TransformPQR)) < 0
  266.        )
  267.         return (code < 0 ? code : e_rangecheck);
  268. #define rRT pcie->RenderTable
  269.     if ( dict_find_string(op, "RenderTable", &pRT) > 0 )
  270.     {    const ref *prte;
  271.         int i;
  272.         uint n2;
  273.         const ref *strings;
  274.         check_read_type(*pRT, t_array);
  275.         prte = pRT->value.const_refs;
  276.         check_type(prte[0], t_integer);
  277.         check_type(prte[1], t_integer);
  278.         check_type(prte[2], t_integer);
  279.         check_read_type(prte[3], t_array);
  280.         check_type(prte[4], t_integer);
  281.         if ( prte[0].value.intval <= 1 ||
  282.              prte[1].value.intval <= 1 ||
  283.              prte[2].value.intval <= 1 ||
  284.              !(prte[4].value.intval == 3 || prte[4].value.intval == 4)
  285.            )
  286.             return_error(e_rangecheck);
  287.         rRT.NA = prte[0].value.intval;
  288.         rRT.NB = prte[1].value.intval;
  289.         rRT.NC = prte[2].value.intval;
  290.         rRT.m = prte[4].value.intval;
  291.         n2 = rRT.m * rRT.NB * rRT.NC;
  292.         if ( r_size(pRT) != rRT.m + 5 || r_size(&prte[3]) != rRT.NA )
  293.             return_error(e_rangecheck);
  294.         strings = prte[3].value.const_refs;
  295.         for ( i = 0; i < rRT.NA; i++ )
  296.         {    const ref *prt2 = strings + i;
  297.             check_read_type(*prt2, t_string);
  298.             if ( r_size(prt2) != n2 )
  299.                 return_error(e_rangecheck);
  300.         }
  301.         prte += 5;
  302.         for ( i = 0; i < rRT.m; i++ )
  303.         {    const ref *prt2 = prte + i;
  304.             check_proc(*prt2);
  305.         }
  306.         /* We allocate the table as a byte array, rather than as */
  307.         /* a byte * array, because it will require special */
  308.         /* relocation after a GC. */
  309.         rRT.table = (const byte **)gs_alloc_byte_array(mem, rRT.NA,
  310.                         sizeof(byte *),
  311.                         "setcolorrendering(table)");
  312.         if ( rRT.table == 0 )
  313.             return_error(e_VMerror);
  314.         for ( i = 0; i < rRT.NA; i++ )
  315.             rRT.table[i] = strings[i].value.bytes;
  316.         make_const_array(&pcprocs->RenderTableT, a_readonly, rRT.m, prte);
  317.     }
  318.     else
  319.     {    rRT.table = 0;
  320.         make_null(&pcprocs->RenderTableT);
  321.     }
  322. #undef rRT
  323.     pcie->EncodeLMN = Encode_default;
  324.     pcie->EncodeABC = Encode_default;
  325.     pcie->TransformPQR = TransformPQR_default;
  326.     pcie->RenderTable.T = RenderTableT_default;
  327.     code = cache_colorrendering(pcie, pcprocs, igs);
  328.     if ( code < 0 )
  329.         return code;
  330.     return gs_setcolorrendering(igs, pcie);
  331. }
  332.  
  333. /* ------ Internal routines ------ */
  334.  
  335. /* Import operators. */
  336. extern int
  337.   zfor(P1(os_ptr)),
  338.   zcvx(P1(os_ptr));
  339. extern int zexec(P1(os_ptr));
  340. /* Import accessors. */
  341. extern gx_cie_joint_caches *gx_currentciecaches(P1(gs_state *));
  342.  
  343. /* Forward declarations */
  344. private int
  345.   cie_cache_finish(P1(os_ptr)),
  346.   cie_exec_tpqr(P1(os_ptr)),
  347.   cie_tpqr_finish(P1(os_ptr));
  348.  
  349. /* Transform a set of ranges. */
  350. private void
  351. cie_transform_range(const gs_range3 *in, const gs_vector3 *col,
  352.   gs_range *out)
  353. {    float umin = col->u * in->u.rmin, umax = col->u * in->u.rmax;
  354.     float vmin = col->v * in->v.rmin, vmax = col->v * in->v.rmax;
  355.     float wmin = col->w * in->w.rmin, wmax = col->w * in->w.rmax;
  356.     float temp;
  357. #define swap(x, y) temp = x, x = y, y = temp
  358.     if ( umin > umax ) swap(umin, umax);
  359.     if ( vmin > vmax ) swap(vmin, vmax);
  360.     if ( wmin > wmax ) swap(wmin, wmax);
  361.     out->rmin = umin + vmin + wmin;
  362.     out->rmax = umax + vmax + wmax;
  363. #undef swap
  364. }
  365. private void
  366. cie_transform_range3(const gs_range3 *in, const gs_matrix3 *mat,
  367.   gs_range3 *out)
  368. {    cie_transform_range(in, &mat->cu, &out->u);
  369.     cie_transform_range(in, &mat->cv, &out->v);
  370.     cie_transform_range(in, &mat->cw, &out->w);
  371. }
  372.  
  373. /* Prepare to cache the values for one or more procedures. */
  374. private int
  375. cie_prepare_caches(const gs_range *domain, const ref *proc,
  376.   gx_cie_cache *pcache, int n)
  377. {    check_estack(n * 8);
  378.     for ( ; --n >= 0; domain++, proc++, pcache++, esp += 8 )
  379.     {    float diff = domain->rmax - domain->rmin;
  380.         float delta = diff / (gx_cie_cache_size - 1);
  381.         register es_ptr ep = esp;
  382.         /* Zero out the cache, since the gs level will try to */
  383.         /* access it before it has been filled. */
  384.         {    register float *pcv = &pcache->values[0];
  385.             register int i;
  386.             for ( i = 0; i < gx_cie_cache_size; i++, pcv++ )
  387.                 *pcv = 0.0;
  388.         }
  389.         pcache->base = domain->rmin - delta / 2;    /* so lookup will round */
  390.         pcache->factor = (delta == 0 ? 0 : 1 / delta);
  391.         pcache->is_identity = r_size(proc) == 0;
  392.         if_debug3('c', "[c]cache 0x%lx base=%g, factor=%g\n",
  393.               (ulong)pcache, pcache->base, pcache->factor);
  394.         make_real(ep + 8, domain->rmin);
  395.         make_real(ep + 7, delta);
  396.         make_real(ep + 6, domain->rmax + delta / 2);
  397.         ep[5] = *proc;
  398.         r_clear_attrs(ep + 5, a_executable);
  399.         make_op_estack(ep + 4, zcvx);
  400.         make_op_estack(ep + 3, zfor);
  401.         make_op_estack(ep + 2, cie_cache_finish);
  402.         make_string(ep + 1, 0, sizeof(*pcache), (byte *)pcache);
  403.     }
  404.     return o_push_estack;
  405. }
  406. /* Prepare to cache the values for 3 procedures. */
  407. #define cie_prepare_cache3(d3,p3,c3)\
  408.   cie_prepare_caches((const gs_range *)(d3), p3, c3, 3)
  409.  
  410. /* Store the result of caching one procedure. */
  411. private int
  412. cie_cache_finish(os_ptr op)
  413. {    gx_cie_cache *pcache;
  414.     int code;
  415.     check_esp(1);
  416.     /* The following should be
  417.         pcache = r_ptr(esp, gx_cie_cache);
  418.      * but we can't do this right now, because the caches are
  419.      * embedded in the middle of another structure.
  420.      */
  421.     pcache = (gx_cie_cache *)esp->value.bytes;
  422.     code = num_params(op, gx_cie_cache_size, &pcache->values[0]);
  423.     if_debug3('c', "[c]cache 0x%lx base=%g, factor=%g:\n",
  424.           (ulong)pcache, pcache->base, pcache->factor);
  425.     if ( code < 0 )
  426.       {    /* We might have underflowed the current stack block. */
  427.         /* Handle the parameters one-by-one. */
  428.         uint i;
  429.         for ( i = 0; i < gx_cie_cache_size; i++ )
  430.           {    code = real_param(ref_stack_index(&o_stack,
  431.                         gx_cie_cache_size - 1 - i),
  432.                       &pcache->values[i]);
  433.             if ( code < 0 )
  434.               return code;
  435.           }
  436.       }
  437. #ifdef DEBUG
  438.     if ( gs_debug_c('c') )
  439.     {    int i;
  440.         for ( i = 0; i < gx_cie_cache_size; i++ )
  441.           dprintf2("[c]cache[%3d]=%g\n", i, pcache->values[i]);
  442.     }
  443. #endif
  444.     ref_stack_pop(&o_stack, gx_cie_cache_size);
  445.     esp--;                /* pop pointer to cache */
  446.     return o_pop_estack;
  447. }
  448.  
  449. /* Install a CIE-based color space. */
  450.  
  451. private int
  452. cs_install_zCIEBasedABC(gs_color_space *pcs, gs_state *pgs)
  453. {    es_ptr ep = esp;
  454.     gs_cie_abc *pcie = pcs->params.abc;
  455.     const int_gstate *pigs = gs_int_gstate(pgs);
  456.     const ref_cie_procs *pcprocs = &pigs->colorspace.procs.cie;
  457.     int code = gx_install_CIEBasedABC(pcs, pgs);    /* former routine */
  458.     if ( code < 0 ) return code;
  459.     code = cie_prepare_cache3(&pcie->RangeABC, pcprocs->Decode.ABC.value.const_refs, &pcie->caches.DecodeABC[0]);
  460.     if ( code < 0 ||
  461.          (code = cache_common(&pcie->common, pcprocs, &pigs->colorrendering.procs, pgs)) < 0
  462.        )
  463.     {    esp = ep;
  464.         return code;
  465.     }
  466.     return o_push_estack;
  467. }
  468.  
  469. private int
  470. cs_install_zCIEBasedA(gs_color_space *pcs, gs_state *pgs)
  471. {    es_ptr ep = esp;
  472.     gs_cie_a *pcie = pcs->params.a;
  473.     const int_gstate *pigs = gs_int_gstate(pgs);
  474.     const ref_cie_procs *pcprocs = &pigs->colorspace.procs.cie;
  475.     int code = gx_install_CIEBasedA(pcs, pgs);    /* former routine */
  476.     if ( code < 0 ) return code;
  477.     code = cie_prepare_caches(&pcie->RangeA, &pcprocs->Decode.A, &pcie->caches.DecodeA, 1);
  478.     if ( code < 0 ||
  479.          (code = cache_common(&pcie->common, pcprocs, &pigs->colorrendering.procs, pgs)) < 0
  480.        )
  481.     {    esp = ep;
  482.         return code;
  483.     }
  484.     return o_push_estack;
  485. }
  486.  
  487. /* Cache the results of the color rendering procedures. */
  488. private int
  489. cache_colorrendering(gs_cie_render *pcie,
  490.   const ref_cie_render_procs *pcrprocs, gs_state *pgs)
  491. {    gs_range3 DomainLMN;
  492.     gs_range3 DomainABC;
  493.     es_ptr ep = esp;
  494.     static const gs_range ranges_01[4] =
  495.         { {0,1}, {0,1}, {0,1}, {0,1} };
  496.     int code = gs_cie_render_init(pcie);    /* sets PQR'*LMN */
  497.     if ( code < 0 ) return code;
  498.     cie_transform_range3(&pcie->RangePQR, &pcie->MatrixPQR_inverse_LMN, &DomainLMN);
  499.     cie_transform_range3(&pcie->RangeLMN, &pcie->MatrixABC, &DomainABC);
  500.     if ( (code = cie_prepare_cache3(&DomainLMN, pcrprocs->EncodeLMN.value.const_refs, &pcie->caches.EncodeLMN[0])) < 0 ||
  501.          (code = cie_prepare_cache3(&DomainABC, pcrprocs->EncodeABC.value.const_refs, &pcie->caches.EncodeABC[0])) < 0 ||
  502.          (pcie->RenderTable.table != 0 &&
  503.           (code = cie_prepare_caches(ranges_01, pcrprocs->RenderTableT.value.const_refs, &pcie->caches.RenderTableT[0], pcie->RenderTable.m)) < 0)
  504.        )
  505.     {    esp = ep;
  506.         return code;
  507.     }
  508.     /* gs_setcolorrendering reinstalls the color space, */
  509.     /* which reloads the joint caches if needed. */
  510.     return o_push_estack;
  511. }
  512.  
  513. /* Common cache code */
  514. private int
  515. cache_common(gs_cie_common *pcie, const ref_cie_procs *pcprocs,
  516.   const ref_cie_render_procs *pcrprocs, gs_state *pgs)
  517. {    int code = cie_prepare_cache3(&pcie->RangeLMN,
  518.                       pcprocs->DecodeLMN.value.const_refs,
  519.                       &pcie->caches.DecodeLMN[0]);
  520.     const gs_cie_render *pcier = gs_currentcolorrendering(pgs);
  521.     /* The former installation procedures have allocated */
  522.     /* the joint caches and filled in points_sd. */
  523.     gx_cie_joint_caches *pjc = gx_currentciecaches(pgs);
  524.     ref pqr_procs;
  525. #define pqr_refs pqr_procs.value.refs
  526.     int i;
  527.     if ( code < 0 ) return code;
  528.     if ( pcier == 0 ) return 0;    /* cache is not used */
  529.     check_estack(2);
  530.     code = ialloc_ref_array(&pqr_procs, a_readonly, 3*(1+3+4*6), "cie_cache_common");
  531.     if ( code < 0 ) return code;
  532.     /* Make sure we deallocate the procs when we're done. */
  533.     push_op_estack(cie_tpqr_finish);
  534.     *++esp = pqr_procs;
  535.     for ( i = 0; i < 3; i++ )
  536.     {    ref *p = pqr_refs + 3 + (3+4*6) * i;
  537.         const float *ppt = (float *)&pjc->points_sd;
  538.         int j;
  539.         make_array(pqr_refs + i, a_readonly + a_executable, 3, p);
  540.         make_array(p, a_readonly, 4*6, p + 3);
  541.         p[1] = pcrprocs->TransformPQR.value.refs[i];
  542.         make_oper(p + 2, 0, cie_exec_tpqr);
  543.         for ( j = 0, p += 3; j < 4*6; j++, p++, ppt++ )
  544.             make_real(p, *ppt);
  545.     }
  546.     return cie_prepare_cache3(&pcier->RangePQR,
  547.                   pqr_procs.value.const_refs,
  548.                   &pjc->TransformPQR[0]);
  549. }
  550.  
  551. /* Private operator to shuffle arguments for the TransformPQR procedure: */
  552. /* v [ws wd bs bd] proc -> ws wd bs bd v proc + exec */
  553. private int
  554. cie_exec_tpqr(register os_ptr op)
  555. {    const ref *ppt = op[-1].value.const_refs;
  556.     int i;
  557.     check_op(3);
  558.     push(3);
  559.     *op = op[-3];        /* proc */
  560.     op[-1] = op[-5];    /* v */
  561.     for ( i = 0; i < 4; i++ )
  562.         make_const_array(op - 5 + i, a_readonly, 6, ppt + i * 6);
  563.     return zexec(op);
  564. }
  565.  
  566. /* Private operator to free procs array. */
  567. private int
  568. cie_tpqr_finish(register os_ptr op)
  569. {    ifree_ref_array(op, "cie_tpqr_finish");
  570.     return 0;
  571. }
  572.  
  573. /* ------ Initialization procedure ------ */
  574.  
  575. op_def zcie_l2_op_defs[] = {
  576.         op_def_begin_level2(),
  577.     {"0currentcolorrendering", zcurrentcolorrendering},
  578.     {"1setcolorrendering", zsetcolorrendering},
  579.         /* Internal operators */
  580.     {"0%cie_cache_finish", cie_cache_finish},
  581.     {"3%cie_exec_tpqr", cie_exec_tpqr},
  582.     {"1%cie_tpqr_finish", cie_tpqr_finish},
  583.     op_def_end(zcie_init)
  584. };
  585.