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.7 by greg, Wed Oct 30 12:30:32 1991 UTC vs.
Revision 2.8 by greg, Mon Mar 10 17:26:26 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   */
# 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 50 | Line 49 | dfname(il, c)                  /* return data file name */
49   struct illum_args  *il;
50   int  c;
51   {
53 –        extern char  *getpath(), *strcpy();
52          char  fname[MAXSTR];
53          register char  *s;
54  
# Line 93 | Line 91 | float  *da;
91   int  n, m;
92   FVECT  u, v, w;
93   {
94 +        float  *Ninv;
95          FILE  *dfp;
96          int  i;
97  
98 +        if ((Ninv = (float *)malloc(3*m*sizeof(float))) == NULL)
99 +                error(SYSTEM, "out of memory in flatout");
100 +        compinv(Ninv, da, m);
101          if (il->flags & IL_COLDST) {
102                  printf("\n%s %s %s%s", VOIDID, ofun[PAT_CDATA].funame,
103                                  il->matname, DSTSUF);
# Line 103 | Line 105 | FVECT  u, v, w;
105                  for (i = 0; i < 3; i++) {
106                          dfp = dfopen(il, DATORD[i]);
107                          fprintf(dfp, "2\n%f %f %d\n%f %f %d\n",
108 <                                        1.-.5/n, .5/n, n,
108 >                                        1.+.5/n, .5/n, n+1,
109                                          0., 2.*PI, m+1);
110 +                        colorout(i, Ninv, 1, m, 1./il->nsamps/il->col[i], dfp);
111                          colorout(i, da, n, m, 1./il->nsamps/il->col[i], dfp);
112 +                        fputeol(dfp);
113                          fclose(dfp);
114                          printf(" %s", dfname(il, DATORD[i]));
115                  }
116          } else {
117                  printf("\n%s %s %s%s", VOIDID, ofun[PAT_BDATA].funame,
118                                  il->matname, DSTSUF);
119 <                printf("\n5 noop");
119 >                printf("\n5 noneg");
120                  dfp = dfopen(il, 0);
121 <                fprintf(dfp, "2\n%f %f %d\n%f %f %d\n", 1.-.5/n, .5/n, n,
121 >                fprintf(dfp, "2\n%f %f %d\n%f %f %d\n",
122 >                                1.+.5/n, .5/n, n+1,
123                                  0., 2.*PI, m+1);
124 +                brightout(Ninv, 1, m, 1./il->nsamps/brt(il->col), dfp);
125                  brightout(da, n, m, 1./il->nsamps/brt(il->col), dfp);
126 +                fputeol(dfp);
127                  fclose(dfp);
128                  printf(" %s", dfname(il, 0));
129          }
# Line 126 | Line 133 | FVECT  u, v, w;
133          printf("\t%f\t%f\t%f\n", v[0], v[1], v[2]);
134          printf("\t%f\t%f\t%f\n", w[0], w[1], w[2]);
135          il->dfnum++;
136 +        free((void *)Ninv);
137   }
138  
139  
# Line 134 | Line 142 | struct illum_args  *il;
142   float  *da;
143   int  n, m;
144   {
145 +        float  *Ninv, *Sinv;
146          FILE  *dfp;
147          int  i;
148  
149 +        if ((Ninv = (float *)malloc(3*m*sizeof(float))) == NULL ||
150 +                        (Sinv = (float *)malloc(3*m*sizeof(float))) == NULL)
151 +                error(SYSTEM, "out of memory in roundout");
152 +        compinv(Ninv, da, m);
153 +        compinv(Sinv, da+3*m*(n-1), m);
154          if (il->flags & IL_COLDST) {
155                  printf("\n%s %s %s%s", VOIDID, ofun[PAT_CDATA].funame,
156                                  il->matname, DSTSUF);
# Line 144 | Line 158 | int  n, m;
158                  for (i = 0; i < 3; i++) {
159                          dfp = dfopen(il, DATORD[i]);
160                          fprintf(dfp, "2\n%f %f %d\n%f %f %d\n",
161 <                                        1.-1./n, -1.+1./n, n,
161 >                                        1.+1./n, -1.-1./n, n+2,
162                                          0., 2.*PI, m+1);
163 +                        colorout(i, Ninv, 1, m, 1./il->nsamps/il->col[i], dfp);
164                          colorout(i, da, n, m, 1./il->nsamps/il->col[i], dfp);
165 +                        colorout(i, Sinv, 1, m, 1./il->nsamps/il->col[i], dfp);
166 +                        fputeol(dfp);
167                          fclose(dfp);
168                          printf(" %s", dfname(il, DATORD[i]));
169                  }
170          } else {
171                  printf("\n%s %s %s%s", VOIDID, ofun[PAT_BDATA].funame,
172                                  il->matname, DSTSUF);
173 <                printf("\n5 noop");
173 >                printf("\n5 noneg");
174                  dfp = dfopen(il, 0);
175 <                fprintf(dfp, "2\n%f %f %d\n%f %f %d\n", 1.-1./n, -1.+1./n, n,
175 >                fprintf(dfp, "2\n%f %f %d\n%f %f %d\n",
176 >                                1.+1./n, -1.-1./n, n+2,
177                                  0., 2.*PI, m+1);
178 +                brightout(Ninv, 1, m, 1./il->nsamps/brt(il->col), dfp);
179                  brightout(da, n, m, 1./il->nsamps/brt(il->col), dfp);
180 +                brightout(Sinv, 1, m, 1./il->nsamps/brt(il->col), dfp);
181 +                fputeol(dfp);
182                  fclose(dfp);
183                  printf(" %s", dfname(il, 0));
184          }
185          printf("\n\t%s il_alt il_azi", FNCFNM);
186          printf("\n0\n0\n");
187          il->dfnum++;
188 +        free((void *)Ninv);
189 +        free((void *)Sinv);
190   }
191  
192  
# Line 196 | Line 219 | OBJREC  *ob;
219   }
220  
221  
222 < average(il, da, n)              /* compute average value for distribution */
223 < register struct illum_args  *il;
222 > compavg(col, da, n)             /* compute average for set of data values */
223 > float  col[3];
224   register float  *da;
225   int  n;
226   {
227          register int  i;
228  
229 <        il->col[0] = il->col[1] = il->col[2] = 0.;
229 >        col[0] = col[1] = col[2] = 0.;
230          i = n;
231          while (i-- > 0) {
232 <                il->col[0] += *da++;
233 <                il->col[1] += *da++;
234 <                il->col[2] += *da++;
232 >                col[0] += *da++;
233 >                col[1] += *da++;
234 >                col[2] += *da++;
235          }
236          for (i = 0; i < 3; i++)
237 <                il->col[i] /= (double)n*il->nsamps;
237 >                col[i] /= (double)n;
238 > }
239  
240 +
241 + compinv(rinv, rp, m)            /* compute other side of row average */
242 + register float  *rinv, *rp;
243 + int  m;
244 + {
245 +        float  avg[3];
246 +
247 +        compavg(avg, rp, m);            /* row average */
248 +        while (m-- > 0) {
249 +                *rinv++ = 2.*avg[0] - *rp++;
250 +                *rinv++ = 2.*avg[1] - *rp++;
251 +                *rinv++ = 2.*avg[2] - *rp++;
252 +        }
253 + }
254 +
255 +
256 + average(il, da, n)              /* evaluate average value for distribution */
257 + register struct illum_args  *il;
258 + float  *da;
259 + int  n;
260 + {
261 +        compavg(il->col, da, n);        /* average */
262 +        if (il->nsamps > 1) {
263 +                il->col[0] /= (double)il->nsamps;
264 +                il->col[1] /= (double)il->nsamps;
265 +                il->col[2] /= (double)il->nsamps;
266 +        }
267 +                                        /* brighter than minimum? */
268          return(brt(il->col) > il->minbrt+FTINY);
269   }
270  
# Line 223 | Line 275 | fputnum(d, fp)                 /* put out a number to fp */
275   double  d;
276   FILE  *fp;
277   {
278 <        if (colmcnt++ % 6 == 0)
278 >        if (colmcnt++ % 5 == 0)
279                  putc('\n', fp);
280          fprintf(fp, " %11e", d);
281   }
# Line 253 | Line 305 | FILE  *fp;
305                  }
306                  fputnum(mult*da[p-3*m], fp);    /* wrap phi */
307          }
256 –        fputeol(fp);
308   }
309  
310  
# Line 272 | Line 323 | FILE  *fp;
323                  }
324                  fputnum(mult*brt(da-3*m), fp);  /* wrap phi */
325          }
275 –        fputeol(fp);
326   }

Diff Legend

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