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 1.10 by greg, Tue Jan 30 11:37:41 1990 UTC vs.
Revision 2.32 by greg, Sat Dec 12 19:01:00 2009 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1987 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   *  rview.c - routines and variables for interactive view generation.
6   *
7 < *     3/24/87
7 > *  External symbols declared in rpaint.h
8   */
9  
10 < #include  "standard.h"
10 > #include "copyright.h"
11  
15 – #include  "color.h"
16 –
17 – #include  "rpaint.h"
18 –
12   #include  <signal.h>
20 –
13   #include  <ctype.h>
14  
15 < VIEW  ourview = STDVIEW;                /* viewing parameters */
16 < int  hresolu, vresolu;                  /* image resolution */
15 > #include  "ray.h"
16 > #include  "rpaint.h"
17  
18 < int  psample = 8;                       /* pixel sample size */
27 < double  maxdiff = .15;                  /* max. sample difference */
18 > #define  CTRL(c)        ((c)-'@')
19  
29 – double  exposure = 1.0;                 /* exposure for scene */
20  
21 < double  dstrsrc = 0.0;                  /* square source distribution */
32 < double  shadthresh = .1;                /* shadow threshold */
33 < double  shadcert = .25;                 /* shadow certainty */
34 <
35 < int  maxdepth = 4;                      /* maximum recursion depth */
36 < double  minweight = 1e-2;               /* minimum ray weight */
37 <
38 < COLOR  ambval = BLKCOLOR;               /* ambient value */
39 < double  ambacc = 0.2;                   /* ambient accuracy */
40 < int  ambres = 8;                        /* ambient resolution */
41 < int  ambdiv = 32;                       /* ambient divisions */
42 < int  ambssamp = 0;                      /* ambient super-samples */
43 < int  ambounce = 0;                      /* ambient bounces */
44 < char  *amblist[128];                    /* ambient include/exclude list */
45 < int  ambincl = -1;                      /* include == 1, exclude == 0 */
46 <
47 < int  greyscale = 0;                     /* map colors to brightness? */
48 < char  *devname = "X";                   /* output device name */
49 <
50 < struct driver  *dev = NULL;             /* driver functions */
51 <
52 < VIEW  oldview;                          /* previous view parameters */
53 <
54 < PNODE  ptrunk;                          /* the base of our image */
55 < RECT  pframe;                           /* current frame boundaries */
56 < int  pdepth;                            /* image depth in current frame */
57 <
58 < #define  CTRL(c)        ('c'-'@')
59 <
60 <
21 > void
22   quit(code)                      /* quit program */
23   int  code;
24   {
25 + #ifdef MSTATS
26 +        if (code == 2 && errno == ENOMEM)
27 +                printmemstats(stderr);
28 + #endif
29 +        if (ray_pnprocs > 0)    /* close children if any */
30 +                ray_pclose(0);          
31          devclose();
32          exit(code);
33   }
34  
35  
36 < devopen(dname)                          /* open device driver */
37 < char  *dname;
36 > void
37 > devopen(                                /* open device driver */
38 >        char  *dname
39 > )
40   {
41          extern char  *progname, *octname;
42          char  *id;
43 <        register int  i;
43 >        int  i;
44  
45          id = octname!=NULL ? octname : progname;
46                                                  /* check device table */
47          for (i = 0; devtable[i].name; i++)
48 <                if (!strcmp(dname, devtable[i].name))
48 >                if (!strcmp(dname, devtable[i].name)) {
49                          if ((dev = (*devtable[i].init)(dname, id)) == NULL) {
50                                  sprintf(errmsg, "cannot initialize %s", dname);
51                                  error(USER, errmsg);
52                          } else
53                                  return;
54 +                }
55                                                  /* not there, try exec */
56          if ((dev = comm_init(dname, id)) == NULL) {
57                  sprintf(errmsg, "cannot start device \"%s\"", dname);
# Line 90 | Line 60 | char  *dname;
60   }
61  
62  
63 < devclose()                              /* close our device */
63 > void
64 > devclose(void)                          /* close our device */
65   {
66          if (dev != NULL)
67                  (*dev->close)();
# Line 98 | Line 69 | devclose()                             /* close our device */
69   }
70  
71  
72 < printdevices()                          /* print list of output devices */
72 > void
73 > printdevices(void)                      /* print list of output devices */
74   {
75 <        register int  i;
75 >        int  i;
76  
77          for (i = 0; devtable[i].name; i++)
78                  printf("%-16s # %s\n", devtable[i].name, devtable[i].descrip);
79   }
80  
81  
82 < rview()                         /* do a view */
82 > void
83 > rview(void)                             /* do a view */
84   {
85          char  buf[32];
86  
87 <        devopen(devname);               /* open device */
88 <        newimage();                     /* set up image */
87 >        devopen(dvcname);               /* open device */
88 >        newimage(NULL);                 /* start image */
89  
90          for ( ; ; ) {                   /* quit in command() */
91 <                while (hresolu <= 1<<pdepth &&
119 <                                vresolu <= 1<<pdepth)
91 >                while (hresolu <= 1<<pdepth && vresolu <= 1<<pdepth)
92                          command("done: ");
93 <
94 <                if (hresolu <= psample<<pdepth &&
123 <                                vresolu <= psample<<pdepth) {
93 >                errno = 0;
94 >                if (hresolu <= psample<<pdepth && vresolu <= psample<<pdepth) {
95                          sprintf(buf, "%d sampling...\n", 1<<pdepth);
96                          (*dev->comout)(buf);
97                          rsample();
98                  } else {
99                          sprintf(buf, "%d refining...\n", 1<<pdepth);
100                          (*dev->comout)(buf);
101 <                        refine(&ptrunk, 0, 0, hresolu, vresolu, pdepth+1);
101 >                        refine(&ptrunk, pdepth+1);
102                  }
103 <                if (dev->inpready)
103 >                if (dev->inpready)              /* noticed some input */
104                          command(": ");
105 <                else
105 >                else                            /* finished this depth */
106                          pdepth++;
107          }
108   }
109  
110  
111 < command(prompt)                 /* get/execute command */
112 < char  *prompt;
111 > void
112 > command(                        /* get/execute command */
113 >        char  *prompt
114 > )
115   {
116 < #define  badcom(s)      strncmp(s, inpbuf, args-inpbuf-1)
144 <        double  atof();
116 > #define  badcom(s)      strncmp(s, inpbuf, args-inpbuf-1)
117          char  inpbuf[256];
118          char  *args;
119   again:
# Line 150 | Line 122 | again:
122                  ;
123          if (*args) *args++ = '\0';
124          else *++args = '\0';
125 +
126 +        if (waitrays() < 0)                     /* clear ray queue */
127 +                quit(1);
128          
129          switch (inpbuf[0]) {
130 <        case 'f':                               /* new frame */
131 <                if (badcom("frame"))
132 <                        goto commerr;
130 >        case 'f':                               /* new frame (|focus|free) */
131 >                if (badcom("frame")) {
132 >                        if (badcom("focus")) {
133 >                                if (badcom("free"))
134 >                                        goto commerr;
135 >                                free_objmem();
136 >                                break;
137 >                        }
138 >                        getfocus(args);
139 >                        break;
140 >                }
141                  getframe(args);
142                  break;
143          case 'v':                               /* view */
# Line 167 | Line 150 | again:
150                          goto commerr;
151                  lastview(args);
152                  break;
153 +        case 'V':                               /* save view */
154 +                if (badcom("V"))
155 +                        goto commerr;
156 +                saveview(args);
157 +                break;
158 +        case 'L':                               /* load view */
159 +                if (badcom("L"))
160 +                        goto commerr;
161 +                loadview(args);
162 +                break;
163          case 'e':                               /* exposure */
164                  if (badcom("exposure"))
165                          goto commerr;
166                  getexposure(args);
167                  break;
168          case 's':                               /* set a parameter */
169 <                if (badcom("set"))
169 >                if (badcom("set")) {
170 > #ifdef  SIGTSTP
171 >                        if (!badcom("stop"))
172 >                                goto dostop;
173 > #endif
174                          goto commerr;
175 +                }
176                  setparam(args);
177                  break;
178          case 'n':                               /* new picture */
179                  if (badcom("new"))
180                          goto commerr;
181 <                newimage();
181 >                newimage(args);
182                  break;
183          case 't':                               /* trace a ray */
184                  if (badcom("trace"))
# Line 192 | Line 190 | again:
190                          goto commerr;
191                  getaim(args);
192                  break;
193 <        case 'm':                               /* move camera */
193 >        case 'm':                               /* move camera (or memstats) */
194                  if (badcom("move"))
195 + #ifdef  MSTATS
196 +                {
197 +                        if (badcom("memory"))
198 +                                goto commerr;
199 +                        printmemstats(stderr);
200 +                        break;
201 +                }
202 + #else
203                          goto commerr;
204 + #endif
205                  getmove(args);
206                  break;
207          case 'r':                               /* rotate/repaint */
208                  if (badcom("rotate")) {
209 <                        if (badcom("repaint"))
210 <                                goto commerr;
209 >                        if (badcom("repaint")) {
210 >                                if (badcom("redraw"))
211 >                                        goto commerr;
212 >                                redraw();
213 >                                break;
214 >                        }
215                          getrepaint(args);
216                          break;
217                  }
218                  getrotate(args);
219                  break;
220          case 'p':                               /* pivot view */
221 <                if (badcom("pivot"))
222 <                        goto commerr;
221 >                if (badcom("pivot")) {
222 >                        if (badcom("pause"))
223 >                                goto commerr;
224 >                        goto again;
225 >                }
226                  getpivot(args);
227                  break;
228 <        case CTRL(R):                           /* redraw */
228 >        case CTRL('R'):                         /* redraw */
229                  redraw();
230                  break;
231          case 'w':                               /* write */
# Line 223 | Line 237 | again:
237                  if (badcom("quit"))
238                          goto commerr;
239                  quit(0);
240 <        case CTRL(C):                           /* interrupt */
240 >        case CTRL('C'):                         /* interrupt */
241                  goto again;
242 < #ifdef  SIGTSTP
243 <        case CTRL(Z):                           /* stop */
242 > #ifdef  SIGTSTP
243 >        case CTRL('Z'):;                        /* stop */
244 > dostop:
245                  devclose();
246                  kill(0, SIGTSTP);
247                  /* pc stops here */
248 <                devopen(devname);
248 >                devopen(dvcname);
249                  redraw();
250                  break;
251   #endif
# Line 246 | Line 261 | commerr:
261                  error(COMMAND, errmsg);
262                  break;
263          }
264 < #undef  badcom
264 > #undef  badcom
265   }
266  
267  
268 < rsample()                       /* sample the image */
268 > void
269 > rsample(void)                   /* sample the image */
270   {
271          int  xsiz, ysiz, y;
256 –        RECT  r;
272          PNODE  *p;
273 <        register RECT  *rl;
274 <        register PNODE  **pl;
260 <        register int  x;
273 >        PNODE   **pl;
274 >        int  x;
275          /*
276           *     We initialize the bottom row in the image at our current
277 <         * resolution.  During sampling, we check super-pixels to the
277 >         * resolution.  During sampling, we check super-pixels to the
278           * right and above by calling bigdiff().  If there is a significant
279           * difference, we subsample the super-pixels.  The testing process
280           * includes initialization of the next row.
281           */
282 <        xsiz = (((pframe.r-pframe.l)<<pdepth)+hresolu-1) / hresolu;
283 <        ysiz = (((pframe.u-pframe.d)<<pdepth)+vresolu-1) / vresolu;
270 <        rl = (RECT *)malloc(xsiz*sizeof(RECT));
282 >        xsiz = (((long)(pframe.r-pframe.l)<<pdepth)+hresolu-1) / hresolu;
283 >        ysiz = (((long)(pframe.u-pframe.d)<<pdepth)+vresolu-1) / vresolu;
284          pl = (PNODE **)malloc(xsiz*sizeof(PNODE *));
285 <        if (rl == NULL || pl == NULL)
286 <                error(SYSTEM, "out of memory in rsample");
285 >        if (pl == NULL)
286 >                return;
287          /*
288           * Initialize the bottom row.
289           */
290 <        rl[0].l = rl[0].d = 0;
278 <        rl[0].r = hresolu; rl[0].u = vresolu;
279 <        pl[0] = findrect(pframe.l, pframe.d, &ptrunk, rl, pdepth);
290 >        pl[0] = findrect(pframe.l, pframe.d, &ptrunk, pdepth);
291          for (x = 1; x < xsiz; x++) {
281 –                rl[x].l = rl[x].d = 0;
282 –                rl[x].r = hresolu; rl[x].u = vresolu;
292                  pl[x] = findrect(pframe.l+((x*hresolu)>>pdepth),
293 <                                pframe.d, &ptrunk, rl+x, pdepth);
293 >                                pframe.d, &ptrunk, pdepth);
294          }
295                                                  /* sample the image */
296          for (y = 0; /* y < ysiz */ ; y++) {
297                  for (x = 0; x < xsiz-1; x++) {
298 <                        if (dev->inpready)
298 >                        if (dev->inpready || errno == ENOMEM)
299                                  goto escape;
300                          /*
301                           * Test super-pixel to the right.
302                           */
303                          if (pl[x] != pl[x+1] && bigdiff(pl[x]->v,
304                                          pl[x+1]->v, maxdiff)) {
305 <                                refine(pl[x], rl[x].l, rl[x].d,
306 <                                                rl[x].r, rl[x].u, 1);
298 <                                refine(pl[x+1], rl[x+1].l, rl[x+1].d,
299 <                                                rl[x+1].r, rl[x+1].u, 1);
305 >                                refine(pl[x], 1);
306 >                                refine(pl[x+1], 1);
307                          }
308                  }
309                  if (y >= ysiz-1)
310                          break;
311                  for (x = 0; x < xsiz; x++) {
312 <                        if (dev->inpready)
312 >                        if (dev->inpready || errno == ENOMEM)
313                                  goto escape;
314                          /*
315                           * Find super-pixel at this position in next row.
316                           */
310 –                        r.l = r.d = 0;
311 –                        r.r = hresolu; r.u = vresolu;
317                          p = findrect(pframe.l+((x*hresolu)>>pdepth),
318                                  pframe.d+(((y+1)*vresolu)>>pdepth),
319 <                                        &ptrunk, &r, pdepth);
319 >                                        &ptrunk, pdepth);
320                          /*
321                           * Test super-pixel in next row.
322                           */
323                          if (pl[x] != p && bigdiff(pl[x]->v, p->v, maxdiff)) {
324 <                                refine(pl[x], rl[x].l, rl[x].d,
325 <                                                rl[x].r, rl[x].u, 1);
321 <                                refine(p, r.l, r.d, r.r, r.u, 1);
324 >                                refine(pl[x], 1);
325 >                                refine(p, 1);
326                          }
327                          /*
328                           * Copy into super-pixel array.
329                           */
326 –                        rl[x].l = r.l; rl[x].d = r.d;
327 –                        rl[x].r = r.r; rl[x].u = r.u;
330                          pl[x] = p;
331                  }
332          }
333   escape:
334 <        free((char *)rl);
333 <        free((char *)pl);
334 >        free((void *)pl);
335   }
336  
337  
338   int
339 < refine(p, xmin, ymin, xmax, ymax, pd)           /* refine a node */
340 < register PNODE  *p;
341 < int  xmin, ymin, xmax, ymax;
342 < int  pd;
339 > refine(                         /* refine a node */
340 >        PNODE   *p,
341 >        int  pd
342 > )
343   {
344          int  growth;
345          int  mx, my;
# Line 350 | Line 351 | int  pd;
351          if (pd <= 0)                            /* depth limit */
352                  return(0);
353  
354 <        mx = (xmin + xmax) >> 1;
355 <        my = (ymin + ymax) >> 1;
354 >        mx = (p->xmin + p->xmax) >> 1;
355 >        my = (p->ymin + p->ymax) >> 1;
356          growth = 0;
357  
358          if (p->kid == NULL) {                   /* subdivide */
359  
360                  if ((p->kid = newptree()) == NULL)
361 <                        error(SYSTEM, "out of memory in refine");
361 >                        return(0);
362 >
363 >                p->kid[UR].xmin = mx;
364 >                p->kid[UR].ymin = my;
365 >                p->kid[UR].xmax = p->xmax;
366 >                p->kid[UR].ymax = p->ymax;
367 >                p->kid[UL].xmin = p->xmin;
368 >                p->kid[UL].ymin = my;
369 >                p->kid[UL].xmax = mx;
370 >                p->kid[UL].ymax = p->ymax;
371 >                p->kid[DR].xmin = mx;
372 >                p->kid[DR].ymin = p->ymin;
373 >                p->kid[DR].xmax = p->xmax;
374 >                p->kid[DR].ymax = my;
375 >                p->kid[DL].xmin = p->xmin;
376 >                p->kid[DL].ymin = p->ymin;
377 >                p->kid[DL].xmax = mx;
378 >                p->kid[DL].ymax = my;
379                  /*
380                   *  The following paint order can leave a black pixel
381 <                 *  when redraw() is called in (*dev->paintr)().
381 >                 *  if redraw() is called in (*dev->paintr)().
382                   */
383                  if (p->x >= mx && p->y >= my)
384                          pcopy(p, p->kid+UR);
385 <                else
386 <                        paint(p->kid+UR, mx, my, xmax, ymax);
385 >                else if (paint(p->kid+UR) < 0)
386 >                        quit(1);
387                  if (p->x < mx && p->y >= my)
388                          pcopy(p, p->kid+UL);
389 <                else
390 <                        paint(p->kid+UL, xmin, my, mx, ymax);
389 >                else if (paint(p->kid+UL) < 0)
390 >                        quit(1);
391                  if (p->x >= mx && p->y < my)
392                          pcopy(p, p->kid+DR);
393 <                else
394 <                        paint(p->kid+DR, mx, ymin, xmax, my);
393 >                else if (paint(p->kid+DR) < 0)
394 >                        quit(1);
395                  if (p->x < mx && p->y < my)
396                          pcopy(p, p->kid+DL);
397 <                else
398 <                        paint(p->kid+DL, xmin, ymin, mx, my);
397 >                else if (paint(p->kid+DL) < 0)
398 >                        quit(1);
399  
400                  growth++;
401          }
402                                                  /* do children */
403          if (mx > pframe.l) {
404                  if (my > pframe.d)
405 <                        growth += refine(p->kid+DL, xmin, ymin, mx, my, pd-1);
405 >                        growth += refine(p->kid+DL, pd-1);
406                  if (my < pframe.u)
407 <                        growth += refine(p->kid+UL, xmin, my, mx, ymax, pd-1);
407 >                        growth += refine(p->kid+UL, pd-1);
408          }
409          if (mx < pframe.r) {
410                  if (my > pframe.d)
411 <                        growth += refine(p->kid+DR, mx, ymin, xmax, my, pd-1);
411 >                        growth += refine(p->kid+DR, pd-1);
412                  if (my < pframe.u)
413 <                        growth += refine(p->kid+UR, mx, my, xmax, ymax, pd-1);
396 <        }
397 <                                                /* recompute sum */
398 <        if (growth) {
399 <                setcolor(p->v, 0.0, 0.0, 0.0);
400 <                for (i = 0; i < 4; i++)
401 <                        addcolor(p->v, p->kid[i].v);
402 <                scalecolor(p->v, 0.25);
413 >                        growth += refine(p->kid+UR, pd-1);
414          }
415          return(growth);
416   }

Diff Legend

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