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

Diff Legend

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