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

Comparing ray/src/util/findglare.c (file contents):
Revision 2.5 by greg, Wed Feb 3 17:55:34 1993 UTC vs.
Revision 2.10 by schorsch, Mon Jul 21 22:30:19 2003 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1991 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   * Find glare sources in a scene or image.
6   *
# Line 16 | Line 13 | static char SCCSid[] = "$SunId$ LBL";
13   #define VEQ(v1,v2)      (FEQ((v1)[0],(v2)[0])&&FEQ((v1)[1],(v2)[1]) \
14                                  &&FEQ((v1)[2],(v2)[2]))
15  
16 < char    *rtargv[32] = {"rtrace", "-h-", "-ov", "-fff"};
17 < int     rtargc = 4;
16 > char    *rtargv[64] = {"rtrace", "-h-", "-ov", "-fff", "-ld-", "-i-", "-I-"};
17 > int     rtargc = 7;
18  
19   VIEW    ourview = STDVIEW;              /* our view */
20   VIEW    pictview = STDVIEW;             /* picture view */
# Line 55 | Line 52 | char   *argv[];
52          progname = argv[0];
53                                          /* process options */
54          for (i = 1; i < argc && argv[i][0] == '-'; i++) {
55 +                                                /* expand arguments */
56 +                while ((rval = expandarg(&argc, &argv, i)) > 0)
57 +                        ;
58 +                if (rval < 0) {
59 +                        fprintf(stderr, "%s: cannot expand '%s'",
60 +                                        argv[0], argv[i]);
61 +                        exit(1);
62 +                }
63                  rval = getviewopt(&ourview, argc-i, argv+i);
64                  if (rval >= 0) {
65                          i += rval;
# Line 105 | Line 110 | char   *argv[];
110                          combine = !combine;
111                          break;
112                  case 'd':
113 <                        if (argv[i][2] == 'v') {
114 <                                rtargv[rtargc++] = argv[i];
113 >                        rtargv[rtargc++] = argv[i];
114 >                        if (argv[i][2] != 'v')
115 >                                rtargv[rtargc++] = argv[++i];
116 >                        break;
117 >                case 'l':
118 >                        if (argv[i][2] == 'd')
119                                  break;
111 –                        }
120                          /* FALL THROUGH */
113 –                case 'l':
121                  case 's':
122                  case 'P':
123                          rtargv[rtargc++] = argv[i];
# Line 167 | Line 174 | char   *argv[];
174                                          progname);
175                          exit(1);
176                  }
177 <                copystruct(&ourview, &pictview);
177 >                ourview = pictview;
178          } else if (picture != NULL && !VEQ(ourview.vp, pictview.vp)) {
179                  fprintf(stderr, "%s: picture must have same viewpoint\n",
180                                  progname);
# Line 196 | Line 203 | char   *argv[];
203                  absorb_specks();                /* eliminate tiny sources */
204          cleanup();                              /* tidy up */
205                                                  /* print header */
206 +        newheader("RADIANCE", stdout);
207          printargs(argc, argv, stdout);
208          fputs(VIEWSTR, stdout);
209          fprintview(&ourview, stdout);
# Line 261 | Line 269 | init()                         /* initialize global variables */
269          if (indirect == NULL)
270                  memerr("indirect illuminances");
271          npixinvw = npixmiss = 0L;
272 <        copystruct(&leftview, &ourview);
273 <        copystruct(&rightview, &ourview);
272 >        leftview = ourview;
273 >        rightview = ourview;
274          spinvector(leftview.vdir, ourview.vdir, ourview.vup, maxtheta);
275          spinvector(rightview.vdir, ourview.vdir, ourview.vup, -maxtheta);
276          setview(&leftview);
# Line 364 | Line 372 | int    x, y;
372                  xo = x-hl;
373          else
374                  xo = 0;
375 <        disc = 1. - (double)(xo*xo + y*y)/(sampdens*sampdens);
376 <        if (disc <= FTINY)
375 >        disc = 1. - (double)((long)xo*xo + (long)y*y)/((long)sampdens*sampdens);
376 >        if (disc <= FTINY*FTINY)
377                  return(0.);
378          return(1./(sampdens*sampdens*sqrt(disc)));
379   }

Diff Legend

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