/*
   Copyright (C) 1996-1998  <Uli Wortmann uliw@erdw.ethz.ch>
   Portions Copyright (C) 1999 Aladdin Enterprises.  All rights reserved.

   This program is free software; you can redistribute it and/or modify
   it under the terms of the GNU General Public License as published by
   the Free Software Foundation; either version 2 of the License, or
   (at your option) any later version.

   This program is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   GNU General Public License for more details.

   You should have received a copy of the GNU General Public License
   along with this program; if not, write to the Free Software
   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111 U.S.A.

   This program may also be distributed as part of Aladdin Ghostscript,
   under the terms of the Aladdin Free Public License (the "License").

   Every copy of Aladdin Ghostscript must include a copy of the
   License, normally in a plain ASCII text file named PUBLIC.  The
   License grants you the right to copy, modify and redistribute
   Aladdin Ghostscript, but only under certain conditions described in
   the License.  Among other things, the License requires that the
   copyright notice and this notice be preserved on all copies.
 */

/*$Id: gdevcd8.c $*/

/*
   A printer driver for the HP670, HP690, HP850, HP855
   HP870, HP890, HP1100 and HP1600 color printers. 

   To be used with the Ghostscript printing system.
   Please report all problems to uliw@erdw.ethz.ch 

   CREDITS: Much of the driver is based on ideas derived
            from the cdj550 driver of George Cameron.
	    
	    The support for the hp670, hp690, hp890
	    and hp1600 was added by Martin Gerbershagen.
*/

/* Note: Depending on how you transfered the files, 
   you might need to remove some CR-codes used on intel-based machines:

   simply type:  unzip -a hp850.zip

   to compile with gs5.x, simply add 

   DEVICE_DEVS4=cdj850.dev cdj670.dev cdj890.dev cdj1600.dev

   to your makefile. 

   BTW, it is always a good idea to read Make.htm found in the
   gs-distrib before attempting to recompile.....

 */

/* 1999-01-07 edited by L. Peter Deutsch <ghost@aladdin.com> to eliminate
   non-const statics and otherwise bring up to date with Ghostscript coding
   style guidelines. */

/* 01.06.98   Version 1.3  Due to the most welcome contribution
   of Martin Gerbershagen (ger@ulm.temic.de),
   support for the hp670, hp690 and hp890
   and hp1600 has been added. Martin has also
   resolved all known bugs.

   Problems  :  Dark colors are still pale.


   The driver no longer needs special switches to be invoked
   except -sDEVICE=cdj850, or -sDEVICE=CDJ890, or sDEVICE=CDJ670
   or -sDEVICE=CDJ1600

   The hp690 is supported through the hp670 device, the hp855, hp870 
   and the hp1100 through the hp850 device.

   The driver implements the following switches:

   -dPapertype= 0  plain paper [default]
   1  bond paper
   2  special paper
   3  glossy film
   4  transparency film

   Note, currently the lookuptables are not suited
   for printing on special paper or transperencies.
   Please revert to the gamma functions in this case.

   -dQuality=  -1 draft     
   0 normal       [default]
   1 presentation 

   -dRetStatus= 0 C-RET off
   1 C-RET on [default]

   -dMasterGamma= 3.0 [default = 1.0]
   __Note__: To take advantage of the calibrated color-transfer
   functions, be sure not to have any Gamma-Statements
   left! If you need to (i.e. overhead sheets), 
   you still can use the gamma-functions, but they will 
   override the built-in calibration. To use gamma in the 
   traditional way, set MasterGamma to any value greater
   1.0 and less 10.0. To adjust individual gamma-values,
   you have to additionally set MasterGamma to a value
   greater 1.0 and less 10.0

   With the next release, gamma functions will be dropped.

   When using the driver, be aware that printing in 600dpi involves
   processing of large amounts of data (> 188MB !). Therefore, the
   driver is not what you would expect to be a fast driver ;-)
   This is no problem when printing a full sized color page (because
   printing itself is slow), but it's really annoying if yoy print only
   text pages. Maybe I can optimize the code for text-only pages in a
   later release. Right now, it is recommended to use the highest
   possible optimisation level your compiler offers....
   For the time beeing, use the cdj550 device with -sBitsPerPixel=3
   for fast proof-prints. If you simply want to print 600dpi b/w data,
   use the cdj550 device with -sBitsPerPixel=8 (or 1).

   Since the printer itself is slow, it may help to set the
   process-priority of the gs-process to regular or even less. On a
   486/100MHZ this is still sufficient to maintain a continuos
   data-flow.
   Note to OS/2 users: Simply put the gs-window into the background,
   or minimize it. Also make sure, that print01.sys is invoked without
   the /irq switch (great speed improvement under warp4).

   The printer default settings compensate for dot-gain by a
   calibrated color-transfer function. If this appears to be to light
   for your business-graphs, or for overhead-sheets, feel free to set
   -dMasterGamma=1.7.

   Furthermore, you may tweak the gammavalues independently by setting
   -dGammaValC, -dGammaValM, -dGammaValY or -dGammaValK (if not set,
   the values default to MasterGamma). This will only work, when
   -dMasterGamma is set to a value greater than 1.0.

   If you want to learn more about gamma, see:

   http://www.erdw.ethz.ch/~bonk/ftp/misc/gammafaq.pdf
   http://www.erdw.ethz.ch/~bonk/ftp/misc/colorfaq.pdf

   Further information, bugs, tips etc, can be found
   at my website.

   Have fun!

   Uli

   uliw@erdw.ethz.ch
   http://www.erdw.ethz.ch/~bonk/bonk.html

 */

/* 25.08.97  Version 1.2. Resolved all but one of the 
   known bugs, introduced a couple
   of perfomance improvements. Complete
   new color-transfer-function handling.
   (see gamma). */

/* 04.05.97  Version 1.1. For added features, */
/* resolved bugs and so forth, please see   */
/* http://bonk.ethz.ch                      */

/* 11.11.96. Initial release of the driver */

#include "math_.h"
#include <stdlib.h>		/* for rand() */
#include <assert.h>
#include "gdevprn.h"
#include "gdevpcl.h"
#include "gsparam.h"

/* Conversion stuff. */
#include "gxlum.h"

/* this holds the initialisation data of the hp850 */
typedef struct hp850_cmyk_init_s {
    byte a[26];
} hp850_cmyk_init_t;
private const hp850_cmyk_init_t hp850_cmyk_init =
{
    {
	0x02,			/* format */
	0x04,			/* number of components */
      /* black */
	0x01,			/* MSB x resolution */
	0x2c,			/* LSB x resolution */
	0x01,			/* MSB y resolution */
	0x2c,			/* LSB y resolution */
	0x00,			/* MSB intensity levels */
	0x02,			/* LSB intensity levels */

      /* cyan */
	0x01,			/* MSB x resolution */
	0x2c,			/* LSB x resolution */
	0x01,			/* MSB y resolution */
	0x2c,			/* LSB y resolution */
	0x00,			/* MSB intensity levels */
	0x02,			/* LSB intensity levels */

      /* magenta */
	0x01,			/* MSB x resolution */
	0x2c,			/* LSB x resolution */
	0x01,			/* MSB y resolution */
	0x2c,			/* LSB y resolution */
	0x00,			/* MSB intensity levels */
	0x02,			/* LSB intensity levels */

      /* yellow */
	0x01,			/* MSB x resolution */
	0x2c,			/* LSB x resolution */
	0x01,			/* MSB y resolution */
	0x2c,			/* LSB y resolution */
	0x00,			/* MSB intensity levels */
	0x02			/* LSB intensity levels */
    }
};

/* this holds the color lookuptable data of the hp850 */
typedef struct {
    byte c[256];		/* Lookuptable for cyan */
    byte m[256];		/* dito for magenta */
    byte y[256];		/* dito for yellow */
    byte k[256];		/* dito for black  */
    int correct[256];		/* potential undercolor black correction */
} Gamma;

private const Gamma gammat850 =
{
  /* Lookup values for cyan */
    {0, 0, 0, 2, 2, 2, 3, 3, 3, 5, 5, 5, 7, 7, 6, 7, 7, 6, 7, 7, 7, 8, 8,
     8, 8, 8, 8, 9, 9, 9, 9, 9, 10, 10, 10, 10, 10, 10, 11, 11, 12, 12, 12,
     12, 12, 12, 13, 13, 14, 14, 14, 15, 15, 16, 16, 15, 16, 16, 17, 17,
     17, 17, 17, 18, 18, 18, 19, 19, 20, 20, 20, 20, 20, 21, 21, 21, 22,
     22, 23, 23, 23, 23, 23, 24, 24, 25, 25, 26, 26, 26, 26, 26, 27, 27,
     27, 27, 28, 28, 29, 28, 28, 29, 29, 30, 30, 31, 31, 32, 32, 33, 34,
     35, 35, 36, 36, 37, 37, 38, 38, 39, 39, 40, 40, 41, 41, 42, 42, 42,
     43, 43, 43, 44, 45, 45, 46, 46, 47, 47, 48, 48, 49, 50, 50, 51, 51,
     52, 52, 53, 54, 54, 54, 55, 55, 56, 57, 58, 58, 59, 60, 60, 61, 62,
     62, 63, 65, 65, 66, 67, 67, 68, 69, 69, 70, 72, 73, 73, 74, 75, 75,
     76, 77, 79, 79, 80, 81, 82, 83, 83, 84, 86, 87, 88, 88, 89, 90, 91,
     92, 93, 94, 95, 96, 97, 97, 99, 100, 101, 102, 103, 104, 105, 106,
     108, 109, 110, 111, 112, 114, 115, 117, 119, 120, 122, 124, 125, 127,
     129, 131, 132, 135, 136, 138, 140, 142, 144, 146, 147, 150, 152, 154,
     157, 159, 162, 164, 166, 168, 171, 174, 176, 180, 182, 187, 192, 197,
     204, 215, 255},
  /* Lookup values for magenta */
    {0, 0, 0, 1, 1, 1, 3, 3, 3, 4, 4, 4, 5, 5, 5, 6, 6, 6, 7, 7, 7, 7, 7,
     7, 8, 8, 8, 9, 9, 10, 10, 9, 10, 10, 10, 11, 11, 11, 11, 11, 12, 12,
     12, 13, 13, 13, 14, 14, 15, 15, 15, 16, 16, 16, 16, 16, 17, 17, 17,
     17, 17, 18, 18, 19, 19, 19, 19, 19, 20, 20, 20, 21, 21, 22, 22, 22,
     23, 23, 24, 24, 25, 25, 25, 26, 26, 27, 27, 28, 29, 29, 29, 29, 30,
     30, 31, 30, 31, 31, 32, 31, 31, 32, 32, 33, 33, 34, 34, 35, 35, 36,
     36, 37, 37, 38, 38, 39, 39, 40, 40, 41, 41, 42, 42, 43, 43, 44, 44,
     45, 45, 46, 46, 47, 48, 48, 49, 49, 50, 50, 51, 51, 52, 53, 53, 54,
     54, 55, 55, 56, 57, 57, 58, 58, 59, 60, 60, 61, 61, 62, 63, 64, 65,
     66, 66, 67, 68, 68, 70, 71, 71, 72, 73, 73, 74, 76, 77, 77, 78, 79,
     79, 80, 81, 82, 83, 84, 85, 86, 87, 87, 88, 89, 90, 91, 91, 92, 93,
     94, 95, 96, 97, 98, 99, 100, 100, 101, 102, 103, 105, 106, 107, 108,
     109, 112, 113, 114, 115, 116, 118, 119, 121, 123, 124, 125, 128, 129,
     130, 133, 134, 135, 138, 139, 142, 144, 145, 148, 150, 152, 154, 157,
     159, 162, 164, 168, 169, 170, 172, 175, 177, 179, 182, 185, 189, 193,
     198, 204, 215, 255},
  /* Lookup values for yellow */
    {0, 0, 0, 2, 2, 2, 3, 3, 3, 5, 5, 5, 7, 7, 6, 7, 7, 6, 7, 7, 7, 8, 8,
     8, 8, 8, 8, 9, 9, 9, 9, 9, 10, 9, 9, 10, 10, 10, 10, 10, 11, 11, 11,
     12, 12, 13, 13, 14, 14, 15, 15, 16, 16, 16, 16, 16, 17, 17, 18, 18,
     18, 19, 18, 19, 19, 19, 20, 20, 21, 21, 21, 22, 22, 22, 22, 22, 23,
     23, 24, 24, 25, 25, 25, 26, 27, 28, 28, 29, 29, 29, 30, 30, 30, 30,
     31, 31, 32, 32, 33, 33, 32, 33, 33, 34, 34, 35, 35, 36, 36, 37, 37,
     38, 38, 39, 39, 40, 40, 41, 41, 42, 42, 43, 43, 44, 44, 45, 45, 45,
     45, 46, 46, 47, 48, 48, 49, 49, 50, 50, 51, 51, 52, 53, 53, 54, 54,
     55, 55, 56, 57, 58, 59, 59, 60, 61, 61, 62, 62, 63, 64, 65, 66, 67,
     67, 68, 69, 69, 70, 71, 72, 73, 74, 74, 75, 76, 77, 77, 78, 79, 79,
     80, 81, 82, 83, 84, 85, 86, 87, 87, 88, 89, 90, 91, 91, 93, 94, 95,
     96, 97, 98, 100, 101, 102, 102, 103, 104, 106, 107, 108, 109, 110,
     111, 113, 114, 115, 116, 117, 118, 119, 121, 123, 124, 126, 128, 130,
     131, 134, 135, 137, 139, 140, 143, 145, 146, 148, 150, 152, 154, 156,
     158, 160, 163, 166, 167, 169, 171, 173, 176, 178, 181, 184, 188, 192,
     198, 204, 215, 255},
  /* Lookup values for black */
    {0, 0, 0, 1, 1, 1, 1, 1, 1, 3, 3, 3, 3, 3, 2, 4, 3, 3, 3, 3, 3, 4, 4,
     4, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 8, 8, 8, 9, 9, 8,
     8, 8, 9, 9, 9, 10, 10, 10, 10, 10, 11, 11, 11, 11, 12, 12, 12, 13, 13,
     12, 12, 12, 13, 13, 13, 13, 13, 14, 14, 14, 14, 14, 15, 15, 16, 16,
     16, 17, 17, 17, 17, 18, 18, 18, 19, 19, 20, 20, 20, 20, 20, 21, 21,
     21, 21, 22, 22, 22, 22, 23, 22, 23, 23, 24, 24, 24, 24, 25, 25, 26,
     26, 26, 26, 27, 27, 28, 28, 28, 28, 29, 29, 30, 30, 31, 31, 31, 32,
     32, 33, 33, 34, 34, 35, 36, 36, 36, 37, 37, 37, 38, 38, 40, 40, 40,
     41, 41, 42, 43, 43, 43, 43, 44, 45, 45, 46, 47, 47, 48, 49, 49, 50,
     52, 52, 53, 54, 54, 56, 56, 57, 58, 59, 60, 60, 61, 62, 63, 63, 64,
     65, 66, 67, 68, 69, 70, 71, 72, 72, 73, 75, 75, 76, 77, 78, 80, 81,
     82, 82, 83, 84, 85, 86, 88, 89, 90, 91, 94, 95, 96, 98, 99, 100, 101,
     103, 105, 106, 107, 110, 111, 112, 115, 116, 118, 120, 121, 124, 126,
     127, 131, 133, 134, 138, 140, 141, 146, 148, 151, 154, 156, 160, 163,
     166, 169, 174, 177, 182, 187, 194, 203, 215, 255}
};

