/* Read in BSDF and interpolate as Klems matrix representation */
int
main(int argc, char *argv[])
{
int dofwd = 0, dobwd = 1;
char buf[2048];
char *cp;
int i, na;
progname = argv[0];
esupport |= E_VARIABLE|E_FUNCTION|E_RCONST;
esupport &= ~(E_INCHAN|E_OUTCHAN);
scompile("PI:3.14159265358979323846", NULL, 0);
biggerlib();
for (i = 1; i < argc && (argv[i][0] == '-') | (argv[i][0] == '+'); i++)
switch (argv[i][1]) { /* get options */
case 'n':
npsamps = atoi(argv[++i]);
if (npsamps <= 0)
goto userr;
break;
case 'e':
scompile(argv[++i], NULL, 0);
single_plane_incident = 0;
break;
case 'f':
if (!argv[i][2]) {
if (strchr(argv[++i], '=') != NULL) {
add_wbsdf("-f", 1);
add_wbsdf(argv[i], 1);
} else {
char *fpath = getpath(argv[i],
getrlibpath(), 0);
if (fpath == NULL) {
fprintf(stderr,
"%s: cannot find file '%s'\n",
argv[0], argv[i]);
return(1);
}
fcompile(fpath);
single_plane_incident = 0;
}
} else
dofwd = (argv[i][0] == '+');
break;
case 'b':
dobwd = (argv[i][0] == '+');
break;
case 'h':
kbasis = klems_half;
add_wbsdf("-a", 1);
add_wbsdf("kh", 1);
break;
case 'q':
kbasis = klems_quarter;
add_wbsdf("-a", 1);
add_wbsdf("kq", 1);
break;
case 'l':
lobe_lim = atoi(argv[++i]);
break;
case 'p':
do_prog = atoi(argv[i]+2);
break;
case 'C':
add_wbsdf(argv[i], 1);
add_wbsdf(argv[++i], 1);
break;
default:
goto userr;
}
if (kbasis == klems_full) { /* default (full) basis? */
add_wbsdf("-a", 1);
add_wbsdf("kf", 1);
}
strcpy(buf, "File produced by: ");
if (convert_commandline(buf+18, sizeof(buf)-18, argv) != NULL) {
add_wbsdf("-C", 1); add_wbsdf(buf, 0);
}
if (single_plane_incident >= 0) { /* function-based BSDF? */
if (i != argc-1 || fundefined(argv[i]) < 3) {
fprintf(stderr,
"%s: need single function with 6 arguments: bsdf(ix,iy,iz,ox,oy,oz)\n",
progname);
fprintf(stderr, "\tor 3 arguments using Dx,Dy,Dz: bsdf(ix,iy,iz)\n");
goto userr;
}
++eclock;
if (dofwd) {
input_orient = -1;
output_orient = -1;
prog_start("Evaluating outside reflectance");
eval_function(argv[i]);
output_orient = 1;
prog_start("Evaluating outside->inside transmission");
eval_function(argv[i]);
}
if (dobwd) {
input_orient = 1;
output_orient = 1;
prog_start("Evaluating inside reflectance");
eval_function(argv[i]);
output_orient = -1;
prog_start("Evaluating inside->outside transmission");
eval_function(argv[i]);
}
return(wrap_up());
}
/* XML input? */
if (i == argc-1 && (cp = argv[i]+strlen(argv[i])-4) > argv[i] &&
!strcasecmp(cp, ".xml")) {
eval_bsdf(argv[i]); /* load & resample BSDF */
return(wrap_up());
}
if (i < argc) { /* open input files if given */
int nbsdf = 0;
for ( ; i < argc; i++) { /* interpolate each component */
char pbuf[256];
FILE *fpin = fopen(argv[i], "rb");
if (fpin == NULL) {
fprintf(stderr, "%s: cannot open BSDF interpolant '%s'\n",
progname, argv[i]);
return(1);
}
if (!load_bsdf_rep(fpin))
return(1);
fclose(fpin);
sprintf(pbuf, "Interpolating component '%s'", argv[i]);
prog_start(pbuf);
eval_rbf();
}
return(wrap_up());
}
SET_FILE_BINARY(stdin); /* load from stdin */
if (!load_bsdf_rep(stdin))
return(1);
prog_start("Interpolating from standard input");
eval_rbf(); /* resample dist. */
return(wrap_up());
userr:
fprintf(stderr,
"Usage: %s [-n spp][-h|-q][-l maxlobes] [bsdf.sir ..] > bsdf.xml\n", progname);
fprintf(stderr,
" or: %s [-n spp][-h|-q] bsdf_in.xml > bsdf_out.xml\n", progname);
fprintf(stderr,
" or: %s [-n spp][-h|-q][{+|-}for[ward]][{+|-}b[ackward]][-e expr][-f file] bsdf_func > bsdf.xml\n",
progname);
return(1);
}
/* Read in BSDF and interpolate as Klems matrix representation */
int
main(int argc, char *argv[])
{
int dofwd = 0, dobwd = 1;
char *cp;
int i, na;
progname = argv[0];
esupport |= E_VARIABLE|E_FUNCTION|E_RCONST;
esupport &= ~(E_INCHAN|E_OUTCHAN);
scompile("PI:3.14159265358979323846", NULL, 0);
biggerlib();
for (i = 1; i < argc && (argv[i][0] == '-') | (argv[i][0] == '+'); i++)
switch (argv[i][1]) { /* get options */
case 'n':
npsamps = atoi(argv[++i]);
if (npsamps <= 0)
goto userr;
break;
case 'e':
scompile(argv[++i], NULL, 0);
single_plane_incident = 0;
break;
case 'f':
if (!argv[i][2]) {
fcompile(argv[++i]);
single_plane_incident = 0;
} else
dofwd = (argv[i][0] == '+');
break;
case 'b':
dobwd = (argv[i][0] == '+');
break;
case 'h':
kbasis = "LBNL/Klems Half";
break;
case 'q':
kbasis = "LBNL/Klems Quarter";
break;
default:
goto userr;
}
if (single_plane_incident >= 0) { /* function-based BSDF? */
if (i != argc-1 || fundefined(argv[i]) != 6) {
fprintf(stderr,
"%s: need single function with 6 arguments: bsdf(ix,iy,iz,ox,oy,oz)\n",
progname);
fprintf(stderr, "\tor 3 arguments using Dx,Dy,Dz: bsdf(ix,iy,iz)\n",
progname);
goto userr;
}
++eclock;
xml_header(argc, argv); /* start XML output */
xml_prologue(NULL);
if (dofwd) {
input_orient = -1;
output_orient = -1;
eval_function(argv[i]); /* outside reflectance */
output_orient = 1;
eval_function(argv[i]); /* outside -> inside */
}
if (dobwd) {
input_orient = 1;
output_orient = 1;
eval_function(argv[i]); /* inside reflectance */
output_orient = -1;
eval_function(argv[i]); /* inside -> outside */
}
xml_epilogue(); /* finish XML output & exit */
return(0);
}
/* XML input? */
if (i == argc-1 && (cp = argv[i]+strlen(argv[i])-4) > argv[i] &&
!strcasecmp(cp, ".xml")) {
xml_header(argc, argv); /* start XML output */
eval_bsdf(argv[i]); /* load & resample BSDF */
xml_epilogue(); /* finish XML output & exit */
return(0);
}
if (i < argc) { /* open input files if given */
int nbsdf = 0;
for ( ; i < argc; i++) { /* interpolate each component */
FILE *fpin = fopen(argv[i], "rb");
if (fpin == NULL) {
fprintf(stderr, "%s: cannot open BSDF interpolant '%s'\n",
progname, argv[i]);
return(1);
}
if (!load_bsdf_rep(fpin))
return(1);
fclose(fpin);
if (!nbsdf++) { /* start XML on first dist. */
xml_header(argc, argv);
xml_prologue(NULL);
}
eval_rbf();
}
xml_epilogue(); /* finish XML output & exit */
return(0);
}
SET_FILE_BINARY(stdin); /* load from stdin */
if (!load_bsdf_rep(stdin))
return(1);
xml_header(argc, argv); /* start XML output */
xml_prologue(NULL);
eval_rbf(); /* resample dist. */
xml_epilogue(); /* finish XML output & exit */
return(0);
userr:
fprintf(stderr,
"Usage: %s [-n spp][-h|-q][bsdf.sir ..] > bsdf.xml\n", progname);
fprintf(stderr,
" or: %s [-n spp][-h|-q] bsdf_in.xml > bsdf_out.xml\n", progname);
fprintf(stderr,
" or: %s [-n spp][-h|-q][{+|-}for[ward]][{+|-}b[ackward]][-e expr][-f file] bsdf_func > bsdf.xml\n",
progname);
return(1);
}