ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/rt/rv2.c
(Generate patch)

Comparing ray/src/rt/rv2.c (file contents):
Revision 2.52 by greg, Thu Jun 16 01:55:44 2005 UTC vs.
Revision 2.65 by greg, Fri Apr 11 20:27:23 2014 UTC

# Line 13 | Line 13 | static const char      RCSid[] = "$Id$";
13   #include  <string.h>
14  
15   #include  "platform.h"
16 + #include  "rtprocess.h" /* win_popen() */
17 + #include  "paths.h"
18   #include  "ray.h"
19   #include  "source.h"
20   #include  "ambient.h"
# Line 24 | Line 26 | extern double  maxdiff;                        /* max. sample difference */
26  
27   #define  CTRL(c)        ((c)-'@')
28  
29 < #ifdef  SMLFLT
28 < #define  sscanvec(s,v)  (sscanf(s,"%f %f %f",v,v+1,v+2)==3)
29 < #else
30 < #define  sscanvec(s,v)  (sscanf(s,"%lf %lf %lf",v,v+1,v+2)==3)
31 < #endif
29 > #define  sscanvec(s,v)  (sscanf(s,FVFORMAT,v,v+1,v+2)==3)
30  
31   extern char  rifname[128];              /* rad input file name */
32  
# Line 36 | Line 34 | extern char  *progname;
34   extern char  *octname;
35  
36  
37 < extern void
37 > void
38   getframe(                               /* get a new frame */
39          char  *s
40   )
# Line 47 | Line 45 | getframe(                              /* get a new frame */
45   }
46  
47  
48 < extern void
48 > void
49   getrepaint(                             /* get area and repaint */
50          char  *s
51   )
# Line 56 | Line 54 | getrepaint(                            /* get area and repaint */
54  
55          if (getrect(s, &box) < 0)
56                  return;
57 <        paintrect(&ptrunk, 0, 0, hresolu, vresolu, &box);
57 >        paintrect(&ptrunk, &box);
58   }
59  
60  
61 < extern void
62 < getview(                                /* get/show view parameters */
61 > void
62 > getview(                                /* get/show/save view parameters */
63          char  *s
64   )
65   {
# Line 80 | Line 78 | getview(                               /* get/show view parameters */
78                          error(COMMAND, "bad view option(s)");
79                  return;
80          }
81 <        if (sscanf(s, "%s", buf) == 1) {        /* write parameters to a file */
81 >        if (nextword(buf, sizeof(buf), s) != NULL) {    /* write to a file */
82                  if ((fname = getpath(buf, NULL, 0)) == NULL ||
83                                  (fp = fopen(fname, "a")) == NULL) {
84                          sprintf(errmsg, "cannot open \"%s\"", buf);
# Line 153 | Line 151 | getview(                               /* get/show view parameters */
151   }
152  
153  
154 < extern void
154 > void
155   lastview(                               /* return to a previous view */
156          char  *s
157   )
# Line 162 | Line 160 | lastview(                              /* return to a previous view */
160          char  *fname;
161          int  success;
162          VIEW  nv;
163 <
164 <        if (sscanf(s, "%s", buf) == 1) {        /* get parameters from a file */
163 >                                        /* get parameters from a file */
164 >        if (nextword(buf, sizeof(buf), s) != NULL) {
165                  nv = stdview;
166                  if ((fname = getpath(buf, "", R_OK)) == NULL ||
167                                  (success = viewfile(fname, &nv, NULL)) == -1) {
# Line 184 | Line 182 | lastview(                              /* return to a previous view */
182          nv = ourview;
183          ourview = oldview;
184          oldview = nv;
185 <        newimage();
185 >        newimage(NULL);
186   }
187  
188  
189 < extern void
189 > void
190   saveview(                               /* save view to rad file */
191          char  *s
192   )
# Line 204 | Line 202 | saveview(                              /* save view to rad file */
202                  }
203                  s = sskip(s);
204          }
205 <        while (isspace(*s))
208 <                s++;
209 <        if (*s)
210 <                atos(rifname, sizeof(rifname), s);
211 <        else if (rifname[0] == '\0') {
205 >        if (nextword(rifname, sizeof(rifname), s) == NULL && !rifname[0]) {
206                  error(COMMAND, "no previous rad file");
207                  return;
208          }
# Line 226 | Line 220 | saveview(                              /* save view to rad file */
220   }
221  
222  
223 < extern void
223 > void
224   loadview(                               /* load view from rad file */
225          char  *s
226   )
# Line 241 | Line 235 | loadview(                              /* load view from rad file */
235                  s = sskip(s);
236          else
237                  strcat(buf, "1");
238 <        if (*s)
245 <                atos(rifname, sizeof(rifname), s);
246 <        else if (rifname[0] == '\0') {
238 >        if (nextword(rifname, sizeof(rifname), s) == NULL && !rifname[0]) {
239                  error(COMMAND, "no previous rad file");
240                  return;
241          }
# Line 269 | Line 261 | loadview(                              /* load view from rad file */
261   }
262  
263  
264 < extern void
264 > void
265   getaim(                         /* aim camera */
266          char  *s
267   )
# Line 284 | Line 276 | getaim(                                /* aim camera */
276   }
277  
278  
279 < extern void
279 > void
280   getfocus(                               /* set focus distance */
281          char *s
282   )
283   {
284 <        FVECT  vc;
284 >        char  buf[64];
285          double  dist;
286  
287          if (sscanf(s, "%lf", &dist) < 1) {
# Line 316 | Line 308 | getfocus(                              /* set focus distance */
308                  return;
309          }
310          ourview.vdist = dist;
311 +        sprintf(buf, "Focus distance set to %f\n", dist);
312 +        (*dev->comout)(buf);
313   }
314  
315  
316 < extern void
316 > void
317   getmove(                                /* move camera */
318          char  *s
319   )
# Line 333 | Line 327 | getmove(                               /* move camera */
327   }
328  
329  
330 < extern void
330 > void
331   getrotate(                              /* rotate camera */
332          char  *s
333   )
334   {
335          VIEW  nv = ourview;
342 –        FVECT  v1;
336          double  angle, elev, zfact;
337          
338          elev = 0.0; zfact = 1.0;
# Line 348 | Line 341 | getrotate(                             /* rotate camera */
341                  return;
342          }
343          spinvector(nv.vdir, ourview.vdir, ourview.vup, angle*(PI/180.));
344 <        if (elev != 0.0) {
345 <                fcross(v1, nv.vdir, ourview.vup);
346 <                normalize(v1);
354 <                spinvector(nv.vdir, nv.vdir, v1, elev*(PI/180.));
355 <        }
344 >        if (elev != 0.0)
345 >                geodesic(nv.vdir, nv.vdir, nv.vup, elev*(PI/180.), GEOD_RAD);
346 >        
347          zoomview(&nv, zfact);
348          newview(&nv);
349   }
350  
351  
352 < extern void
352 > void
353   getpivot(                               /* pivot viewpoint */
354 <        register char  *s
354 >        char  *s
355   )
356   {
357          FVECT  vc;
# Line 377 | Line 368 | getpivot(                              /* pivot viewpoint */
368   }
369  
370  
371 < extern void
371 > void
372 > getorigin(                              /* origin viewpoint */
373 >        char  *s
374 > )
375 > {
376 >        VIEW    nv = ourview;
377 >        double  d;
378 >                                        /* get new view origin */
379 >        if (!sscanvec(s, nv.vp)) {
380 >                int     x, y;
381 >                RAY     thisray;
382 >                if (dev->getcur == NULL)
383 >                        return;
384 >                (*dev->comout)("Pick point on surface for new origin\n");
385 >                if ((*dev->getcur)(&x, &y) == ABORT)
386 >                        return;
387 >                if ((thisray.rmax = viewray(thisray.rorg, thisray.rdir,
388 >                        &ourview, (x+.5)/hresolu, (y+.5)/vresolu)) < -FTINY) {
389 >                        error(COMMAND, "not on image");
390 >                        return;
391 >                }
392 >                rayorigin(&thisray, PRIMARY, NULL, NULL);
393 >                if (!localhit(&thisray, &thescene)) {
394 >                        error(COMMAND, "not a local object");
395 >                        return;
396 >                }
397 >                if (thisray.rod < 0.0)  /* don't look through other side */
398 >                        flipsurface(&thisray);
399 >                VSUM(nv.vp, thisray.rop, thisray.ron, 2.0*FTINY);
400 >                VCOPY(nv.vdir, thisray.ron);
401 >        } else if (!sscanvec(sskip2(s,3), nv.vdir) || normalize(nv.vdir) == 0.0)
402 >                VCOPY(nv.vdir, ourview.vdir);
403 >
404 >        d = DOT(nv.vdir, nv.vup);       /* need different up vector? */
405 >        if (d*d >= 1.-2.*FTINY) {
406 >                int     i;
407 >                nv.vup[0] = nv.vup[1] = nv.vup[2] = 0.0;
408 >                for (i = 3; i--; )
409 >                        if (nv.vdir[i]*nv.vdir[i] < 0.34)
410 >                                break;
411 >                nv.vup[i] = 1.;
412 >        }
413 >        newview(&nv);
414 > }
415 >
416 >
417 > void
418   getexposure(                            /* get new exposure */
419          char  *s
420   )
421   {
422          char  buf[128];
423 <        register char  *cp;
387 <        RECT  r;
423 >        char  *cp;
424          int  x, y;
425 <        register PNODE  *p = &ptrunk;
425 >        PNODE  *p = &ptrunk;
426          int  adapt = 0;
427          double  e = 1.0;
428  
# Line 403 | Line 439 | getexposure(                           /* get new exposure */
439                  (*dev->comout)("Pick point for exposure\n");
440                  if ((*dev->getcur)(&x, &y) == ABORT)
441                          return;
442 <                r.l = r.d = 0;
407 <                r.r = hresolu; r.u = vresolu;
408 <                p = findrect(x, y, &ptrunk, &r, -1);
442 >                p = findrect(x, y, &ptrunk, -1);
443          } else {
444                  if (*cp == '=') {       /* absolute setting */
445                          p = NULL;
# Line 428 | Line 462 | getexposure(                           /* get new exposure */
462                          e *= atof(cp);
463          }
464          if (p != NULL) {                /* relative setting */
465 +                compavg(p);
466                  if (bright(p->v) < 1e-15) {
467                          error(COMMAND, "cannot normalize to zero");
468                          return;
# Line 445 | Line 480 | getexposure(                           /* get new exposure */
480   }
481  
482   typedef union {int i; double d; COLOR C;}       *MyUptr;
483 + #define  FEQ(x,y)     (fabs((x)-(y)) <= FTINY)
484  
485 < extern int
485 > int
486   getparam(               /* get variable from user */
487          char  *str,
488          char  *dsc,
# Line 454 | Line 490 | getparam(              /* get variable from user */
490          void  *p
491   )
492   {
493 <        register MyUptr  ptr = (MyUptr)p;
493 >        MyUptr  ptr = (MyUptr)p;
494          int  i0;
495          double  d0, d1, d2;
496          char  buf[48];
# Line 469 | Line 505 | getparam(              /* get variable from user */
505                          if (sscanf(buf, "%d", &i0) != 1)
506                                  return(0);
507                  }
508 +                if (ptr->i == i0)
509 +                        return(0);
510                  ptr->i = i0;
511 <                return(1);
511 >                break;
512          case 'r':                       /* real */
513                  if (sscanf(str, "%lf", &d0) != 1) {
514                          (*dev->comout)(dsc);
# Line 480 | Line 518 | getparam(              /* get variable from user */
518                          if (sscanf(buf, "%lf", &d0) != 1)
519                                  return(0);
520                  }
521 +                if (FEQ(ptr->d, d0))
522 +                        return(0);
523                  ptr->d = d0;
524 <                return(1);
524 >                break;
525          case 'b':                       /* boolean */
526                  if (sscanf(str, "%1s", buf) != 1) {
527                          (*dev->comout)(dsc);
528                          sprintf(buf, "? (%c): ", ptr->i ? 'y' : 'n');
529                          (*dev->comout)(buf);
530                          (*dev->comin)(buf, NULL);
531 <                        if (buf[0] == '\0' ||
492 <                                        strchr("yY+1tTnN-0fF", buf[0]) == NULL)
531 >                        if (buf[0] == '\0')
532                                  return(0);
533                  }
534 <                ptr->i = strchr("yY+1tT", buf[0]) != NULL;
535 <                return(1);
534 >                if (strchr("yY+1tTnN-0fF", buf[0]) == NULL)
535 >                        return(0);
536 >                i0 = strchr("yY+1tT", buf[0]) != NULL;
537 >                if (ptr->i == i0)
538 >                        return(0);
539 >                ptr->i = i0;
540 >                break;
541          case 'C':                       /* color */
542                  if (sscanf(str, "%lf %lf %lf", &d0, &d1, &d2) != 3) {
543                          (*dev->comout)(dsc);
# Line 506 | Line 550 | getparam(              /* get variable from user */
550                          if (sscanf(buf, "%lf %lf %lf", &d0, &d1, &d2) != 3)
551                                  return(0);
552                  }
553 +                if (FEQ(colval(ptr->C,RED), d0) &&
554 +                                FEQ(colval(ptr->C,GRN), d1) &&
555 +                                FEQ(colval(ptr->C,BLU), d2))
556 +                        return(0);
557                  setcolor(ptr->C, d0, d1, d2);
558 <                return(1);
558 >                break;
559 >        default:
560 >                return(0);              /* shouldn't happen */
561          }
562 <        return 0; /* nothing matched */
562 >        newparam++;
563 >        return(1);
564   }
565  
566  
567 < extern void
567 > void
568   setparam(                               /* get/set program parameter */
569 <        register char  *s
569 >        char  *s
570   )
571   {
572 +        int  prev_newp = newparam;
573          char  buf[128];
574          
575          if (s[0] == '\0') {
576                  (*dev->comout)(
577 <                "aa ab ad ar as av aw b bv dc dv dj ds dt i lr lw me ma mg ms ps pt sj st u: ");
577 >                "aa ab ad ar as av aw b bv dc dv dj ds dt i lr lw me ma mg ms ps pt ss st u: ");
578                  (*dev->comin)(buf, NULL);
579                  s = buf;
580          }
# Line 583 | Line 635 | setparam(                              /* get/set program parameter */
635                  case 'n': case 'N': case 'f': case 'F': case '0': case '-':
636                          getparam(s+1, "black and white", 'b',
637                                          (void *)&greyscale);
638 +                        newparam = prev_newp;
639                          break;
640                  default:
641                          goto badparam;
# Line 665 | Line 718 | setparam(                              /* get/set program parameter */
718                  default:
719                          goto badparam;
720                  }
721 +                newparam = prev_newp;
722                  break;
723          case 's':                       /* specular */
724                  switch (s[1]) {
725 <                case 'j':                       /* jitter */
726 <                        getparam(s+2, "specular jitter", 'r',
725 >                case 's':                       /* sampling */
726 >                        getparam(s+2, "specular sampling", 'r',
727                                          (void *)&specjitter);
728                          break;
729                  case 't':                       /* threshold */
# Line 692 | Line 746 | badparam:
746   }
747  
748  
749 < extern void
750 < traceray(s)                             /* trace a single ray */
751 < char  *s;
749 > void
750 > traceray(                               /* trace a single ray */
751 >        char  *s
752 > )
753   {
754 <        char  buf[128];
755 <        int  x, y;
701 <        OBJREC  *ino;
702 <        RAY  thisray;
754 >        RAY     thisray;
755 >        char    buf[512];
756  
757          thisray.rmax = 0.0;
758  
759          if (!sscanvec(s, thisray.rorg) ||
760                          !sscanvec(sskip2(s,3), thisray.rdir)) {
761 +                int  x, y;
762  
763                  if (dev->getcur == NULL)
764                          return;
# Line 723 | Line 777 | char  *s;
777                  return;
778          }
779  
780 <        rayorigin(&thisray, PRIMARY, NULL, NULL);
727 <        
728 <        rayvalue(&thisray);
780 >        ray_trace(&thisray);
781  
782          if (thisray.ro == NULL)
783                  (*dev->comout)("ray hit nothing");
# Line 733 | Line 785 | char  *s;
785                  OBJREC  *mat = NULL;
786                  OBJREC  *mod = NULL;
787                  char    matspec[256];
788 +                OBJREC  *ino;
789 +
790                  matspec[0] = '\0';
791                  if (thisray.ro->omod != OVOID) {
792                          mod = objptr(thisray.ro->omod);
# Line 775 | Line 829 | char  *s;
829   }
830  
831  
832 < extern void
833 < writepict(s)                            /* write the picture to a file */
834 < char  *s;
832 > void
833 > writepict(                              /* write the picture to a file */
834 >        char  *s
835 > )
836   {
837          static char  buf[128];
838          char  *fname;
839          FILE  *fp;
840          COLR  *scanline;
841          int  y;
842 <
843 <        while (isspace(*s))
789 <                s++;
790 <        if (*s)
791 <                atos(buf, sizeof(buf), s);
792 <        else if (buf[0] == '\0') {
842 >                                /* XXX relies on words.c 2.11 behavior */
843 >        if (nextword(buf, sizeof(buf), s) == NULL && !buf[0]) {
844                  error(COMMAND, "no file");
845                  return;
846          }

Diff Legend

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