--- ray/src/cv/bsdf2ttree.c 2012/11/10 19:47:42 2.7 +++ ray/src/cv/bsdf2ttree.c 2013/03/23 02:21:14 2.13 @@ -1,5 +1,5 @@ #ifndef lint -static const char RCSid[] = "$Id: bsdf2ttree.c,v 2.7 2012/11/10 19:47:42 greg Exp $"; +static const char RCSid[] = "$Id: bsdf2ttree.c,v 2.13 2013/03/23 02:21:14 greg Exp $"; #endif /* * Load measured BSDF interpolant and write out as XML file with tensor tree. @@ -12,6 +12,7 @@ static const char RCSid[] = "$Id: bsdf2ttree.c,v 2.7 2 #include #include #include "platform.h" +#include "calcomp.h" #include "bsdfrep.h" /* global argv[0] */ char *progname; @@ -20,22 +21,92 @@ double pctcull = 90.; /* sampling order */ int samp_order = 6; +/* Output XML prologue to stdout */ +static void +xml_prologue(int ac, char *av[]) +{ + puts(""); + puts(""); + fputs(""); + puts("System"); + puts("BSDF"); + puts(""); + puts(""); + puts("\t"); + puts("\t\tName"); + puts("\t\tManufacturer"); + puts("\t\tOther"); + puts("\t"); + puts("\t"); + printf("\t\tTensorTree%c\n", + single_plane_incident ? '3' : '4'); + puts("\t"); +} + +/* Output XML data prologue to stdout */ +static void +data_prologue() +{ + static const char *bsdf_type[4] = { + "Reflection Front", + "Transmission Front", + "Transmission Back", + "Reflection Back" + }; + + puts("\t"); + puts("\t\tSystem"); + puts("\t\tVisible"); + puts("\t\tCIE Illuminant D65 1nm.ssp"); + puts("\t\tASTM E308 1931 Y.dsp"); + puts("\t\t"); + printf("\t\t\t%s\n", + bsdf_type[(input_orient>0)<<1 | (output_orient>0)]); + puts("\t\t\tLBNL/Shirley-Chiu"); + puts("\t\t\tBTDF"); + puts("\t\t\t"); +} + +/* Output XML data epilogue to stdout */ +static void +data_epilogue(void) +{ + puts("\t\t\t"); + puts("\t\t"); + puts("\t"); +} + +/* Output XML epilogue to stdout */ +static void +xml_epilogue(void) +{ + puts(""); + puts(""); + puts(""); +} + /* Interpolate and output isotropic BSDF data */ static void -interp_isotropic() +eval_isotropic(char *funame) { const int sqres = 1<= 0) fprintf(stderr, "data with %f%% culling\n", pctcull); + if (pctcull >= 0) fprintf(stderr, "data with %.1f%% culling\n", pctcull); else fputs("raw data\n", stderr); #endif - if (pctcull >= 0) { /* begin output */ + data_prologue(); /* begin output */ + if (pctcull >= 0) { sprintf(cmd, "rttree_reduce -h -a -ff -r 3 -t %f -g %d", pctcull, samp_order); fflush(stdout); @@ -50,17 +121,22 @@ interp_isotropic() fputs("{\n", stdout); /* run through directions */ for (ix = 0; ix < sqres/2; ix++) { - RBFNODE *rbf; - SDsquare2disk(ivec, (ix+.5)/sqres, .5); - ivec[2] = input_orient * - sqrt(1. - ivec[0]*ivec[0] - ivec[1]*ivec[1]); - rbf = advect_rbf(ivec); + RBFNODE *rbf = NULL; + iovec[0] = (ix+.5)/sqres - 1.; + iovec[1] = .0; + iovec[2] = input_orient * sqrt(1. - iovec[0]*iovec[0]); + if (funame == NULL) + rbf = advect_rbf(iovec); for (ox = 0; ox < sqres; ox++) for (oy = 0; oy < sqres; oy++) { - SDsquare2disk(ovec, (ox+.5)/sqres, (oy+.5)/sqres); - ovec[2] = output_orient * - sqrt(1. - ovec[0]*ovec[0] - ovec[1]*ovec[1]); - bsdf = eval_rbfrep(rbf, ovec) * output_orient/ovec[2]; + SDsquare2disk(iovec+3, (ox+.5)/sqres, (oy+.5)/sqres); + iovec[5] = output_orient * + sqrt(1. - iovec[3]*iovec[3] - iovec[4]*iovec[4]); + if (funame == NULL) + bsdf = eval_rbfrep(rbf, iovec+3) * + output_orient/iovec[5]; + else + bsdf = funvalue(funame, 6, iovec); if (pctcull >= 0) fwrite(&bsdf, sizeof(bsdf), 1, ofp); else @@ -80,24 +156,26 @@ interp_isotropic() fputs("\t0\n", stdout); fputs("}\n", stdout); } + data_epilogue(); } /* Interpolate and output anisotropic BSDF data */ static void -interp_anisotropic() +eval_anisotropic(char *funame) { const int sqres = 1<= 0) fprintf(stderr, "data with %f%% culling\n", pctcull); + if (pctcull >= 0) fprintf(stderr, "data with %.1f%% culling\n", pctcull); else fputs("raw data\n", stderr); #endif - if (pctcull >= 0) { /* begin output */ + data_prologue(); /* begin output */ + if (pctcull >= 0) { sprintf(cmd, "rttree_reduce -h -a -ff -r 4 -t %f -g %d", pctcull, samp_order); fflush(stdout); @@ -112,17 +190,23 @@ interp_anisotropic() /* run through directions */ for (ix = 0; ix < sqres; ix++) for (iy = 0; iy < sqres; iy++) { - RBFNODE *rbf; - SDsquare2disk(ivec, (ix+.5)/sqres, (iy+.5)/sqres); - ivec[2] = input_orient * - sqrt(1. - ivec[0]*ivec[0] - ivec[1]*ivec[1]); - rbf = advect_rbf(ivec); + RBFNODE *rbf = NULL; /* Klems reversal */ + SDsquare2disk(iovec, (ix+.5)/sqres, (iy+.5)/sqres); + iovec[0] = -iovec[0]; iovec[1] = -iovec[1]; + iovec[2] = input_orient * + sqrt(1. - iovec[0]*iovec[0] - iovec[1]*iovec[1]); + if (funame == NULL) + rbf = advect_rbf(iovec); for (ox = 0; ox < sqres; ox++) for (oy = 0; oy < sqres; oy++) { - SDsquare2disk(ovec, (ox+.5)/sqres, (oy+.5)/sqres); - ovec[2] = output_orient * - sqrt(1. - ovec[0]*ovec[0] - ovec[1]*ovec[1]); - bsdf = eval_rbfrep(rbf, ovec) * output_orient/ovec[2]; + SDsquare2disk(iovec+3, (ox+.5)/sqres, (oy+.5)/sqres); + iovec[5] = output_orient * + sqrt(1. - iovec[3]*iovec[3] - iovec[4]*iovec[4]); + if (funame == NULL) + bsdf = eval_rbfrep(rbf, iovec+3) * + output_orient/iovec[5]; + else + bsdf = funvalue(funame, 6, iovec); if (pctcull >= 0) fwrite(&bsdf, sizeof(bsdf), 1, ofp); else @@ -139,77 +223,49 @@ interp_anisotropic() } } else fputs("}\n", stdout); + data_epilogue(); } -/* Output XML prologue to stdout */ -static void -xml_prologue(int ac, char *av[]) -{ - static const char *bsdf_type[4] = { - "Reflection Front", - "Transmission Front", - "Transmission Back", - "Reflection Back" - }; - - puts(""); - puts(""); - fputs(""); - puts("System"); - puts("BSDF"); - puts(""); - puts(""); - puts("\t"); - puts("\t\tName"); - puts("\t\tManufacturer"); - puts("\t\tOther"); - puts("\t"); - puts("\t"); - printf("\t\tTensorTree%c\n", - single_plane_incident ? '3' : '4'); - puts("\t"); - puts("\t"); - puts("\t\tSystem"); - puts("\t\tVisible"); - puts("\t\tCIE Illuminant D65 1nm.ssp"); - puts("\t\tASTM E308 1931 Y.dsp"); - puts("\t\t"); - printf("\t\t\t%s\n", - bsdf_type[(input_orient>0)<<1 | (output_orient>0)]); - puts("\t\t\tLBNL/Shirley-Chiu"); - puts("\t\t\tBTDF"); - puts("\t\t\t"); -} - -/* Output XML epilogue to stdout */ -static void -xml_epilogue(void) -{ - puts("\t\t\t"); - puts("\t\t"); - puts("\t"); - puts(""); - puts(""); - puts(""); -} - /* Read in BSDF and interpolate as tensor tree representation */ int main(int argc, char *argv[]) { - FILE *fpin = stdin; - int i; + int dofwd = 0, dobwd = 1; + int i, na; progname = argv[0]; - for (i = 1; i < argc-1 && argv[i][0] == '-'; i++) - switch (argv[i][1]) { /* get option */ + esupport |= E_VARIABLE|E_FUNCTION|E_RCONST; + esupport &= ~(E_INCHAN|E_OUTCHAN); + scompile("PI:3.14159265358979323846", NULL, 0); + biggerlib(); + for (i = 1; i < argc-1 && (argv[i][0] == '-') | (argv[i][0] == '+'); i++) + switch (argv[i][1]) { /* get options */ + case 'e': + scompile(argv[++i], NULL, 0); + break; + case 'f': + if (!argv[i][2]) + fcompile(argv[++i]); + else + dofwd = (argv[i][0] == '+'); + break; + case 'b': + dobwd = (argv[i][0] == '+'); + break; case 't': - pctcull = atof(argv[++i]); + switch (argv[i][2]) { + case '3': + single_plane_incident = 1; + break; + case '4': + single_plane_incident = 0; + break; + case '\0': + pctcull = atof(argv[++i]); + break; + default: + goto userr; + } break; case 'g': samp_order = atoi(argv[++i]); @@ -217,30 +273,71 @@ main(int argc, char *argv[]) default: goto userr; } - - if (i == argc-1) { /* open input if given */ - fpin = fopen(argv[i], "r"); - if (fpin == NULL) { - fprintf(stderr, "%s: cannot open BSDF interpolant '%s'\n", - progname, argv[1]); - return(1); + if (single_plane_incident >= 0) { /* function-based BSDF? */ + void (*evf)(char *s) = single_plane_incident ? + &eval_isotropic : &eval_anisotropic; + 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); + goto userr; } - } else if (i < argc-1) - goto userr; - SET_FILE_BINARY(fpin); /* load BSDF interpolant */ - if (!load_bsdf_rep(fpin)) + xml_prologue(argc, argv); /* start XML output */ + if (dofwd) { + input_orient = -1; + output_orient = -1; + (*evf)(argv[i]); /* outside reflectance */ + output_orient = 1; + (*evf)(argv[i]); /* outside -> inside */ + } + if (dobwd) { + input_orient = 1; + output_orient = 1; + (*evf)(argv[i]); /* inside reflectance */ + output_orient = -1; + (*evf)(argv[i]); /* inside -> outside */ + } + 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_prologue(argc, argv); + if (single_plane_incident) + eval_isotropic(NULL); + else + eval_anisotropic(NULL); + } + xml_epilogue(); /* finish XML output & exit */ + return(0); + } + SET_FILE_BINARY(stdin); /* load from stdin */ + if (!load_bsdf_rep(stdin)) return(1); - fclose(fpin); xml_prologue(argc, argv); /* start XML output */ if (single_plane_incident) /* resample dist. */ - interp_isotropic(); + eval_isotropic(NULL); else - interp_anisotropic(); - xml_epilogue(); /* finish XML output */ + eval_anisotropic(NULL); + xml_epilogue(); /* finish XML output & exit */ return(0); userr: fprintf(stderr, - "Usage: %s [-t pctcull][-g log2grid] [bsdf.sir] > bsdf.xml\n", + "Usage: %s [-g Nlog2][-t pctcull] [bsdf.sir ..] > bsdf.xml\n", + progname); + fprintf(stderr, + " or: %s -t{3|4} [-g Nlog2][-t pctcull][{+|-}for[ward]][{+|-}b[ackward]][-e expr][-f file] bsdf_func > bsdf.xml\n", progname); return(1); }