ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/common/color.c
(Generate patch)

Comparing ray/src/common/color.c (file contents):
Revision 1.13 by greg, Fri Oct 19 11:13:02 1990 UTC vs.
Revision 2.14 by greg, Tue Dec 9 15:51:42 2003 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1986 Regents of the University of California */
2 –
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ LBL";
2 > static const char       RCSid[] = "$Id$";
3   #endif
6 –
4   /*
5   *  color.c - routines for color calculations.
6   *
7 < *     10/10/85
7 > *  Externals declared in color.h
8   */
9  
10 + #include "copyright.h"
11 +
12   #include  <stdio.h>
13  
14 + #include  <stdlib.h>
15 +
16 + #include  <math.h>
17 +
18   #include  "color.h"
19  
20 + #ifdef getc_unlocked            /* avoid horrendous overhead of flockfile */
21 + #define getc    getc_unlocked
22 + #define putc    putc_unlocked
23 + #endif
24  
25 + #define  MINELEN        8       /* minimum scanline length for encoding */
26 + #define  MAXELEN        0x7fff  /* maximum scanline length for encoding */
27 + #define  MINRUN         4       /* minimum run length */
28 +
29 +
30 + char *
31 + tempbuffer(len)                 /* get a temporary buffer */
32 + unsigned int  len;
33 + {
34 +        static char  *tempbuf = NULL;
35 +        static unsigned  tempbuflen = 0;
36 +
37 +        if (len > tempbuflen) {
38 +                if (tempbuflen > 0)
39 +                        tempbuf = (char *)realloc((void *)tempbuf, len);
40 +                else
41 +                        tempbuf = (char *)malloc(len);
42 +                tempbuflen = tempbuf==NULL ? 0 : len;
43 +        }
44 +        return(tempbuf);
45 + }
46 +
47 +
48 + int
49   fwritecolrs(scanline, len, fp)          /* write out a colr scanline */
50   register COLR  *scanline;
51   int  len;
52   register FILE  *fp;
53   {
54 <        COLR  lastcolr;
55 <        int  rept;
54 >        register int  i, j, beg, cnt = 1;
55 >        int  c2;
56          
57 <        lastcolr[RED] = lastcolr[GRN] = lastcolr[BLU] = 1;
58 <        lastcolr[EXP] = 0;
59 <        rept = 0;
60 <        
61 <        while (len > 0) {
62 <                if (scanline[0][EXP] == lastcolr[EXP] &&
63 <                                scanline[0][RED] == lastcolr[RED] &&
64 <                                scanline[0][GRN] == lastcolr[GRN] &&
65 <                                scanline[0][BLU] == lastcolr[BLU])
66 <                        rept++;
67 <                else {
68 <                        while (rept) {          /* write out count */
69 <                                putc(1, fp);
70 <                                putc(1, fp);
71 <                                putc(1, fp);
72 <                                putc(rept & 255, fp);
73 <                                rept >>= 8;
57 >        if ((len < MINELEN) | (len > MAXELEN))  /* OOBs, write out flat */
58 >                return(fwrite((char *)scanline,sizeof(COLR),len,fp) - len);
59 >                                        /* put magic header */
60 >        putc(2, fp);
61 >        putc(2, fp);
62 >        putc(len>>8, fp);
63 >        putc(len&255, fp);
64 >                                        /* put components seperately */
65 >        for (i = 0; i < 4; i++) {
66 >            for (j = 0; j < len; j += cnt) {    /* find next run */
67 >                for (beg = j; beg < len; beg += cnt) {
68 >                    for (cnt = 1; cnt < 127 && beg+cnt < len &&
69 >                            scanline[beg+cnt][i] == scanline[beg][i]; cnt++)
70 >                        ;
71 >                    if (cnt >= MINRUN)
72 >                        break;                  /* long enough */
73 >                }
74 >                if (beg-j > 1 && beg-j < MINRUN) {
75 >                    c2 = j+1;
76 >                    while (scanline[c2++][i] == scanline[j][i])
77 >                        if (c2 == beg) {        /* short run */
78 >                            putc(128+beg-j, fp);
79 >                            putc(scanline[j][i], fp);
80 >                            j = beg;
81 >                            break;
82                          }
44 –                        putc(scanline[0][RED], fp);     /* new color */
45 –                        putc(scanline[0][GRN], fp);
46 –                        putc(scanline[0][BLU], fp);
47 –                        putc(scanline[0][EXP], fp);
48 –                        copycolr(lastcolr, scanline[0]);
49 –                        rept = 0;
83                  }
84 <                scanline++;
85 <                len--;
84 >                while (j < beg) {               /* write out non-run */
85 >                    if ((c2 = beg-j) > 128) c2 = 128;
86 >                    putc(c2, fp);
87 >                    while (c2--)
88 >                        putc(scanline[j++][i], fp);
89 >                }
90 >                if (cnt >= MINRUN) {            /* write out run */
91 >                    putc(128+cnt, fp);
92 >                    putc(scanline[beg][i], fp);
93 >                } else
94 >                    cnt = 0;
95 >            }
96          }
54 –        while (rept) {          /* write out count */
55 –                putc(1, fp);
56 –                putc(1, fp);
57 –                putc(1, fp);
58 –                putc(rept & 255, fp);
59 –                rept >>= 8;
60 –        }
97          return(ferror(fp) ? -1 : 0);
98   }
99  
100  
101 < freadcolrs(scanline, len, fp)           /* read in a colr scanline */
101 > static int
102 > oldreadcolrs(scanline, len, fp)         /* read in an old colr scanline */
103   register COLR  *scanline;
104   int  len;
105   register FILE  *fp;
# Line 98 | Line 135 | register FILE  *fp;
135   }
136  
137  
138 + int
139 + freadcolrs(scanline, len, fp)           /* read in an encoded colr scanline */
140 + register COLR  *scanline;
141 + int  len;
142 + register FILE  *fp;
143 + {
144 +        register int  i, j;
145 +        int  code, val;
146 +                                        /* determine scanline type */
147 +        if ((len < MINELEN) | (len > MAXELEN))
148 +                return(oldreadcolrs(scanline, len, fp));
149 +        if ((i = getc(fp)) == EOF)
150 +                return(-1);
151 +        if (i != 2) {
152 +                ungetc(i, fp);
153 +                return(oldreadcolrs(scanline, len, fp));
154 +        }
155 +        scanline[0][GRN] = getc(fp);
156 +        scanline[0][BLU] = getc(fp);
157 +        if ((i = getc(fp)) == EOF)
158 +                return(-1);
159 +        if (scanline[0][GRN] != 2 || scanline[0][BLU] & 128) {
160 +                scanline[0][RED] = 2;
161 +                scanline[0][EXP] = i;
162 +                return(oldreadcolrs(scanline+1, len-1, fp));
163 +        }
164 +        if ((scanline[0][BLU]<<8 | i) != len)
165 +                return(-1);             /* length mismatch! */
166 +                                        /* read each component */
167 +        for (i = 0; i < 4; i++)
168 +            for (j = 0; j < len; ) {
169 +                if ((code = getc(fp)) == EOF)
170 +                    return(-1);
171 +                if (code > 128) {       /* run */
172 +                    code &= 127;
173 +                    if ((val = getc(fp)) == EOF)
174 +                        return -1;
175 +                    while (code--)
176 +                        scanline[j++][i] = val;
177 +                } else                  /* non-run */
178 +                    while (code--) {
179 +                        if ((val = getc(fp)) == EOF)
180 +                            return -1;
181 +                        scanline[j++][i] = val;
182 +                    }
183 +            }
184 +        return(0);
185 + }
186 +
187 +
188 + int
189   fwritescan(scanline, len, fp)           /* write out a scanline */
190   register COLOR  *scanline;
191   int  len;
192 < register FILE  *fp;
192 > FILE  *fp;
193   {
194 <        COLR  lastcolr, thiscolr;
195 <        int  rept;
196 <        
197 <        lastcolr[RED] = lastcolr[GRN] = lastcolr[BLU] = 1;
198 <        lastcolr[EXP] = 0;
199 <        rept = 0;
200 <        
201 <        while (len > 0) {
202 <                setcolr(thiscolr, scanline[0][RED],
194 >        COLR  *clrscan;
195 >        int  n;
196 >        register COLR  *sp;
197 >                                        /* get scanline buffer */
198 >        if ((sp = (COLR *)tempbuffer(len*sizeof(COLR))) == NULL)
199 >                return(-1);
200 >        clrscan = sp;
201 >                                        /* convert scanline */
202 >        n = len;
203 >        while (n-- > 0) {
204 >                setcolr(sp[0], scanline[0][RED],
205                                    scanline[0][GRN],
206                                    scanline[0][BLU]);
117 –                if (thiscolr[EXP] == lastcolr[EXP] &&
118 –                                thiscolr[RED] == lastcolr[RED] &&
119 –                                thiscolr[GRN] == lastcolr[GRN] &&
120 –                                thiscolr[BLU] == lastcolr[BLU])
121 –                        rept++;
122 –                else {
123 –                        while (rept) {          /* write out count */
124 –                                putc(1, fp);
125 –                                putc(1, fp);
126 –                                putc(1, fp);
127 –                                putc(rept & 255, fp);
128 –                                rept >>= 8;
129 –                        }
130 –                        putc(thiscolr[RED], fp);        /* new color */
131 –                        putc(thiscolr[GRN], fp);
132 –                        putc(thiscolr[BLU], fp);
133 –                        putc(thiscolr[EXP], fp);
134 –                        copycolr(lastcolr, thiscolr);
135 –                        rept = 0;
136 –                }
207                  scanline++;
208 <                len--;
208 >                sp++;
209          }
210 <        while (rept) {          /* write out count */
141 <                putc(1, fp);
142 <                putc(1, fp);
143 <                putc(1, fp);
144 <                putc(rept & 255, fp);
145 <                rept >>= 8;
146 <        }
147 <        return(ferror(fp) ? -1 : 0);
210 >        return(fwritecolrs(clrscan, len, fp));
211   }
212  
213  
214 + int
215   freadscan(scanline, len, fp)            /* read in a scanline */
216   register COLOR  *scanline;
217   int  len;
218 < register FILE  *fp;
218 > FILE  *fp;
219   {
220 <        COLR  thiscolr;
221 <        int  rshift;
222 <        register int  i;
223 <        
224 <        rshift = 0;
225 <        
226 <        while (len > 0) {
227 <                thiscolr[RED] = getc(fp);
228 <                thiscolr[GRN] = getc(fp);
229 <                thiscolr[BLU] = getc(fp);
230 <                thiscolr[EXP] = getc(fp);
231 <                if (feof(fp) || ferror(fp))
232 <                        return(-1);
233 <                if (thiscolr[RED] == 1 &&
234 <                                thiscolr[GRN] == 1 &&
235 <                                thiscolr[BLU] == 1) {
236 <                        for (i = thiscolr[EXP] << rshift; i > 0; i--) {
173 <                                copycolor(scanline[0], scanline[-1]);
174 <                                scanline++;
175 <                                len--;
176 <                        }
177 <                        rshift += 8;
178 <                } else {
179 <                        colr_color(scanline[0], thiscolr);
180 <                        scanline++;
181 <                        len--;
182 <                        rshift = 0;
183 <                }
220 >        register COLR  *clrscan;
221 >
222 >        if ((clrscan = (COLR *)tempbuffer(len*sizeof(COLR))) == NULL)
223 >                return(-1);
224 >        if (freadcolrs(clrscan, len, fp) < 0)
225 >                return(-1);
226 >                                        /* convert scanline */
227 >        colr_color(scanline[0], clrscan[0]);
228 >        while (--len > 0) {
229 >                scanline++; clrscan++;
230 >                if (clrscan[0][RED] == clrscan[-1][RED] &&
231 >                            clrscan[0][GRN] == clrscan[-1][GRN] &&
232 >                            clrscan[0][BLU] == clrscan[-1][BLU] &&
233 >                            clrscan[0][EXP] == clrscan[-1][EXP])
234 >                        copycolor(scanline[0], scanline[-1]);
235 >                else
236 >                        colr_color(scanline[0], clrscan[0]);
237          }
238          return(0);
239   }
240  
241  
242 + void
243   setcolr(clr, r, g, b)           /* assign a short color value */
244   register COLR  clr;
245   double  r, g, b;
246   {
193 –        double  frexp();
247          double  d;
248          int  e;
249          
# Line 203 | Line 256 | double  r, g, b;
256                  return;
257          }
258  
259 <        d = frexp(d, &e) * 256.0 / d;
259 >        d = frexp(d, &e) * 255.9999 / d;
260  
261 <        clr[RED] = r * d;
262 <        clr[GRN] = g * d;
263 <        clr[BLU] = b * d;
261 >        if (r > 0.0)
262 >                clr[RED] = r * d;
263 >        else
264 >                clr[RED] = 0;
265 >        if (g > 0.0)
266 >                clr[GRN] = g * d;
267 >        else
268 >                clr[GRN] = 0;
269 >        if (b > 0.0)
270 >                clr[BLU] = b * d;
271 >        else
272 >                clr[BLU] = 0;
273 >
274          clr[EXP] = e + COLXS;
275   }
276  
277  
278 + void
279   colr_color(col, clr)            /* convert short to float color */
280   register COLOR  col;
281   register COLR  clr;
# Line 229 | Line 293 | register COLR  clr;
293   }
294  
295  
296 + int
297   bigdiff(c1, c2, md)                     /* c1 delta c2 > md? */
298   register COLOR  c1, c2;
299   double  md;

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines