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 2.18 by greg, Sat Feb 22 02:07:29 2003 UTC vs.
Revision 2.37 by greg, Fri Oct 18 17:04:13 2013 UTC

# Line 7 | Line 7 | static const char      RCSid[] = "$Id$";
7   *  External symbols declared in rpaint.h
8   */
9  
10 < /* ====================================================================
11 < * The Radiance Software License, Version 1.0
12 < *
13 < * Copyright (c) 1990 - 2002 The Regents of the University of California,
14 < * through Lawrence Berkeley National Laboratory.   All rights reserved.
15 < *
16 < * Redistribution and use in source and binary forms, with or without
17 < * modification, are permitted provided that the following conditions
18 < * are met:
19 < *
20 < * 1. Redistributions of source code must retain the above copyright
21 < *         notice, this list of conditions and the following disclaimer.
22 < *
23 < * 2. Redistributions in binary form must reproduce the above copyright
24 < *       notice, this list of conditions and the following disclaimer in
25 < *       the documentation and/or other materials provided with the
26 < *       distribution.
27 < *
28 < * 3. The end-user documentation included with the redistribution,
29 < *           if any, must include the following acknowledgment:
30 < *             "This product includes Radiance software
31 < *                 (http://radsite.lbl.gov/)
32 < *                 developed by the Lawrence Berkeley National Laboratory
33 < *               (http://www.lbl.gov/)."
34 < *       Alternately, this acknowledgment may appear in the software itself,
35 < *       if and wherever such third-party acknowledgments normally appear.
36 < *
37 < * 4. The names "Radiance," "Lawrence Berkeley National Laboratory"
38 < *       and "The Regents of the University of California" must
39 < *       not be used to endorse or promote products derived from this
40 < *       software without prior written permission. For written
41 < *       permission, please contact [email protected].
42 < *
43 < * 5. Products derived from this software may not be called "Radiance",
44 < *       nor may "Radiance" appear in their name, without prior written
45 < *       permission of Lawrence Berkeley National Laboratory.
46 < *
47 < * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESSED OR IMPLIED
48 < * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
49 < * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
50 < * DISCLAIMED.   IN NO EVENT SHALL Lawrence Berkeley National Laboratory OR
51 < * ITS CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
52 < * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
53 < * LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF
54 < * USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
55 < * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
56 < * OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT
57 < * OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
58 < * SUCH DAMAGE.
59 < * ====================================================================
60 < *
61 < * This software consists of voluntary contributions made by many
62 < * individuals on behalf of Lawrence Berkeley National Laboratory.   For more
63 < * information on Lawrence Berkeley National Laboratory, please see
64 < * <http://www.lbl.gov/>.
65 < */
10 > #include "copyright.h"
11  
67 – #include  "ray.h"
68 –
69 – #include  "rpaint.h"
70 –
12   #include  <signal.h>
72 –
13   #include  <ctype.h>
14  
15 < CUBE  thescene;                         /* our scene */
16 < OBJECT  nsceneobjs;                     /* number of objects in our scene */
15 > #include  "ray.h"
16 > #include  "rpaint.h"
17  
78 – int  dimlist[MAXDIM];                   /* sampling dimensions */
79 – int  ndims = 0;                         /* number of sampling dimensions */
80 – int  samplendx = 0;                     /* index for this sample */
81 –
82 – extern void  ambnotify();
83 – void  (*addobjnotify[])() = {ambnotify, NULL};
84 –
85 – VIEW  ourview = STDVIEW;                /* viewing parameters */
86 – int  hresolu, vresolu;                  /* image resolution */
87 –
88 – void  (*trace)() = NULL;                /* trace call */
89 –
90 – int  do_irrad = 0;                      /* compute irradiance? */
91 –
92 – int  psample = 8;                       /* pixel sample size */
93 – double  maxdiff = .15;                  /* max. sample difference */
94 –
95 – double  exposure = 1.0;                 /* exposure for scene */
96 –
97 – double  dstrsrc = 0.0;                  /* square source distribution */
98 – double  shadthresh = .1;                /* shadow threshold */
99 – double  shadcert = .25;                 /* shadow certainty */
100 – int  directrelay = 0;                   /* number of source relays */
101 – int  vspretest = 128;                   /* virtual source pretest density */
102 – int  directvis = 1;                     /* sources visible? */
103 – double  srcsizerat = 0.;                /* maximum ratio source size/dist. */
104 –
105 – COLOR  cextinction = BLKCOLOR;          /* global extinction coefficient */
106 – COLOR  salbedo = BLKCOLOR;              /* global scattering albedo */
107 – double  seccg = 0.;                     /* global scattering eccentricity */
108 – double  ssampdist = 0.;                 /* scatter sampling distance */
109 –
110 – double  specthresh = .3;                /* specular sampling threshold */
111 – double  specjitter = 1.;                /* specular sampling jitter */
112 –
113 – int  backvis = 1;                       /* back face visibility */
114 –
115 – int  maxdepth = 4;                      /* maximum recursion depth */
116 – double  minweight = 1e-2;               /* minimum ray weight */
117 –
118 – char  *ambfile = NULL;                  /* ambient file name */
119 – COLOR  ambval = BLKCOLOR;               /* ambient value */
120 – int  ambvwt = 0;                        /* initial weight for ambient value */
121 – double  ambacc = 0.2;                   /* ambient accuracy */
122 – int  ambres = 8;                        /* ambient resolution */
123 – int  ambdiv = 32;                       /* ambient divisions */
124 – int  ambssamp = 0;                      /* ambient super-samples */
125 – int  ambounce = 0;                      /* ambient bounces */
126 – char  *amblist[128];                    /* ambient include/exclude list */
127 – int  ambincl = -1;                      /* include == 1, exclude == 0 */
128 –
129 – int  greyscale = 0;                     /* map colors to brightness? */
130 – char  *dvcname = dev_default;           /* output device name */
131 –
132 – struct driver  *dev = NULL;             /* driver functions */
133 –
134 – char  rifname[128];                     /* rad input file name */
135 –
136 – VIEW  oldview;                          /* previous view parameters */
137 –
138 – PNODE  ptrunk;                          /* the base of our image */
139 – RECT  pframe;                           /* current frame boundaries */
140 – int  pdepth;                            /* image depth in current frame */
141 –
142 – static char  *reserve_mem = NULL;       /* pre-allocated reserve memory */
143 –
144 – #define RESERVE_AMT     32768           /* amount of memory to reserve */
145 –
18   #define  CTRL(c)        ((c)-'@')
19  
20  
# Line 150 | Line 22 | void
22   quit(code)                      /* quit program */
23   int  code;
24   {
25 < #ifdef MSTATS
26 <        if (code == 2 && errno == ENOMEM)
27 <                printmemstats(stderr);
28 < #endif
157 <        devclose();
25 >        if (ray_pnprocs > 0)    /* close children if any */
26 >                ray_pclose(0);
27 >        if (!ray_pnprocs)       /* in parent */
28 >                devclose();
29          exit(code);
30   }
31  
32  
33   void
34 < devopen(dname)                          /* open device driver */
35 < char  *dname;
34 > devopen(                                /* open device driver */
35 >        char  *dname
36 > )
37   {
38          extern char  *progname, *octname;
39          char  *id;
40 <        register int  i;
40 >        int  i;
41  
42          id = octname!=NULL ? octname : progname;
43                                                  /* check device table */
44          for (i = 0; devtable[i].name; i++)
45 <                if (!strcmp(dname, devtable[i].name))
45 >                if (!strcmp(dname, devtable[i].name)) {
46                          if ((dev = (*devtable[i].init)(dname, id)) == NULL) {
47                                  sprintf(errmsg, "cannot initialize %s", dname);
48                                  error(USER, errmsg);
49                          } else
50                                  return;
51 < #ifndef NIX                                                
51 >                }
52                                                  /* not there, try exec */
53          if ((dev = comm_init(dname, id)) == NULL) {
54                  sprintf(errmsg, "cannot start device \"%s\"", dname);
55                  error(USER, errmsg);
56          }
185 – #endif
57   }
58  
59  
60   void
61 < devclose()                              /* close our device */
61 > devclose(void)                          /* close our device */
62   {
63          if (dev != NULL)
64                  (*dev->close)();
# Line 196 | Line 67 | devclose()                             /* close our device */
67  
68  
69   void
70 < printdevices()                          /* print list of output devices */
70 > printdevices(void)                      /* print list of output devices */
71   {
72 <        register int  i;
72 >        int  i;
73  
74          for (i = 0; devtable[i].name; i++)
75                  printf("%-16s # %s\n", devtable[i].name, devtable[i].descrip);
# Line 206 | Line 77 | printdevices()                         /* print list of output devices */
77  
78  
79   void
80 < rview()                         /* do a view */
80 > rview(void)                             /* do a view */
81   {
82          char  buf[32];
83  
84          devopen(dvcname);               /* open device */
85 <        newimage();                     /* start image (calls fillreserves) */
85 >        newimage(NULL);                 /* start image */
86  
87          for ( ; ; ) {                   /* quit in command() */
88                  while (hresolu <= 1<<pdepth && vresolu <= 1<<pdepth)
89                          command("done: ");
219 –                while (reserve_mem == NULL)
220 –                        command("out of memory: ");
90                  errno = 0;
91                  if (hresolu <= psample<<pdepth && vresolu <= psample<<pdepth) {
92                          sprintf(buf, "%d sampling...\n", 1<<pdepth);
# Line 226 | Line 95 | rview()                                /* do a view */
95                  } else {
96                          sprintf(buf, "%d refining...\n", 1<<pdepth);
97                          (*dev->comout)(buf);
98 <                        refine(&ptrunk, 0, 0, hresolu, vresolu, pdepth+1);
98 >                        refine(&ptrunk, pdepth+1);
99                  }
100 <                if (errno == ENOMEM)            /* ran out of memory */
232 <                        freereserves();
233 <                else if (dev->inpready)         /* noticed some input */
100 >                if (dev->inpready)              /* noticed some input */
101                          command(": ");
102                  else                            /* finished this depth */
103                          pdepth++;
# Line 239 | Line 106 | rview()                                /* do a view */
106  
107  
108   void
109 < fillreserves()                  /* fill memory reserves */
109 > command(                        /* get/execute command */
110 >        char  *prompt
111 > )
112   {
244 –        if (reserve_mem != NULL)
245 –                return;
246 –        reserve_mem = (char *)malloc(RESERVE_AMT);
247 – }
248 –
249 –
250 – void
251 – freereserves()                  /* free memory reserves */
252 – {
253 –        if (reserve_mem == NULL)
254 –                return;
255 –        free(reserve_mem);
256 –        reserve_mem = NULL;
257 – }
258 –
259 –
260 – void
261 – command(prompt)                 /* get/execute command */
262 – char  *prompt;
263 – {
113   #define  badcom(s)      strncmp(s, inpbuf, args-inpbuf-1)
114          char  inpbuf[256];
115          char  *args;
# Line 270 | Line 119 | again:
119                  ;
120          if (*args) *args++ = '\0';
121          else *++args = '\0';
122 +
123 +        if (waitrays() < 0)                     /* clear ray queue */
124 +                quit(1);
125          
126          switch (inpbuf[0]) {
127 <        case 'f':                               /* new frame (or free mem.) */
127 >        case 'f':                               /* new frame (|focus|free) */
128                  if (badcom("frame")) {
129 <                        if (badcom("free"))
130 <                                goto commerr;
131 <                        free_objmem();
129 >                        if (badcom("focus")) {
130 >                                if (badcom("free"))
131 >                                        goto commerr;
132 >                                free_objmem();
133 >                                break;
134 >                        }
135 >                        getfocus(args);
136                          break;
137                  }
138                  getframe(args);
# Line 319 | Line 175 | again:
175          case 'n':                               /* new picture */
176                  if (badcom("new"))
177                          goto commerr;
178 <                newimage();
178 >                newimage(args);
179                  break;
180          case 't':                               /* trace a ray */
181                  if (badcom("trace"))
# Line 333 | Line 189 | again:
189                  break;
190          case 'm':                               /* move camera (or memstats) */
191                  if (badcom("move"))
336 – #ifdef  MSTATS
337 –                {
338 –                        if (badcom("memory"))
339 –                                goto commerr;
340 –                        printmemstats(stderr);
341 –                        break;
342 –                }
343 – #else
192                          goto commerr;
345 – #endif
193                  getmove(args);
194                  break;
195          case 'r':                               /* rotate/repaint */
# Line 366 | Line 213 | again:
213                  }
214                  getpivot(args);
215                  break;
216 +        case 'o':                               /* origin view */
217 +                if (badcom("origin"))
218 +                        goto commerr;
219 +                getorigin(args);
220 +                break;
221          case CTRL('R'):                         /* redraw */
222                  redraw();
223                  break;
# Line 402 | Line 254 | commerr:
254                  error(COMMAND, errmsg);
255                  break;
256          }
257 +        if (newparam && ray_pnprocs)            /* drop into immediate mode */
258 +                ray_pclose(0);
259   #undef  badcom
260   }
261  
262  
263   void
264 < rsample()                       /* sample the image */
264 > rsample(void)                   /* sample the image */
265   {
266          int  xsiz, ysiz, y;
413 –        RECT  r;
267          PNODE  *p;
268 <        register RECT  *rl;
269 <        register PNODE  **pl;
417 <        register int  x;
268 >        PNODE   **pl;
269 >        int  x;
270          /*
271           *     We initialize the bottom row in the image at our current
272           * resolution.  During sampling, we check super-pixels to the
# Line 424 | Line 276 | rsample()                      /* sample the image */
276           */
277          xsiz = (((long)(pframe.r-pframe.l)<<pdepth)+hresolu-1) / hresolu;
278          ysiz = (((long)(pframe.u-pframe.d)<<pdepth)+vresolu-1) / vresolu;
427 –        rl = (RECT *)malloc(xsiz*sizeof(RECT));
428 –        if (rl == NULL)
429 –                return;
279          pl = (PNODE **)malloc(xsiz*sizeof(PNODE *));
280 <        if (pl == NULL) {
432 <                free((void *)rl);
280 >        if (pl == NULL)
281                  return;
434 –        }
282          /*
283           * Initialize the bottom row.
284           */
285 <        rl[0].l = rl[0].d = 0;
439 <        rl[0].r = hresolu; rl[0].u = vresolu;
440 <        pl[0] = findrect(pframe.l, pframe.d, &ptrunk, rl, pdepth);
285 >        pl[0] = findrect(pframe.l, pframe.d, &ptrunk, pdepth);
286          for (x = 1; x < xsiz; x++) {
442 –                rl[x].l = rl[x].d = 0;
443 –                rl[x].r = hresolu; rl[x].u = vresolu;
287                  pl[x] = findrect(pframe.l+((x*hresolu)>>pdepth),
288 <                                pframe.d, &ptrunk, rl+x, pdepth);
288 >                                pframe.d, &ptrunk, pdepth);
289          }
290                                                  /* sample the image */
291          for (y = 0; /* y < ysiz */ ; y++) {
292                  for (x = 0; x < xsiz-1; x++) {
293 <                        if (dev->inpready || errno == ENOMEM)
293 >                        if (dev->inpready)
294                                  goto escape;
295                          /*
296                           * Test super-pixel to the right.
297                           */
298                          if (pl[x] != pl[x+1] && bigdiff(pl[x]->v,
299                                          pl[x+1]->v, maxdiff)) {
300 <                                refine(pl[x], rl[x].l, rl[x].d,
301 <                                                rl[x].r, rl[x].u, 1);
459 <                                refine(pl[x+1], rl[x+1].l, rl[x+1].d,
460 <                                                rl[x+1].r, rl[x+1].u, 1);
300 >                                refine(pl[x], 1);
301 >                                refine(pl[x+1], 1);
302                          }
303                  }
304                  if (y >= ysiz-1)
305                          break;
306                  for (x = 0; x < xsiz; x++) {
307 <                        if (dev->inpready || errno == ENOMEM)
307 >                        if (dev->inpready)
308                                  goto escape;
309                          /*
310                           * Find super-pixel at this position in next row.
311                           */
471 –                        r.l = r.d = 0;
472 –                        r.r = hresolu; r.u = vresolu;
312                          p = findrect(pframe.l+((x*hresolu)>>pdepth),
313                                  pframe.d+(((y+1)*vresolu)>>pdepth),
314 <                                        &ptrunk, &r, pdepth);
314 >                                        &ptrunk, pdepth);
315                          /*
316                           * Test super-pixel in next row.
317                           */
318                          if (pl[x] != p && bigdiff(pl[x]->v, p->v, maxdiff)) {
319 <                                refine(pl[x], rl[x].l, rl[x].d,
320 <                                                rl[x].r, rl[x].u, 1);
482 <                                refine(p, r.l, r.d, r.r, r.u, 1);
319 >                                refine(pl[x], 1);
320 >                                refine(p, 1);
321                          }
322                          /*
323                           * Copy into super-pixel array.
324                           */
487 –                        rl[x].l = r.l; rl[x].d = r.d;
488 –                        rl[x].r = r.r; rl[x].u = r.u;
325                          pl[x] = p;
326                  }
327          }
328   escape:
493 –        free((void *)rl);
329          free((void *)pl);
330   }
331  
332  
333   int
334 < refine(p, xmin, ymin, xmax, ymax, pd)           /* refine a node */
335 < register PNODE  *p;
336 < int  xmin, ymin, xmax, ymax;
337 < int  pd;
334 > refine(                         /* refine a node */
335 >        PNODE   *p,
336 >        int  pd
337 > )
338   {
339          int  growth;
340          int  mx, my;
506 –        int  i;
341  
342          if (dev->inpready)                      /* quit for input */
343                  return(0);
# Line 511 | Line 345 | int  pd;
345          if (pd <= 0)                            /* depth limit */
346                  return(0);
347  
348 <        mx = (xmin + xmax) >> 1;
349 <        my = (ymin + ymax) >> 1;
348 >        mx = (p->xmin + p->xmax) >> 1;
349 >        my = (p->ymin + p->ymax) >> 1;
350          growth = 0;
351  
352          if (p->kid == NULL) {                   /* subdivide */
353  
354                  if ((p->kid = newptree()) == NULL)
355                          return(0);
356 +
357 +                p->kid[UR].xmin = mx;
358 +                p->kid[UR].ymin = my;
359 +                p->kid[UR].xmax = p->xmax;
360 +                p->kid[UR].ymax = p->ymax;
361 +                p->kid[UL].xmin = p->xmin;
362 +                p->kid[UL].ymin = my;
363 +                p->kid[UL].xmax = mx;
364 +                p->kid[UL].ymax = p->ymax;
365 +                p->kid[DR].xmin = mx;
366 +                p->kid[DR].ymin = p->ymin;
367 +                p->kid[DR].xmax = p->xmax;
368 +                p->kid[DR].ymax = my;
369 +                p->kid[DL].xmin = p->xmin;
370 +                p->kid[DL].ymin = p->ymin;
371 +                p->kid[DL].xmax = mx;
372 +                p->kid[DL].ymax = my;
373                  /*
374                   *  The following paint order can leave a black pixel
375                   *  if redraw() is called in (*dev->paintr)().
376                   */
377                  if (p->x >= mx && p->y >= my)
378                          pcopy(p, p->kid+UR);
379 <                else
380 <                        paint(p->kid+UR, mx, my, xmax, ymax);
379 >                else if (paint(p->kid+UR) < 0)
380 >                        quit(1);
381                  if (p->x < mx && p->y >= my)
382                          pcopy(p, p->kid+UL);
383 <                else
384 <                        paint(p->kid+UL, xmin, my, mx, ymax);
383 >                else if (paint(p->kid+UL) < 0)
384 >                        quit(1);
385                  if (p->x >= mx && p->y < my)
386                          pcopy(p, p->kid+DR);
387 <                else
388 <                        paint(p->kid+DR, mx, ymin, xmax, my);
387 >                else if (paint(p->kid+DR) < 0)
388 >                        quit(1);
389                  if (p->x < mx && p->y < my)
390                          pcopy(p, p->kid+DL);
391 <                else
392 <                        paint(p->kid+DL, xmin, ymin, mx, my);
391 >                else if (paint(p->kid+DL) < 0)
392 >                        quit(1);
393  
394                  growth++;
395          }
396                                                  /* do children */
397          if (mx > pframe.l) {
398                  if (my > pframe.d)
399 <                        growth += refine(p->kid+DL, xmin, ymin, mx, my, pd-1);
399 >                        growth += refine(p->kid+DL, pd-1);
400                  if (my < pframe.u)
401 <                        growth += refine(p->kid+UL, xmin, my, mx, ymax, pd-1);
401 >                        growth += refine(p->kid+UL, pd-1);
402          }
403          if (mx < pframe.r) {
404                  if (my > pframe.d)
405 <                        growth += refine(p->kid+DR, mx, ymin, xmax, my, pd-1);
405 >                        growth += refine(p->kid+DR, pd-1);
406                  if (my < pframe.u)
407 <                        growth += refine(p->kid+UR, mx, my, xmax, ymax, pd-1);
557 <        }
558 <                                                /* recompute sum */
559 <        if (growth) {
560 <                setcolor(p->v, 0.0, 0.0, 0.0);
561 <                for (i = 0; i < 4; i++)
562 <                        addcolor(p->v, p->kid[i].v);
563 <                scalecolor(p->v, 0.25);
407 >                        growth += refine(p->kid+UR, pd-1);
408          }
409          return(growth);
410   }

Diff Legend

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