1 % Copyright (C) 1995 Aladdin Enterprises. All rights reserved.
3 % This software is provided AS-IS with no warranty, either express or
6 % This software is distributed under license and may not be copied,
7 % modified or distributed except as expressly authorized under the terms
8 % of the license contained in the file LICENSE in this distribution.
10 % For more information about licensing, please refer to
11 % http://www.ghostscript.com/licensing/. For information on
12 % commercial licensing, go to http://www.artifex.com/licensing/ or
13 % contact Artifex Software, Inc., 101 Lucas Valley Road #110,
14 % San Rafael, CA 94903, U.S.A., +1(415)492-9861.
16 % $Id: docie.ps,v 1.5 2003/04/12 18:08:18 ray Exp $
18 % Emulate CIE algorithms in PostScript.
20 % ---------------- Auxiliary procedures ---------------- %
23 /r3default [0 1 0 1 0 1] def
25 /apply3 % <u> <v> <w> [<pu> <pv> <pw>] apply3 <u'> <v'> <w'>
26 { { 4 -1 roll exch exec } forall
29 /restrict % <u> <min> <max> restrict <u'>
33 /restrict3 % <u> <v> <w> [<minu> ... <maxw>] restrict3 <u'> <v'> <w'>
35 7 -1 roll 3 1 roll restrict 7 1 roll
36 5 -1 roll 3 1 roll restrict 5 1 roll
40 /rescale % <u> <min> <max> rescale <u'>
41 { 1 index sub 3 1 roll sub exch div 0 .max 1 .min
44 /rescale3 % <u> <v> <w> [<minu> ... <maxw>] rescale3 <u'> <v'> <w'>
46 7 -1 roll 3 1 roll rescale 7 1 roll
47 5 -1 roll 3 1 roll rescale 5 1 roll
51 /mmult3 % <u> <v> <w> [<uu> <uv> ... <wv> <ww>] mmult3
53 { 4 -1 roll dup dup 6 -1 roll dup dup 8 -1 roll dup dup
54 10 -1 roll { 10 -1 roll mul } forall
55 % Stack: u1 v1 w1 u2 v2 w2 u3 v3 w3
56 4 -1 roll add 6 -1 roll add
57 % Stack: u1 v1 u2 v2 u3 v3 w'
58 7 1 roll 3 -1 roll add 4 -1 roll add
59 % Stack: w' u1 u2 u3 v'
60 5 1 roll add add 3 1 roll
63 /minvert3 % [<uu> <uv> ... <wv> <ww>] minvert3
64 % [<uu'> <uv'> ... <wv'> <ww'>]
66 aload pop { I H G F E D C B A } { exch def } forall
67 /coa E I mul F H mul sub def
68 /cob F G mul D I mul sub def
69 /coc D H mul E G mul sub def
70 /det A coa mul B cob mul add C coc mul add def
72 C H mul B I mul sub det div
73 B F mul C E mul sub det div
75 A I mul C G mul sub det div
76 C D mul A F mul sub det div
78 B G mul A H mul sub det div
79 A E mul B D mul sub det div
89 { print 3 array astore dup == aload pop
92 % ---------------- Mapping to XYZ ---------------- %
94 /csmap % <csdict> <l> <m> <n> csmap <csdict> <x> <y> <z>
95 { 3 index /RangeLMN .knownget not { r3default } if restrict3
96 DOCIEDEBUG { (After RangeLMN Decode: ) print3 } if
97 3 index /DecodeLMN .knownget { apply3 } if
98 DOCIEDEBUG { (After DecodeLMN Decode: ) print3 } if
99 3 index /MatrixLMN .knownget { mmult3 } if
100 DOCIEDEBUG { (After MatrixLMN Decode: ) print3 } if
103 /csciea % <csdict> <a> csciea <csdict> <x> <y> <z>
104 { 1 index /RangeA .knownget not { r1default } if restrict
105 DOCIEDEBUG { (After RangeA Decode: ) print1 } if
106 1 index /DecodeA .knownget { exec } if
107 DOCIEDEBUG { (After DecodeA Decode: ) print1 } if
108 1 index /MatrixA .knownget
109 { { 1 index mul exch } forall pop }
112 DOCIEDEBUG { (After MatrixA Decode: ) print3 } if
116 /cscieabc % <csdict> <a> <b> <c> cscieabc <csdict> <x> <y> <z>
117 { 3 index /RangeABC .knownget not { r3default } if restrict3
118 DOCIEDEBUG { (After RangeABC Decode: ) print3 } if
119 3 index /DecodeABC .knownget { apply3 } if
120 DOCIEDEBUG { (After DecodeABC Decode: ) print3 } if
121 3 index /MatrixABC .knownget { mmult3 } if
122 DOCIEDEBUG { (After MatrixABC Decode: ) print3 } if
126 % ---------------- Rendering from XYZ ---------------- %
128 /lookup3 % <rtable> <a[0..1]> <b[0..1]> <c[0..1]> lookup3
130 { 3 -1 roll 3 index 0 get 1 sub mul
131 3 -1 roll 3 index 1 get 1 sub mul
132 3 -1 roll 3 index 2 get 1 sub mul
133 % Stack: rtable ia ib ic
134 DOCIEDEBUG { (RenderTable indices: ) print3 mark 5 1 roll } if
135 3 -1 roll round cvi 3 index 3 get exch get
136 % Stack: rtable ib ic string
137 3 -1 roll round cvi 3 index 2 get mul
138 % Stack: rtable ic string ib*nc
139 3 -1 roll round cvi add 2 index 4 get mul
140 % Stack: rtable string index
141 2 index 4 get getinterval
142 % Stack: rtable bytes
143 DOCIEDEBUG { (RenderTable values: ) print (<) print (%stdout) (w) file 1 index writehexstring (>) = } if
146 /bpdefault [0 0 0] def
148 /crmap % <csdict> <crdict> <x> <y> <z> crmap <v1> ...
150 DOCIEDEBUG { (CIE XYZ = ) print3 } if
151 3 index /MatrixPQR .knownget { mmult3 } if
152 DOCIEDEBUG { (After MatrixPQR: ) print3 } if
153 4 index /WhitePoint get
154 5 index /BlackPoint .knownget not { bpdefault } if
155 5 index /WhitePoint get
156 6 index /BlackPoint .knownget not { bpdefault } if
158 { 4 -1 roll aload pop
159 % Stack: csdict crdict x y z pt pt pt px py pz
160 3 copy 12 index /MatrixPQR .knownget { mmult3 } if 6 array astore
163 % Stack: csdict crdict x y z wps+ bps+ wpd+ bpd+
164 9 -1 roll pop % get rid of csdict
166 7 index /TransformPQR get
167 { % Stack: crdict wps+ bps+ wpd+ bpd+ u v w proc
168 8 copy exch pop exch pop
169 exec exch pop 4 -1 roll pop
172 7 3 roll pop pop pop pop % get rid of White/BlackPoints
173 DOCIEDEBUG { (After TransformPQR: ) print3 } if
174 3 index /MatrixPQR .knownget { minvert3 mmult3 } if
175 DOCIEDEBUG { (After MatrixPQR': ) print3 } if
176 3 index /MatrixLMN .knownget { mmult3 } if
177 DOCIEDEBUG { (After MatrixLMN Encode: ) print3 } if
178 3 index /EncodeLMN .knownget { apply3 } if
179 DOCIEDEBUG { (After EncodeLMN Encode: ) print3 } if
180 3 index /RangeLMN .knownget not { r3default } if restrict3
181 DOCIEDEBUG { (After RangeLMN Encode: ) print3 } if
182 3 index /MatrixABC .knownget { mmult3 } if
183 DOCIEDEBUG { (After MatrixABC Encode: ) print3 } if
184 3 index /EncodeABC .knownget { apply3 } if
185 DOCIEDEBUG { (After EncodeABC Encode: ) print3 } if
186 3 index /RangeABC .knownget not { r3default } if
187 5 -1 roll /RenderTable .knownget
188 { % Stack: u v w ranges rtable
190 DOCIEDEBUG { (Rescaled ABC: ) print3 } if
191 % Stack: rtable a b c
193 % Stack: rtable bytes
194 0 1 3 index 4 get 1 sub
195 { % Stack: values rtable bytes c
197 % Stack: values rtable bytes c v
198 3 index 3 -1 roll 5 add get exec 3 1 roll
201 DOCIEDEBUG { (After RenderTableT: ) print ] dup == aload pop } if
204 DOCIEDEBUG { (After RangeABC Encode: ) print3 } if
209 % ---------------- Top level control ---------------- %
212 /CIEBasedA { 1 get exch csciea currentcolorrendering 4 1 roll crmap } bind
213 /CIEBasedABC { 1 get 4 1 roll cscieabc currentcolorrendering 4 1 roll crmap } bind
215 /mapcie % <a> mapcie <v1> ...
216 % <a> <b> <c> mapcie <v1> ...
217 { currentcolorspace dup 0 get //mapdict exch get exec