private const Gamma gammat890 =
{
/* Lookup values for cyan */
    {0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15,
     16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31,
     32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47,
     48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63,
     64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79,
     80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95,
     96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111,
     112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125,
     126, 127,
     128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141,
     142, 143,
     144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157,
     158, 159,
     160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173,
     174, 175,
     176, 177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 187, 188, 189,
     190, 191,
     192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204, 205,
     206, 207,
     208, 209, 210, 211, 212, 213, 214, 215, 216, 217, 218, 219, 220, 221,
     222, 223,
     224, 225, 226, 227, 228, 229, 230, 231, 232, 233, 234, 235, 236, 237,
     238, 239,
     240, 241, 242, 243, 244, 245, 246, 247, 248, 249, 250, 251, 252, 253,
     254, 255},

/* Lookup values for magenta */
    {0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15,
     16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31,
     32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47,
     48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63,
     64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79,
     80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95,
     96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111,
     112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125,
     126, 127,
     128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141,
     142, 143,
     144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157,
     158, 159,
     160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173,
     174, 175,
     176, 177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 187, 188, 189,
     190, 191,
     192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204, 205,
     206, 207,
     208, 209, 210, 211, 212, 213, 214, 215, 216, 217, 218, 219, 220, 221,
     222, 223,
     224, 225, 226, 227, 228, 229, 230, 231, 232, 233, 234, 235, 236, 237,
     238, 239,
     240, 241, 242, 243, 244, 245, 246, 247, 248, 249, 250, 251, 252, 253,
     254, 255},

/* Lookup values for yellow */
    {0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15,
     16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31,
     32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47,
     48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63,
     64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79,
     80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95,
     96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111,
     112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125,
     126, 127,
     128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141,
     142, 143,
     144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157,
     158, 159,
     160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173,
     174, 175,
     176, 177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 187, 188, 189,
     190, 191,
     192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204, 205,
     206, 207,
     208, 209, 210, 211, 212, 213, 214, 215, 216, 217, 218, 219, 220, 221,
     222, 223,
     224, 225, 226, 227, 228, 229, 230, 231, 232, 233, 234, 235, 236, 237,
     238, 239,
     240, 241, 242, 243, 244, 245, 246, 247, 248, 249, 250, 251, 252, 253,
     254, 255},

 /* Lookup values for black */
    {0, 0, 0, 1, 1, 1, 1, 1, 1, 3, 3, 3, 3, 3, 2, 4, 3, 3, 3, 3, 3, 4, 4,
     4, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 8, 8, 8, 9, 9, 8,
     8, 8, 9, 9, 9, 10, 10, 10, 10, 10, 11, 11, 11, 11, 12, 12, 12, 13, 13,
     12, 12, 12, 13, 13, 13, 13, 13, 14, 14, 14, 14, 14, 15, 15, 16, 16,
     16, 17, 17, 17, 17, 18, 18, 18, 19, 19, 20, 20, 20, 20, 20, 21, 21,
     21, 21, 22, 22, 22, 22, 23, 22, 23, 23, 24, 24, 24, 24, 25, 25, 26,
     26, 26, 26, 27, 27, 28, 28, 28, 28, 29, 29, 30, 30, 31, 31, 31, 32,
     32, 33, 33, 34, 34, 35, 36, 36, 36, 37, 37, 37, 38, 38, 40, 40, 40,
     41, 41, 42, 43, 43, 43, 43, 44, 45, 45, 46, 47, 47, 48, 49, 49, 50,
     52, 52, 53, 54, 54, 56, 56, 57, 58, 59, 60, 60, 61, 62, 63, 63, 64,
     65, 66, 67, 68, 69, 70, 71, 72, 72, 73, 75, 75, 76, 77, 78, 80, 81,
     82, 82, 83, 84, 85, 86, 88, 89, 90, 91, 94, 95, 96, 98, 99, 100, 101,
     103, 105, 106, 107, 110, 111, 112, 115, 116, 118, 120, 121, 124, 126,
     127, 131, 133, 134, 138, 140, 141, 146, 148, 151, 154, 156, 160, 163,
     166, 169, 174, 177, 182, 187, 194, 203, 215, 255}
};

private const Gamma * const gammat[] =
{
    &gammat850,			/* CDJ670 */
    &gammat850,			/* CDJ850 */
    &gammat890,			/* CDJ890 */
    &gammat850			/* CDJ1600 */
};

private int
    rescale_byte_wise1x1(P4(int bytecount, const byte * inbytea,
			    const byte * inbyteb, byte * outbyte));
private int
    rescale_byte_wise2x1(P4(int bytecount, const byte * inbytea,
			    const byte * inbyteb, byte * outbyte));
private int
    rescale_byte_wise1x2(P4(int bytecount, const byte * inbytea,
			    const byte * inbyteb, byte * outbyte));
private int
    rescale_byte_wise2x2(P4(int bytecount, const byte * inbytea,
			    const byte * inbyteb, byte * outbyte));

private int (* const rescale_color_plane[2][2]) (P4(int, const byte *, const byte *, byte *)) = {
    {
	rescale_byte_wise1x1, rescale_byte_wise1x2
    },
    {
	rescale_byte_wise2x1, rescale_byte_wise2x2
    }
};

/*
 * Drivers stuff.
 *
 */
#define DESKJET_PRINT_LIMIT  0.04	/* 'real' top margin? */
/* Margins are left, bottom, right, top. */
#define DESKJET_MARGINS_LETTER   0.25, 0.50, 0.25, 0.167
#define DESKJET_MARGINS_A4       0.13, 0.46, 0.13, 0.04
/* Define bits-per-pixel - default is 32-bit cmyk-mode */
#ifndef BITSPERPIXEL
#  define BITSPERPIXEL 32
#endif
#define DOFFSET (dev_t_margin(pdev) - DESKJET_PRINT_LIMIT)	/* Print position */


#define W sizeof(word)
#define I sizeof(int)

/* paper types */
typedef enum {
    PLAIN_PAPER, BOND_PAPER, SPECIAL_PAPER, GLOSSY_FILM, TRANSPARENCY_FILM
} cdj_paper_type_t;

/* quality */
typedef enum {
    DRAFT = -1, NORMAL = 0, PRESENTATION = 1
} cdj_quality_t;

/* Printer types */
typedef enum {
    DJ670C, DJ850C, DJ890C, DJ1600C
} cdj_printer_type_t;

/* No. of ink jets (used to minimise head movements) */
#define HEAD_ROWS_MONO 50
#define HEAD_ROWS_COLOUR 16

/* Colour mapping procedures */
private dev_proc_map_cmyk_color(gdev_cmyk_map_cmyk_color);
private dev_proc_map_rgb_color(gdev_cmyk_map_rgb_color);
private dev_proc_map_color_rgb(gdev_cmyk_map_color_rgb);

private dev_proc_map_rgb_color(gdev_pcl_map_rgb_color);
private dev_proc_map_color_rgb(gdev_pcl_map_color_rgb);

/* Print-page, parameters and miscellaneous procedures */
private dev_proc_open_device(hp_colour_open);

private dev_proc_get_params(cdj850_get_params);
private dev_proc_put_params(cdj850_put_params);

private dev_proc_print_page(cdj850_print_page);

/* The device descriptors */

/* The basic structure for all printers. Note the presence of the cmyk, depth
   and correct fields even if some are not used by all printers. */

#define prn_colour_device_body(dtype, procs, dname, w10, h10, xdpi, ydpi, lm, bm, rm, tm, ncomp, depth, mg, mc, dg, dc, print_page, cmyk, correct)\
    prn_device_body(dtype, procs, dname, w10, h10, xdpi, ydpi, lm, bm, rm, tm, ncomp, depth, mg, mc, dg, dc, print_page), cmyk, depth /* default */, correct



#define gx_prn_colour_device_common \
    gx_prn_device_common; \
    int cmyk;	  	/* 0: not CMYK-capable, > 0: printing CMYK, */ \
		  	/* < 0 : CMYK-capable, not printing CMYK */ \
    uint default_depth;	/* Used only for CMYK-capable printers now. */ \
    uint correction


/* some definitions needed later */
struct error_val_field {
    int c;		/* Current value of Cyan error during dithering */
    int m;		/* Current value of Magenta error during dithering */
    int y;		/* Current value of Yellow error during dithering */
    int k;		/* Current value of Black error during dithering */
};

/* this structure holds all the pointers to the different values
   in all those data fields */
 /*
    * The principal data pointers are stored as pairs of values, with
    * the selection being made by the 'scan' variable. The function of the
    * scan variable is overloaded, as it controls both the alternating
    * raster scan direction used in the Floyd-Steinberg dithering and also
    * the buffer alternation required for line-difference compression.
    *
    * Thus, the number of pointers required is as follows:
  */

struct ptr_arrays {
    byte *data[4];		/* 4 600dpi data, scan direction and alternating buffers */
    byte *data_c[4];		/* 4 300dpi data, as above, */
    byte *plane_data[4][4];	/*4 b/w-planes, scan direction and alternating buffers */
    byte *plane_data_c[4][8];	/* as above, but for 8 planes */
    byte *out_data;		/* output buffer for the b/w data, one 600dpi plane */
    byte *test_data[4];		/* holds a copy of the last plane */
    int *errors[2];		/* 2 b/w dithering erros (scan direction only) */
    int *errors_c[2];		/* 2 color dithering errors (scan direction only) */
    word *storage;		/* pointer to the beginning of the b/w-buffer */
    word *storage_start;	/* used for debugging */
    word *storage_end;		/* used for debugging */
    word *storage_size;		/* used for debugging */
};

/* Some miscellaneous variables */
struct misc_struct {
    int line_size;		/* size of scan_line */
    int line_size_c;		/* size of rescaled scan_line */
    int line_size_words;	/* size of scan_line in words */
    int paper_size;		/* size of paper */
    int num_comps;		/* number of color components (1 - 4) */
    int bits_per_pixel;		/* bits per pixel 1,4,8,16,24,32 */
    int storage_bpp;		/* = bits_per_pixel */
    int expanded_bpp;		/* = bits_per_pixel */
    int plane_size;		/* size of b/w bit plane */
    int plane_size_c;		/* size of color bit plane */
    int databuff_size;		/* size of databuffer for b/w data */
    int databuff_size_c;	/* size of databuffer for color data */
    int errbuff_size;		/* size of error buffer b/w -data */
    int errbuff_size_c;		/* size of error buffer color -data */
    int outbuff_size;		/* size of output buffer for b/w data */
    int scan;			/* scan-line variable [0,1] */
    int cscan;			/* dito for the color-planes */
    int is_two_pass;		/* checks if b/w data has already been printed */
    int zero_row_count;		/* How many empty lines */
    uint storage_size_words;	/* size of storage in words for b/w data */
    uint storage_size_words_c;	/* size of storage in words for c-data */
    int is_color_data;		/* indicates whether there is color data */
};

    /* function pointer typedefs for device driver struct */
typedef void (*StartRasterMode) (P3(gx_device_printer * pdev, int paper_size,
				    FILE * prn_stream));
typedef void (*PrintNonBlankLines) (P6(gx_device_printer * pdev,
				       struct ptr_arrays *data_ptrs,
				       struct misc_struct *misc_vars,
				       struct error_val_field *error_values,
				       const Gamma *gamma,
				       FILE * prn_stream));

typedef void (*TerminatePage) (P2(gx_device_printer * pdev, FILE * prn_stream));

