--- ray/src/px/clrtab.c	1993/02/05 09:45:40	2.4
+++ ray/src/px/clrtab.c	2003/05/13 17:58:33	2.14
@@ -1,9 +1,6 @@
-/* Copyright (c) 1993 Regents of the University of California */
-
 #ifndef lint
-static char SCCSid[] = "$SunId$ LBL";
+static const char	RCSid[] = "$Id: clrtab.c,v 2.14 2003/05/13 17:58:33 greg Exp $";
 #endif
-
 /*
  * Simple median-cut color quantization based on colortab.c
  */
@@ -23,7 +20,7 @@ static char SCCSid[] = "$SunId$ LBL";
 #define part(cn)	((cn)>>2)
 #define prim(cn)	((cn)&3)
 				/* our color table (global) */
-BYTE	clrtab[256][3];
+extern BYTE	clrtab[256][3];
 				/* histogram of colors / color assignments */
 static unsigned	histo[NRED][NGRN][NBLU];
 #define cndx(c)		histo[((c)[RED]*NRED)>>8][((c)[GRN]*NGRN)>>8][((c)[BLU]*NBLU)>>8]
@@ -40,10 +37,16 @@ static int	CLRCUBE[3][2] = {0,NRED,0,NGRN,0,NBLU};
 #endif
 #endif
 
+static	cut(), mktabent(), closest(), addneigh(), setclosest();
+static int	split();
+static unsigned	dist();
 
-new_histo()		/* clear our histogram */
+
+new_histo(n)		/* clear our histogram */
+int	n;
 {
-	bzero((char *)histo, sizeof(histo));
+	bzero((void *)histo, sizeof(histo));
+	return(0);
 }
 
 
@@ -77,6 +80,8 @@ int	ncolors;
 #ifdef CLOSEST
 	closest(ncolors);	/* ensure colors picked are closest */
 #endif
+				/* reset dithering function */
+	dith_colrs((BYTE *)NULL, (COLR *)NULL, 0);
 				/* return new color table size */
 	return(ncolors);
 }
@@ -107,16 +112,17 @@ register BYTE	*bs;
 register COLR	*cs;
 int	n;
 {
-	static short	(*cerr)[3];
+	static short	(*cerr)[3] = NULL;
 	static int	N = 0;
 	int	err[3], errp[3];
 	register int	x, i;
 
 	if (n != N) {		/* get error propogation array */
-		if (N)
-			cerr = (short (*)[3])realloc((char *)cerr,
-					3*n*sizeof(short));
-		else
+		if (N) {
+			free((void *)cerr);
+			cerr = NULL;
+		}
+		if (n)
 			cerr = (short (*)[3])malloc(3*n*sizeof(short));
 		if (cerr == NULL) {
 			N = 0;
@@ -124,7 +130,7 @@ int	n;
 			return;
 		}
 		N = n;
-		bzero((char *)cerr, 3*N*sizeof(short));
+		bzero((void *)cerr, 3*N*sizeof(short));
 	}
 	err[0] = err[1] = err[2] = 0;
 	for (x = 0; x < n; x++) {
@@ -163,7 +169,7 @@ int	c0, c1;
 	}
 					/* split box */
 	branch = split(box);
-	bcopy((char *)box, (char *)kb, sizeof(kb));
+	bcopy((void *)box, (void *)kb, sizeof(kb));
 						/* do left (lesser) branch */
 	kb[prim(branch)][1] = part(branch);
 	cut(kb, c0, (c0+c1)>>1);
@@ -181,7 +187,7 @@ register int	box[3][2];
 #define c0	r
 	register int	r, g, b;
 	int	pri;
-	int	t[HMAX], med;
+	long	t[HMAX], med;
 					/* find dominant axis */
 	pri = RED;
 	if (box[GRN][1]-box[GRN][0] > box[pri][1]-box[pri][0])
@@ -237,9 +243,10 @@ mktabent(p, box)	/* compute average color for box and 
 int	p;
 register int	box[3][2];
 {
-	long	sum[3];
-	int	r, g, n;
-	register int	b, c;
+	unsigned long	sum[3];
+	unsigned	r, g;
+	unsigned long	n;
+	register unsigned	b, c;
 						/* sum pixels in box */
 	n = 0;
 	sum[RED] = sum[GRN] = sum[BLU] = 0;
@@ -254,6 +261,12 @@ register int	box[3][2];
 		    }
 		    histo[r][g][b] = p;		/* assign pixel */
 		}
+	if (n >= (1L<<23)/HMAX) {		/* avoid overflow */
+		sum[RED] /= n;
+		sum[GRN] /= n;
+		sum[BLU] /= n;
+		n = 1;
+	}
 	if (n) {				/* compute average */
 		clrtab[p][RED] = sum[RED]*256/NRED/n;
 		clrtab[p][GRN] = sum[GRN]*256/NGRN/n;
@@ -324,7 +337,8 @@ int	j;
 		if (nl[i][t] == i) {		/* add to list */
 			nl[i][t++] = j;
 			if (t % NBSIZ == 0) {	/* enlarge list */
-				if ((nnl = realloc(nl[i], t+NBSIZ)) == NULL)
+				if ((nnl = realloc((void *)nl[i],
+						t+NBSIZ)) == NULL)
 					t--;
 				else
 					nl[i] = (BYTE *)nnl;