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.13 by greg, Fri May 3 15:43:32 1991 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  <setjmp.h>
19 > CUBE  thescene;                         /* our scene */
20 > OBJECT  nsceneobjs;                     /* number of objects in our scene */
21  
22 < #include  <ctype.h>
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 = 8;                        /* 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 = dev_default;           /* 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  
60 – static jmp_buf  mainloop;               /* longjmp back to main loop */
86   static char  *reserve_mem = NULL;       /* pre-allocated reserve memory */
87  
88 < #define RESERVE_AMT     8192            /* amount of memory to reserve */
88 > #define RESERVE_AMT     32768           /* amount of memory to reserve */
89  
90 < #define  CTRL(c)        ('c'-'@')
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   {
# Line 83 | Line 114 | char  *dname;
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))
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          if ((dev = comm_init(dname, id)) == NULL) {
126                  sprintf(errmsg, "cannot start device \"%s\"", dname);
# Line 97 | Line 129 | char  *dname;
129   }
130  
131  
132 + void
133   devclose()                              /* close our device */
134   {
135          if (dev != NULL)
# Line 105 | Line 138 | devclose()                             /* close our device */
138   }
139  
140  
141 + void
142   printdevices()                          /* print list of output devices */
143   {
144          register int  i;
# Line 114 | 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 */
158 <        setjmp(mainloop);
156 >        devopen(dvcname);               /* open device */
157 >        newimage();                     /* start image (calls fillreserves) */
158 >
159          for ( ; ; ) {                   /* quit in command() */
160 <                while (hresolu <= 1<<pdepth &&
126 <                                vresolu <= 1<<pdepth)
160 >                while (hresolu <= 1<<pdepth && vresolu <= 1<<pdepth)
161                          command("done: ");
162 <
163 <                if (hresolu <= psample<<pdepth &&
164 <                                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();
# Line 136 | Line 171 | rview()                                /* do a view */
171                          (*dev->comout)(buf);
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 < memreserve()                    /* fill memory reserves */
184 > void
185 > fillreserves()                  /* fill memory reserves */
186   {
187          if (reserve_mem != NULL)
188 <                return;                 /* got some already */
189 <        reserve_mem = malloc(RESERVE_AMT);
188 >                return;
189 >        reserve_mem = (char *)malloc(RESERVE_AMT);
190   }
191  
192  
193 < memerror(detail)                /* try and rescue a memory error */
194 < char    *detail;
193 > void
194 > freereserves()                  /* free memory reserves */
195   {
196 <        if (reserve_mem == NULL) {
197 <                sprintf(errmsg, "out of memory %s", detail);
160 <                error(SYSTEM, errmsg);
161 <        }
196 >        if (reserve_mem == NULL)
197 >                return;
198          free(reserve_mem);
199          reserve_mem = NULL;
164 –        do
165 –                command("out of memory: ");
166 –        while (reserve_mem == NULL);
167 –        longjmp(mainloop, 1);
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)
175 <        double  atof();
207 > #define  badcom(s)      strncmp(s, inpbuf, args-inpbuf-1)
208          char  inpbuf[256];
209          char  *args;
210   again:
# Line 183 | Line 215 | again:
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 198 | 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 223 | 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/repaint */
292                  if (badcom("rotate")) {
293 <                        if (badcom("repaint"))
294 <                                goto commerr;
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 254 | 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 277 | 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)+hresolu-1) / hresolu;
369 <        ysiz = (((pframe.u-pframe.d)<<pdepth)+vresolu-1) / 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 <                memerror("in rsample");
372 >                return;
373          pl = (PNODE **)malloc(xsiz*sizeof(PNODE *));
374 <        if (pl == NULL)
375 <                memerror("in rsample");
374 >        if (pl == NULL) {
375 >                free((void *)rl);
376 >                return;
377 >        }
378          /*
379           * Initialize the bottom row.
380           */
# Line 319 | Line 390 | rsample()                      /* sample the image */
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 335 | 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.
# Line 362 | 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 390 | Line 461 | int  pd;
461          if (p->kid == NULL) {                   /* subdivide */
462  
463                  if ((p->kid = newptree()) == NULL)
464 <                        memerror("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