typedef struct gx_device_cdj850_s {
    gx_device_common;
    gx_prn_colour_device_common;
    int /*cdj_quality_t*/ quality;  /* -1 draft, 0 normal, 1 best */
    int /*cdj_paper_type_t*/ papertype;  /* papertype [0,4] */
    int intensities;		/* intensity values per pixel [2,4] */
    int xscal;			/* boolean to indicate x scaling by 2 */
    int yscal;			/* boolean to indicate y scaling by 2 */
    int /*cdj_printer_type_t*/ ptype;  /* printer type, one of DJ670C, DJ850C, DJ890C, DJ1600C */
    int compression;		/* compression level */
    float mastergamma;		/* Gammavalue applied to all colors */
    float gammavalc;		/* range to which gamma-correction is
				   applied to bw values */
    float gammavalm;		/* amount of gamma correction for bw */
    float gammavaly;		/* range to which gamma-correction i
				   applied to color values */
    float gammavalk;		/* amount of gamma correction for color */
    float blackcorrect;		/* amount of gamma correction for color */
    StartRasterMode start_raster_mode;	/* output function to start raster mode */
    PrintNonBlankLines print_non_blank_lines;	/* output function to print a non blank line */
    TerminatePage terminate_page;	/* page termination output function */
} gx_device_cdj850;

typedef struct {
    gx_device_common;
    gx_prn_colour_device_common;
} gx_device_colour_prn;


/* Use the cprn_device macro to access generic fields (like cmyk,
   default_depth and correction), and specific macros for specific
   devices. */

#define cprn_device     ((gx_device_colour_prn*) pdev)
#define cdj850    ((gx_device_cdj850 *)pdev)

#define prn_cmyk_colour_device(dtype, procs, dev_name, x_dpi, y_dpi, bpp, print_page, correct)\
    prn_colour_device_body(dtype, procs, dev_name,\
    DEFAULT_WIDTH_10THS, DEFAULT_HEIGHT_10THS, x_dpi, y_dpi, 0, 0, 0, 0,\
    ((bpp == 1 || bpp == 4) ? 1 : 4), bpp,\
    (bpp > 8 ? 255 : 1), (1 << (bpp >> 2)) - 1, /* max_gray, max_color */\
    (bpp > 8 ? 5 : 2), (bpp > 8 ? 5 : bpp > 1 ? 2 : 0),\
    print_page, 1 /* cmyk */, correct)


#define prn_cmy_colour_device(dtype, procs, dev_name, x_dpi, y_dpi, bpp, print_page, correct)\
    prn_colour_device_body(dtype, procs, dev_name,\
    DEFAULT_WIDTH_10THS, DEFAULT_HEIGHT_10THS, x_dpi, y_dpi, 0, 0, 0, 0,\
    ((bpp == 1 || bpp == 4) ? 1 : 3), bpp,\
    (bpp > 8 ? 255 : 1), (bpp > 8 ? 255 : 1), /* max_gray, max_color */\
    (bpp > 8 ? 5 : 2), (bpp > 8 ? 5 : bpp > 1 ? 2 : 0),\
    print_page, -1 /* cmyk */, correct)


#define cdj_850_device(procs, dev_name, x_dpi, y_dpi, bpp, print_page, correction, quality, papertype, intensities,ptype,compression,mastergamma,gammavalc,gammavalm,gammavaly,gammavalk,blackcorrect,start_raster_mode,print_non_blank_line,terminate_page)\
{ prn_cmyk_colour_device(gx_device_cdj850, procs, dev_name, x_dpi, y_dpi, bpp, print_page, correction),\
    quality,\
    papertype,\
    intensities,\
    0, 0, /* xscal, yscal */\
    ptype,\
    compression,\
    mastergamma,\
    gammavalc,\
    gammavalm,\
    gammavaly,\
    gammavalk,\
    blackcorrect,\
    start_raster_mode,\
    print_non_blank_line,\
    terminate_page\
}

#define cdj_1600_device(procs, dev_name, x_dpi, y_dpi, bpp, print_page, correction, quality, papertype, intensities,ptype,compression,mastergamma,gammavalc,gammavalm,gammavaly,gammavalk,blackcorrect,start_raster_mode,print_non_blank_line,terminate_page)\
{ prn_cmy_colour_device(gx_device_cdj850, procs, dev_name, x_dpi, y_dpi, bpp, print_page, correction),\
    quality,\
    papertype,\
    intensities,\
    0, 0, /* xscal, yscal */\
    ptype,\
    compression,\
    mastergamma,\
    gammavalc,\
    gammavalm,\
    gammavaly,\
    gammavalk,\
    blackcorrect,\
    start_raster_mode,\
    print_non_blank_line,\
    terminate_page\
}

#define cmyk_colour_procs(proc_colour_open, proc_get_params, proc_put_params, \
			  map_rgb_color, map_color_rgb, map_cmyk_color) {\
	proc_colour_open,\
	gx_default_get_initial_matrix,\
	gx_default_sync_output,\
	gdev_prn_output_page,\
	gdev_prn_close,\
	map_rgb_color,\
	map_color_rgb,\
	NULL /* fill_rectangle */,\
	NULL /* tile_rectangle */,\
	NULL /* copy_mono */,\
	NULL /* copy_color */,\
	NULL /* draw_line */,\
	gx_default_get_bits,\
	proc_get_params,\
	proc_put_params,\
        map_cmyk_color\
}

private void
     cdj850_start_raster_mode(P3(gx_device_printer * pdev,
				 int papersize, FILE * prn_stream));

private void
     cdj850_print_non_blank_lines(P6(gx_device_printer * pdev,
				     struct ptr_arrays *data_ptrs,
				     struct misc_struct *misc_vars,
				     struct error_val_field *error_values,
				     const Gamma *gamma,
				     FILE * prn_stream));
private void
     cdj850_terminate_page(P2(gx_device_printer * pdev, FILE * prn_stream));

private void
     cdj1600_start_raster_mode(P3(gx_device_printer * pdev,
				  int papersize, FILE * prn_stream));
private void
     cdj1600_print_non_blank_lines(P6(gx_device_printer * pdev,
				      struct ptr_arrays *data_ptrs,
				      struct misc_struct *misc_vars,
				      struct error_val_field *error_values,
				      const Gamma *gamma,
				      FILE * prn_stream));
private void
     cdj1600_terminate_page(P2(gx_device_printer * pdev, FILE * prn_stream));



private const gx_device_procs cdj670_procs =
cmyk_colour_procs(hp_colour_open, cdj850_get_params, cdj850_put_params,
		  NULL, gdev_cmyk_map_color_rgb, gdev_cmyk_map_cmyk_color);

private const gx_device_procs cdj850_procs =
cmyk_colour_procs(hp_colour_open, cdj850_get_params, cdj850_put_params,
		  NULL, gdev_cmyk_map_color_rgb, gdev_cmyk_map_cmyk_color);

private const gx_device_procs cdj890_procs =
cmyk_colour_procs(hp_colour_open, cdj850_get_params, cdj850_put_params,
		  NULL, gdev_cmyk_map_color_rgb, gdev_cmyk_map_cmyk_color);

private const gx_device_procs cdj1600_procs =
cmyk_colour_procs(hp_colour_open, cdj850_get_params, cdj850_put_params,
		  gdev_pcl_map_rgb_color, gdev_pcl_map_color_rgb, NULL);

const gx_device_cdj850 gs_cdj670_device =
cdj_850_device(cdj670_procs, "cdj670", 600, 600, 32, cdj850_print_page, 0,
	       PRESENTATION, PLAIN_PAPER, 2, DJ670C, 9,
	       1.0, 0.0, 0.0, 0.0, 0.0, 1.0,
	       cdj850_start_raster_mode, cdj850_print_non_blank_lines,
	       cdj850_terminate_page);

const gx_device_cdj850 gs_cdj850_device =
cdj_850_device(cdj850_procs, "cdj850", 600, 600, 32, cdj850_print_page, 0,
	       PRESENTATION, PLAIN_PAPER, 4, DJ850C, 9,
	       1.0, 0.0, 0.0, 0.0, 0.0, 1.0,
	       cdj850_start_raster_mode, cdj850_print_non_blank_lines,
	       cdj850_terminate_page);

const gx_device_cdj850 gs_cdj890_device =
cdj_850_device(cdj890_procs, "cdj890", 600, 600, 32, cdj850_print_page, 0,
	       PRESENTATION, PLAIN_PAPER, 4, DJ890C, 9,
	       1.0, 0.0, 0.0, 0.0, 0.0, 1.0,
	       cdj850_start_raster_mode, cdj850_print_non_blank_lines,
	       cdj850_terminate_page);

const gx_device_cdj850 gs_cdj1600_device =
cdj_1600_device(cdj1600_procs, "cdj1600", 300, 300, 24, cdj850_print_page, 0,
		PRESENTATION, PLAIN_PAPER, 2, DJ1600C, 3,
		1.0, 0.0, 0.0, 0.0, 0.0, 1.0,
		cdj1600_start_raster_mode, cdj1600_print_non_blank_lines,
		cdj1600_terminate_page);

/* Forward references */
private int cdj_put_param_int(P6(gs_param_list *, gs_param_name,
				 int *, int, int, int));
private int cdj_put_param_float(P6(gs_param_list *, gs_param_name, float
				   *, float, float, int));
private int cdj_put_param_bpp(P5(gx_device *, gs_param_list *, int, int, int));
private int cdj_set_bpp(P3(gx_device *, int, int));


/* Open the printer and set up the margins. */
private int
hp_colour_open(gx_device * pdev)
{				/* Change the margins if necessary. */
    static const float dj_a4[4] = {
	DESKJET_MARGINS_A4
    };
    static const float dj_letter[4] = {
	DESKJET_MARGINS_LETTER
    };

    /* margins for DJ1600C from manual */
    static const float m_cdj1600[4] = {
	0.25, 0.5, 0.25, 0.5
    };

    const float *m = (float *)0;

    /* Set up colour params if put_params has not already done so */
    if (pdev->color_info.num_components == 0) {
	int code = cdj_set_bpp(pdev, pdev->color_info.depth,
			       pdev->color_info.num_components);

	if (code < 0)
	    return code;
    }
    /* assign printer type and set resolution dependent on printer type */
    switch (cdj850->ptype) {
    case DJ670C:
	if (cdj850->papertype <= SPECIAL_PAPER) {	/* paper */
	    if (cdj850->quality == DRAFT) {
		gx_device_set_resolution(pdev, 300.0, 300.0);
		cdj850->xscal = 0;
		cdj850->yscal = 0;
	    } else if (cdj850->quality == NORMAL) {
		gx_device_set_resolution(pdev, 600.0, 300.0);
		cdj850->xscal = 1;
		cdj850->yscal = 0;
	    } else {		/* quality == PRESENTATION */
		gx_device_set_resolution(pdev, 600.0, 600.0);
		cdj850->xscal = 1;
		cdj850->yscal = 1;
	    }
	} else {		/* film */
	    gx_device_set_resolution(pdev, 600.0, 300.0);
	    cdj850->xscal = 0;
	    cdj850->yscal = 0;
	}
	m = (gdev_pcl_paper_size(pdev) == PAPER_SIZE_A4 ? dj_a4 : dj_letter);
	break;
    case DJ850C:
    case DJ890C:
	if (cdj850->quality == DRAFT) {
	    gx_device_set_resolution(pdev, 300.0, 300.0);
	    cdj850->xscal = 0;
	    cdj850->yscal = 0;
	    cdj850->intensities = 2;
	} else if (cdj850->quality == NORMAL) {
	    gx_device_set_resolution(pdev, 600.0, 300.0);
	    cdj850->xscal = 1;
	    cdj850->yscal = 0;
	    /* only 3 intensities for normal paper */
	    if (cdj850->papertype <= PLAIN_PAPER) {
		cdj850->intensities = 3;
	    }			/* else cdj850->intensities = 4 from initialization */
	} else {		/* quality == PRESENTATION */
	    gx_device_set_resolution(pdev, 600.0, 600.0);
	    cdj850->xscal = 1;
	    cdj850->yscal = 1;
	    /* intensities = 4 from initialization */
	}
	m = (gdev_pcl_paper_size(pdev) == PAPER_SIZE_A4 ? dj_a4 : dj_letter);
	break;
    case DJ1600C:
	gx_device_set_resolution(pdev, 300.0, 300.0);
	m = m_cdj1600;
	break;
    default:
	assert(0);
    }
    gx_device_set_margins(pdev, m, true);
    return gdev_prn_open(pdev);
}

/* Added parameters for DeskJet 850C */
private int
cdj850_get_params(gx_device * pdev, gs_param_list * plist)
{
    int code = gdev_prn_get_params(pdev, plist);

    if (code < 0 ||
	(code = param_write_int(plist, "Quality", &cdj850->quality)) < 0 ||
	(code = param_write_int(plist, "Papertype", &cdj850->papertype)) < 0 ||
	(code = param_write_float(plist, "MasterGamma", &cdj850->gammavalc))
	< 0 ||
	(code = param_write_float(plist, "GammaValC", &cdj850->gammavalc)) <
	0 ||
	(code = param_write_float(plist, "GammaValM", &cdj850->gammavalm)) <
	0 ||
	(code = param_write_float(plist, "GammaValY", &cdj850->gammavaly)) <
	0 ||
	(code = param_write_float(plist, "GammaValK", &cdj850->gammavalk)) <
	0 ||
	(code = param_write_float(plist, "BlackCorrect",
				  &cdj850->blackcorrect)) < 0
	)
	return code;

    return code;
}

