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

Comparing ray/src/rt/rview.c (file contents):
Revision 2.28 by greg, Mon Jun 13 20:07:56 2005 UTC vs.
Revision 2.29 by greg, Thu Aug 21 07:05:59 2008 UTC

# Line 15 | Line 15 | static const char      RCSid[] = "$Id$";
15   #include  "ray.h"
16   #include  "rpaint.h"
17  
18 –
19 – CUBE  thescene;                         /* our scene */
20 – OBJECT  nsceneobjs;                     /* number of objects in our scene */
21 –
22 – int  dimlist[MAXDIM];                   /* sampling dimensions */
23 – int  ndims = 0;                         /* number of sampling dimensions */
24 – int  samplendx = 0;                     /* index for this sample */
25 –
26 – extern void  ambnotify();
27 – void  (*addobjnotify[])() = {ambnotify, NULL};
28 –
29 – VIEW  ourview = STDVIEW;                /* viewing parameters */
30 – int  hresolu, vresolu;                  /* image resolution */
31 –
32 – void  (*trace)() = NULL;                /* trace call */
33 –
34 – int  do_irrad = 0;                      /* compute irradiance? */
35 –
36 – int  rand_samp = 0;                     /* pure Monte Carlo sampling? */
37 –
38 – int  psample = 8;                       /* pixel sample size */
39 – double  maxdiff = .15;                  /* max. sample difference */
40 –
41 – double  exposure = 1.0;                 /* exposure for scene */
42 –
43 – double  dstrsrc = 0.0;                  /* square source distribution */
44 – double  shadthresh = .1;                /* shadow threshold */
45 – double  shadcert = .25;                 /* shadow certainty */
46 – int  directrelay = 0;                   /* number of source relays */
47 – int  vspretest = 128;                   /* virtual source pretest density */
48 – int  directvis = 1;                     /* sources visible? */
49 – double  srcsizerat = 0.;                /* maximum ratio source size/dist. */
50 –
51 – COLOR  cextinction = BLKCOLOR;          /* global extinction coefficient */
52 – COLOR  salbedo = BLKCOLOR;              /* global scattering albedo */
53 – double  seccg = 0.;                     /* global scattering eccentricity */
54 – double  ssampdist = 0.;                 /* scatter sampling distance */
55 –
56 – double  specthresh = .3;                /* specular sampling threshold */
57 – double  specjitter = 1.;                /* specular sampling jitter */
58 –
59 – int  backvis = 1;                       /* back face visibility */
60 –
61 – int  maxdepth = 6;                      /* maximum recursion depth */
62 – double  minweight = 1e-2;               /* minimum ray weight */
63 –
64 – char  *ambfile = NULL;                  /* ambient file name */
65 – COLOR  ambval = BLKCOLOR;               /* ambient value */
66 – int  ambvwt = 0;                        /* initial weight for ambient value */
67 – double  ambacc = 0.3;                   /* ambient accuracy */
68 – int  ambres = 32;                       /* ambient resolution */
69 – int  ambdiv = 256;                      /* ambient divisions */
70 – int  ambssamp = 64;                     /* ambient super-samples */
71 – int  ambounce = 0;                      /* ambient bounces */
72 – char  *amblist[AMBLLEN];                /* ambient include/exclude list */
73 – int  ambincl = -1;                      /* include == 1, exclude == 0 */
74 –
75 – int  greyscale = 0;                     /* map colors to brightness? */
76 – char  *dvcname = dev_default;           /* output device name */
77 –
78 – struct driver  *dev = NULL;             /* driver functions */
79 –
80 – char  rifname[128];                     /* rad input file name */
81 –
82 – VIEW  oldview;                          /* previous view parameters */
83 –
84 – PNODE  ptrunk;                          /* the base of our image */
85 – RECT  pframe;                           /* current frame boundaries */
86 – int  pdepth;                            /* image depth in current frame */
87 –
88 – static char  *reserve_mem = NULL;       /* pre-allocated reserve memory */
89 –
90 – #define RESERVE_AMT     32768           /* amount of memory to reserve */
91 –
18   #define  CTRL(c)        ((c)-'@')
19  
20  
21   void
22 < quit(code)                      /* quit program */
23 < int  code;
22 > devopen(                                /* open device driver */
23 >        char  *dname
24 > )
25   {
99 – #ifdef MSTATS
100 –        if (code == 2 && errno == ENOMEM)
101 –                printmemstats(stderr);
102 – #endif
103 –        devclose();
104 –        exit(code);
105 – }
106 –
107 –
108 – void
109 – devopen(dname)                          /* open device driver */
110 – char  *dname;
111 – {
26          extern char  *progname, *octname;
27          char  *id;
28 <        register int  i;
28 >        int  i;
29  
30          id = octname!=NULL ? octname : progname;
31                                                  /* check device table */
# Line 132 | Line 46 | char  *dname;
46  
47  
48   void
49 < devclose()                              /* close our device */
49 > devclose(void)                          /* close our device */
50   {
51          if (dev != NULL)
52                  (*dev->close)();
# Line 141 | Line 55 | devclose()                             /* close our device */
55  
56  
57   void
58 < printdevices()                          /* print list of output devices */
58 > printdevices(void)                      /* print list of output devices */
59   {
60 <        register int  i;
60 >        int  i;
61  
62          for (i = 0; devtable[i].name; i++)
63                  printf("%-16s # %s\n", devtable[i].name, devtable[i].descrip);
# Line 151 | Line 65 | printdevices()                         /* print list of output devices */
65  
66  
67   void
68 < rview()                         /* do a view */
68 > rview(void)                             /* do a view */
69   {
70          char  buf[32];
71  
72          devopen(dvcname);               /* open device */
73 <        newimage();                     /* start image (calls fillreserves) */
73 >        newimage(NULL);                 /* start image */
74  
75          for ( ; ; ) {                   /* quit in command() */
76                  while (hresolu <= 1<<pdepth && vresolu <= 1<<pdepth)
77                          command("done: ");
164 –                while (reserve_mem == NULL)
165 –                        command("out of memory: ");
78                  errno = 0;
79                  if (hresolu <= psample<<pdepth && vresolu <= psample<<pdepth) {
80                          sprintf(buf, "%d sampling...\n", 1<<pdepth);
# Line 171 | Line 83 | rview()                                /* do a view */
83                  } else {
84                          sprintf(buf, "%d refining...\n", 1<<pdepth);
85                          (*dev->comout)(buf);
86 <                        refine(&ptrunk, 0, 0, hresolu, vresolu, pdepth+1);
86 >                        refine(&ptrunk, pdepth+1);
87                  }
88 <                if (errno == ENOMEM)            /* ran out of memory */
89 <                        freereserves();
90 <                else if (dev->inpready)         /* noticed some input */
88 >                if (waitrays() < 0)
89 >                        quit(1);
90 >                if (dev->inpready)              /* noticed some input */
91                          command(": ");
92                  else                            /* finished this depth */
93                          pdepth++;
# Line 184 | Line 96 | rview()                                /* do a view */
96  
97  
98   void
99 < fillreserves()                  /* fill memory reserves */
99 > command(                        /* get/execute command */
100 >        char  *prompt
101 > )
102   {
189 –        if (reserve_mem != NULL)
190 –                return;
191 –        reserve_mem = (char *)malloc(RESERVE_AMT);
192 – }
193 –
194 –
195 – void
196 – freereserves()                  /* free memory reserves */
197 – {
198 –        if (reserve_mem == NULL)
199 –                return;
200 –        free(reserve_mem);
201 –        reserve_mem = NULL;
202 – }
203 –
204 –
205 – void
206 – command(prompt)                 /* get/execute command */
207 – char  *prompt;
208 – {
103   #define  badcom(s)      strncmp(s, inpbuf, args-inpbuf-1)
104          char  inpbuf[256];
105          char  *args;
# Line 268 | Line 162 | again:
162          case 'n':                               /* new picture */
163                  if (badcom("new"))
164                          goto commerr;
165 <                newimage();
165 >                newimage(args);
166                  break;
167          case 't':                               /* trace a ray */
168                  if (badcom("trace"))
# Line 356 | Line 250 | commerr:
250  
251  
252   void
253 < rsample()                       /* sample the image */
253 > rsample(void)                   /* sample the image */
254   {
255          int  xsiz, ysiz, y;
362 –        RECT  r;
256          PNODE  *p;
257 <        register RECT  *rl;
258 <        register PNODE  **pl;
366 <        register int  x;
257 >        PNODE   **pl;
258 >        int  x;
259          /*
260           *     We initialize the bottom row in the image at our current
261           * resolution.  During sampling, we check super-pixels to the
# Line 373 | Line 265 | rsample()                      /* sample the image */
265           */
266          xsiz = (((long)(pframe.r-pframe.l)<<pdepth)+hresolu-1) / hresolu;
267          ysiz = (((long)(pframe.u-pframe.d)<<pdepth)+vresolu-1) / vresolu;
376 –        rl = (RECT *)malloc(xsiz*sizeof(RECT));
377 –        if (rl == NULL)
378 –                return;
268          pl = (PNODE **)malloc(xsiz*sizeof(PNODE *));
269 <        if (pl == NULL) {
381 <                free((void *)rl);
269 >        if (pl == NULL)
270                  return;
383 –        }
271          /*
272           * Initialize the bottom row.
273           */
274 <        rl[0].l = rl[0].d = 0;
388 <        rl[0].r = hresolu; rl[0].u = vresolu;
389 <        pl[0] = findrect(pframe.l, pframe.d, &ptrunk, rl, pdepth);
274 >        pl[0] = findrect(pframe.l, pframe.d, &ptrunk, pdepth);
275          for (x = 1; x < xsiz; x++) {
391 –                rl[x].l = rl[x].d = 0;
392 –                rl[x].r = hresolu; rl[x].u = vresolu;
276                  pl[x] = findrect(pframe.l+((x*hresolu)>>pdepth),
277 <                                pframe.d, &ptrunk, rl+x, pdepth);
277 >                                pframe.d, &ptrunk, pdepth);
278          }
279                                                  /* sample the image */
280          for (y = 0; /* y < ysiz */ ; y++) {
# Line 403 | Line 286 | rsample()                      /* sample the image */
286                           */
287                          if (pl[x] != pl[x+1] && bigdiff(pl[x]->v,
288                                          pl[x+1]->v, maxdiff)) {
289 <                                refine(pl[x], rl[x].l, rl[x].d,
290 <                                                rl[x].r, rl[x].u, 1);
408 <                                refine(pl[x+1], rl[x+1].l, rl[x+1].d,
409 <                                                rl[x+1].r, rl[x+1].u, 1);
289 >                                refine(pl[x], 1);
290 >                                refine(pl[x+1], 1);
291                          }
292                  }
293                  if (y >= ysiz-1)
# Line 417 | Line 298 | rsample()                      /* sample the image */
298                          /*
299                           * Find super-pixel at this position in next row.
300                           */
420 –                        r.l = r.d = 0;
421 –                        r.r = hresolu; r.u = vresolu;
301                          p = findrect(pframe.l+((x*hresolu)>>pdepth),
302                                  pframe.d+(((y+1)*vresolu)>>pdepth),
303 <                                        &ptrunk, &r, pdepth);
303 >                                        &ptrunk, pdepth);
304                          /*
305                           * Test super-pixel in next row.
306                           */
307                          if (pl[x] != p && bigdiff(pl[x]->v, p->v, maxdiff)) {
308 <                                refine(pl[x], rl[x].l, rl[x].d,
309 <                                                rl[x].r, rl[x].u, 1);
431 <                                refine(p, r.l, r.d, r.r, r.u, 1);
308 >                                refine(pl[x], 1);
309 >                                refine(p, 1);
310                          }
311                          /*
312                           * Copy into super-pixel array.
313                           */
436 –                        rl[x].l = r.l; rl[x].d = r.d;
437 –                        rl[x].r = r.r; rl[x].u = r.u;
314                          pl[x] = p;
315                  }
316          }
317   escape:
442 –        free((void *)rl);
318          free((void *)pl);
319   }
320  
321  
322   int
323 < refine(p, xmin, ymin, xmax, ymax, pd)           /* refine a node */
324 < register PNODE  *p;
325 < int  xmin, ymin, xmax, ymax;
326 < int  pd;
323 > refine(                         /* refine a node */
324 >        PNODE   *p,
325 >        int  pd
326 > )
327   {
328          int  growth;
329          int  mx, my;
# Line 460 | Line 335 | int  pd;
335          if (pd <= 0)                            /* depth limit */
336                  return(0);
337  
338 <        mx = (xmin + xmax) >> 1;
339 <        my = (ymin + ymax) >> 1;
338 >        mx = (p->xmin + p->xmax) >> 1;
339 >        my = (p->ymin + p->ymax) >> 1;
340          growth = 0;
341  
342          if (p->kid == NULL) {                   /* subdivide */
343  
344                  if ((p->kid = newptree()) == NULL)
345                          return(0);
346 +
347 +                p->kid[UR].xmin = mx;
348 +                p->kid[UR].ymin = my;
349 +                p->kid[UR].xmax = p->xmax;
350 +                p->kid[UR].ymax = p->ymax;
351 +                p->kid[UL].xmin = p->xmin;
352 +                p->kid[UL].ymin = my;
353 +                p->kid[UL].xmax = mx;
354 +                p->kid[UL].ymax = p->ymax;
355 +                p->kid[DR].xmin = mx;
356 +                p->kid[DR].ymin = p->ymin;
357 +                p->kid[DR].xmax = p->xmax;
358 +                p->kid[DR].ymax = my;
359 +                p->kid[DL].xmin = p->xmin;
360 +                p->kid[DL].ymin = p->ymin;
361 +                p->kid[DL].xmax = mx;
362 +                p->kid[DL].ymax = my;
363                  /*
364                   *  The following paint order can leave a black pixel
365                   *  if redraw() is called in (*dev->paintr)().
366                   */
367                  if (p->x >= mx && p->y >= my)
368                          pcopy(p, p->kid+UR);
369 <                else
370 <                        paint(p->kid+UR, mx, my, xmax, ymax);
369 >                else if (paint(p->kid+UR) < 0)
370 >                        quit(1);
371                  if (p->x < mx && p->y >= my)
372                          pcopy(p, p->kid+UL);
373 <                else
374 <                        paint(p->kid+UL, xmin, my, mx, ymax);
373 >                else if (paint(p->kid+UL) < 0)
374 >                        quit(1);
375                  if (p->x >= mx && p->y < my)
376                          pcopy(p, p->kid+DR);
377 <                else
378 <                        paint(p->kid+DR, mx, ymin, xmax, my);
377 >                else if (paint(p->kid+DR) < 0)
378 >                        quit(1);
379                  if (p->x < mx && p->y < my)
380                          pcopy(p, p->kid+DL);
381 <                else
382 <                        paint(p->kid+DL, xmin, ymin, mx, my);
381 >                else if (paint(p->kid+DL) < 0)
382 >                        quit(1);
383  
384                  growth++;
385          }
386                                                  /* do children */
387          if (mx > pframe.l) {
388                  if (my > pframe.d)
389 <                        growth += refine(p->kid+DL, xmin, ymin, mx, my, pd-1);
389 >                        growth += refine(p->kid+DL, pd-1);
390                  if (my < pframe.u)
391 <                        growth += refine(p->kid+UL, xmin, my, mx, ymax, pd-1);
391 >                        growth += refine(p->kid+UL, pd-1);
392          }
393          if (mx < pframe.r) {
394                  if (my > pframe.d)
395 <                        growth += refine(p->kid+DR, mx, ymin, xmax, my, pd-1);
395 >                        growth += refine(p->kid+DR, pd-1);
396                  if (my < pframe.u)
397 <                        growth += refine(p->kid+UR, mx, my, xmax, ymax, pd-1);
397 >                        growth += refine(p->kid+UR, pd-1);
398          }
399                                                  /* recompute sum */
400          if (growth) {

Diff Legend

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