ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/cv/pabopto2bsdf.c
(Generate patch)

Comparing ray/src/cv/pabopto2bsdf.c (file contents):
Revision 2.11 by greg, Thu Jan 30 18:19:43 2014 UTC vs.
Revision 2.18 by greg, Wed Mar 19 00:19:11 2014 UTC

# Line 14 | Line 14 | static const char RCSid[] = "$Id$";
14   #include <math.h>
15   #include "platform.h"
16   #include "bsdfrep.h"
17 + #include "resolu.h"
18                                  /* global argv[0] */
19   char                    *progname;
20  
# Line 24 | Line 25 | typedef struct {
25          long            dstart;         /* data start offset in file */
26   } PGINPUT;
27  
28 + double          angle_eps = 0;  /* epsilon for angle comparisons */
29 +
30   PGINPUT         *inpfile;       /* input files sorted by incidence */
31   int             ninpfiles;      /* number of input files */
32  
# Line 34 | Line 37 | cmp_inang(const void *p1, const void *p2)
37          const PGINPUT   *inp1 = (const PGINPUT *)p1;
38          const PGINPUT   *inp2 = (const PGINPUT *)p2;
39          
40 <        if (inp1->theta > inp2->theta+FTINY)
40 >        if (inp1->theta > inp2->theta+angle_eps)
41                  return(1);
42 <        if (inp1->theta < inp2->theta-FTINY)
42 >        if (inp1->theta < inp2->theta-angle_eps)
43                  return(-1);
44 <        if (inp1->phi > inp2->phi+FTINY)
44 >        if (inp1->phi > inp2->phi+angle_eps)
45                  return(1);
46 <        if (inp1->phi < inp2->phi-FTINY)
46 >        if (inp1->phi < inp2->phi-angle_eps)
47                  return(-1);
48          return(0);
49   }
# Line 121 | Line 124 | add_pabopto_inp(const int i)
124                  return(0);
125          }
126                                          /* prepare input grid */
127 +        angle_eps = 180./2./GRIDRES;
128          if (!i || cmp_inang(&inpfile[i-1], &inpfile[i])) {
129                  if (i)                  /* need to process previous incidence */
130                          make_rbfrep();
# Line 180 | Line 184 | main(int argc, char *argv[])
184          for (i = 0; i < ninpfiles; i++)
185                  if (!init_pabopto_inp(i, argv[i+1]))
186                          return(1);
187 +        angle_eps = 0;
188          qsort(inpfile, ninpfiles, sizeof(PGINPUT), &cmp_inang);
189                                                  /* compile measurements */
190          for (i = 0; i < ninpfiles; i++)
# Line 219 | Line 224 | main(int argc, char *argv[])
224   #ifdef DEBUG
225          fprintf(stderr, "Minimum BSDF = %.4f\n", bsdf_min);
226   #endif
227 +        min_log = log(bsdf_min*.5 + 1e-5);
228   #if 1                                           /* produce spheres at meas. */
229          puts("void plastic yellow\n0\n0\n5 .6 .4 .01 .04 .08\n");
224 –        min_log = log(bsdf_min*.5);
230          n = 0;
231          for (i = 0; i < GRIDRES; i++)
232              for (j = 0; j < GRIDRES; j++)
233 <                if (dsf_grid[i][j].nval > 0) {
233 >                if (dsf_grid[i][j].sum.n > 0) {
234                          ovec_from_pos(dir, i, j);
235 <                        bsdf = dsf_grid[i][j].vsum /
236 <                                (dsf_grid[i][j].nval*output_orient*dir[2]);
235 >                        bsdf = dsf_grid[i][j].sum.v /
236 >                                (dsf_grid[i][j].sum.n*output_orient*dir[2]);
237                          if (bsdf <= bsdf_min*.6)
238                                  continue;
239 <                        bsdf = log(bsdf) - min_log;
239 >                        bsdf = log(bsdf + 1e-5) - min_log;
240                          ovec_from_pos(dir, i, j);
241                          printf("yellow sphere s%04d\n0\n0\n", ++n);
242                          printf("4 %.6g %.6g %.6g %.6g\n\n",
# Line 245 | Line 250 | main(int argc, char *argv[])
250                  RBFVAL  *rbf = &dsf_list->rbfa[n];
251                  ovec_from_pos(dir, rbf->gx, rbf->gy);
252                  bsdf = eval_rbfrep(dsf_list, dir) / (output_orient*dir[2]);
253 <                bsdf = log(bsdf) - min_log;
253 >                bsdf = log(bsdf + 1e-5) - min_log;
254                  printf("red sphere p%04d\n0\n0\n", ++n);
255                  printf("4 %.6g %.6g %.6g %.6g\n\n",
256                                  dir[0]*bsdf, dir[1]*bsdf, dir[2]*bsdf,
# Line 265 | Line 270 | main(int argc, char *argv[])
270              for (j = 0; j < GRIDRES; j++) {
271                  ovec_from_pos(dir, i, j);
272                  bsdf = eval_rbfrep(dsf_list, dir) / (output_orient*dir[2]);
273 <                bsdf = log(bsdf) - min_log;
273 >                bsdf = log(bsdf + 1e-5) - min_log;
274                  fprintf(pfp, "%.8e %.8e %.8e\n",
275                                  dir[0]*bsdf, dir[1]*bsdf, dir[2]*bsdf);
276              }

Diff Legend

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