private int
cdj850_put_params(gx_device * pdev, gs_param_list * plist)
{
    int quality = cdj850->quality;
    int papertype = cdj850->papertype;
    float mastergamma = cdj850->mastergamma;
    float gammavalc = cdj850->gammavalc;
    float gammavalm = cdj850->gammavalm;
    float gammavaly = cdj850->gammavaly;
    float gammavalk = cdj850->gammavalk;
    float blackcorrect = cdj850->blackcorrect;
    int bpp = 0;
    int code = 0;

    code = cdj_put_param_int(plist, "BitsPerPixel", &bpp, 1, 32, code);
    code = cdj_put_param_int(plist, "Quality", &quality, 0, 2, code);
    code = cdj_put_param_int(plist, "Papertype", &papertype, 0, 4, code);
    code = cdj_put_param_float(plist, "MasterGamma", &mastergamma, 0.1, 9.0, code);
    code = cdj_put_param_float(plist, "GammaValC", &gammavalc, 0.0, 9.0, code);
    code = cdj_put_param_float(plist, "GammaValM", &gammavalm, 0.0, 9.0, code);
    code = cdj_put_param_float(plist, "GammaValY", &gammavaly, 0.0, 9.0, code);
    code = cdj_put_param_float(plist, "GammaValK", &gammavalk, 0.0, 9.0, code);
    code = cdj_put_param_float(plist, "BlackCorrect", &blackcorrect, 0.0,
			       9.0, code);


    if (code < 0)
	return code;
    code = cdj_put_param_bpp(pdev, plist, bpp, bpp, 0);
    if (code < 0)
	return code;

    cdj850->quality = quality;
    cdj850->papertype = papertype;
    cdj850->mastergamma = mastergamma;
    cdj850->gammavalc = gammavalc;
    cdj850->gammavalm = gammavalm;
    cdj850->gammavaly = gammavaly;
    cdj850->gammavalk = gammavalk;
    cdj850->blackcorrect = blackcorrect;
    return 0;
}

/* ------ Internal routines ------ */
/* The DeskJet850C can compress (mode 9) */


/* Some convenient shorthand .. */
#define x_dpi        (pdev->x_pixels_per_inch)
#define y_dpi        (pdev->y_pixels_per_inch)

/* To calculate buffer size as next greater multiple of both parameter and W */
#define calc_buffsize(a, b) (((((a) + ((b) * W) - 1) / ((b) * W))) * W)

/* internal functions */
private void
     FSDlinebw(P7(int scan, int plane_size,
		  struct error_val_field *error_values,
		  byte * kP,
		  int n, int *ep, byte * dp));
private void
     FSDlinec2(P9(int scan, int plane_size,
		  struct error_val_field *error_values,
		  byte * cPa, byte * mPa, byte * yPa, int n,
		  byte * dp, int *ep));
private void
     FSDlinec3(P12(int scan, int plane_size,
		   struct error_val_field *error_values,
		   byte * cPa, byte * mPa, byte * yPa,
		   byte * cPb, byte * mPb, byte * yPb,
		   int n, byte * dp, int *ep));
private void
     FSDlinec4(P12(int scan, int plane_size,
		   struct error_val_field *error_values,
		   byte * cPa, byte * mPa, byte * yPa,
		   byte * cPb, byte * mPb, byte * yPb,
		   int n, byte * dp, int *ep));
private void
     init_error_buffer(struct misc_struct *misc_vars,
		       struct ptr_arrays *data_ptrs);
private void
     do_floyd_steinberg(P8(int scan, int cscan, int plane_size,
			   int plane_size_c, int n,
			   struct ptr_arrays *data_ptrs,
			   gx_device_printer * pdev,
			   struct error_val_field *error_values));
private int
    do_gcr(P7(int bytecount, byte * inbyte, const byte * kvalues,
	      const byte * cvalues, const byte * mvalues,
	      const byte * yvalues, const int *kcorrect));

/* UNUSED
 *private int
 *test_scan (P4(int size, 
 *            byte * current,
 *            byte * last, 
 *            byte * control));
 *private void
 *save_color_data(P3(int size,
 *                 byte * current,
 *                 byte * saved));
 *
 */
private void
     send_scan_lines(P6(gx_device_printer * pdev,
			struct ptr_arrays *data_ptrs,
			struct misc_struct *misc_vars,
			struct error_val_field *error_values,
			const Gamma *gamma,
			FILE * prn_stream));
private void
     do_gamma(P3(float mastergamma, float gammaval, byte * values));
private void
     do_black_correction(P2(float kvalue, int *kcorrect));

private void
     init_data_structure(P3(gx_device_printer * pdev,
			    struct ptr_arrays *data_ptrs,
			    struct misc_struct *misc_vars));
private void
     calculate_memory_size(P2(gx_device_printer * pdev,
			      struct misc_struct *misc_vars));


private void
assign_dpi(int dpi, byte * msb)
{
    if (dpi == 600) {
	msb[0] = 0x02;
	msb[1] = 0x58;
    } else {
	msb[0] = 0x01;
	msb[1] = 0x2c;
    }
}

private void
cdj850_terminate_page(gx_device_printer * pdev, FILE * prn_stream)
{
    fputs("0M", prn_stream);	/* Reset compression */
    fputs("\033*rC\033E", prn_stream);	/* End Graphics, Reset */
    fputs("\033&l0H", prn_stream);	/* eject page */
}

/* Here comes the hp850 output routine -------------------- */
private int
cdj850_print_page(gx_device_printer * pdev, FILE * prn_stream)
{

    struct error_val_field error_values;
    struct ptr_arrays data_ptrs;
    struct misc_struct misc_vars;

    Gamma gamma;

    /* make a local writable copy of the Gamma tables */
    memcpy(&gamma, gammat[cdj850->ptype], sizeof(Gamma));
    /* if mastergamma, don't use the built in functions */
    if (cdj850->mastergamma > 1.0) {
	/* prepare the bw lookup table */
	do_gamma(cdj850->mastergamma, cdj850->gammavalk, gamma.k);
	/* prepare the color lookup table */
	do_gamma(cdj850->mastergamma, cdj850->gammavalc, gamma.c);
	do_gamma(cdj850->mastergamma, cdj850->gammavalm, gamma.m);
	do_gamma(cdj850->mastergamma, cdj850->gammavaly, gamma.y);
    }
    /* prepare the black correction table for the unbunt mask */
    do_black_correction(cdj850->blackcorrect, gamma.correct);

    /* Calculate the needed memory */
    calculate_memory_size(pdev, &misc_vars);

    /* and allocate the memory */

    /* Since we need 600 and 300 dpi, we set up several buffers:
       storage contains the data as copied from gs, as well as the
       plane-data and the out_row buffer.
       storagec will contain the rescaled color data. It also contains the
       plane_data for the color-planes - these are needed by the
       compression routine, but would be overwritten by the
       b/w-dithering. The color planes allow for overwriting the
       color-data by the error-data. Since we might use the
       2bpp feature of the hp850 someday, it is sized like storage.
       storagee contains the errors from b/w fs-ditherng */

    data_ptrs.storage = (ulong *) gs_malloc(misc_vars.storage_size_words, W,
					    "cdj850_print_page");

    /* if we can't allocate working area */
    if (data_ptrs.storage == 0) {
	return_error(gs_error_VMerror);
    }
    /* Initialise the needed pointers */
    init_data_structure(pdev, &data_ptrs, &misc_vars);

    /* Start Raster mode */
    (*cdj850->start_raster_mode) (pdev, misc_vars.paper_size, prn_stream);

    /* Send each scan line in turn */
    send_scan_lines(pdev, &data_ptrs, &misc_vars,
		    &error_values, &gamma, prn_stream);

    /* terminate page and eject paper */
    (*cdj850->terminate_page) (pdev, prn_stream);

    /* Free Memory */
    gs_free((char *)data_ptrs.storage, misc_vars.storage_size_words, W,
	    "hp850_print_page");

    return 0;
}

#define odd(i) ((i & 01) != 0)

private int
GetScanLine(gx_device_printer * pdev, int *lnum,
	    struct ptr_arrays *data_ptrs,
	    struct misc_struct *misc_vars,
	    word rmask)
{
    word *data_words = (word *) data_ptrs->data[misc_vars->scan];
    register word *end_data = data_words + misc_vars->line_size_words;

    ++(*lnum);
    gdev_prn_copy_scan_lines(pdev, *lnum, (byte *) data_words, misc_vars->line_size);

    misc_vars->scan = 1 - misc_vars->scan;	/* toggle scan direction */
    misc_vars->is_two_pass = odd(*lnum);	/* color output for odd lines */

    /* Mask off 1-bits beyond the line width. */
    end_data[-1] &= rmask;

    /* Remove trailing 0s. */
    while (end_data > data_words && end_data[-1] == 0)
	end_data--;

    return end_data - data_words;
}

/* Send the scan lines to the printer */
private void
send_scan_lines(gx_device_printer * pdev,
		struct ptr_arrays *data_ptrs,
		struct misc_struct *misc_vars,
		struct error_val_field *error_values,
		const Gamma *gamma,
		FILE * prn_stream)
{
    int lnum, lend, llen;
    int num_blank_lines = 0;

    word rmask =
    ~(word) 0 << ((-pdev->width * misc_vars->storage_bpp) & (W * 8 - 1));

    lend = pdev->height - (dev_t_margin(pdev) + dev_b_margin(pdev)) * y_dpi;

    error_values->c = error_values->m = error_values->y =
	error_values->k = 0;

    /* init the error buffer */
    init_error_buffer(misc_vars, data_ptrs);

    misc_vars->zero_row_count = 0;
    lnum = -1;
    llen = GetScanLine(pdev, &lnum, data_ptrs, misc_vars, rmask);
    while (lnum < lend) {
	num_blank_lines = 0;
	while (lnum < lend && llen == 0) {
	    ++num_blank_lines;
	    llen = GetScanLine(pdev, &lnum, data_ptrs, misc_vars, rmask);
	}
	if (lnum >= lend) {
	    break;
	}
	/* Skip blank lines if any */
	if (num_blank_lines > 0) {
	    fprintf(prn_stream, "%dy", num_blank_lines / (cdj850->yscal + 1));
	    memset(data_ptrs->plane_data[0][0], 0,
		   (misc_vars->plane_size * 2 * misc_vars->num_comps));
	    memset(data_ptrs->plane_data_c[0][0], 0,
		   (misc_vars->plane_size_c * 2 * misc_vars->num_comps));

	}
	/* all blank lines printed, now for the non-blank lines */
	if (cdj850->yscal && odd(lnum)) {
	    /* output a blank black plane for odd lines */
	    putc('v', prn_stream);
	}
	/* now output all non blank lines */
	while (lnum < lend && llen != 0) {
	    misc_vars->is_color_data = 0;	/* maybe we have color ? */
	    (*cdj850->print_non_blank_lines) (pdev, data_ptrs, misc_vars,
					      error_values, gamma, prn_stream);
	    llen = GetScanLine(pdev, &lnum, data_ptrs, misc_vars, rmask);
	}
	if (cdj850->yscal && odd(lnum)) {	/* output empty line for odd lines */
	    (*cdj850->print_non_blank_lines) (pdev, data_ptrs, misc_vars,
					      error_values, gamma, prn_stream);
	}
	/* the current line is empty => run the next iteration */
    }
}

/* print_line compresses (mode 9) and outputs one plane */
private void
print_c9plane(FILE * prn_stream, char plane_code, int plane_size,
	      const byte * curr, const byte * prev, byte * out_data)
{
    /* Compress the output data */
    int out_count = gdev_pcl_mode9compress(plane_size, curr, prev, out_data);

    /* and output the data */
    if (out_count > 0) {
	fprintf(prn_stream, "%d%c", out_count, plane_code);
	fwrite(out_data, sizeof(byte), out_count, prn_stream);
    } else {
	putc(plane_code, prn_stream);
    }
}

/* Printing non-blank lines */
private void
cdj850_print_non_blank_lines(gx_device_printer * pdev,
			     struct ptr_arrays *data_ptrs,
			     struct misc_struct *misc_vars,
			     struct error_val_field *error_values,
			     const Gamma *gamma,
			     FILE * prn_stream)
{
    static const char *const plane_code[2] =
    {"wvvv", "vvvv"};

    int i;
    byte *kP = data_ptrs->plane_data[misc_vars->scan + 2][3];
    byte *dp = data_ptrs->data[misc_vars->scan + 2];
    int *ep = data_ptrs->errors[misc_vars->scan];

    /* we need cmyk color separation befor all the rest, since
       black may be contained in the color fields. This needs to
       be done on all pixel-rows, since even unused color-bytes
       might generate black */

    misc_vars->is_color_data =
	do_gcr(misc_vars->databuff_size, data_ptrs->data[misc_vars->scan],
	       gamma->k, gamma->c, gamma->m, gamma->y, gamma->correct);

    /* dithering the black-plane */
    FSDlinebw(misc_vars->scan, misc_vars->plane_size,
	      error_values, kP, misc_vars->num_comps, ep, dp);

    /* output the black plane */
    print_c9plane(prn_stream, 'v', misc_vars->plane_size,
		  data_ptrs->plane_data[misc_vars->scan][3],
		  data_ptrs->plane_data[1 - misc_vars->scan][3],
		  data_ptrs->out_data);

    /* since color resolution is only half of the b/w-resolution,
       we only output every second row */
    if (!cdj850->yscal || misc_vars->is_two_pass) {

	int plane_size_c = (*rescale_color_plane[cdj850->xscal][cdj850->yscal])
	(misc_vars->databuff_size,
	 data_ptrs->data[misc_vars->scan],
	 data_ptrs->data[1 - misc_vars->scan],
	 data_ptrs->data_c[misc_vars->cscan]) / misc_vars->storage_bpp;

	/* dither the color planes */
	do_floyd_steinberg(misc_vars->scan, misc_vars->cscan,
			   misc_vars->plane_size, plane_size_c,
			   misc_vars->num_comps, data_ptrs, pdev, error_values);

	/* Transfer raster graphics in the order C, M, Y, that is
	   planes 2,1,0 */
	for (i = misc_vars->num_comps - 2; i >= 0; i--) {

	    /* output the lower color planes */
	    print_c9plane(prn_stream, plane_code[cdj850->intensities > 2][i],
			  plane_size_c,
			  data_ptrs->plane_data_c[misc_vars->cscan][i],
			  data_ptrs->plane_data_c[1 - misc_vars->cscan][i],
			  data_ptrs->out_data);

	    /* output the upper color planes */
	    if (cdj850->intensities > 2) {
		print_c9plane(prn_stream, plane_code[0][i], plane_size_c,
			      data_ptrs->plane_data_c[misc_vars->cscan][i + 4],
			      data_ptrs->plane_data_c[1 -
									    misc_vars->cscan][i
									    + 4],
			      data_ptrs->out_data);
	    }			/* end cdj850->intensities > 2 */
	}			/* End For i = num_comps */
	misc_vars->cscan = 1 - misc_vars->cscan;
    }				/* End of is_two_pass */
    return;
}

