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

Comparing src/common/colrops.c (file contents):
Revision 2.1 by greg, Tue Nov 12 16:55:43 1991 UTC vs.
Revision 2.3 by greg, Tue Sep 8 10:04:33 1992 UTC

# Line 1 | Line 1
1 < /* Copyright (c) 1990 Regents of the University of California */
1 > /* Copyright (c) 1992 Regents of the University of California */
2  
3   #ifndef lint
4   static char SCCSid[] = "$SunId$ LBL";
# Line 10 | Line 10 | static char SCCSid[] = "$SunId$ LBL";
10  
11   #include "color.h"
12  
13 < #define MAXGSHIFT       15              /* maximum shift for gamma table */
13 > #define NULL            0
14  
15 < static BYTE     g_mant[256], g_nexp[256];
15 > extern char     *bmalloc();
16  
17 < static BYTE     g_bval[MAXGSHIFT+1][256];
17 > #define MAXGSHIFT       31              /* maximum shift for gamma table */
18  
19 + static BYTE     *g_mant = NULL, *g_nexp = NULL;
20  
21 < setcolrgam(g)                   /* set gamma conversion */
22 < double  g;
21 > static BYTE     (*g_bval)[256] = NULL;
22 >
23 > extern double   pow();
24 >
25 >
26 > setcolrcor(f, a2)               /* set brightness correction */
27 > double  (*f)();
28 > double  a2;
29   {
23 –        extern double   pow();
30          double  mult;
31          register int    i, j;
32 +                                        /* allocate tables */
33 +        if (g_bval == NULL && (g_bval =
34 +                        (BYTE (*)[256])bmalloc((MAXGSHIFT+1)*256)) == NULL)
35 +                return(-1);
36                                          /* compute colr -> gamb mapping */
37 +        mult = 1.0/256.0;
38          for (i = 0; i <= MAXGSHIFT; i++) {
28 –                mult = pow(0.5, (double)(i+8));
39                  for (j = 0; j < 256; j++)
40 <                        g_bval[i][j] = 256.0 * pow((j+.5)*mult, 1.0/g);
40 >                        g_bval[i][j] = 256.0 * (*f)((j+.5)*mult, a2);
41 >                mult *= 0.5;
42          }
43 +        return(0);
44 + }
45 +
46 +
47 + setcolrinv(f, a2)               /* set inverse brightness correction */
48 + double  (*f)();
49 + double  a2;
50 + {
51 +        double  mult;
52 +        register int    i, j;
53 +                                        /* allocate tables */
54 +        if (g_mant == NULL && (g_mant = (BYTE *)bmalloc(256)) == NULL)
55 +                return(-1);
56 +        if (g_nexp == NULL && (g_nexp = (BYTE *)bmalloc(256)) == NULL)
57 +                return(-1);
58                                          /* compute gamb -> colr mapping */
59          i = 0;
60          mult = 256.0;
61          for (j = 255; j > 0; j--) {
62 <                while ((g_mant[j] = mult * pow(j/256.0, g)) < 128) {
62 >                while ((g_mant[j] = mult * (*f)(j/256.0, a2)) < 128) {
63                          i++;
64                          mult *= 2.0;
65                  }
# Line 41 | Line 67 | double g;
67          }
68          g_mant[0] = 0;
69          g_nexp[0] = COLXS;
70 +        return(0);
71   }
72  
73  
74 + setcolrgam(g)                   /* set gamma conversion */
75 + double  g;
76 + {
77 +        if (setcolrcor(pow, 1.0/g) < 0)
78 +                return(-1);
79 +        return(setcolrinv(pow, g));
80 + }
81 +
82 +
83   colrs_gambs(scan, len)          /* convert scanline of colrs to gamma bytes */
84   register COLR   *scan;
85   int     len;
86   {
87          register int    i, expo;
88  
89 +        if (g_bval == NULL)
90 +                return(-1);
91          while (len-- > 0) {
92                  expo = scan[0][EXP] - COLXS;
93                  if (expo < -MAXGSHIFT) {
# Line 87 | Line 125 | int    len;
125                  scan[0][EXP] = COLXS;
126                  scan++;
127          }
128 +        return(0);
129   }
130  
131  
# Line 96 | Line 135 | int    len;
135   {
136          register int    nexpo;
137  
138 +        if (g_mant == NULL | g_nexp == NULL)
139 +                return(-1);
140          while (len-- > 0) {
141                  nexpo = g_nexp[scan[0][RED]];
142                  if (g_nexp[scan[0][GRN]] < nexpo)
# Line 120 | Line 161 | int    len;
161                  scan[0][EXP] = COLXS - nexpo;
162                  scan++;
163          }
164 +        return(0);
165   }
166  
167  
# Line 128 | Line 170 | register COLR  *scan;
170   register int    len;
171   register int    adjust;
172   {
173 +        int     minexp;
174 +
175 +        if (adjust == 0)
176 +                return;
177 +        minexp = adjust < 0 ? -adjust : 0;
178          while (len-- > 0) {
179 <                scan[0][EXP] += adjust;
179 >                if (scan[0][EXP] <= minexp)
180 >                        scan[0][RED] = scan[0][GRN] = scan[0][BLU] =
181 >                        scan[0][EXP] = 0;
182 >                else
183 >                        scan[0][EXP] += adjust;
184                  scan++;
185          }
186   }

Diff Legend

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