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 1.15 by greg, Fri Apr 5 14:26:47 1991 UTC vs.
Revision 1.19 by greg, Fri Apr 19 17:40:59 1991 UTC

# Line 36 | Line 36 | ANGLE  glarang[180] = {AEND};          /* glare calculation ang
36   int     nglarangs = 0;
37   double  maxtheta;                       /* maximum angle (in radians) */
38   int     hsize;                          /* horizontal size */
39 – int     hlim;                           /* central limit of horizontal */
39  
40   struct illum    *indirect;              /* array of indirect illuminances */
41  
# Line 48 | Line 47 | main(argc, argv)
47   int     argc;
48   char    *argv[];
49   {
50 +        int     combine = 1;
51          int     gotview = 0;
52          int     rval, i;
53          char    *err;
# Line 101 | Line 101 | char   *argv[];
101                  case 'p':
102                          picture = argv[++i];
103                          break;
104 +                case 'c':
105 +                        combine = !combine;
106 +                        break;
107                  case 'd':
108                  case 'l':
109                          rtargv[rtargc++] = argv[i];
# Line 179 | Line 182 | char   *argv[];
182          if (threshold <= FTINY)
183                  comp_thresh();                  /* compute glare threshold */
184          analyze();                              /* analyze view */
185 +        if (combine)
186 +                absorb_specks();                /* eliminate tiny sources */
187          cleanup();                              /* tidy up */
188                                                  /* print header */
189          printargs(argc, argv, stdout);
190          fputs(VIEWSTR, stdout);
191          fprintview(&ourview, stdout);
192 <        printf("\n\n");
192 >        printf("\n");
193 >        fputformat("ASCII", stdout);
194 >        printf("\n");
195          printsources();                         /* print glare sources */
196          printillum();                           /* print illuminances */
197          exit(0);
# Line 219 | Line 226 | init()                         /* initialize global variables */
226                  maxtheta = (PI/180.)*glarang[nglarangs-1];
227          else
228                  maxtheta = 0.0;
229 <        hlim = sampdens*maxtheta;
223 <        hsize = hlim + sampdens - 1;
229 >        hsize = hlim(0) + sampdens - 1;
230          if (hsize > (int)(PI*sampdens))
231                  hsize = PI*sampdens;
232          indirect = (struct illum *)calloc(nglardirs, sizeof(struct illum));
# Line 275 | Line 281 | init()                         /* initialize global variables */
281   cleanup()                               /* close files, wait for children */
282   {
283          if (verbose)
284 <                fprintf(stderr, "%s: cleaning up...\n", progname);
284 >                fprintf(stderr, "%s: cleaning up...        \n", progname);
285          if (picture != NULL)
286                  close_pict();
287          if (octree != NULL)
# Line 290 | Line 296 | compdir(vd, x, y)                      /* compute direction for x,y */
296   FVECT   vd;
297   int     x, y;
298   {
299 <        static int      cury = 10000;
294 <        static double   err, cmpval;
295 <        long    t;
299 >        int     hl;
300          FVECT   org;                    /* dummy variable */
301  
302 <        if (x <= -hlim)                 /* left region */
302 >        hl = hlim(y);
303 >        if (x <= -hl)                   /* left region */
304                  return(viewray(org, vd, &leftview,
305 <                                (double)(x+hlim)/(2*sampdens)+.5,
305 >                                (double)(x+hl)/(2*sampdens)+.5,
306                                  (double)y/(2*sampdens)+.5));
307 <        if (x >= hlim)                  /* right region */
307 >        if (x >= hl)                    /* right region */
308                  return(viewray(org, vd, &rightview,
309 <                                (double)(x-hlim)/(2*sampdens)+.5,
309 >                                (double)(x-hl)/(2*sampdens)+.5,
310                                  (double)y/(2*sampdens)+.5));
311                                          /* central region */
307 –                                /* avoid over-counting of poles */
308 –        if (cury != y) {
309 –                err = 0.0;
310 –                cmpval = sqrt(1.0 - (double)((long)y*y)/((long)vsize*vsize));
311 –                cury = y;
312 –        }
313 –        err += cmpval;
314 –        if (err <= 0.5)
315 –                return(-1);
316 –        err -= 1.0;
312          if (viewray(org, vd, &ourview, .5, (double)y/(2*sampdens)+.5) < 0)
313                  return(-1);
314 <        spinvector(vd, vd, ourview.vup, h_theta(x));
314 >        spinvector(vd, vd, ourview.vup, h_theta(x,y));
315          return(0);
316 + }
317 +
318 +
319 + double
320 + pixsize(x, y)           /* return the solid angle of pixel at (x,y) */
321 + int     x, y;
322 + {
323 +        register int    hl, xo;
324 +        double  disc;
325 +
326 +        hl = hlim(y);
327 +        if (x < -hl)
328 +                xo = x+hl;
329 +        else if (x > hl)
330 +                xo = x-hl;
331 +        else
332 +                xo = 0;
333 +        disc = 1. - (xo*xo + y*y)/(sampdens*sampdens);
334 +        if (disc <= FTINY)
335 +                return(0.);
336 +        return(1./(sampdens*sampdens*sqrt(disc)));
337   }
338  
339  

Diff Legend

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