/* moved that code into his own subroutine, otherwise things get
   somewhat clumsy */
private void
do_floyd_steinberg(int scan, int cscan, int plane_size,
		   int plane_size_c, int n,
		   struct ptr_arrays *data_ptrs,
		   gx_device_printer * pdev,
		   struct error_val_field *error_values)
{
    /* the color pointers */
    byte *cPa, *mPa, *yPa, *cPb, *mPb, *yPb;
    byte *dpc;
    int *epc;

    /* the b/w pointers */
    byte *kP, *dp;
    int *ep;

    /* the color pointers, lower byte */
    cPa = data_ptrs->plane_data_c[cscan + 2][2];
    mPa = data_ptrs->plane_data_c[cscan + 2][1];
    yPa = data_ptrs->plane_data_c[cscan + 2][0];
    /* upper byte */
    cPb = data_ptrs->plane_data_c[cscan + 2][6];
    mPb = data_ptrs->plane_data_c[cscan + 2][5];
    yPb = data_ptrs->plane_data_c[cscan + 2][4];
    /* data and error */
    dpc = data_ptrs->data_c[cscan + 2];
    epc = data_ptrs->errors_c[cscan];
    /* the b/w pointers */
    kP = data_ptrs->plane_data[scan + 2][3];
    dp = data_ptrs->data[scan + 2];
    ep = data_ptrs->errors[scan];

    switch (cdj850->intensities) {
	case 2:
	FSDlinec2(cscan, plane_size_c, error_values,
		  cPa, mPa, yPa, n, dpc, epc);
	break;
	case 3:
	FSDlinec3(cscan, plane_size_c, error_values,
		  cPa, mPa, yPa, cPb, mPb, yPb, n, dpc, epc);
	break;
	case 4:
	FSDlinec4(cscan, plane_size_c, error_values,
		  cPa, mPa, yPa, cPb, mPb, yPb, n, dpc, epc);
	break;
	default:
	assert(0);
    }
    return;
}

/* here we do our own gamma-correction */
private void
do_gamma(float mastergamma, float gammaval, byte values[256])
{
    int i;
    float gamma;

    if (gammaval > 0.0) {
	gamma = gammaval;
    } else {
	gamma = mastergamma;
    }

    for (i = 0; i < 256; i++) {
	values[i] = (byte) (255.0 *
			    (1.0 - pow(((double)(255.0 - (float)i) / 255.0),
			     (double)(1.0 / gamma))));
    }

    return;
}

/* here we calculate a lookup-table which is used to compensate the
   relativ loss of color due to undercolor-removal */
private void
do_black_correction(float kvalue, int kcorrect[256])
{
    int i;

    for (i = 0; i < 256; i++) {
	kcorrect[i] = (int)
	    (100.0 * kvalue * (
				  pow(10.0,
				      pow((i / 255.0), 3.0)
				  )
				  - 1.0
	     )
	    );
    }

    return;
}

/* For Better Performance we use a macro here */
#define DOUCR(col1, col2, col3, col4)\
{\
  /* determine how far we are from the grey axis. This is  */\
  /* traditionally done by computing MAX(CMY)-MIN(CMY).    */\
  /* However, if two colors are very similar, we could     */\
  /* as either CMYRGB and K. Therefore we calculate the    */\
  /* the distance col1-col2 and col2-col3, and use the     */\
  /* smaller one.                                          */\
  a = *col1 - *col2;\
  b = *col2 - *col3;\
  if (a >= b) {\
    grey_distance = 1.0 - (b/255.0);\
  } else {\
    grey_distance = 1.0 - (a/255.0);\
  }\
  ucr   = (byte) (*col3 * grey_distance); \
  *col4 = *col4 + ucr;  /* add removed black to black */\
  /* remove only as much color as black is surviving the   */\
  /* gamma correction */\
  ucr   = *(kvalues + ucr);\
  *col1 = *col1 - ucr ;\
  *col2 = *col2 - ucr ;\
  *col3 = *col3 - ucr ;\
}

/* For Better Performance we use a macro here */
#define DOGCR(col1, col2, col3, col4)\
{\
  ucr = (int) *col3;\
  *col1 -= ucr ;\
  *col2 -= ucr ;\
  *col3 -= ucr ;\
  *col4 += ucr;  /* add removed black to black */\
  kadd  = ucr + *(kcorrect + ucr);\
  uca_fac = 1.0 + (kadd/255.0);\
  *col1 *= uca_fac;\
  *col2 *= uca_fac;\
}

/* Since resolution can be different on different planes, we need to
   do real color separation, here we try a real grey component
   replacement */
private int
do_gcr(int bytecount, byte * inbyte, const byte kvalues[256],
       const byte cvalues[256], const byte mvalues[256],
       const byte yvalues[256], const int kcorrect[256])
{
    int i, ucr, kadd, is_color = 0;
    byte *black, *cyan, *magenta, *yellow;
    float uca_fac;

    /* Grey component replacement */
    for (i = 0; i < bytecount; i += 4) {
	black = inbyte++;	/* Assign to black the current address of  inbyte */
	cyan = inbyte++;
	magenta = inbyte++;
	yellow = inbyte++;

	if (*magenta + *yellow + *cyan > 0) {	/* if any color at all */

	    is_color = 1;

	    if ((*cyan >= *magenta)
		&& (*magenta >= *yellow)
		&& (*yellow > 0)) {	/* if any grey component */
		DOGCR(cyan, magenta, yellow, black);
	    } else if ((*cyan >= *yellow)
		       && (*yellow >= *magenta)
		       && (*magenta > 0)) {
		DOGCR(cyan, yellow, magenta, black);
	    } else if ((*yellow >= *magenta)
		       && (*magenta >= *cyan)
		       && (*cyan > 0)) {
		DOGCR(yellow, magenta, cyan, black);
	    } else if ((*yellow >= *cyan)
		       && (*cyan >= *magenta)
		       && (*magenta > 0)) {
		DOGCR(yellow, cyan, magenta, black);
	    } else if ((*magenta >= *yellow)
		       && (*yellow >= *cyan)
		       && (*cyan > 0)) {
		DOGCR(magenta, yellow, cyan, black);
	    } else if ((*magenta >= *cyan)
		       && (*cyan >= *yellow)
		       && (*yellow > 0)) {
		DOGCR(magenta, cyan, yellow, black);
	    } else {		/* do gamma only if no black */
	    }
	    *cyan = *(cvalues + *cyan);
	    *magenta = *(mvalues + *magenta);
	    *yellow = *(yvalues + *yellow);
	}			/* end of if c+m+y > 0 */
	*black = *(kvalues + *black);
    }				/* end of for bytecount */
    return is_color;
}

/* Since resolution can be different on different planes, we need to
   rescale the data byte by byte */
private int
rescale_byte_wise2x2(int bytecount, const byte * inbytea, const byte * inbyteb,
		     byte * outbyte)
{
    register int i, j;
    int max = bytecount / 2;

    for (i = 0; i < max; i += 4) {
	j = 2 * i;
	/* cyan */
	outbyte[i + 1] = (inbytea[j + 1] + inbytea[j + 5] + inbyteb[j + 1] +
			  inbyteb[j + 5]) / 4;
	/* magenta */
	outbyte[i + 2] = (inbytea[j + 2] + inbytea[j + 6] + inbyteb[j + 2] +
			  inbyteb[j + 6]) / 4;
	/* yellow */
	outbyte[i + 3] = (inbytea[j + 3] + inbytea[j + 7] + inbyteb[j + 3] +
			  inbyteb[j + 7]) / 4;
    }
    return max;
}

/* Since resolution can be different on different planes, we need to
   rescale the data byte by byte */
private int
rescale_byte_wise2x1(int bytecount, const byte * inbytea, const byte * inbyteb,
		     byte * outbyte)
{
    register int i, j;
    int max = bytecount / 2;

    for (i = 0; i < max; i += 4) {
	j = 2 * i;
	/* cyan */
	outbyte[i + 1] = (inbytea[j + 1] + inbytea[j + 5]) / 2;
	/* magenta */
	outbyte[i + 2] = (inbytea[j + 2] + inbytea[j + 6]) / 2;
	/* yellow */
	outbyte[i + 3] = (inbytea[j + 3] + inbytea[j + 7]) / 2;
    }
    return max;
}

/* Since resolution can be different on different planes, we need to
   rescale the data byte by byte */
private int
rescale_byte_wise1x2(int bytecount, const byte * inbytea, const byte * inbyteb,
		     byte * outbyte)
{
    register int i;

    for (i = 0; i < bytecount; i += 4) {
	/* cyan */
	outbyte[i + 1] = (inbytea[i + 1] + inbyteb[i + 1]) / 2;
	/* magenta */
	outbyte[i + 2] = (inbytea[i + 2] + inbyteb[i + 2]) / 2;
	/* yellow */
	outbyte[i + 3] = (inbytea[i + 3] + inbyteb[i + 3]) / 2;
    }
    return bytecount;
}

/* Since resolution can be different on different planes, we need to
   rescale the data byte by byte */
private int
rescale_byte_wise1x1(int bytecount, const byte * inbytea, const byte * inbyteb,
		     byte * outbyte)
{
    register int i;

    for (i = 0; i < bytecount; i += 4) {
	/* cyan */
	outbyte[i + 1] = inbytea[i + 1];
	/* magenta */
	outbyte[i + 2] = inbytea[i + 2];
	/* yellow */
	outbyte[i + 3] = inbytea[i + 3];
    }
    return bytecount;
}

/* MACROS FOR DITHERING (we use macros for compact source and faster code) */
/* Floyd-Steinberg dithering. Often results in a dramatic improvement in
 * subjective image quality, but can also produce dramatic increases in
 * amount of printer data generated and actual printing time!! Mode 9 2D
 * compression is still useful for fairly flat colour or blank areas but its
 * compression is much less effective in areas where the dithering has
 * effectively randomised the dot distribution. */

#define RSHIFT ((I * 8) - 16)
#define SHIFT ((I * 8) - 13)
#define MAXVALUE  (255 << SHIFT)
#define RANDOM (((rand() << RSHIFT) % (MAXVALUE / 2))  - MAXVALUE /4);
#define MINVALUE  0
#define C 8

#define THRESHOLD (128 << SHIFT)

/* --- needed for the hp850 -- */
#define SHIFTS ((I * 8) - 14)
#define SHIFTM ((I * 8) - 13)
#define SHIFTL ((I * 8) - 12)

#define MAXVALUES  (160 << SHIFTM)
#define MAXVALUEM  (226 << SHIFTM)
#define MAXVALUEL  (255 << SHIFTM)

#define THRESHOLDS (128 << SHIFTM)
#define THRESHOLDM (192 << SHIFTM)
#define THRESHOLDL (226 << SHIFTM)
/* --------------------------- */

/* initialise the error_buffer */
private void
init_error_buffer(struct misc_struct *misc_vars,
		  struct ptr_arrays *data_ptrs)
{
    int i;
    int *ep;
    int *epc;

    ep = data_ptrs->errors[0];
    epc = data_ptrs->errors_c[0];

    if (misc_vars->bits_per_pixel > 4) {	/* Randomly seed initial error
						   buffer */
	/* Otherwise, the first dithered rows would look rather uniform */
	for (i = 0; i < misc_vars->databuff_size; i++) {	/* 600dpi planes */
	    *ep++ = RANDOM;
	}

	/* Now for the 2 * 300dpi color planes */
	for (i = 0; i < misc_vars->databuff_size_c; i++) {
	    *epc++ = RANDOM;
	}
    }
    return;
}

#define FSdither(inP, out, errP, Err, Bit, Offset, Element)\
{\
   oldErr = Err;\
	Err = (*(errP + Element)\
	       + ((Err * 7 + C) >> 4)\
	       + ((int)*(inP + Element) << SHIFT));\
	if (Err > THRESHOLD) {\
	  out |= Bit;\
	  Err -= MAXVALUE;\
        }\
	*(errP + (Element + Offset)) += ((Err * 3 + C) >> 4);\
	*(errP + Element) = ((Err * 5 + oldErr + C) >> 4);\
}

/* The hp850c has 600dpi black and 300 dpi color. Therefore, we need
   an adapted dither algorythm */
private void
FSDlinebw(int scan, int plane_size,
	  struct error_val_field *error_values,
	  byte * kP, int n, int *ep, byte * dp)
{
    if (scan == 0) {		/* going_up */
	byte k, bitmask;	/* k = outbyte byte, whereas bitmask defines the

				   bit to be set within k */
	int oldErr, i;

	for (i = 0; i < plane_size; i++) {
	    bitmask = 0x80;
	    for (k = 0; bitmask != 0; bitmask >>= 1) {
		/* dp points to the first word of the input data which is in
		   kcmy-format */
		/* k points to the beginning of the first outbut byte, which
		   is filled up, bit by bit while looping over bytemask */
		/* ep points to the first word of the error-plane which
		   contains the errors kcmy format */
		/* err_values->k tempararily holds the error-value */
		/* bitmask selects the bit to be set in the outbyte */
		/* n gives the offset for the byte selection within
		   words. With simple cmyk-printing, this should be 4 */
		/* 0 points to the active color within the input-word, i.e. 0
		   = black, 1 = cyan, 2 = yellow, 3 = magenta */

		FSdither(dp, k, ep, error_values->k, bitmask, -n, 0);
		dp += n, ep += n;	/* increment the input and error pointer one
					   word (=4 byte) further, in order to
					   convert the next word into an bit */
	    }
	    *kP++ = k;		/* fill the output-plane byte with the computet byte
				   and increment the output plane pointer  one byte */
	}

    } else {			/* going_down */
	byte k, bitmask;
	int oldErr, i;

	for (i = 0; i < plane_size; i++) {
	    bitmask = 0x01;
	    for (k = 0; bitmask != 0; bitmask <<= 1) {
		dp -= n, ep -= n;
		FSdither(dp, k, ep, error_values->k, bitmask, n, 0);
	    }
	    *--kP = k;
	}
    }
    return;
}

