home *** CD-ROM | disk | FTP | other *** search
/ Acorn User 10 / AU_CD10.iso / Updates / GhostScript / !GhostScr / 6_01 / lib / docie.ps < prev    next >
Text File  |  2000-03-09  |  7KB  |  221 lines
  1. %    Copyright (C) 1995 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. % $Id: docie.ps,v 1.1 2000/03/09 08:40:39 lpd Exp $
  19. % docie.ps
  20. % Emulate CIE algorithms in PostScript.
  21.  
  22. % ---------------- Auxiliary procedures ---------------- %
  23.  
  24. /r1default [0 1] def
  25. /r3default [0 1 0 1 0 1] def
  26.  
  27. /apply3        % <u> <v> <w> [<pu> <pv> <pw>] apply3 <u'> <v'> <w'>
  28.  {  { 4 -1 roll exch exec } forall
  29.  } bind def
  30.  
  31. /restrict    % <u> <min> <max> restrict <u'>
  32.  { 3 1 roll .max .min
  33.  } bind def
  34.  
  35. /restrict3    % <u> <v> <w> [<minu> ... <maxw>] restrict3 <u'> <v'> <w'>
  36.  { aload pop
  37.    7 -1 roll 3 1 roll restrict 7 1 roll
  38.    5 -1 roll 3 1 roll restrict 5 1 roll
  39.    restrict 3 1 roll
  40.  } bind def
  41.  
  42. /rescale    % <u> <min> <max> rescale <u'>
  43.  { 1 index sub 3 1 roll sub exch div 0 .max 1 .min
  44.  } bind def
  45.  
  46. /rescale3    % <u> <v> <w> [<minu> ... <maxw>] rescale3 <u'> <v'> <w'>
  47.  { aload pop
  48.    7 -1 roll 3 1 roll rescale 7 1 roll
  49.    5 -1 roll 3 1 roll rescale 5 1 roll
  50.    rescale 3 1 roll
  51.  } bind def
  52.  
  53. /mmult3        % <u> <v> <w> [<uu> <uv> ... <wv> <ww>] mmult3
  54.         %   <u'> <v'> <w'>
  55.  { 4 -1 roll dup dup 6 -1 roll dup dup 8 -1 roll dup dup
  56.    10 -1 roll { 10 -1 roll mul } forall
  57.         % Stack: u1 v1 w1 u2 v2 w2 u3 v3 w3
  58.    4 -1 roll add 6 -1 roll add
  59.         % Stack: u1 v1 u2 v2 u3 v3 w'
  60.    7 1 roll 3 -1 roll add 4 -1 roll add
  61.         % Stack: w' u1 u2 u3 v'
  62.    5 1 roll add add 3 1 roll
  63.  } bind def
  64.  
  65. /minvert3    % [<uu> <uv> ... <wv> <ww>] minvert3
  66.         %   [<uu'> <uv'> ... <wv'> <ww'>]
  67.  { 16 dict begin
  68.    aload pop { I H G F E D C B A } { exch def } forall
  69.    /coa E I mul F H mul sub def
  70.    /cob F G mul D I mul sub def
  71.    /coc D H mul E G mul sub def
  72.    /det A coa mul B cob mul add C coc mul add def
  73.    [ coa det div
  74.      C H mul B I mul sub det div
  75.      B F mul C E mul sub det div
  76.      cob det div
  77.      A I mul C G mul sub det div
  78.      C D mul A F mul sub det div
  79.      coc det div
  80.      B G mul A H mul sub det div
  81.      A E mul B D mul sub det div
  82.    ]
  83.    end
  84.  } bind def
  85.  
  86. /print1
  87.  { print dup ==
  88.  } bind def
  89.  
  90. /print3
  91.  { print 3 array astore dup == aload pop
  92.  } bind def
  93.  
  94. % ---------------- Mapping to XYZ ---------------- %
  95.  
  96. /csmap        % <csdict> <l> <m> <n> csmap <csdict> <x> <y> <z>
  97.  { 3 index /RangeLMN .knownget not { r3default } if restrict3
  98. DEBUG { (After RangeLMN Decode: ) print3 } if
  99.    3 index /DecodeLMN .knownget { apply3 } if
  100. DEBUG { (After DecodeLMN Decode: ) print3 } if
  101.    3 index /MatrixLMN .knownget { mmult3 } if
  102. DEBUG { (After MatrixLMN Decode: ) print3 } if
  103.  } bind def
  104.  
  105. /csciea        % <csdict> <a> csciea <csdict> <x> <y> <z>
  106.  { 1 index /RangeA .knownget not { r1default } if restrict
  107. DEBUG { (After RangeA Decode: ) print1 } if
  108.    1 index /DecodeA .knownget { exec } if
  109. DEBUG { (After DecodeA Decode: ) print1 } if
  110.    1 index /MatrixA .knownget
  111.     { { 1 index mul exch } forall pop }
  112.     { dup dup }
  113.    ifelse
  114. DEBUG { (After MatrixA Decode: ) print3 } if
  115.    csmap
  116.  } bind def
  117.  
  118. /cscieabc    % <csdict> <a> <b> <c> cscieabc <csdict> <x> <y> <z>
  119.  { 3 index /RangeABC .knownget not { r3default } if restrict3
  120. DEBUG { (After RangeABC Decode: ) print3 } if
  121.    3 index /DecodeABC .knownget { apply3 } if
  122. DEBUG { (After DecodeABC Decode: ) print3 } if
  123.    3 index /MatrixABC .knownget { mmult3 } if
  124. DEBUG { (After MatrixABC Decode: ) print3 } if
  125.    csmap
  126.  } bind def
  127.  
  128. % ---------------- Rendering from XYZ ---------------- %
  129.  
  130. /lookup3    % <rtable> <a[0..1]> <b[0..1]> <c[0..1]> lookup3
  131.         %   <rtable> <bytes>
  132.  { 3 -1 roll 3 index 0 get 1 sub mul
  133.    3 -1 roll 3 index 1 get 1 sub mul
  134.    3 -1 roll 3 index 2 get 1 sub mul
  135.         % Stack: rtable ia ib ic
  136. DEBUG { (RenderTable indices: ) print3 mark 5 1 roll } if
  137.    3 -1 roll round cvi 3 index 3 get exch get
  138.         % Stack: rtable ib ic string
  139.    3 -1 roll round cvi 3 index 2 get mul
  140.         % Stack: rtable ic string ib*nc
  141.    3 -1 roll round cvi add 2 index 4 get mul
  142.         % Stack: rtable string index
  143.    2 index 4 get getinterval
  144.         % Stack: rtable bytes
  145. DEBUG { (RenderTable values: ) print (<) print (%stdout) (w) file 1 index writehexstring (>) = } if
  146.  } bind def
  147.  
  148. /bpdefault [0 0 0] def
  149.  
  150. /crmap        % <csdict> <crdict> <x> <y> <z> crmap <v1> ...
  151.  {
  152. DEBUG { (CIE XYZ = ) print3 } if
  153.    3 index /MatrixPQR .knownget { mmult3 } if
  154. DEBUG { (After MatrixPQR: ) print3 } if
  155.    4 index /WhitePoint get
  156.    5 index /BlackPoint .knownget not { bpdefault } if
  157.    5 index /WhitePoint get
  158.    6 index /BlackPoint .knownget not { bpdefault } if
  159.    4
  160.      { 4 -1 roll aload pop
  161.             % Stack: csdict crdict x y z pt pt pt px py pz
  162.        3 copy 12 index /MatrixPQR .knownget { mmult3 } if 6 array astore
  163.      }
  164.    repeat
  165.             % Stack: csdict crdict x y z wps+ bps+ wpd+ bpd+
  166.    9 -1 roll pop        % get rid of csdict
  167.    7 4 roll
  168.    7 index /TransformPQR get
  169.     {            % Stack: crdict wps+ bps+ wpd+ bpd+ u v w proc
  170.       8 copy exch pop exch pop
  171.       exec exch pop 4 -1 roll pop
  172.     }
  173.    forall
  174.    7 3 roll pop pop pop pop    % get rid of White/BlackPoints
  175. DEBUG { (After TransformPQR: ) print3 } if
  176.    3 index /MatrixPQR .knownget { minvert3 mmult3 } if
  177. DEBUG { (After MatrixPQR': ) print3 } if
  178.    3 index /MatrixLMN .knownget { mmult3 } if
  179. DEBUG { (After MatrixLMN Encode: ) print3 } if
  180.    3 index /EncodeLMN .knownget { apply3 } if
  181. DEBUG { (After EncodeLMN Encode: ) print3 } if
  182.    3 index /RangeLMN .knownget not { r3default } if restrict3
  183. DEBUG { (After RangeLMN Encode: ) print3 } if
  184.    3 index /MatrixABC .knownget { mmult3 } if
  185. DEBUG { (After MatrixABC Encode: ) print3 } if
  186.    3 index /EncodeABC .knownget { apply3 } if
  187. DEBUG { (After EncodeABC Encode: ) print3 } if
  188.    3 index /RangeABC .knownget not { r3default } if
  189.    5 -1 roll /RenderTable .knownget
  190.     {        % Stack: u v w ranges rtable
  191.       5 1 roll rescale3
  192. DEBUG { (Rescaled ABC: ) print3 } if
  193.         % Stack: rtable a b c
  194.       lookup3
  195.         % Stack: rtable bytes
  196.       0 1 3 index 4 get 1 sub
  197.        {    % Stack: values rtable bytes c
  198.      2 copy get 255 div
  199.         % Stack: values rtable bytes c v
  200.      3 index 3 -1 roll 5 add get exec 3 1 roll
  201.        }
  202.       for pop pop
  203. DEBUG { (After RenderTableT: ) print ] dup == aload pop } if
  204.     }
  205.     { restrict3
  206. DEBUG { (After RangeABC Encode: ) print3 } if
  207.     }
  208.    ifelse
  209.  } bind def
  210.  
  211. % ---------------- Top level control ---------------- %
  212.  
  213. /mapdict mark
  214.   /CIEBasedA { 1 get exch csciea currentcolorrendering 4 1 roll crmap } bind
  215.   /CIEBasedABC { 1 get 4 1 roll cscieabc currentcolorrendering 4 1 roll crmap } bind
  216. .dicttomark def
  217. /mapcie        % <a> mapcie <v1> ...
  218.         % <a> <b> <c> mapcie <v1> ...
  219.  { currentcolorspace dup 0 get //mapdict exch get exec
  220.  } bind def
  221.