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

Diff Legend

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