/* Since bw has already been dithered for the hp850c, we need
   an adapted dither algorythm */
private void
FSDlinec2(int scan, int plane_size,
	  struct error_val_field *error_values,
	  byte * cPa, byte * mPa, byte * yPa, int n,
	  byte * dp, int *ep)
{
    if (scan == 0) {		/* going_up */
	int oldErr, i;
	byte ca, ya, ma, bitmask;

	for (i = 0; i < plane_size; i++) {
	    bitmask = 0x80;
	    ca = ya = ma = 0;
	    for (ca = 0; bitmask != 0; bitmask >>= 1) {
		FSdither(dp, ca, ep, error_values->c, bitmask, -n, n - 3);
		FSdither(dp, ma, ep, error_values->m, bitmask, -n, n - 2);
		FSdither(dp, ya, ep, error_values->y, bitmask, -n, n - 1);
		dp += n, ep += n;
	    }
	    *cPa++ = ca;
	    *mPa++ = ma;
	    *yPa++ = ya;
	}

    } else {			/* going_down */
	byte ca, ya, ma, bitmask;
	int oldErr, i;

	for (i = 0; i < plane_size; i++) {
	    bitmask = 0x01;
	    ca = ya = ma = 0;
	    for (ca = 0; bitmask != 0; bitmask <<= 1) {
		dp -= n, ep -= n;
		FSdither(dp, ya, ep, error_values->y, bitmask, n, n - 1);
		FSdither(dp, ma, ep, error_values->m, bitmask, n, n - 2);
		FSdither(dp, ca, ep, error_values->c, bitmask, n, n - 3);
	    }
	    *--yPa = ya;
	    *--mPa = ma;
	    *--cPa = ca;
	}
    }
    return;
}

/* while printing on paper, we only use 3 -intensities */
#define FSdither8503(inP, outa, outb, errP, Err, Bit, Offset, Element)\
{\
	oldErr = Err;\
	Err = (*(errP + Element)\
	       + ((Err * 7 + C) >> 4)\
	       + ((int) *(inP + Element) << SHIFT));\
	if ((Err > THRESHOLDS) && (Err <= THRESHOLDM)) {\
	  outa |= Bit;\
	  Err -= MAXVALUES;\
	}\
	if (Err > THRESHOLDM) {\
	  outb |= Bit;\
	  Err -= MAXVALUEM;\
	}\
	*(errP + (Element + Offset)) += ((Err * 3 + C) >> 4);\
	*(errP + Element) = ((Err * 5 + oldErr + C) >> 4);\
}

/* On ordinary paper, we'll only use 3 intensities with the hp850  */
private void
FSDlinec3(int scan, int plane_size,
	  struct error_val_field *error_values,
	  byte * cPa, byte * mPa, byte * yPa,
	  byte * cPb, byte * mPb, byte * yPb,
	  int n, byte * dp, int *ep)
{
    if (scan == 0) {		/* going_up */
	byte ca, ya, ma, cb, yb, mb, bitmask;
	int oldErr, i;

	for (i = 0; i < plane_size; i++) {
	    bitmask = 0x80;
	    ca = ya = ma = cb = yb = mb = 0;
	    for (ca = 0; bitmask != 0; bitmask >>= 1) {
		FSdither8503(dp, ca, cb, ep, error_values->c, bitmask, -n, n
			     - 3);
		FSdither8503(dp, ma, mb, ep, error_values->m, bitmask, -n, n
			     - 2);
		FSdither8503(dp, ya, yb, ep, error_values->y, bitmask, -n, n
			     - 1);
		dp += n, ep += n;
	    }
	    *cPa++ = ca;
	    *mPa++ = ma;
	    *yPa++ = ya;
	    *cPb++ = cb;
	    *mPb++ = mb;
	    *yPb++ = yb;
	}
    } else {			/* going_down */
	byte ca, ya, ma, cb, yb, mb, bitmask;
	int oldErr, i;

	for (i = 0; i < plane_size; i++) {
	    bitmask = 0x01;
	    ca = ya = ma = cb = yb = mb = 0;
	    for (ca = 0; bitmask != 0; bitmask <<= 1) {
		dp -= n, ep -= n;
		FSdither8503(dp, ya, yb, ep, error_values->y, bitmask, n, n
			     - 1);
		FSdither8503(dp, ma, mb, ep, error_values->m, bitmask, n, n
			     - 2);
		FSdither8503(dp, ca, cb, ep, error_values->c, bitmask, n, n
			     - 3);
	    }
	    *--yPa = ya;
	    *--mPa = ma;
	    *--cPa = ca;
	    *--yPb = yb;
	    *--mPb = mb;
	    *--cPb = cb;
	}
    }
    return;
}


/* the hp850 knows about 4 different color intensities per color */
#define FSdither8504(inP, outa, outb, errP, Err, Bit, Offset, Element)\
{\
	oldErr = Err;\
	Err = (*(errP + Element)\
	       + ((Err * 7 + C) >> 4)\
	       + ((int) *(inP + Element) << SHIFT));\
	if ((Err > THRESHOLDS) && (Err <= THRESHOLDM)) {\
	  outa |= Bit;\
	  Err -= MAXVALUES;\
	}\
	if ((Err > THRESHOLDM) && (Err <= THRESHOLDL)) {\
	  outb |= Bit;\
	  Err -= MAXVALUEM;\
	}\
	if (Err > THRESHOLDL) {\
          outa |= Bit;\
	  outb |= Bit;\
	  Err -= MAXVALUEL;\
	}\
	*(errP + (Element + Offset)) += ((Err * 3 + C) >> 4);\
	*(errP + Element) = ((Err * 5 + oldErr + C) >> 4);\
}

/* The hp850c knows about 4 intensity levels per color. Once more, we need
   an adapted dither algorythm */
private void
FSDlinec4(int scan, int plane_size,
	  struct error_val_field *error_values,
	  byte * cPa, byte * mPa, byte * yPa,
	  byte * cPb, byte * mPb, byte * yPb,
	  int n, byte * dp, int *ep)
{
    if (scan == 0) {		/* going_up */
	byte ca, ya, ma, cb, yb, mb, bitmask;
	int oldErr, i;

	for (i = 0; i < plane_size; i++) {
	    bitmask = 0x80;
	    ca = ya = ma = cb = yb = mb = 0;
	    for (ca = 0; bitmask != 0; bitmask >>= 1) {
		FSdither8504(dp, ca, cb, ep, error_values->c, bitmask, -n, n
			     - 3);
		FSdither8504(dp, ma, mb, ep, error_values->m, bitmask, -n, n
			     - 2);
		FSdither8504(dp, ya, yb, ep, error_values->y, bitmask, -n, n
			     - 1);
		dp += n, ep += n;
	    }
	    *cPa++ = ca;
	    *mPa++ = ma;
	    *yPa++ = ya;
	    *cPb++ = cb;
	    *mPb++ = mb;
	    *yPb++ = yb;
	}
    } else {			/* going_down */
	byte ca, ya, ma, cb, yb, mb, bitmask;
	int oldErr, i;

	for (i = 0; i < plane_size; i++) {
	    bitmask = 0x01;
	    ca = ya = ma = cb = yb = mb = 0;
	    for (ca = 0; bitmask != 0; bitmask <<= 1) {
		dp -= n, ep -= n;
		FSdither8504(dp, ya, yb, ep, error_values->y, bitmask, n, n
			     - 1);
		FSdither8504(dp, ma, mb, ep, error_values->m, bitmask, n, n
			     - 2);
		FSdither8504(dp, ca, cb, ep, error_values->c, bitmask, n, n
			     - 3);
	    }
	    *--yPa = ya;
	    *--mPa = ma;
	    *--cPa = ca;
	    *--yPb = yb;
	    *--mPb = mb;
	    *--cPb = cb;
	}
    }
    return;
}


/* calculate the needed memory */
private void
calculate_memory_size(gx_device_printer * pdev,
		      struct misc_struct *misc_vars)
{
    int xfac = cdj850->xscal ? 2 : 1;

    misc_vars->line_size = gdev_prn_raster(pdev);
    misc_vars->line_size_c = misc_vars->line_size / xfac;
    misc_vars->line_size_words = (misc_vars->line_size + W - 1) / W;
    misc_vars->paper_size = gdev_pcl_paper_size((gx_device *) pdev);
    misc_vars->num_comps = pdev->color_info.num_components;
    misc_vars->bits_per_pixel = pdev->color_info.depth;
    misc_vars->storage_bpp = misc_vars->num_comps * 8;
    misc_vars->expanded_bpp = misc_vars->num_comps * 8;
    misc_vars->errbuff_size = 0;
    misc_vars->errbuff_size_c = 0;

    misc_vars->plane_size = calc_buffsize(misc_vars->line_size, misc_vars->storage_bpp);

    /* plane_size_c is dependedend on the bits used for
       dithering. Currently 2 bits are sufficient  */
    misc_vars->plane_size_c = 2 * misc_vars->plane_size / xfac;

    /* 4n extra values for line ends */
    /* might be wrong, see gdevcdj.c */
    misc_vars->errbuff_size =
	calc_buffsize((misc_vars->plane_size * misc_vars->expanded_bpp +
		       misc_vars->num_comps * 4) * I, 1);

    /* 4n extra values for line ends */
    misc_vars->errbuff_size_c =
	calc_buffsize((misc_vars->plane_size_c / 2 * misc_vars->expanded_bpp
		       + misc_vars->num_comps * 4) * I, 1);

    misc_vars->databuff_size =
	misc_vars->plane_size * misc_vars->storage_bpp;

    misc_vars->databuff_size_c =
	misc_vars->plane_size_c / 2 * misc_vars->storage_bpp;


    misc_vars->outbuff_size = misc_vars->plane_size * 4;

    misc_vars->storage_size_words = (((misc_vars->plane_size)
				      * 2
				      * misc_vars->num_comps)
				     + misc_vars->databuff_size
				     + misc_vars->errbuff_size
				     + misc_vars->outbuff_size
				     + ((misc_vars->plane_size_c)
					* 2
					* misc_vars->num_comps)
				     + misc_vars->databuff_size_c
				     + misc_vars->errbuff_size_c
				     + (4 * misc_vars->plane_size_c))
	/ W;

    return;
}


/* Initialise the needed pointers */
private void
init_data_structure(gx_device_printer * pdev,
		    struct ptr_arrays *data_ptrs,
		    struct misc_struct *misc_vars)
{
    int i;
    byte *p = (byte *) data_ptrs->storage;

    misc_vars->scan = 0;
    misc_vars->cscan = 0;
    misc_vars->is_two_pass = 1;

    /* the b/w pointer */
    data_ptrs->data[0] = data_ptrs->data[1] = data_ptrs->data[2] = p;
    data_ptrs->data[3] = p + misc_vars->databuff_size;
    /* Note: The output data will overwrite part of the input-data */

    if (misc_vars->bits_per_pixel > 1) {
	p += misc_vars->databuff_size;
    }
    if (misc_vars->bits_per_pixel > 4) {
	data_ptrs->errors[0] = (int *)p + misc_vars->num_comps * 2;
	data_ptrs->errors[1] = data_ptrs->errors[0] + misc_vars->databuff_size;
	p += misc_vars->errbuff_size;
    }
    for (i = 0; i < misc_vars->num_comps; i++) {
	data_ptrs->plane_data[0][i] = data_ptrs->plane_data[2][i] = p;
	p += misc_vars->plane_size;
    }
    for (i = 0; i < misc_vars->num_comps; i++) {
	data_ptrs->plane_data[1][i] = p;
	data_ptrs->plane_data[3][i] = p + misc_vars->plane_size;
	p += misc_vars->plane_size;
    }
    data_ptrs->out_data = p;
    p += misc_vars->outbuff_size;

    /* ---------------------------------------------------------
       now for the color pointers 
       --------------------------------------------------------- */

    data_ptrs->data_c[0] = data_ptrs->data_c[1] = data_ptrs->data_c[2] = p;
    data_ptrs->data_c[3] = p + misc_vars->databuff_size_c;
    /* Note: The output data will overwrite part of the input-data */

    if (misc_vars->bits_per_pixel > 1) {
	p += misc_vars->databuff_size_c;
    }
    if (misc_vars->bits_per_pixel > 4) {
	data_ptrs->errors_c[0] = (int *)p + misc_vars->num_comps * 2;
	data_ptrs->errors_c[1] = data_ptrs->errors_c[0] + misc_vars->databuff_size_c;
	p += misc_vars->errbuff_size_c;
    }
    /* pointer for the lower bits of the output data */
    for (i = 0; i < misc_vars->num_comps; i++) {
	data_ptrs->plane_data_c[0][i] = data_ptrs->plane_data_c[2][i] = p;
	p += misc_vars->plane_size_c / 2;
    }
    for (i = 0; i < misc_vars->num_comps; i++) {
	data_ptrs->plane_data_c[1][i] = p;
	data_ptrs->plane_data_c[3][i] = p + misc_vars->plane_size_c / 2;
	p += misc_vars->plane_size_c / 2;
    }

