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.1 by greg, Thu Feb 2 10:34:29 1989 UTC vs.
Revision 1.12 by greg, Sat Sep 22 10:45:42 1990 UTC

# Line 14 | Line 14 | static char SCCSid[] = "$SunId$ LBL";
14  
15   #include  "color.h"
16  
17 – #ifdef  SPEC_RGB
18 – /*
19 – *      The following table contains the CIE tristimulus integrals
20 – *  for X, Y, and Z.  The table is cumulative, so that
21 – *  each color coordinate integrates to 1.
22 – */
17  
24 – #define  STARTWL        380             /* starting wavelength (nanometers) */
25 – #define  INCWL          10              /* wavelength increment */
26 – #define  NINC           40              /* # of values */
27 –
28 – static BYTE  chroma[3][NINC] = {
29 –        {                                                       /* X */
30 –                0,   0,   0,   2,   6,   13,  22,  30,  36,  41,
31 –                42,  43,  43,  44,  46,  52,  60,  71,  87,  106,
32 –                128, 153, 178, 200, 219, 233, 243, 249, 252, 254,
33 –                255, 255, 255, 255, 255, 255, 255, 255, 255, 255
34 –        }, {                                                    /* Y */
35 –                0,   0,   0,   0,   0,   1,   2,   4,   7,   11,
36 –                17,  24,  34,  48,  64,  84,  105, 127, 148, 169,
37 –                188, 205, 220, 232, 240, 246, 250, 253, 254, 255,
38 –                255, 255, 255, 255, 255, 255, 255, 255, 255, 255
39 –        }, {                                                    /* Z */
40 –                0,   0,   2,   10,  32,  66,  118, 153, 191, 220,
41 –                237, 246, 251, 253, 254, 255, 255, 255, 255, 255,
42 –                255, 255, 255, 255, 255, 255, 255, 255, 255, 255,
43 –                255, 255, 255, 255, 255, 255, 255, 255, 255, 255
44 –        }
45 – };
46 –
47 –
48 – spec_rgb(col, s, e)             /* comput RGB color from spectral range */
49 – COLOR  col;
50 – int  s, e;
51 – {
52 –        COLOR  ciecolor;
53 –
54 –        spec_cie(ciecolor, s, e);
55 –        cie_rgb(col, ciecolor);
56 – }
57 –
58 –
59 – spec_cie(col, s, e)             /* compute a color from a spectral range */
60 – COLOR  col;             /* returned color */
61 – int  s, e;              /* starting and ending wavelengths */
62 – {
63 –        register int  i, d, r;
64 –        
65 –        s -= STARTWL;
66 –        if (s < 0)
67 –                s = 0;
68 –
69 –        e -= STARTWL;
70 –        if (e >= INCWL*(NINC - 1))
71 –                e = INCWL*(NINC - 1) - 1;
72 –
73 –        d = e / INCWL;                  /* interpolate values */
74 –        r = e % INCWL;
75 –        for (i = 0; i < 3; i++)
76 –                col[i] = chroma[i][d]*(INCWL - r) + chroma[i][d + 1]*r;
77 –
78 –        d = s / INCWL;
79 –        r = s % INCWL;
80 –        for (i = 0; i < 3; i++)
81 –                col[i] -= chroma[i][d]*(INCWL - r) - chroma[i][d + 1]*r;
82 –
83 –        col[RED] = (col[RED] + 0.5) / (256*INCWL);
84 –        col[GRN] = (col[GRN] + 0.5) / (256*INCWL);
85 –        col[BLU] = (col[BLU] + 0.5) / (256*INCWL);
86 – }
87 –
88 –
89 – cie_rgb(rgbcolor, ciecolor)             /* convert CIE to RGB (NTSC) */
90 – register COLOR  rgbcolor, ciecolor;
91 – {
92 –        static float  cmat[3][3] = {
93 –                1.73, -.48, -.26,
94 –                -.81, 1.65, -.02,
95 –                 .08, -.17, 1.28,
96 –        };
97 –        register int  i;
98 –
99 –        for (i = 0; i < 3; i++) {
100 –                rgbcolor[i] =   cmat[i][0]*ciecolor[0] +
101 –                                cmat[i][1]*ciecolor[1] +
102 –                                cmat[i][2]*ciecolor[2] ;
103 –                if (rgbcolor[i] < 0.0)
104 –                        rgbcolor[i] = 0.0;
105 –        }
106 – }
107 – #endif
108 –
109 –
18   fwritecolrs(scanline, len, fp)          /* write out a colr scanline */
19   register COLR  *scanline;
20   int  len;
# Line 289 | Line 197 | double  r, g, b;
197          d = r > g ? r : g;
198          if (b > d) d = b;
199  
200 <        if (d <= 0.0) {
200 >        if (d <= 1e-32) {
201                  clr[RED] = clr[GRN] = clr[BLU] = 0;
202                  clr[EXP] = 0;
203                  return;
# Line 308 | Line 216 | colr_color(col, clr)           /* convert short to float color
216   register COLOR  col;
217   register COLR  clr;
218   {
219 <        double  ldexp(), f;
219 >        double  f;
220          
221          if (clr[EXP] == 0)
222                  col[RED] = col[GRN] = col[BLU] = 0.0;
223          else {
224 <                f = ldexp(1.0, clr[EXP]-(COLXS+8));
224 >                f = ldexp(1.0, (int)clr[EXP]-(COLXS+8));
225                  col[RED] = (clr[RED] + 0.5)*f;
226                  col[GRN] = (clr[GRN] + 0.5)*f;
227                  col[BLU] = (clr[BLU] + 0.5)*f;
# Line 321 | Line 229 | register COLR  clr;
229   }
230  
231  
232 < #ifdef  FREXP
233 < double
234 < frexp(x, ip)            /* call it paranoia, I've seen the lib version */
235 < register double  x;
328 < int  *ip;
232 > normcolrs(scan, len, adjust)    /* normalize a scanline of colrs */
233 > register COLR  *scan;
234 > int  len;
235 > int  adjust;
236   {
237 <        int  neg;
238 <        register int  i;
237 >        register int  c;
238 >        register int  shift;
239  
240 <        if (neg = (x < 0.0))
241 <                x = -x;
242 <        else if (x == 0.0) {
243 <                *ip = 0;
244 <                return(0.0);
240 >        while (len-- > 0) {
241 >                shift = scan[0][EXP] + adjust - COLXS;
242 >                if (shift > 0) {
243 >                        if (shift > 8) {
244 >                                scan[0][RED] =
245 >                                scan[0][GRN] =
246 >                                scan[0][BLU] = 255;
247 >                        } else {
248 >                                shift--;
249 >                                c = (scan[0][RED]<<1 | 1) << shift;
250 >                                scan[0][RED] = c > 255 ? 255 : c;
251 >                                c = (scan[0][GRN]<<1 | 1) << shift;
252 >                                scan[0][GRN] = c > 255 ? 255 : c;
253 >                                c = (scan[0][BLU]<<1 | 1) << shift;
254 >                                scan[0][BLU] = c > 255 ? 255 : c;
255 >                        }
256 >                } else if (shift < 0) {
257 >                        if (shift < -8) {
258 >                                scan[0][RED] =
259 >                                scan[0][GRN] =
260 >                                scan[0][BLU] = 0;
261 >                        } else {
262 >                                shift = -1-shift;
263 >                                scan[0][RED] = ((scan[0][RED]>>shift)+1)>>1;
264 >                                scan[0][GRN] = ((scan[0][GRN]>>shift)+1)>>1;
265 >                                scan[0][BLU] = ((scan[0][BLU]>>shift)+1)>>1;
266 >                        }
267 >                }
268 >                scan[0][EXP] = COLXS - adjust;
269 >                scan++;
270          }
339 –        if (x < 0.5)
340 –                for (i = 0; x < 0.5; i--)
341 –                        x *= 2.0;
342 –        else
343 –                for (i = 0; x >= 1.0; i++)
344 –                        x /= 2.0;
345 –        *ip = i;
346 –        if (neg)
347 –                return(-x);
348 –        else
349 –                return(x);
271   }
272 < #endif
272 >
273 >
274 > bigdiff(c1, c2, md)                     /* c1 delta c2 > md? */
275 > register COLOR  c1, c2;
276 > double  md;
277 > {
278 >        register int  i;
279 >
280 >        for (i = 0; i < 3; i++)
281 >                if (colval(c1,i)-colval(c2,i) > md*colval(c2,i) ||
282 >                        colval(c2,i)-colval(c1,i) > md*colval(c1,i))
283 >                        return(1);
284 >        return(0);
285 > }

Diff Legend

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