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.8 by greg, Tue Nov 12 10:52:37 1991 UTC vs.
Revision 2.6 by greg, Wed Aug 12 14:24:08 1992 UTC

# Line 24 | Line 24 | register OBJREC  *obj;
24   {
25          register int  i;
26  
27 +        if (issurface(obj->otype) && !strcmp(mod, VOIDID))
28 +                return;         /* don't print void surfaces */
29          printf("\n%s %s %s", mod, ofun[obj->otype].funame, obj->oname);
30          printf("\n%d", obj->oargs.nsargs);
31          for (i = 0; i < obj->oargs.nsargs; i++)
# Line 93 | Line 95 | float  *da;
95   int  n, m;
96   FVECT  u, v, w;
97   {
98 +        float  *Ninv;
99          FILE  *dfp;
100          int  i;
101  
102 +        if ((Ninv = (float *)malloc(3*m*sizeof(float))) == NULL)
103 +                error(SYSTEM, "out of memory in flatout");
104 +        compinv(Ninv, da, m);
105          if (il->flags & IL_COLDST) {
106                  printf("\n%s %s %s%s", VOIDID, ofun[PAT_CDATA].funame,
107                                  il->matname, DSTSUF);
# Line 103 | Line 109 | FVECT  u, v, w;
109                  for (i = 0; i < 3; i++) {
110                          dfp = dfopen(il, DATORD[i]);
111                          fprintf(dfp, "2\n%f %f %d\n%f %f %d\n",
112 <                                        1.-.5/n, .5/n, n,
112 >                                        1.+.5/n, .5/n, n+1,
113                                          0., 2.*PI, m+1);
114 +                        colorout(i, Ninv, 1, m, 1./il->nsamps/il->col[i], dfp);
115                          colorout(i, da, n, m, 1./il->nsamps/il->col[i], dfp);
116 +                        fputeol(dfp);
117                          fclose(dfp);
118                          printf(" %s", dfname(il, DATORD[i]));
119                  }
# Line 114 | Line 122 | FVECT  u, v, w;
122                                  il->matname, DSTSUF);
123                  printf("\n5 noneg");
124                  dfp = dfopen(il, 0);
125 <                fprintf(dfp, "2\n%f %f %d\n%f %f %d\n", 1.-.5/n, .5/n, n,
125 >                fprintf(dfp, "2\n%f %f %d\n%f %f %d\n",
126 >                                1.+.5/n, .5/n, n+1,
127                                  0., 2.*PI, m+1);
128 +                brightout(Ninv, 1, m, 1./il->nsamps/brt(il->col), dfp);
129                  brightout(da, n, m, 1./il->nsamps/brt(il->col), dfp);
130 +                fputeol(dfp);
131                  fclose(dfp);
132                  printf(" %s", dfname(il, 0));
133          }
# Line 126 | Line 137 | FVECT  u, v, w;
137          printf("\t%f\t%f\t%f\n", v[0], v[1], v[2]);
138          printf("\t%f\t%f\t%f\n", w[0], w[1], w[2]);
139          il->dfnum++;
140 +        free((char *)Ninv);
141   }
142  
143  
# Line 134 | Line 146 | struct illum_args  *il;
146   float  *da;
147   int  n, m;
148   {
149 +        float  *Ninv, *Sinv;
150          FILE  *dfp;
151          int  i;
152  
153 +        if ((Ninv = (float *)malloc(3*m*sizeof(float))) == NULL ||
154 +                        (Sinv = (float *)malloc(3*m*sizeof(float))) == NULL)
155 +                error(SYSTEM, "out of memory in roundout");
156 +        compinv(Ninv, da, m);
157 +        compinv(Sinv, da+3*m*(n-1), m);
158          if (il->flags & IL_COLDST) {
159                  printf("\n%s %s %s%s", VOIDID, ofun[PAT_CDATA].funame,
160                                  il->matname, DSTSUF);
# Line 144 | Line 162 | int  n, m;
162                  for (i = 0; i < 3; i++) {
163                          dfp = dfopen(il, DATORD[i]);
164                          fprintf(dfp, "2\n%f %f %d\n%f %f %d\n",
165 <                                        1.-1./n, -1.+1./n, n,
165 >                                        1.+1./n, -1.-1./n, n+2,
166                                          0., 2.*PI, m+1);
167 +                        colorout(i, Ninv, 1, m, 1./il->nsamps/il->col[i], dfp);
168                          colorout(i, da, n, m, 1./il->nsamps/il->col[i], dfp);
169 +                        colorout(i, Sinv, 1, m, 1./il->nsamps/il->col[i], dfp);
170 +                        fputeol(dfp);
171                          fclose(dfp);
172                          printf(" %s", dfname(il, DATORD[i]));
173                  }
# Line 155 | Line 176 | int  n, m;
176                                  il->matname, DSTSUF);
177                  printf("\n5 noneg");
178                  dfp = dfopen(il, 0);
179 <                fprintf(dfp, "2\n%f %f %d\n%f %f %d\n", 1.-1./n, -1.+1./n, n,
179 >                fprintf(dfp, "2\n%f %f %d\n%f %f %d\n",
180 >                                1.+1./n, -1.-1./n, n+2,
181                                  0., 2.*PI, m+1);
182 +                brightout(Ninv, 1, m, 1./il->nsamps/brt(il->col), dfp);
183                  brightout(da, n, m, 1./il->nsamps/brt(il->col), dfp);
184 +                brightout(Sinv, 1, m, 1./il->nsamps/brt(il->col), dfp);
185 +                fputeol(dfp);
186                  fclose(dfp);
187                  printf(" %s", dfname(il, 0));
188          }
189          printf("\n\t%s il_alt il_azi", FNCFNM);
190          printf("\n0\n0\n");
191          il->dfnum++;
192 +        free((char *)Ninv);
193 +        free((char *)Sinv);
194   }
195  
196  
# Line 196 | Line 223 | OBJREC  *ob;
223   }
224  
225  
226 < average(il, da, n)              /* compute average value for distribution */
227 < register struct illum_args  *il;
226 > compavg(col, da, n)             /* compute average for set of data values */
227 > float  col[3];
228   register float  *da;
229   int  n;
230   {
231          register int  i;
232  
233 <        il->col[0] = il->col[1] = il->col[2] = 0.;
233 >        col[0] = col[1] = col[2] = 0.;
234          i = n;
235          while (i-- > 0) {
236 <                il->col[0] += *da++;
237 <                il->col[1] += *da++;
238 <                il->col[2] += *da++;
236 >                col[0] += *da++;
237 >                col[1] += *da++;
238 >                col[2] += *da++;
239          }
240          for (i = 0; i < 3; i++)
241 <                il->col[i] /= (double)n*il->nsamps;
241 >                col[i] /= (double)n;
242 > }
243  
244 +
245 + compinv(rinv, rp, m)            /* compute other side of row average */
246 + register float  *rinv, *rp;
247 + int  m;
248 + {
249 +        float  avg[3];
250 +
251 +        compavg(avg, rp, m);            /* row average */
252 +        while (m-- > 0) {
253 +                *rinv++ = 2.*avg[0] - *rp++;
254 +                *rinv++ = 2.*avg[1] - *rp++;
255 +                *rinv++ = 2.*avg[2] - *rp++;
256 +        }
257 + }
258 +
259 +
260 + average(il, da, n)              /* evaluate average value for distribution */
261 + register struct illum_args  *il;
262 + float  *da;
263 + int  n;
264 + {
265 +        compavg(il->col, da, n);        /* average */
266 +        if (il->nsamps > 1) {
267 +                il->col[0] /= (double)il->nsamps;
268 +                il->col[1] /= (double)il->nsamps;
269 +                il->col[2] /= (double)il->nsamps;
270 +        }
271 +                                        /* brighter than minimum? */
272          return(brt(il->col) > il->minbrt+FTINY);
273   }
274  
# Line 223 | Line 279 | fputnum(d, fp)                 /* put out a number to fp */
279   double  d;
280   FILE  *fp;
281   {
282 <        if (colmcnt++ % 6 == 0)
282 >        if (colmcnt++ % 5 == 0)
283                  putc('\n', fp);
284          fprintf(fp, " %11e", d);
285   }
# Line 253 | Line 309 | FILE  *fp;
309                  }
310                  fputnum(mult*da[p-3*m], fp);    /* wrap phi */
311          }
256 –        fputeol(fp);
312   }
313  
314  
# Line 272 | Line 327 | FILE  *fp;
327                  }
328                  fputnum(mult*brt(da-3*m), fp);  /* wrap phi */
329          }
275 –        fputeol(fp);
330   }

Diff Legend

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