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.2 by greg, Wed Jun 7 08:35:29 1989 UTC vs.
Revision 2.14 by greg, Wed Dec 21 09:52:07 1994 UTC

# Line 1 | Line 1
1 < /* Copyright (c) 1987 Regents of the University of California */
1 > /* Copyright (c) 1992 Regents of the University of California */
2  
3   #ifndef lint
4   static char SCCSid[] = "$SunId$ LBL";
# Line 10 | Line 10 | static char SCCSid[] = "$SunId$ LBL";
10   *     3/24/87
11   */
12  
13 < #include  "standard.h"
13 > #include  "ray.h"
14  
15 – #include  "color.h"
16 –
15   #include  "rpaint.h"
16  
17   #include  <signal.h>
18  
19   #include  <ctype.h>
20  
21 < VIEW  ourview = STDVIEW(470);           /* viewing parameters */
21 > VIEW  ourview = STDVIEW;                /* viewing parameters */
22 > int  hresolu, vresolu;                  /* image resolution */
23  
24 + int  dimlist[MAXDIM];                   /* sampling dimensions */
25 + int  ndims = 0;                         /* number of sampling dimensions */
26 + int  samplendx = 0;                     /* index for this sample */
27 +
28   int  psample = 8;                       /* pixel sample size */
29 < double  maxdiff = .15;                  /* max. sample difference */
29 > double  maxdiff = .15;                  /* max. sample difference */
30  
31 < double  exposure = 1.0;                 /* exposure for scene */
31 > double  exposure = 1.0;                 /* exposure for scene */
32  
33 < double  dstrsrc = 0.0;                  /* square source distribution */
34 < double  shadthresh = .05;               /* shadow threshold */
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 = 128;                   /* virtual source pretest density */
38 > int  directvis = 1;                     /* sources visible? */
39 > double  srcsizerat = 0.;                /* maximum ratio source size/dist. */
40  
41 + double  specthresh = .3;                /* specular sampling threshold */
42 + double  specjitter = 1.;                /* specular sampling jitter */
43 +
44 + int  backvis = 1;                       /* back face visibility */
45 +
46   int  maxdepth = 4;                      /* maximum recursion depth */
47 < double  minweight = 1e-2;               /* minimum ray weight */
47 > double  minweight = 1e-2;               /* minimum ray weight */
48  
49   COLOR  ambval = BLKCOLOR;               /* ambient value */
50 < double  ambacc = 0.2;                   /* ambient accuracy */
51 < int  ambres = 64;                       /* ambient resolution */
50 > double  ambacc = 0.2;                   /* ambient accuracy */
51 > int  ambres = 8;                        /* ambient resolution */
52   int  ambdiv = 32;                       /* ambient divisions */
53   int  ambssamp = 0;                      /* ambient super-samples */
54   int  ambounce = 0;                      /* ambient bounces */
# Line 43 | Line 56 | char  *amblist[128];                   /* ambient include/exclude list
56   int  ambincl = -1;                      /* include == 1, exclude == 0 */
57  
58   int  greyscale = 0;                     /* map colors to brightness? */
59 < char  *devname = "X";                   /* output device name */
59 > char  *devname = dev_default;           /* output device name */
60  
61   struct driver  *dev = NULL;             /* driver functions */
62  
63 + char  rifname[128];                     /* rad input file name */
64 +
65   VIEW  oldview;                          /* previous view parameters */
66  
67   PNODE  ptrunk;                          /* the base of our image */
68   RECT  pframe;                           /* current frame boundaries */
69   int  pdepth;                            /* image depth in current frame */
70  
71 < #define  CTRL(c)        ('c'-'@')
71 > static char  *reserve_mem = NULL;       /* pre-allocated reserve memory */
72  
73 + #define RESERVE_AMT     32768           /* amount of memory to reserve */
74  
75 + #define  CTRL(c)        ((c)-'@')
76 +
77 +
78   quit(code)                      /* quit program */
79   int  code;
80   {
81 + #ifdef MSTATS
82 +        if (code == 2 && errno == ENOMEM)
83 +                printmemstats(stderr);
84 + #endif
85          devclose();
86          exit(code);
87   }
# Line 67 | Line 90 | int  code;
90   devopen(dname)                          /* open device driver */
91   char  *dname;
92   {
93 <        extern char  *progname;
94 <        char  *devargv[3];
93 >        extern char  *progname, *octname;
94 >        char  *id;
95          register int  i;
96 +
97 +        id = octname!=NULL ? octname : progname;
98                                                  /* check device table */
99          for (i = 0; devtable[i].name; i++)
100                  if (!strcmp(dname, devtable[i].name))
101 <                        if ((dev = (*devtable[i].init)(progname)) == NULL) {
101 >                        if ((dev = (*devtable[i].init)(dname, id)) == NULL) {
102                                  sprintf(errmsg, "cannot initialize %s", dname);
103                                  error(USER, errmsg);
104                          } else
105                                  return;
106 + #ifndef NIX                                                
107                                                  /* not there, try exec */
108 <        devargv[0] = dname;
83 <        devargv[1] = progname;
84 <        devargv[2] = NULL;
85 <        if ((dev = comm_init(devargv)) == NULL) {
108 >        if ((dev = comm_init(dname, id)) == NULL) {
109                  sprintf(errmsg, "cannot start device \"%s\"", dname);
110                  error(USER, errmsg);
111          }
112 + #endif
113   }
114  
115  
# Line 111 | Line 135 | rview()                                /* do a view */
135          char  buf[32];
136  
137          devopen(devname);               /* open device */
138 <        newimage();                     /* set up image */
138 >        newimage();                     /* start image (calls fillreserves) */
139  
140          for ( ; ; ) {                   /* quit in command() */
141 <                while (ourview.hresolu <= 1<<pdepth &&
118 <                                ourview.vresolu <= 1<<pdepth)
141 >                while (hresolu <= 1<<pdepth && vresolu <= 1<<pdepth)
142                          command("done: ");
143 <
144 <                if (ourview.hresolu <= psample<<pdepth &&
145 <                                ourview.vresolu <= psample<<pdepth) {
143 >                while (reserve_mem == NULL)
144 >                        command("out of memory: ");
145 >                errno = 0;
146 >                if (hresolu <= psample<<pdepth && vresolu <= psample<<pdepth) {
147                          sprintf(buf, "%d sampling...\n", 1<<pdepth);
148                          (*dev->comout)(buf);
149                          rsample();
150                  } else {
151                          sprintf(buf, "%d refining...\n", 1<<pdepth);
152                          (*dev->comout)(buf);
153 <                        refine(&ptrunk, 0, 0, ourview.hresolu,
130 <                                        ourview.vresolu, pdepth+1);
153 >                        refine(&ptrunk, 0, 0, hresolu, vresolu, pdepth+1);
154                  }
155 <                if (dev->inpready)
155 >                if (errno == ENOMEM)            /* ran out of memory */
156 >                        freereserves();
157 >                else if (dev->inpready)         /* noticed some input */
158                          command(": ");
159 <                else
159 >                else                            /* finished this depth */
160                          pdepth++;
161          }
162   }
163  
164  
165 + fillreserves()                  /* fill memory reserves */
166 + {
167 +        if (reserve_mem != NULL)
168 +                return;
169 +        reserve_mem = malloc(RESERVE_AMT);
170 + }
171 +
172 +
173 + freereserves()                  /* free memory reserves */
174 + {
175 +        if (reserve_mem == NULL)
176 +                return;
177 +        free(reserve_mem);
178 +        reserve_mem = NULL;
179 + }
180 +
181 +
182   command(prompt)                 /* get/execute command */
183   char  *prompt;
184   {
185 < #define  badcom(s)      strncmp(s, inpbuf, args-inpbuf-1)
144 <        double  atof();
185 > #define  badcom(s)      strncmp(s, inpbuf, args-inpbuf-1)
186          char  inpbuf[256];
187          char  *args;
188   again:
189 <        (*dev->comout)(prompt);                 /* get command + arguments */
149 <        (*dev->comin)(inpbuf);
189 >        (*dev->comin)(inpbuf, prompt);          /* get command + arguments */
190          for (args = inpbuf; *args && *args != ' '; args++)
191                  ;
192          if (*args) *args++ = '\0';
193          else *++args = '\0';
194          
195          switch (inpbuf[0]) {
196 <        case 'f':                               /* new frame */
197 <                if (badcom("frame"))
198 <                        goto commerr;
196 >        case 'f':                               /* new frame (or free mem.) */
197 >                if (badcom("frame")) {
198 >                        if (badcom("free"))
199 >                                goto commerr;
200 >                        free_objmem();
201 >                        break;
202 >                }
203                  getframe(args);
204                  break;
205          case 'v':                               /* view */
# Line 168 | Line 212 | again:
212                          goto commerr;
213                  lastview(args);
214                  break;
215 +        case 'V':                               /* save view */
216 +                if (badcom("V"))
217 +                        goto commerr;
218 +                saveview(args);
219 +                break;
220 +        case 'L':                               /* load view */
221 +                if (badcom("L"))
222 +                        goto commerr;
223 +                loadview(args);
224 +                break;
225          case 'e':                               /* exposure */
226                  if (badcom("exposure"))
227                          goto commerr;
228                  getexposure(args);
229                  break;
230          case 's':                               /* set a parameter */
231 <                if (badcom("set"))
231 >                if (badcom("set")) {
232 > #ifdef  SIGTSTP
233 >                        if (!badcom("stop"))
234 >                                goto dostop;
235 > #endif
236                          goto commerr;
237 +                }
238                  setparam(args);
239                  break;
240          case 'n':                               /* new picture */
# Line 193 | Line 252 | again:
252                          goto commerr;
253                  getaim(args);
254                  break;
255 <        case 'm':                               /* move camera */
255 >        case 'm':                               /* move camera (or memstats) */
256                  if (badcom("move"))
257 + #ifdef  MSTATS
258 +                {
259 +                        if (badcom("memory"))
260 +                                goto commerr;
261 +                        printmemstats(stderr);
262 +                        break;
263 +                }
264 + #else
265                          goto commerr;
266 + #endif
267                  getmove(args);
268                  break;
269 <        case 'r':                               /* rotate camera */
270 <                if (badcom("rotate"))
271 <                        goto commerr;
269 >        case 'r':                               /* rotate/repaint */
270 >                if (badcom("rotate")) {
271 >                        if (badcom("repaint")) {
272 >                                if (badcom("redraw"))
273 >                                        goto commerr;
274 >                                redraw();
275 >                                break;
276 >                        }
277 >                        getrepaint(args);
278 >                        break;
279 >                }
280                  getrotate(args);
281                  break;
282          case 'p':                               /* pivot view */
283 <                if (badcom("pivot"))
284 <                        goto commerr;
283 >                if (badcom("pivot")) {
284 >                        if (badcom("pause"))
285 >                                goto commerr;
286 >                        goto again;
287 >                }
288                  getpivot(args);
289                  break;
290 <        case CTRL(R):                           /* redraw */
290 >        case CTRL('R'):                         /* redraw */
291                  redraw();
292                  break;
293          case 'w':                               /* write */
# Line 220 | Line 299 | again:
299                  if (badcom("quit"))
300                          goto commerr;
301                  quit(0);
302 <        case CTRL(C):                           /* interrupt */
302 >        case CTRL('C'):                         /* interrupt */
303                  goto again;
304 < #ifdef  SIGTSTP
305 <        case CTRL(Z):                           /* stop */
304 > #ifdef  SIGTSTP
305 >        case CTRL('Z'):;                        /* stop */
306 > dostop:
307                  devclose();
308                  kill(0, SIGTSTP);
309                  /* pc stops here */
# Line 243 | Line 323 | commerr:
323                  error(COMMAND, errmsg);
324                  break;
325          }
326 < #undef  badcom
326 > #undef  badcom
327   }
328  
329  
# Line 253 | Line 333 | rsample()                      /* sample the image */
333          RECT  r;
334          PNODE  *p;
335          register RECT  *rl;
336 <        register PNODE  **pl;
336 >        register PNODE  **pl;
337          register int  x;
338          /*
339           *     We initialize the bottom row in the image at our current
340 <         * resolution.  During sampling, we check super-pixels to the
340 >         * resolution.  During sampling, we check super-pixels to the
341           * right and above by calling bigdiff().  If there is a significant
342           * difference, we subsample the super-pixels.  The testing process
343           * includes initialization of the next row.
344           */
345 <        xsiz = (((pframe.r-pframe.l)<<pdepth)+ourview.hresolu-1) /
346 <                        ourview.hresolu;
267 <        ysiz = (((pframe.u-pframe.d)<<pdepth)+ourview.vresolu-1) /
268 <                        ourview.vresolu;
345 >        xsiz = (((long)(pframe.r-pframe.l)<<pdepth)+hresolu-1) / hresolu;
346 >        ysiz = (((long)(pframe.u-pframe.d)<<pdepth)+vresolu-1) / vresolu;
347          rl = (RECT *)malloc(xsiz*sizeof(RECT));
348 +        if (rl == NULL)
349 +                return;
350          pl = (PNODE **)malloc(xsiz*sizeof(PNODE *));
351 <        if (rl == NULL || pl == NULL)
352 <                error(SYSTEM, "out of memory in rsample");
351 >        if (pl == NULL) {
352 >                free((char *)rl);
353 >                return;
354 >        }
355          /*
356           * Initialize the bottom row.
357           */
358          rl[0].l = rl[0].d = 0;
359 <        rl[0].r = ourview.hresolu; rl[0].u = ourview.vresolu;
359 >        rl[0].r = hresolu; rl[0].u = vresolu;
360          pl[0] = findrect(pframe.l, pframe.d, &ptrunk, rl, pdepth);
361          for (x = 1; x < xsiz; x++) {
362                  rl[x].l = rl[x].d = 0;
363 <                rl[x].r = ourview.hresolu; rl[x].u = ourview.vresolu;
364 <                pl[x] = findrect(pframe.l+((x*ourview.hresolu)>>pdepth),
363 >                rl[x].r = hresolu; rl[x].u = vresolu;
364 >                pl[x] = findrect(pframe.l+((x*hresolu)>>pdepth),
365                                  pframe.d, &ptrunk, rl+x, pdepth);
366          }
367                                                  /* sample the image */
368          for (y = 0; /* y < ysiz */ ; y++) {
369                  for (x = 0; x < xsiz-1; x++) {
370 <                        if (dev->inpready)
370 >                        if (dev->inpready || errno == ENOMEM)
371                                  goto escape;
372                          /*
373                           * Test super-pixel to the right.
# Line 301 | Line 383 | rsample()                      /* sample the image */
383                  if (y >= ysiz-1)
384                          break;
385                  for (x = 0; x < xsiz; x++) {
386 <                        if (dev->inpready)
386 >                        if (dev->inpready || errno == ENOMEM)
387                                  goto escape;
388                          /*
389                           * Find super-pixel at this position in next row.
390                           */
391                          r.l = r.d = 0;
392 <                        r.r = ourview.hresolu; r.u = ourview.vresolu;
393 <                        p = findrect(pframe.l+((x*ourview.hresolu)>>pdepth),
394 <                                pframe.d+(((y+1)*ourview.vresolu)>>pdepth),
392 >                        r.r = hresolu; r.u = vresolu;
393 >                        p = findrect(pframe.l+((x*hresolu)>>pdepth),
394 >                                pframe.d+(((y+1)*vresolu)>>pdepth),
395                                          &ptrunk, &r, pdepth);
396                          /*
397                           * Test super-pixel in next row.
# Line 335 | Line 417 | escape:
417  
418   int
419   refine(p, xmin, ymin, xmax, ymax, pd)           /* refine a node */
420 < register PNODE  *p;
420 > register PNODE  *p;
421   int  xmin, ymin, xmax, ymax;
422   int  pd;
423   {
# Line 356 | Line 438 | int  pd;
438          if (p->kid == NULL) {                   /* subdivide */
439  
440                  if ((p->kid = newptree()) == NULL)
441 <                        error(SYSTEM, "out of memory in refine");
441 >                        return(0);
442                  /*
443                   *  The following paint order can leave a black pixel
444 <                 *  when redraw() is called in (*dev->paintr)().
444 >                 *  if redraw() is called in (*dev->paintr)().
445                   */
446                  if (p->x >= mx && p->y >= my)
447                          pcopy(p, p->kid+UR);

Diff Legend

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