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

Comparing ray/src/common/tmapcolrs.c (file contents):
Revision 3.9 by gwlarson, Wed Oct 28 09:26:03 1998 UTC vs.
Revision 3.15 by schorsch, Fri Jan 2 11:36:26 2004 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1998 Silicon Graphics, Inc. */
2 –
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ SGI";
2 > static const char       RCSid[] = "$Id$";
3   #endif
6 –
4   /*
5   * Routines for tone mapping on Radiance RGBE and XYZE pictures.
6 < * See tonemap.h for detailed function descriptions.
6 > *
7 > * Externals declared in tonemap.h
8   */
9  
10 + #include "copyright.h"
11 +
12   #include        <stdio.h>
13 + #include        <string.h>
14   #include        <math.h>
15 + #include        <time.h>
16   #include        "tmprivat.h"
17   #include        "resolu.h"
18  
19 + #ifndef TM_PIC_CTRANS
20 + #define TM_PIC_CTRANS   1               /* transform colors? (expensive) */
21 + #endif
22  
18 – extern char     *tempbuffer();
19 –
23   #define GAMTSZ  1024
24  
25   typedef struct {
# Line 25 | Line 28 | typedef struct {
28          TMbright        inpsfb;         /* encoded tm->inpsf */
29   } COLRDATA;
30  
31 < static MEM_PTR  colrInit();
32 < static void     colrNewSpace();
33 < extern void     free();
31 > static MEM_PTR  colrInit(struct tmStruct *);
32 > static void     colrNewSpace(struct tmStruct *);
33 > static gethfunc headline;
34 >
35   static struct tmPackage colrPkg = {     /* our package functions */
36          colrInit, colrNewSpace, free
37   };
# Line 35 | Line 39 | static int     colrReg = -1;           /* our package registration
39  
40   #define LOGISZ  260
41   static TMbright logi[LOGISZ];
38 – static BYTE     photofact[BMESUPPER-BMESLOWER];
42  
43  
44   int
# Line 52 | Line 55 | int    len;
55  
56          if (tmTop == NULL)
57                  returnErr(TM_E_TMINVAL);
58 <        if (ls == NULL | scan == NULL | len < 0)
58 >        if ((ls == NULL) | (scan == NULL) | (len < 0))
59                  returnErr(TM_E_ILLEGAL);
60 + #if TM_PIC_CTRANS
61          if (tmNeedMatrix(tmTop)) {              /* need floating point */
62 + #else
63 +        if (tmTop->inppri == TM_XYZPRIM) {      /* no way around this */
64 + #endif
65                  register COLOR  *newscan;
66                  newscan = (COLOR *)tempbuffer(len*sizeof(COLOR));
67                  if (newscan == NULL)
# Line 70 | Line 77 | int    len;
77                  for (i = 256; i--; )
78                          logi[i] = TM_BRTSCALE*log((i+.5)/256.) - .5;
79                  for (i = 256; i < LOGISZ; i++)
80 <                        logi[i] = logi[255];
81 <                for (i = BMESLOWER; i < BMESUPPER; i++)
75 <                        photofact[i-BMESLOWER] = 256. *
76 <                                        (tmLuminance(i) - LMESLOWER) /
77 <                                        (LMESUPPER - LMESLOWER);
80 >                        logi[i] = 0;
81 >                tmMkMesofact();
82          }
83          if ((cd = (COLRDATA *)tmPkgData(tmTop,colrReg)) == NULL)
84                  returnErr(TM_E_NOMEM);
# Line 84 | Line 88 | int    len;
88                  li =  ( cd->clfb[RED]*cmon[RED] +
89                          cd->clfb[GRN]*cmon[GRN] +
90                          cd->clfb[BLU]*cmon[BLU] ) >> 8;
91 <                bi = BRT2SCALE*(cmon[EXP]-COLXS) +
91 >                bi = BRT2SCALE(cmon[EXP]-COLXS) +
92                                  logi[li] + cd->inpsfb;
93 <                if (bi < MINBRT) {
94 <                        bi = MINBRT-1;                  /* bogus value */
95 <                        li++;                           /* avoid li==0 */
93 >                if (li <= 0) {
94 >                        bi = TM_NOBRT;                  /* bogus value */
95 >                        li = 1;                         /* avoid li==0 */
96                  }
97                  ls[i] = bi;
98                  if (cs == TM_NOCHROM)                   /* no color? */
# Line 101 | Line 105 | int    len;
105                          else {
106                                  if (tmTop->flags & TM_F_BW)
107                                          cmon[RED] = cmon[GRN] = cmon[BLU] = li;
108 <                                pf = photofact[bi-BMESLOWER];
108 >                                pf = tmMesofact[bi-BMESLOWER];
109                                  sli *= 256 - pf;
110                                  cmon[RED] = ( sli + pf*cmon[RED] ) >> 8;
111                                  cmon[GRN] = ( sli + pf*cmon[GRN] ) >> 8;
# Line 110 | Line 114 | int    len;
114                  } else if (tmTop->flags & TM_F_BW) {
115                          cmon[RED] = cmon[GRN] = cmon[BLU] = li;
116                  }
117 <                bi = ( (int4)GAMTSZ*cd->clfb[RED]*cmon[RED]/li ) >> 8;
117 >                bi = ( (int32)GAMTSZ*cd->clfb[RED]*cmon[RED]/li ) >> 8;
118                  cs[3*i  ] = bi>=GAMTSZ ? 255 : cd->gamb[bi];
119 <                bi = ( (int4)GAMTSZ*cd->clfb[GRN]*cmon[GRN]/li ) >> 8;
119 >                bi = ( (int32)GAMTSZ*cd->clfb[GRN]*cmon[GRN]/li ) >> 8;
120                  cs[3*i+1] = bi>=GAMTSZ ? 255 : cd->gamb[bi];
121 <                bi = ( (int4)GAMTSZ*cd->clfb[BLU]*cmon[BLU]/li ) >> 8;
121 >                bi = ( (int32)GAMTSZ*cd->clfb[BLU]*cmon[BLU]/li ) >> 8;
122                  cs[3*i+2] = bi>=GAMTSZ ? 255 : cd->gamb[bi];
123          }
124          returnOK;
# Line 135 | Line 139 | static struct radhead {
139  
140  
141   static int
142 < headline(s, rh)                 /* grok a header line */
143 < register char   *s;
144 < register struct radhead *rh;
142 > headline(                       /* grok a header line */
143 >        register char   *s,
144 >        void    *vrh
145 > )
146   {
147          char    fmt[32];
148 +        register struct radhead *rh = vrh;
149  
150          if (formatval(fmt, s)) {
151                  if (!strcmp(fmt, COLRFMT))
# Line 180 | Line 186 | FILE   *fp;
186                                                  /* check arguments */
187          if (tmTop == NULL)
188                  returnErr(TM_E_TMINVAL);
189 <        if (lpp == NULL | xp == NULL | yp == NULL |
190 <                        (fname == NULL & fp == TM_GETFILE))
189 >        if ((lpp == NULL) | (xp == NULL) | (yp == NULL) |
190 >                        ((fname == NULL) & (fp == TM_GETFILE)))
191                  returnErr(TM_E_ILLEGAL);
192          *xp = *yp = 0;                          /* error precaution */
193          if ((inpf = fp) == TM_GETFILE && (inpf = fopen(fname, "r")) == NULL)
# Line 189 | Line 195 | FILE   *fp;
195          *lpp = NULL;
196          if (cpp != TM_NOCHROMP) *cpp = NULL;
197          info = rhdefault;                       /* get our header */
198 <        getheader(inpf, headline, (MEM_PTR)&info);
199 <        if (info.format == FMTBAD | info.expos <= 0. ||
198 >        getheader(inpf, headline, &info);
199 >        if ((info.format == FMTBAD) | (info.expos <= 0.) ||
200                          fgetresolu(xp, yp, inpf) < 0) {
201                  err = TM_E_BADFILE; goto done;
202          }
# Line 243 | Line 249 | done:                                          /* clean up */
249  
250  
251   #ifdef PCOND
252 < int                                     /* run pcond to map picture */
252 > static int                                      /* run pcond to map picture */
253   dopcond(psp, xp, yp, flags, monpri, gamval, Lddyn, Ldmax, fname)
254   BYTE    **psp;
255   int     *xp, *yp;
# Line 253 | Line 259 | double gamval, Lddyn, Ldmax;
259   char    *fname;
260   {
261          char    *funcName = fname;
262 <        char    cmdbuf[512];
262 >        char    cmdbuf[1024];
263          FILE    *infp;
264          register COLR   *scan;
265          register BYTE   *rp;
# Line 298 | Line 304 | char   *fname;
304                  rp = (BYTE *)malloc(sizeof(BYTE) * *xp * *yp);
305          else
306                  rp = (BYTE *)malloc(3*sizeof(BYTE) * *xp * *yp);
307 <        if ((*psp = rp) == NULL | scan == NULL) {
307 >        if (((*psp = rp) == NULL) | (scan == NULL)) {
308                  pclose(infp);
309                  returnErr(TM_E_NOMEM);
310          }
# Line 344 | Line 350 | FILE   *fp;
350          TMbright        *lp;
351          int     err;
352                                                  /* check arguments */
353 <        if (psp == NULL | xp == NULL | yp == NULL | monpri == NULL |
354 <                        (fname == NULL & fp == TM_GETFILE))
353 >        if ((psp == NULL) | (xp == NULL) | (yp == NULL) | (monpri == NULL) |
354 >                        ((fname == NULL) & (fp == TM_GETFILE)))
355                  returnErr(TM_E_ILLEGAL);
356                                                  /* set defaults */
357          if (gamval < MINGAM) gamval = DEFGAM;

Diff Legend

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