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.4 by greg, Tue Jun 13 10:57:49 1989 UTC vs.
Revision 2.25 by greg, Sun Mar 28 16:31:14 2004 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  
12 < #include  "color.h"
12 > #include  <signal.h>
13 > #include  <ctype.h>
14  
15 + #include  "ray.h"
16   #include  "rpaint.h"
17  
19 – #include  <signal.h>
18  
19 < #include  <ctype.h>
19 > CUBE  thescene;                         /* our scene */
20 > OBJECT  nsceneobjs;                     /* number of objects in our scene */
21  
22 < VIEW  ourview = STDVIEW(470);           /* viewing parameters */
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  psample = 8;                       /* pixel sample size */
37 < double  maxdiff = .15;                  /* max. sample difference */
37 > double  maxdiff = .15;                  /* max. sample difference */
38  
39 < double  exposure = 1.0;                 /* exposure for scene */
39 > double  exposure = 1.0;                 /* exposure for scene */
40  
41 < double  dstrsrc = 0.0;                  /* square source distribution */
42 < double  shadthresh = .1;                /* shadow threshold */
43 < double  shadcert = .25;                 /* shadow certainty */
41 > double  dstrsrc = 0.0;                  /* square source distribution */
42 > double  shadthresh = .1;                /* shadow threshold */
43 > double  shadcert = .25;                 /* shadow certainty */
44 > int  directrelay = 0;                   /* number of source relays */
45 > int  vspretest = 128;                   /* virtual source pretest density */
46 > int  directvis = 1;                     /* sources visible? */
47 > double  srcsizerat = 0.;                /* maximum ratio source size/dist. */
48  
49 < int  maxdepth = 4;                      /* maximum recursion depth */
50 < double  minweight = 1e-2;               /* minimum ray weight */
49 > COLOR  cextinction = BLKCOLOR;          /* global extinction coefficient */
50 > COLOR  salbedo = BLKCOLOR;              /* global scattering albedo */
51 > double  seccg = 0.;                     /* global scattering eccentricity */
52 > double  ssampdist = 0.;                 /* scatter sampling distance */
53  
54 + double  specthresh = .3;                /* specular sampling threshold */
55 + double  specjitter = 1.;                /* specular sampling jitter */
56 +
57 + int  backvis = 1;                       /* back face visibility */
58 +
59 + int  maxdepth = 6;                      /* maximum recursion depth */
60 + double  minweight = 1e-2;               /* minimum ray weight */
61 +
62 + char  *ambfile = NULL;                  /* ambient file name */
63   COLOR  ambval = BLKCOLOR;               /* ambient value */
64 < double  ambacc = 0.2;                   /* ambient accuracy */
65 < int  ambres = 64;                       /* ambient resolution */
66 < int  ambdiv = 32;                       /* ambient divisions */
67 < int  ambssamp = 0;                      /* ambient super-samples */
64 > int  ambvwt = 0;                        /* initial weight for ambient value */
65 > double  ambacc = 0.3;                   /* ambient accuracy */
66 > int  ambres = 32;                       /* ambient resolution */
67 > int  ambdiv = 256;                      /* ambient divisions */
68 > int  ambssamp = 64;                     /* ambient super-samples */
69   int  ambounce = 0;                      /* ambient bounces */
70   char  *amblist[128];                    /* ambient include/exclude list */
71   int  ambincl = -1;                      /* include == 1, exclude == 0 */
72  
73   int  greyscale = 0;                     /* map colors to brightness? */
74 < char  *devname = "X";                   /* output device name */
74 > char  *dvcname = dev_default;           /* output device name */
75  
76   struct driver  *dev = NULL;             /* driver functions */
77  
78 + char  rifname[128];                     /* rad input file name */
79 +
80   VIEW  oldview;                          /* previous view parameters */
81  
82   PNODE  ptrunk;                          /* the base of our image */
83   RECT  pframe;                           /* current frame boundaries */
84   int  pdepth;                            /* image depth in current frame */
85  
86 < #define  CTRL(c)        ('c'-'@')
86 > static char  *reserve_mem = NULL;       /* pre-allocated reserve memory */
87  
88 + #define RESERVE_AMT     32768           /* amount of memory to reserve */
89  
90 + #define  CTRL(c)        ((c)-'@')
91 +
92 +
93 + void
94   quit(code)                      /* quit program */
95   int  code;
96   {
97 + #ifdef MSTATS
98 +        if (code == 2 && errno == ENOMEM)
99 +                printmemstats(stderr);
100 + #endif
101          devclose();
102          exit(code);
103   }
104  
105  
106 + void
107   devopen(dname)                          /* open device driver */
108   char  *dname;
109   {
110 <        extern char  *progname;
111 <        char  *devargv[3];
110 >        extern char  *progname, *octname;
111 >        char  *id;
112          register int  i;
113 +
114 +        id = octname!=NULL ? octname : progname;
115                                                  /* check device table */
116          for (i = 0; devtable[i].name; i++)
117 <                if (!strcmp(dname, devtable[i].name))
118 <                        if ((dev = (*devtable[i].init)(progname)) == NULL) {
117 >                if (!strcmp(dname, devtable[i].name)) {
118 >                        if ((dev = (*devtable[i].init)(dname, id)) == NULL) {
119                                  sprintf(errmsg, "cannot initialize %s", dname);
120                                  error(USER, errmsg);
121                          } else
122                                  return;
123 +                }
124                                                  /* not there, try exec */
125 <        devargv[0] = dname;
84 <        devargv[1] = progname;
85 <        devargv[2] = NULL;
86 <        if ((dev = comm_init(devargv)) == NULL) {
125 >        if ((dev = comm_init(dname, id)) == NULL) {
126                  sprintf(errmsg, "cannot start device \"%s\"", dname);
127                  error(USER, errmsg);
128          }
129   }
130  
131  
132 + void
133   devclose()                              /* close our device */
134   {
135          if (dev != NULL)
# Line 98 | Line 138 | devclose()                             /* close our device */
138   }
139  
140  
141 + void
142   printdevices()                          /* print list of output devices */
143   {
144          register int  i;
# Line 107 | Line 148 | printdevices()                         /* print list of output devices */
148   }
149  
150  
151 + void
152   rview()                         /* do a view */
153   {
154          char  buf[32];
155  
156 <        devopen(devname);               /* open device */
157 <        newimage();                     /* set up image */
156 >        devopen(dvcname);               /* open device */
157 >        newimage();                     /* start image (calls fillreserves) */
158  
159          for ( ; ; ) {                   /* quit in command() */
160 <                while (ourview.hresolu <= 1<<pdepth &&
119 <                                ourview.vresolu <= 1<<pdepth)
160 >                while (hresolu <= 1<<pdepth && vresolu <= 1<<pdepth)
161                          command("done: ");
162 <
163 <                if (ourview.hresolu <= psample<<pdepth &&
164 <                                ourview.vresolu <= psample<<pdepth) {
162 >                while (reserve_mem == NULL)
163 >                        command("out of memory: ");
164 >                errno = 0;
165 >                if (hresolu <= psample<<pdepth && vresolu <= psample<<pdepth) {
166                          sprintf(buf, "%d sampling...\n", 1<<pdepth);
167                          (*dev->comout)(buf);
168                          rsample();
169                  } else {
170                          sprintf(buf, "%d refining...\n", 1<<pdepth);
171                          (*dev->comout)(buf);
172 <                        refine(&ptrunk, 0, 0, ourview.hresolu,
131 <                                        ourview.vresolu, pdepth+1);
172 >                        refine(&ptrunk, 0, 0, hresolu, vresolu, pdepth+1);
173                  }
174 <                if (dev->inpready)
174 >                if (errno == ENOMEM)            /* ran out of memory */
175 >                        freereserves();
176 >                else if (dev->inpready)         /* noticed some input */
177                          command(": ");
178 <                else
178 >                else                            /* finished this depth */
179                          pdepth++;
180          }
181   }
182  
183  
184 + void
185 + fillreserves()                  /* fill memory reserves */
186 + {
187 +        if (reserve_mem != NULL)
188 +                return;
189 +        reserve_mem = (char *)malloc(RESERVE_AMT);
190 + }
191 +
192 +
193 + void
194 + freereserves()                  /* free memory reserves */
195 + {
196 +        if (reserve_mem == NULL)
197 +                return;
198 +        free(reserve_mem);
199 +        reserve_mem = NULL;
200 + }
201 +
202 +
203 + void
204   command(prompt)                 /* get/execute command */
205   char  *prompt;
206   {
207 < #define  badcom(s)      strncmp(s, inpbuf, args-inpbuf-1)
145 <        double  atof();
207 > #define  badcom(s)      strncmp(s, inpbuf, args-inpbuf-1)
208          char  inpbuf[256];
209          char  *args;
210   again:
211 <        (*dev->comout)(prompt);                 /* get command + arguments */
150 <        (*dev->comin)(inpbuf);
211 >        (*dev->comin)(inpbuf, prompt);          /* get command + arguments */
212          for (args = inpbuf; *args && *args != ' '; args++)
213                  ;
214          if (*args) *args++ = '\0';
215          else *++args = '\0';
216          
217          switch (inpbuf[0]) {
218 <        case 'f':                               /* new frame */
219 <                if (badcom("frame"))
220 <                        goto commerr;
218 >        case 'f':                               /* new frame (or free mem.) */
219 >                if (badcom("frame")) {
220 >                        if (badcom("free"))
221 >                                goto commerr;
222 >                        free_objmem();
223 >                        break;
224 >                }
225                  getframe(args);
226                  break;
227          case 'v':                               /* view */
# Line 169 | Line 234 | again:
234                          goto commerr;
235                  lastview(args);
236                  break;
237 +        case 'V':                               /* save view */
238 +                if (badcom("V"))
239 +                        goto commerr;
240 +                saveview(args);
241 +                break;
242 +        case 'L':                               /* load view */
243 +                if (badcom("L"))
244 +                        goto commerr;
245 +                loadview(args);
246 +                break;
247          case 'e':                               /* exposure */
248                  if (badcom("exposure"))
249                          goto commerr;
250                  getexposure(args);
251                  break;
252          case 's':                               /* set a parameter */
253 <                if (badcom("set"))
253 >                if (badcom("set")) {
254 > #ifdef  SIGTSTP
255 >                        if (!badcom("stop"))
256 >                                goto dostop;
257 > #endif
258                          goto commerr;
259 +                }
260                  setparam(args);
261                  break;
262          case 'n':                               /* new picture */
# Line 194 | Line 274 | again:
274                          goto commerr;
275                  getaim(args);
276                  break;
277 <        case 'm':                               /* move camera */
277 >        case 'm':                               /* move camera (or memstats) */
278                  if (badcom("move"))
279 + #ifdef  MSTATS
280 +                {
281 +                        if (badcom("memory"))
282 +                                goto commerr;
283 +                        printmemstats(stderr);
284 +                        break;
285 +                }
286 + #else
287                          goto commerr;
288 + #endif
289                  getmove(args);
290                  break;
291 <        case 'r':                               /* rotate camera */
292 <                if (badcom("rotate"))
293 <                        goto commerr;
291 >        case 'r':                               /* rotate/repaint */
292 >                if (badcom("rotate")) {
293 >                        if (badcom("repaint")) {
294 >                                if (badcom("redraw"))
295 >                                        goto commerr;
296 >                                redraw();
297 >                                break;
298 >                        }
299 >                        getrepaint(args);
300 >                        break;
301 >                }
302                  getrotate(args);
303                  break;
304          case 'p':                               /* pivot view */
305 <                if (badcom("pivot"))
306 <                        goto commerr;
305 >                if (badcom("pivot")) {
306 >                        if (badcom("pause"))
307 >                                goto commerr;
308 >                        goto again;
309 >                }
310                  getpivot(args);
311                  break;
312 <        case CTRL(R):                           /* redraw */
312 >        case CTRL('R'):                         /* redraw */
313                  redraw();
314                  break;
315          case 'w':                               /* write */
# Line 221 | Line 321 | again:
321                  if (badcom("quit"))
322                          goto commerr;
323                  quit(0);
324 <        case CTRL(C):                           /* interrupt */
324 >        case CTRL('C'):                         /* interrupt */
325                  goto again;
326 < #ifdef  SIGTSTP
327 <        case CTRL(Z):                           /* stop */
326 > #ifdef  SIGTSTP
327 >        case CTRL('Z'):;                        /* stop */
328 > dostop:
329                  devclose();
330                  kill(0, SIGTSTP);
331                  /* pc stops here */
332 <                devopen(devname);
332 >                devopen(dvcname);
333                  redraw();
334                  break;
335   #endif
# Line 244 | Line 345 | commerr:
345                  error(COMMAND, errmsg);
346                  break;
347          }
348 < #undef  badcom
348 > #undef  badcom
349   }
350  
351  
352 + void
353   rsample()                       /* sample the image */
354   {
355          int  xsiz, ysiz, y;
356          RECT  r;
357          PNODE  *p;
358          register RECT  *rl;
359 <        register PNODE  **pl;
359 >        register PNODE  **pl;
360          register int  x;
361          /*
362           *     We initialize the bottom row in the image at our current
363 <         * resolution.  During sampling, we check super-pixels to the
363 >         * resolution.  During sampling, we check super-pixels to the
364           * right and above by calling bigdiff().  If there is a significant
365           * difference, we subsample the super-pixels.  The testing process
366           * includes initialization of the next row.
367           */
368 <        xsiz = (((pframe.r-pframe.l)<<pdepth)+ourview.hresolu-1) /
369 <                        ourview.hresolu;
268 <        ysiz = (((pframe.u-pframe.d)<<pdepth)+ourview.vresolu-1) /
269 <                        ourview.vresolu;
368 >        xsiz = (((long)(pframe.r-pframe.l)<<pdepth)+hresolu-1) / hresolu;
369 >        ysiz = (((long)(pframe.u-pframe.d)<<pdepth)+vresolu-1) / vresolu;
370          rl = (RECT *)malloc(xsiz*sizeof(RECT));
371 +        if (rl == NULL)
372 +                return;
373          pl = (PNODE **)malloc(xsiz*sizeof(PNODE *));
374 <        if (rl == NULL || pl == NULL)
375 <                error(SYSTEM, "out of memory in rsample");
374 >        if (pl == NULL) {
375 >                free((void *)rl);
376 >                return;
377 >        }
378          /*
379           * Initialize the bottom row.
380           */
381          rl[0].l = rl[0].d = 0;
382 <        rl[0].r = ourview.hresolu; rl[0].u = ourview.vresolu;
382 >        rl[0].r = hresolu; rl[0].u = vresolu;
383          pl[0] = findrect(pframe.l, pframe.d, &ptrunk, rl, pdepth);
384          for (x = 1; x < xsiz; x++) {
385                  rl[x].l = rl[x].d = 0;
386 <                rl[x].r = ourview.hresolu; rl[x].u = ourview.vresolu;
387 <                pl[x] = findrect(pframe.l+((x*ourview.hresolu)>>pdepth),
386 >                rl[x].r = hresolu; rl[x].u = vresolu;
387 >                pl[x] = findrect(pframe.l+((x*hresolu)>>pdepth),
388                                  pframe.d, &ptrunk, rl+x, pdepth);
389          }
390                                                  /* sample the image */
391          for (y = 0; /* y < ysiz */ ; y++) {
392                  for (x = 0; x < xsiz-1; x++) {
393 <                        if (dev->inpready)
393 >                        if (dev->inpready || errno == ENOMEM)
394                                  goto escape;
395                          /*
396                           * Test super-pixel to the right.
# Line 302 | Line 406 | rsample()                      /* sample the image */
406                  if (y >= ysiz-1)
407                          break;
408                  for (x = 0; x < xsiz; x++) {
409 <                        if (dev->inpready)
409 >                        if (dev->inpready || errno == ENOMEM)
410                                  goto escape;
411                          /*
412                           * Find super-pixel at this position in next row.
413                           */
414                          r.l = r.d = 0;
415 <                        r.r = ourview.hresolu; r.u = ourview.vresolu;
416 <                        p = findrect(pframe.l+((x*ourview.hresolu)>>pdepth),
417 <                                pframe.d+(((y+1)*ourview.vresolu)>>pdepth),
415 >                        r.r = hresolu; r.u = vresolu;
416 >                        p = findrect(pframe.l+((x*hresolu)>>pdepth),
417 >                                pframe.d+(((y+1)*vresolu)>>pdepth),
418                                          &ptrunk, &r, pdepth);
419                          /*
420                           * Test super-pixel in next row.
# Line 329 | Line 433 | rsample()                      /* sample the image */
433                  }
434          }
435   escape:
436 <        free((char *)rl);
437 <        free((char *)pl);
436 >        free((void *)rl);
437 >        free((void *)pl);
438   }
439  
440  
441   int
442   refine(p, xmin, ymin, xmax, ymax, pd)           /* refine a node */
443 < register PNODE  *p;
443 > register PNODE  *p;
444   int  xmin, ymin, xmax, ymax;
445   int  pd;
446   {
# Line 357 | Line 461 | int  pd;
461          if (p->kid == NULL) {                   /* subdivide */
462  
463                  if ((p->kid = newptree()) == NULL)
464 <                        error(SYSTEM, "out of memory in refine");
464 >                        return(0);
465                  /*
466                   *  The following paint order can leave a black pixel
467 <                 *  when redraw() is called in (*dev->paintr)().
467 >                 *  if redraw() is called in (*dev->paintr)().
468                   */
469                  if (p->x >= mx && p->y >= my)
470                          pcopy(p, p->kid+UR);

Diff Legend

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