55 |
|
progname = argv[0]; |
56 |
|
/* process options */ |
57 |
|
for (i = 1; i < argc && argv[i][0] == '-'; i++) { |
58 |
+ |
/* expand arguments */ |
59 |
+ |
while (rval = expandarg(&argc, &argv, i)) |
60 |
+ |
if (rval < 0) { |
61 |
+ |
fprintf(stderr, "%s: cannot expand '%s'", |
62 |
+ |
argv[0], argv[i]); |
63 |
+ |
exit(1); |
64 |
+ |
} |
65 |
|
rval = getviewopt(&ourview, argc-i, argv+i); |
66 |
|
if (rval >= 0) { |
67 |
|
i += rval; |
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); |
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 |
|
} |