    /* pointer for the upper bits of the output data */
    for (i = 0; i < misc_vars->num_comps; i++) {
	data_ptrs->plane_data_c[0][i + 4] = data_ptrs->plane_data_c[2][i +
	    4] = p;
	p += misc_vars->plane_size_c / 2;
    }
    for (i = 0; i < misc_vars->num_comps; i++) {
	data_ptrs->plane_data_c[1][i + 4] = p;
	data_ptrs->plane_data_c[3][i + 4] = p + misc_vars->plane_size_c / 2;
	p += misc_vars->plane_size_c / 2;
    }

    for (i = 0; i < misc_vars->num_comps; i++) {
	data_ptrs->test_data[i] = p;
	p += misc_vars->plane_size_c / 2;
    }

    /* Clear temp storage */
    memset(data_ptrs->storage, 0, misc_vars->storage_size_words * W);

    return;
}				/* end init_data_structure */

/* Configure the printer and start Raster mode */
private void
cdj850_start_raster_mode(gx_device_printer * pdev, int paper_size,
			 FILE * prn_stream)
{
    int xres, yres;		/* x,y resolution for color planes */
    hp850_cmyk_init_t init;

    init = hp850_cmyk_init;
    init.a[13] = cdj850->intensities;	/* Intensity levels cyan */
    init.a[19] = cdj850->intensities;	/* Intensity levels magenta */
    init.a[25] = cdj850->intensities;	/* Intensity levels yellow */

    /* black plane resolution */
    assign_dpi(cdj850->x_pixels_per_inch, init.a + 2);
    assign_dpi(cdj850->y_pixels_per_inch, init.a + 4);
    /* color plane resolution */
    xres = cdj850->x_pixels_per_inch / (cdj850->xscal + 1);
    yres = cdj850->y_pixels_per_inch / (cdj850->yscal + 1);
    /* cyan */
    assign_dpi(xres, init.a + 8);
    assign_dpi(yres, init.a + 10);
    /* magenta */
    assign_dpi(xres, init.a + 14);
    assign_dpi(yres, init.a + 16);
    /* yellow */
    assign_dpi(xres, init.a + 20);
    assign_dpi(yres, init.a + 22);

    fputs("\033*rbC", prn_stream);	/* End raster graphics */
    fputs("\033E", prn_stream);	/* Reset */
    /* Page size, orientation, top margin & perforation skip */
    fprintf(prn_stream, "\033&l%daolE", paper_size);

    /* Print Quality, -1 = draft, 0 = normal, 1 = presentation */
    fprintf(prn_stream, "\033*o%dM", cdj850->quality);
    /* Media Type,0 = plain paper, 1 = bond paper, 2 = special
       paper, 3 = glossy film, 4 = transparency film */
    fprintf(prn_stream, "\033&l%dM", cdj850->papertype);

    /* Move to top left of printed area */
    fprintf(prn_stream, "\033*p%dY", (int)(600 * DOFFSET));

    /* This will start and configure the raster-mode */
    fprintf(prn_stream, "\033*g%dW", (int)sizeof(init.a));	/* The new configure
									   raster data comand */
    fwrite(init.a, sizeof(byte), sizeof(init.a),
	   prn_stream);		/* Transmit config
				   data */
    /* From now on, all escape commands start with \033*b, so we
     * combine them (if the printer supports this). */
    fputs("\033*b", prn_stream);
    /* Set compression if the mode has been defined. */
    if (cdj850->compression)
	fprintf(prn_stream, "%dm", cdj850->compression);

    return;
}				/* end configure raster-mode */

private int
cdj_put_param_int(gs_param_list * plist, gs_param_name pname, int *pvalue,
		  int minval, int maxval, int ecode)
{
    int code, value;

    switch (code = param_read_int(plist, pname, &value)) {
	default:
	return code;
	case 1:
	return ecode;
	case 0:
	if (value < minval || value > maxval)
	    param_signal_error(plist, pname, gs_error_rangecheck);
	*pvalue = value;
	return (ecode < 0 ? ecode : 1);
    }
}

private int
cdj_put_param_float(gs_param_list * plist, gs_param_name pname, float *pvalue,
		    float minval, float maxval, int ecode)
{
    int code;
    float value;

    switch (code = param_read_float(plist, pname, &value)) {
	default:
	return code;
	case 1:
	return ecode;
	case 0:
	if (value < minval || value > maxval)
	    param_signal_error(plist, pname, gs_error_rangecheck);
	*pvalue = value;
	return (ecode < 0 ? ecode : 1);
    }
}

private int
cdj_set_bpp(gx_device * pdev, int bpp, int ccomps)
{
    gx_device_color_info *ci = &pdev->color_info;

    if (ccomps && bpp == 0) {
	if (cprn_device->cmyk) {
	    switch (ccomps) {
		default:
		return gs_error_rangecheck;
		/*NOTREACHED */
		break;

		case 1:
		bpp = 1;
		break;

		case 3:
		bpp = 24;
		break;

		case 4:
		switch (ci->depth) {
		    case 8:
		    case 16:
		    case 24:
		    case 32:
		    break;

		    default:
		    bpp = cprn_device->default_depth;
		    break;
		}
		break;
	    }
	}
    }
    if (bpp == 0) {
	bpp = ci->depth;	/* Use the current setting. */
    }
    if (cprn_device->cmyk < 0) {

	/* Reset procedures because we may have been in another mode. */

	dev_proc(pdev, map_cmyk_color) = gdev_cmyk_map_cmyk_color;
	dev_proc(pdev, map_rgb_color) = NULL;
	dev_proc(pdev, map_color_rgb) = gdev_cmyk_map_color_rgb;

	if (pdev->is_open)
	    gs_closedevice(pdev);
    }
    /* Check for valid bpp values */

    switch (bpp) {
	case 16:
	case 32:
	if (cprn_device->cmyk && ccomps && ccomps != 4)
	    goto bppe;
	break;

	case 24:
	if (!cprn_device->cmyk || ccomps == 0 || ccomps == 4) {
	    break;
	} else if (ccomps == 1) {
	    goto bppe;
	} else {

	    /* 3 components 24 bpp printing for CMYK device. */

	    cprn_device->cmyk = -1;
	}
	break;

	case 8:
	if (cprn_device->cmyk) {
	    if (ccomps) {
		if (ccomps == 3) {
		    cprn_device->cmyk = -1;
		    bpp = 3;
		} else if (ccomps != 1 && ccomps != 4) {
		    goto bppe;
		}
	    }
	    if (ccomps != 1)
		break;
	} else {
	    break;
	}

	case 1:
	if (ccomps != 1)
	    goto bppe;

	if (cprn_device->cmyk && bpp != pdev->color_info.depth) {
	    dev_proc(pdev, map_cmyk_color) = NULL;
	    dev_proc(pdev, map_rgb_color) = gdev_cmyk_map_rgb_color;

	    if (pdev->is_open) {
		gs_closedevice(pdev);
	    }
	}
	break;

	case 3:
	if (!cprn_device->cmyk) {
	    break;
	}
	default:
      bppe:return gs_error_rangecheck;
    }


    if (cprn_device->cmyk == -1) {
	dev_proc(pdev, map_cmyk_color) = NULL;
	dev_proc(pdev, map_rgb_color) = gdev_pcl_map_rgb_color;
	dev_proc(pdev, map_color_rgb) = gdev_pcl_map_color_rgb;

	if (pdev->is_open) {
	    gs_closedevice(pdev);
	}
    }
    switch (ccomps) {
	case 0:
	break;

	case 1:
	if (bpp != 1 && bpp != 8)
	    goto cce;
	break;

	case 4:
	if (cprn_device->cmyk) {
	    if (bpp >= 8)
		break;
	}
	case 3:
	if (bpp == 1 || bpp == 3 || bpp == 8 || bpp == 16
	    || bpp == 24 || bpp == 32) {
	    break;
	}
	cce: default:
	return gs_error_rangecheck;
    }

    if (cprn_device->cmyk) {
	if (cprn_device->cmyk > 0) {
	    ci->num_components = ccomps ? ccomps : (bpp < 8 ? 1 : 4);
	} else {
	    ci->num_components = ccomps ? ccomps : (bpp < 8 ? 1 : 3);
	}
	if (bpp != 1 && ci->num_components == 1) {	/* We do dithered grays. */
	    bpp = bpp < 8 ? 8 : bpp;
	}
	ci->max_color = (1 << (bpp >> 2)) - 1;
	ci->max_gray = (bpp >= 8 ? 255 : 1);

	if (ci->num_components == 1) {
	    ci->dither_grays = (bpp >= 8 ? 5 : 2);
	    ci->dither_colors = (bpp >= 8 ? 5 : bpp > 1 ? 2 : 0);
	} else {
	    ci->dither_grays = (bpp > 8 ? 5 : 2);
	    ci->dither_colors = (bpp > 8 ? 5 : bpp > 1 ? 2 : 0);
	}
    } else {
	ci->num_components = (bpp == 1 || bpp == 8 ? 1 : 3);
	ci->max_color = (bpp >= 8 ? 255 : bpp > 1 ? 1 : 0);
	ci->max_gray = (bpp >= 8 ? 255 : 1);
	ci->dither_grays = (bpp >= 8 ? 5 : 2);
	ci->dither_colors = (bpp >= 8 ? 5 : bpp > 1 ? 2 : 0);
    }

    ci->depth = ((bpp > 1) && (bpp < 8) ? 8 : bpp);

    return 0;
}

/*
 * Map a CMYK color to a color index. We just use depth / 4 bits per color
 * to produce the color index.
 *
 * Important note: CMYK values are stored in the order K, C, M, Y because of
 * the way the HP drivers work.
 *
 */

#define gx_color_value_to_bits(cv, b) \
    ((cv) >> (gx_color_value_bits - (b)))
#define gx_bits_to_color_value(cv, b) \
    ((cv) << (gx_color_value_bits - (b)))

#define gx_cmyk_value_bits(c, m, y, k, b) \
    ((gx_color_value_to_bits((k), (b)) << (3 * (b))) | \
     (gx_color_value_to_bits((c), (b)) << (2 * (b))) | \
     (gx_color_value_to_bits((m), (b)) << (b)) | \
     (gx_color_value_to_bits((y), (b))))

#define gx_value_cmyk_bits(v, c, m, y, k, b) \
    (k) = gx_bits_to_color_value(((v) >> (3 * (b))) & ((1 << (b)) - 1), (b)), \
    (c) = gx_bits_to_color_value(((v) >> (2 * (b))) & ((1 << (b)) - 1), (b)), \
    (m) = gx_bits_to_color_value(((v) >> (b)) & ((1 << (b)) - 1), (b)), \
    (y) = gx_bits_to_color_value((v) & ((1 << (b)) - 1), (b))

private gx_color_index
gdev_cmyk_map_cmyk_color(gx_device * pdev,
			 gx_color_value cyan, gx_color_value magenta,
			 gx_color_value yellow,
			 gx_color_value black)
{

    gx_color_index color;

    switch (pdev->color_info.depth) {
	case 1:
	color = (cyan | magenta | yellow | black) > gx_max_color_value / 2 ?
	    (gx_color_index) 1 : (gx_color_index) 0;
	break;

	default:{
	    int nbits = pdev->color_info.depth;

	    if (cyan == magenta && magenta == yellow) {
		/* Convert CMYK to gray -- Red Book 6.2.2 */
		float bpart = ((float)cyan) * (lum_red_weight / 100.) +
		((float)magenta) * (lum_green_weight / 100.) +
		((float)yellow) * (lum_blue_weight / 100.) +
		(float)black;

		cyan = magenta = yellow = (gx_color_index) 0;
		black = (gx_color_index) (bpart > gx_max_color_value ?
					  gx_max_color_value : bpart);
	    }
	    color = gx_cmyk_value_bits(cyan, magenta, yellow, black,
				       nbits >> 2);
	}
    }

    return color;
}

/* Mapping of RGB colors to gray values. */

private gx_color_index
gdev_cmyk_map_rgb_color(gx_device * pdev, gx_color_value r, gx_color_value
			g, gx_color_value b)
{
    if (gx_color_value_to_byte(r & g & b) == 0xff) {
	return (gx_color_index) 0;	/* White */
    } else {
	gx_color_value c = gx_max_color_value - r;
	gx_color_value m = gx_max_color_value - g;
	gx_color_value y = gx_max_color_value - b;

	switch (pdev->color_info.depth) {
	    case 1:
	    return (c | m | y) > gx_max_color_value / 2 ?
		(gx_color_index) 1 : (gx_color_index) 0;
	    /*NOTREACHED */
	    break;

	    case 8:
	    return ((ulong) c * lum_red_weight * 10
		    + (ulong) m * lum_green_weight * 10
		    + (ulong) y * lum_blue_weight * 10)
		>> (gx_color_value_bits + 2);
	    /*NOTREACHED */
	    break;
	}
    }

    return (gx_color_index) 0;	/* This should never happen. */
}

/* Mapping of CMYK colors. */
private int
gdev_cmyk_map_color_rgb(gx_device * pdev, gx_color_index color,
			gx_color_value prgb[3])
{
    switch (pdev->color_info.depth) {
	case 1:
	prgb[0] = prgb[1] = prgb[2] = gx_max_color_value * (1 - color);
	break;

	case 8:
	if (pdev->color_info.num_components == 1) {
	    gx_color_value value = (gx_color_value) color ^ 0xff;

	    prgb[0] = prgb[1] = prgb[2] = (value << 8) + value;

	    break;
	}
	default:{
	    unsigned long bcyan, bmagenta, byellow, black;
	    int nbits = pdev->color_info.depth;

	    gx_value_cmyk_bits(color, bcyan, bmagenta, byellow, black,
			       nbits >> 2);

#ifdef USE_ADOBE_CMYK_RGB

	    /* R = 1.0 - min(1.0, C + K), etc. */

	    bcyan += black, bmagenta += black, byellow += black;
	    prgb[0] = (bcyan > gx_max_color_value ? (gx_color_value) 0 :
		       gx_max_color_value - bcyan);
	    prgb[1] = (bmagenta > gx_max_color_value ? (gx_color_value) 0 :
		       gx_max_color_value - bmagenta);
	    prgb[2] = (byellow > gx_max_color_value ? (gx_color_value) 0 :
		       gx_max_color_value - byellow);

#else

	    /* R = (1.0 - C) * (1.0 - K), etc. */

	    prgb[0] = (gx_color_value)
		((ulong) (gx_max_color_value - bcyan) *
		 (gx_max_color_value - black) / gx_max_color_value);
	    prgb[1] = (gx_color_value)
		((ulong) (gx_max_color_value - bmagenta) *
		 (gx_max_color_value - black) / gx_max_color_value);
	    prgb[2] = (gx_color_value)
		((ulong) (gx_max_color_value - byellow) *
		 (gx_max_color_value - black) / gx_max_color_value);

#endif

	}
    }

    return 0;
}

