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

Comparing ray/src/gen/mkillum3.c (file contents):
Revision 1.6 by greg, Mon Aug 26 10:16:59 1991 UTC vs.
Revision 2.7 by greg, Sat Feb 22 02:07:24 2003 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1991 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   * Routines to print mkillum objects
6   */
7  
8   #include  "mkillum.h"
9  
10 < #define  brt(col)       (.295*(col)[0] + .635*(col)[1] + .070*(col)[2])
10 > #define  brt(col)       (.263*(col)[0]+.655*(col)[1]+.082*(col)[2])
11  
12   char    DATORD[] = "RGB";               /* data ordering */
13   char    DATSUF[] = ".dat";              /* data file suffix */
# Line 24 | Line 21 | register OBJREC  *obj;
21   {
22          register int  i;
23  
24 +        if (issurface(obj->otype) && !strcmp(mod, VOIDID))
25 +                return;         /* don't print void surfaces */
26          printf("\n%s %s %s", mod, ofun[obj->otype].funame, obj->oname);
27          printf("\n%d", obj->oargs.nsargs);
28          for (i = 0; i < obj->oargs.nsargs; i++)
# Line 93 | Line 92 | float  *da;
92   int  n, m;
93   FVECT  u, v, w;
94   {
95 +        float  *Ninv;
96          FILE  *dfp;
97          int  i;
98  
99 +        if ((Ninv = (float *)malloc(3*m*sizeof(float))) == NULL)
100 +                error(SYSTEM, "out of memory in flatout");
101 +        compinv(Ninv, da, m);
102          if (il->flags & IL_COLDST) {
103                  printf("\n%s %s %s%s", VOIDID, ofun[PAT_CDATA].funame,
104                                  il->matname, DSTSUF);
# Line 103 | Line 106 | FVECT  u, v, w;
106                  for (i = 0; i < 3; i++) {
107                          dfp = dfopen(il, DATORD[i]);
108                          fprintf(dfp, "2\n%f %f %d\n%f %f %d\n",
109 <                                        1.-.5/n, .5/n, n,
109 >                                        1.+.5/n, .5/n, n+1,
110                                          0., 2.*PI, m+1);
111 +                        colorout(i, Ninv, 1, m, 1./il->nsamps/il->col[i], dfp);
112                          colorout(i, da, n, m, 1./il->nsamps/il->col[i], dfp);
113 +                        fputeol(dfp);
114                          fclose(dfp);
115                          printf(" %s", dfname(il, DATORD[i]));
116                  }
117          } else {
118                  printf("\n%s %s %s%s", VOIDID, ofun[PAT_BDATA].funame,
119                                  il->matname, DSTSUF);
120 <                printf("\n5 noop");
120 >                printf("\n5 noneg");
121                  dfp = dfopen(il, 0);
122 <                fprintf(dfp, "2\n%f %f %d\n%f %f %d\n", 1.-.5/n, .5/n, n,
122 >                fprintf(dfp, "2\n%f %f %d\n%f %f %d\n",
123 >                                1.+.5/n, .5/n, n+1,
124                                  0., 2.*PI, m+1);
125 +                brightout(Ninv, 1, m, 1./il->nsamps/brt(il->col), dfp);
126                  brightout(da, n, m, 1./il->nsamps/brt(il->col), dfp);
127 +                fputeol(dfp);
128                  fclose(dfp);
129                  printf(" %s", dfname(il, 0));
130          }
# Line 126 | Line 134 | FVECT  u, v, w;
134          printf("\t%f\t%f\t%f\n", v[0], v[1], v[2]);
135          printf("\t%f\t%f\t%f\n", w[0], w[1], w[2]);
136          il->dfnum++;
137 +        free((void *)Ninv);
138   }
139  
140  
# Line 134 | Line 143 | struct illum_args  *il;
143   float  *da;
144   int  n, m;
145   {
146 +        float  *Ninv, *Sinv;
147          FILE  *dfp;
148          int  i;
149  
150 +        if ((Ninv = (float *)malloc(3*m*sizeof(float))) == NULL ||
151 +                        (Sinv = (float *)malloc(3*m*sizeof(float))) == NULL)
152 +                error(SYSTEM, "out of memory in roundout");
153 +        compinv(Ninv, da, m);
154 +        compinv(Sinv, da+3*m*(n-1), m);
155          if (il->flags & IL_COLDST) {
156                  printf("\n%s %s %s%s", VOIDID, ofun[PAT_CDATA].funame,
157                                  il->matname, DSTSUF);
# Line 144 | Line 159 | int  n, m;
159                  for (i = 0; i < 3; i++) {
160                          dfp = dfopen(il, DATORD[i]);
161                          fprintf(dfp, "2\n%f %f %d\n%f %f %d\n",
162 <                                        1.-1./n, -1.+1./n, n,
162 >                                        1.+1./n, -1.-1./n, n+2,
163                                          0., 2.*PI, m+1);
164 +                        colorout(i, Ninv, 1, m, 1./il->nsamps/il->col[i], dfp);
165                          colorout(i, da, n, m, 1./il->nsamps/il->col[i], dfp);
166 +                        colorout(i, Sinv, 1, m, 1./il->nsamps/il->col[i], dfp);
167 +                        fputeol(dfp);
168                          fclose(dfp);
169                          printf(" %s", dfname(il, DATORD[i]));
170                  }
171          } else {
172                  printf("\n%s %s %s%s", VOIDID, ofun[PAT_BDATA].funame,
173                                  il->matname, DSTSUF);
174 <                printf("\n5 noop");
174 >                printf("\n5 noneg");
175                  dfp = dfopen(il, 0);
176 <                fprintf(dfp, "2\n%f %f %d\n%f %f %d\n", 1.-1./n, -1.+1./n, n,
176 >                fprintf(dfp, "2\n%f %f %d\n%f %f %d\n",
177 >                                1.+1./n, -1.-1./n, n+2,
178                                  0., 2.*PI, m+1);
179 +                brightout(Ninv, 1, m, 1./il->nsamps/brt(il->col), dfp);
180                  brightout(da, n, m, 1./il->nsamps/brt(il->col), dfp);
181 +                brightout(Sinv, 1, m, 1./il->nsamps/brt(il->col), dfp);
182 +                fputeol(dfp);
183                  fclose(dfp);
184                  printf(" %s", dfname(il, 0));
185          }
186          printf("\n\t%s il_alt il_azi", FNCFNM);
187          printf("\n0\n0\n");
188          il->dfnum++;
189 +        free((void *)Ninv);
190 +        free((void *)Sinv);
191   }
192  
193  
# Line 196 | Line 220 | OBJREC  *ob;
220   }
221  
222  
223 < average(il, da, n)              /* compute average value for distribution */
224 < register struct illum_args  *il;
223 > compavg(col, da, n)             /* compute average for set of data values */
224 > float  col[3];
225   register float  *da;
226   int  n;
227   {
228          register int  i;
229  
230 <        il->col[0] = il->col[1] = il->col[2] = 0.;
230 >        col[0] = col[1] = col[2] = 0.;
231          i = n;
232          while (i-- > 0) {
233 <                il->col[0] += *da++;
234 <                il->col[1] += *da++;
235 <                il->col[2] += *da++;
233 >                col[0] += *da++;
234 >                col[1] += *da++;
235 >                col[2] += *da++;
236          }
237          for (i = 0; i < 3; i++)
238 <                il->col[i] /= (double)n*il->nsamps;
238 >                col[i] /= (double)n;
239 > }
240  
241 +
242 + compinv(rinv, rp, m)            /* compute other side of row average */
243 + register float  *rinv, *rp;
244 + int  m;
245 + {
246 +        float  avg[3];
247 +
248 +        compavg(avg, rp, m);            /* row average */
249 +        while (m-- > 0) {
250 +                *rinv++ = 2.*avg[0] - *rp++;
251 +                *rinv++ = 2.*avg[1] - *rp++;
252 +                *rinv++ = 2.*avg[2] - *rp++;
253 +        }
254 + }
255 +
256 +
257 + average(il, da, n)              /* evaluate average value for distribution */
258 + register struct illum_args  *il;
259 + float  *da;
260 + int  n;
261 + {
262 +        compavg(il->col, da, n);        /* average */
263 +        if (il->nsamps > 1) {
264 +                il->col[0] /= (double)il->nsamps;
265 +                il->col[1] /= (double)il->nsamps;
266 +                il->col[2] /= (double)il->nsamps;
267 +        }
268 +                                        /* brighter than minimum? */
269          return(brt(il->col) > il->minbrt+FTINY);
270   }
271  
# Line 223 | Line 276 | fputnum(d, fp)                 /* put out a number to fp */
276   double  d;
277   FILE  *fp;
278   {
279 <        if (colmcnt++ % 6 == 0)
279 >        if (colmcnt++ % 5 == 0)
280                  putc('\n', fp);
281          fprintf(fp, " %11e", d);
282   }
# Line 253 | Line 306 | FILE  *fp;
306                  }
307                  fputnum(mult*da[p-3*m], fp);    /* wrap phi */
308          }
256 –        fputeol(fp);
309   }
310  
311  
# Line 272 | Line 324 | FILE  *fp;
324                  }
325                  fputnum(mult*brt(da-3*m), fp);  /* wrap phi */
326          }
275 –        fputeol(fp);
327   }

Diff Legend

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