14 |
|
#include <math.h> |
15 |
|
#include "rtio.h" |
16 |
|
#include "platform.h" |
17 |
+ |
#include "paths.h" |
18 |
|
#include "bsdfrep.h" |
19 |
|
/* global argv[0] */ |
20 |
|
char *progname; |
33 |
|
|
34 |
|
int rev_orient = 0; /* shall we reverse surface orientation? */ |
35 |
|
|
36 |
+ |
double lim_graze = 0; /* limit scattering near grazing (deg above) */ |
37 |
+ |
|
38 |
|
/* Compare incident angles */ |
39 |
|
static int |
40 |
|
cmp_indir(const void *p1, const void *p2) |
181 |
|
fclose(fp); |
182 |
|
return(1); |
183 |
|
} |
184 |
< |
if (rev_orient) { /* reverse Z-axis to face outside */ |
184 |
> |
/* check if scatter angle is too low */ |
185 |
> |
if (fabs(theta_out - 90.) < lim_graze-FTINY) |
186 |
> |
continue; |
187 |
> |
if (rev_orient) { /* reverse Z-axis to face outside? */ |
188 |
|
theta_out = 180. - theta_out; |
189 |
|
phi_out = 360. - phi_out; |
190 |
|
} |
230 |
|
main(int argc, char *argv[]) |
231 |
|
{ |
232 |
|
extern int nprocs; |
233 |
+ |
static char gval_buf[16]; |
234 |
+ |
char * auto_grazing = NULL; |
235 |
|
const char *symmetry = "0"; |
236 |
|
int ninpfiles, totinc; |
237 |
|
int a, i; |
248 |
|
case 's': |
249 |
|
symmetry = argv[++a]; |
250 |
|
break; |
251 |
+ |
case 'g': |
252 |
+ |
if (toupper(argv[a+1][0]) == 'A') |
253 |
+ |
auto_grazing = argv[a+1] = gval_buf; |
254 |
+ |
else |
255 |
+ |
lim_graze = atof(argv[a+1]); |
256 |
+ |
++a; |
257 |
+ |
break; |
258 |
|
default: |
259 |
|
goto userr; |
260 |
|
} |
266 |
|
inpfile = (PGINPUT *)malloc(sizeof(PGINPUT)*totinc); |
267 |
|
if (inpfile == NULL) |
268 |
|
return(1); |
269 |
< |
for (i = 0; i < ninpfiles; i++) |
269 |
> |
if (auto_grazing) |
270 |
> |
lim_graze = 90.; |
271 |
> |
for (i = 0; i < ninpfiles; i++) { |
272 |
|
if (!init_pabopto_inp(i, argv[a+i])) |
273 |
|
return(1); |
274 |
< |
|
274 |
> |
if (auto_grazing && fabs(inpfile[i].theta - 90.) < lim_graze) |
275 |
> |
lim_graze = fabs(inpfile[i].theta - 90.); |
276 |
> |
} |
277 |
> |
if (auto_grazing) |
278 |
> |
sprintf(auto_grazing, "%.2f", lim_graze); |
279 |
|
for (i = ninpfiles; i < totinc; i++) { /* copy for "up" symmetry */ |
280 |
|
inpfile[i] = inpfile[i-ninpfiles]; |
281 |
|
inpfile[i].phi += 180.; /* invert duplicate data */ |
329 |
|
save_bsdf_rep(stdout); /* complete header + data */ |
330 |
|
return(0); |
331 |
|
userr: |
332 |
< |
fprintf(stderr, "Usage: %s [-t][-n nproc][-s symmetry] meas1.dat meas2.dat .. > bsdf.sir\n", |
332 |
> |
fprintf(stderr, "Usage: %s [-t][-n nproc][-s symmetry][-g angle|'A'] meas1.dat meas2.dat .. > bsdf.sir\n", |
333 |
|
progname); |
334 |
|
return(1); |
335 |
|
} |