private gx_color_index
gdev_pcl_map_rgb_color(gx_device * pdev, gx_color_value r,
		       gx_color_value g, gx_color_value b)
{
    if (gx_color_value_to_byte(r & g & b) == 0xff)
	return (gx_color_index) 0;	/* white */
    else {
	gx_color_value c = gx_max_color_value - r;
	gx_color_value m = gx_max_color_value - g;
	gx_color_value y = gx_max_color_value - b;

	switch (pdev->color_info.depth) {
	    case 1:
	    return ((c | m | y) > gx_max_color_value / 2 ?
		    (gx_color_index) 1 : (gx_color_index) 0);
	    case 8:
	    if (pdev->color_info.num_components >= 3)
#define gx_color_value_to_1bit(cv) ((cv) >> (gx_color_value_bits - 1))
		return (gx_color_value_to_1bit(c) +
			(gx_color_value_to_1bit(m) << 1) +
			(gx_color_value_to_1bit(y) << 2));
	    else
#define red_weight 306
#define green_weight 601
#define blue_weight 117
		return ((((ulong) c * red_weight +
			  (ulong) m * green_weight +
			  (ulong) y * blue_weight)
			 >> (gx_color_value_bits + 2)));
	    case 16:
#define gx_color_value_to_5bits(cv) ((cv) >> (gx_color_value_bits - 5))
#define gx_color_value_to_6bits(cv) ((cv) >> (gx_color_value_bits - 6))
	    return (gx_color_value_to_5bits(y) +
		    (gx_color_value_to_6bits(m) << 5) +
		    (gx_color_value_to_5bits(c) << 11));
	    case 24:
	    return (gx_color_value_to_byte(y) +
		    (gx_color_value_to_byte(m) << 8) +
		    ((ulong) gx_color_value_to_byte(c) << 16));
	    case 32:
	    {
		return ((c == m && c == y) ? ((ulong)
					      gx_color_value_to_byte(c) << 24)
			: (gx_color_value_to_byte(y) +
			   (gx_color_value_to_byte(m) << 8) +
			   ((ulong) gx_color_value_to_byte(c) << 16)));
	    }
	}
    }
    return (gx_color_index) 0;	/* This never happens */
}

/* Map a color index to a r-g-b color. */
private int
gdev_pcl_map_color_rgb(gx_device * pdev, gx_color_index color,
		       gx_color_value prgb[3])
{
    /* For the moment, we simply ignore any black correction */
    switch (pdev->color_info.depth) {
	case 1:
	prgb[0] = prgb[1] = prgb[2] = -((gx_color_value) color ^ 1);
	break;
	case 8:
	if (pdev->color_info.num_components >= 3) {
	    gx_color_value c = (gx_color_value) color ^ 7;

	    prgb[0] = -(c & 1);
	    prgb[1] = -((c >> 1) & 1);
	    prgb[2] = -(c >> 2);
	} else {
	    gx_color_value value = (gx_color_value) color ^ 0xff;

	    prgb[0] = prgb[1] = prgb[2] = (value << 8) + value;
	}
	break;
	case 16:
	{
	    gx_color_value c = (gx_color_value) color ^ 0xffff;
	    ushort value = c >> 11;

	    prgb[0] = ((value << 11) + (value << 6) + (value << 1) +
		       (value >> 4)) >> (16 - gx_color_value_bits);
	    value = (c >> 6) & 0x3f;
	    prgb[1] = ((value << 10) + (value << 4) + (value >> 2))
		>> (16 - gx_color_value_bits);
	    value = c & 0x1f;
	    prgb[2] = ((value << 11) + (value << 6) + (value << 1) +
		       (value >> 4)) >> (16 - gx_color_value_bits);
	}
	break;
	case 24:
	{
	    gx_color_value c = (gx_color_value) color ^ 0xffffff;

	    prgb[0] = gx_color_value_from_byte(c >> 16);
	    prgb[1] = gx_color_value_from_byte((c >> 8) & 0xff);
	    prgb[2] = gx_color_value_from_byte(c & 0xff);
	}
	break;
	case 32:
#define  gx_maxcol gx_color_value_from_byte(gx_color_value_to_byte(gx_max_color_value))
	{
	    gx_color_value w = gx_maxcol - gx_color_value_from_byte(color >> 24);

	    prgb[0] = w - gx_color_value_from_byte((color >> 16) & 0xff);
	    prgb[1] = w - gx_color_value_from_byte((color >> 8) & 0xff);
	    prgb[2] = w - gx_color_value_from_byte(color & 0xff);
	}
	break;
    }
    return 0;
}

/* new_bpp == save_bpp or new_bpp == 0 means don't change bpp.
   ccomps == 0 means don't change number of color comps.
   If new_bpp != 0, it must be the value of the BitsPerPixel element of
   the plist; real_bpp may differ from new_bpp.
 */
private int
cdj_put_param_bpp(gx_device * pdev, gs_param_list * plist, int new_bpp,
		  int real_bpp, int ccomps)
{
    if (new_bpp == 0 && ccomps == 0)
	return gdev_prn_put_params(pdev, plist);
    else {
	gx_device_color_info save_info;
	int save_bpp;
	int code;

	save_info = pdev->color_info;
	save_bpp = save_info.depth;
#define save_ccomps save_info.num_components
	if (save_bpp == 8 && save_ccomps == 3 && !cprn_device->cmyk)
	    save_bpp = 3;
	code = cdj_set_bpp(pdev, real_bpp, ccomps);
	if (code < 0) {
	    param_signal_error(plist, "BitsPerPixel", code);
	    param_signal_error(plist, "ProcessColorModel", code);
	    return code;
	}
	pdev->color_info.depth = new_bpp;	/* cdj_set_bpp maps 3/6 to 8 */
	code = gdev_prn_put_params(pdev, plist);
	if (code < 0) {
	    cdj_set_bpp(pdev, save_bpp, save_ccomps);
	    return code;
	}
	cdj_set_bpp(pdev, real_bpp, ccomps);	/* reset depth if needed */
	if ((cdj850->color_info.depth != save_bpp ||
	     (ccomps != 0 && ccomps != save_ccomps))
	    && pdev->is_open)
	    return gs_closedevice(pdev);
	return 0;
#undef save_ccomps
    }
}

/* the following code was in the original driver but is unused

 * private int
 * x_mul_div (int a, int b, int c)
 * { 
 *   int result;
 *   
 *   result = (int) ((a * b) / c) ;
 *  return result;
 * }
 *   
 * private void 
 * save_color_data(int size,
 *              byte * current,
 *              byte * saved)
 * {
 *   int i;
 *   for (i=0;i<size;i++){
 *     *saved++ = *current++;
 *   }
 *   return;
 * }
 * 
 * private int 
 * test_scan (int size, 
 *         byte * current,
 *         byte * last, 
 *         byte * control)
 *   {
 *   int error = 0;
 *   int i;
 *   
 *   for (i=0;i<size;i++){
 *     if (*control != *last){
 *       error = 1;
 *     }
 *     *control = *current;
 * 
 *     control++;
 *     last++;
 *     current++;
 *   }
 *   return error;
 * }
 * 
 * * Transform from cmy into hsv 
 * private void
 * cmy2hsv(int *c, int *m, int *y, int *h, int *s, int *v)
 * {
 *   int hue;
 *   int r, g, b;
 *   int r1, g1, b1;
 *   int maxValue, minValue, diff;
 *   
 *   r = 255 - *c;
 *   g = 255 - *m;
 *   b = 255 - *y;
 *   
 *   maxValue = max(r, max(g,b));
 *   minValue = min(r,min(g,b));
 *   diff = maxValue - minValue;
 *   *v = maxValue;
 *   
 *   if (maxValue != 0)
 *     *s = x_mul_div(diff,255,maxValue);
 *   else
 *     *s = 0;
 *   
 *   if (*s == 0)
 *     {
 *       hue = 0;
 *     }
 *   else
 *     {
 *       r1 = x_mul_div(maxValue - r,255,diff);
 *       g1 = x_mul_div(maxValue - g,255,diff);
 *       b1 = x_mul_div(maxValue - b,255,diff);
 *       
 *       if (r == maxValue)
 *      hue = b1 - g1;
 *       else if (g == maxValue)
 *      hue = 510 + r1 - b1;
 *       else
 *      hue = 1020 + g1 - r1;
 *       
 *       if (hue < 0)
 *      hue += 1530;
 *     }
 *   
 *   *h = (hue + 3) / 6;
 * 
 *   return;
 * }
 * end of unused code */


/************************ the routines for the cdj1600 printer ***************/

/* Configure the printer and start Raster mode */
private void
cdj1600_start_raster_mode(gx_device_printer * pdev, int paper_size,
			  FILE * prn_stream)
{
    uint raster_width = pdev->width -
    pdev->x_pixels_per_inch * (dev_l_margin(pdev) + dev_r_margin(pdev));

    /* switch to PCL control language */
    fputs("\033%-12345X@PJL enter language = PCL\n", prn_stream);

    fputs("\033*rbC", prn_stream);	/* End raster graphics */
    fputs("\033E", prn_stream);	/* Reset */

    /* resolution */
    fprintf(prn_stream, "\033*t%dR", (int)cdj850->x_pixels_per_inch);

    /* Page size, orientation, top margin & perforation skip */
    fprintf(prn_stream, "\033&l%daolE", paper_size);

    /* no negative motion */
    fputs("\033&a1N", prn_stream);

    /* Print Quality, -1 = draft, 0 = normal, 1 = presentation */
    fprintf(prn_stream, "\033*o%dQ", cdj850->quality);

    /* Media Type,0 = plain paper, 1 = bond paper, 2 = special
       paper, 3 = glossy film, 4 = transparency film */
    fprintf(prn_stream, "\033&l%dM", cdj850->papertype);

    /* Move to top left of printed area */
    fprintf(prn_stream, "\033*p%dY", (int)(300.0 * DOFFSET));

    /* raster width and number of planes */
    fprintf(prn_stream, "\033*r%ds-%du0A",
	    raster_width, pdev->color_info.num_components);

    /* start raster graphics */
    fputs("\033*r1A", prn_stream);

    /* From now on, all escape commands start with \033*b, so we
     * combine them (if the printer supports this). */
    fputs("\033*b", prn_stream);

    /* Set compression if the mode has been defined. */
    if (cdj850->compression)
	fprintf(prn_stream, "%dm", cdj850->compression);

    return;
}				/* end configure raster-mode */

/* print_plane compresses (mode 3) and outputs one plane */
private void
print_c3plane(FILE * prn_stream, char plane_code, int plane_size,
	      const byte * curr, byte * prev, byte * out_data)
{
    /* Compress the output data */
    int out_count = gdev_pcl_mode3compress(plane_size, curr, prev, out_data);

    /* and output the data */
    if (out_count > 0) {
	fprintf(prn_stream, "%d%c", out_count, plane_code);
	fwrite(out_data, sizeof(byte), out_count, prn_stream);
    } else {
	putc(plane_code, prn_stream);
    }
}

private int
copy_color_data(byte * dest, const byte * src, int n)
{
    /* copy word by word */
    register int i = n / 4;
    register word *d = (word *) dest;
    register const word *s = (const word *)src;

    while (i-- > 0) {
	*d++ = *s++;
    }
    return n;
}

/* Printing non-blank lines */
private void
cdj1600_print_non_blank_lines(gx_device_printer * pdev,
			      struct ptr_arrays *data_ptrs,
			      struct misc_struct *misc_vars,
			      struct error_val_field *error_values,
			      const Gamma *gamma,
			      FILE * prn_stream)
{
    int i, plane_size_c;

    /* copy data to data_c in order to make do_floyd_steinberg work */
    plane_size_c = copy_color_data
	(data_ptrs->data_c[misc_vars->cscan],
	 data_ptrs->data[misc_vars->scan],
	 misc_vars->databuff_size) / misc_vars->storage_bpp;

    /* dither the color planes */
    do_floyd_steinberg(misc_vars->scan, misc_vars->cscan,
		       misc_vars->plane_size, plane_size_c,
		       misc_vars->num_comps, data_ptrs, pdev, error_values);

    /* Transfer raster graphics in the order C, M, Y, that is
       planes 2,1,0 */
    for (i = misc_vars->num_comps - 1; i >= 0; i--) {

	/* output the lower color planes */
	print_c3plane(prn_stream, "wvv"[i], plane_size_c,
		      data_ptrs->plane_data_c[misc_vars->cscan][i],
		      data_ptrs->plane_data_c[1 - misc_vars->cscan][i],
		      data_ptrs->out_data);
    }				/* End For i = num_comps */
    misc_vars->cscan = 1 - misc_vars->cscan;
}

private void
cdj1600_terminate_page(gx_device_printer * pdev, FILE * prn_stream)
{
    cdj850_terminate_page(pdev, prn_stream);
    fputs("\033%-12345X", prn_stream);
}
