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"; |
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 |
> |
#ifndef pow |
24 |
> |
extern double pow(); |
25 |
> |
#endif |
26 |
> |
|
27 |
> |
|
28 |
> |
setcolrcor(f, a2) /* set brightness correction */ |
29 |
> |
double (*f)(); |
30 |
> |
double a2; |
31 |
|
{ |
23 |
– |
extern double pow(); |
32 |
|
double mult; |
33 |
|
register int i, j; |
34 |
+ |
/* allocate tables */ |
35 |
+ |
if (g_bval == NULL && (g_bval = |
36 |
+ |
(BYTE (*)[256])bmalloc((MAXGSHIFT+1)*256)) == NULL) |
37 |
+ |
return(-1); |
38 |
|
/* compute colr -> gamb mapping */ |
39 |
+ |
mult = 1.0/256.0; |
40 |
|
for (i = 0; i <= MAXGSHIFT; i++) { |
28 |
– |
mult = pow(0.5, (double)(i+8)); |
41 |
|
for (j = 0; j < 256; j++) |
42 |
< |
g_bval[i][j] = 256.0 * pow((j+.5)*mult, 1.0/g); |
42 |
> |
g_bval[i][j] = 256.0 * (*f)((j+.5)*mult, a2); |
43 |
> |
mult *= 0.5; |
44 |
|
} |
45 |
+ |
return(0); |
46 |
+ |
} |
47 |
+ |
|
48 |
+ |
|
49 |
+ |
setcolrinv(f, a2) /* set inverse brightness correction */ |
50 |
+ |
double (*f)(); |
51 |
+ |
double a2; |
52 |
+ |
{ |
53 |
+ |
double mult; |
54 |
+ |
register int i, j; |
55 |
+ |
/* allocate tables */ |
56 |
+ |
if (g_mant == NULL && (g_mant = (BYTE *)bmalloc(256)) == NULL) |
57 |
+ |
return(-1); |
58 |
+ |
if (g_nexp == NULL && (g_nexp = (BYTE *)bmalloc(256)) == NULL) |
59 |
+ |
return(-1); |
60 |
|
/* compute gamb -> colr mapping */ |
61 |
|
i = 0; |
62 |
|
mult = 256.0; |
63 |
|
for (j = 255; j > 0; j--) { |
64 |
< |
while ((g_mant[j] = mult * pow(j/256.0, g)) < 128) { |
64 |
> |
while ((g_mant[j] = mult * (*f)(j/256.0, a2)) < 128) { |
65 |
|
i++; |
66 |
|
mult *= 2.0; |
67 |
|
} |
69 |
|
} |
70 |
|
g_mant[0] = 0; |
71 |
|
g_nexp[0] = COLXS; |
72 |
+ |
return(0); |
73 |
|
} |
74 |
|
|
75 |
|
|
76 |
+ |
setcolrgam(g) /* set gamma conversion */ |
77 |
+ |
double g; |
78 |
+ |
{ |
79 |
+ |
if (setcolrcor(pow, 1.0/g) < 0) |
80 |
+ |
return(-1); |
81 |
+ |
return(setcolrinv(pow, g)); |
82 |
+ |
} |
83 |
+ |
|
84 |
+ |
|
85 |
|
colrs_gambs(scan, len) /* convert scanline of colrs to gamma bytes */ |
86 |
|
register COLR *scan; |
87 |
|
int len; |
88 |
|
{ |
89 |
|
register int i, expo; |
90 |
|
|
91 |
+ |
if (g_bval == NULL) |
92 |
+ |
return(-1); |
93 |
|
while (len-- > 0) { |
94 |
|
expo = scan[0][EXP] - COLXS; |
95 |
|
if (expo < -MAXGSHIFT) { |
127 |
|
scan[0][EXP] = COLXS; |
128 |
|
scan++; |
129 |
|
} |
130 |
+ |
return(0); |
131 |
|
} |
132 |
|
|
133 |
|
|
137 |
|
{ |
138 |
|
register int nexpo; |
139 |
|
|
140 |
+ |
if (g_mant == NULL | g_nexp == NULL) |
141 |
+ |
return(-1); |
142 |
|
while (len-- > 0) { |
143 |
|
nexpo = g_nexp[scan[0][RED]]; |
144 |
|
if (g_nexp[scan[0][GRN]] < nexpo) |
163 |
|
scan[0][EXP] = COLXS - nexpo; |
164 |
|
scan++; |
165 |
|
} |
166 |
+ |
return(0); |
167 |
|
} |
168 |
|
|
169 |
|
|
172 |
|
register int len; |
173 |
|
register int adjust; |
174 |
|
{ |
175 |
+ |
int minexp; |
176 |
+ |
|
177 |
+ |
if (adjust == 0) |
178 |
+ |
return; |
179 |
+ |
minexp = adjust < 0 ? -adjust : 0; |
180 |
|
while (len-- > 0) { |
181 |
< |
scan[0][EXP] += adjust; |
181 |
> |
if (scan[0][EXP] <= minexp) |
182 |
> |
scan[0][RED] = scan[0][GRN] = scan[0][BLU] = |
183 |
> |
scan[0][EXP] = 0; |
184 |
> |
else |
185 |
> |
scan[0][EXP] += adjust; |
186 |
|
scan++; |
187 |
|
} |
188 |
|
} |