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.8 by greg, Tue Oct 22 04:29:27 2013 UTC vs.
Revision 2.14 by greg, Thu Feb 20 09:21:32 2014 UTC

# Line 70 | Line 70 | init_pabopto_inp(const int i, const char *fname)
70                          fclose(fp);
71                          return(0);
72                  }
73 +                if (sscanf(buf, "sample_name \"%[^\"]\"", bsdf_name) == 1)
74 +                        continue;
75                  if (sscanf(buf, "format: theta phi %s", typ) == 1) {
76                          if (!strcasecmp(typ, "DSF")) {
77                                  inpfile[i].isDSF = 1;
# Line 214 | Line 216 | main(int argc, char *argv[])
216                  return(1);
217                                                  /* reduce data set */
218          make_rbfrep();
219 <                                                /* produce spheres at meas. */
219 > #ifdef DEBUG
220 >        fprintf(stderr, "Minimum BSDF = %.4f\n", bsdf_min);
221 > #endif
222 >        min_log = log(bsdf_min*.5 + 1e-5);
223 > #if 1                                           /* produce spheres at meas. */
224          puts("void plastic yellow\n0\n0\n5 .6 .4 .01 .04 .08\n");
219 –        min_log = log(bsdf_min*.5);
225          n = 0;
226          for (i = 0; i < GRIDRES; i++)
227              for (j = 0; j < GRIDRES; j++)
# Line 226 | Line 231 | main(int argc, char *argv[])
231                                  (dsf_grid[i][j].nval*output_orient*dir[2]);
232                          if (bsdf <= bsdf_min*.6)
233                                  continue;
234 <                        bsdf = log(bsdf) - min_log;
234 >                        bsdf = log(bsdf + 1e-5) - min_log;
235                          ovec_from_pos(dir, i, j);
236                          printf("yellow sphere s%04d\n0\n0\n", ++n);
237                          printf("4 %.6g %.6g %.6g %.6g\n\n",
238                                          dir[0]*bsdf, dir[1]*bsdf, dir[2]*bsdf,
239                                          .007*bsdf);
240                  }
241 <                                                /* output continuous surface */
241 > #endif
242 > #if 1                                           /* spheres at RBF peaks */
243 >        puts("void plastic red\n0\n0\n5 .8 .01 .01 .04 .08\n");
244 >        for (n = 0; n < dsf_list->nrbf; n++) {
245 >                RBFVAL  *rbf = &dsf_list->rbfa[n];
246 >                ovec_from_pos(dir, rbf->gx, rbf->gy);
247 >                bsdf = eval_rbfrep(dsf_list, dir) / (output_orient*dir[2]);
248 >                bsdf = log(bsdf + 1e-5) - min_log;
249 >                printf("red sphere p%04d\n0\n0\n", ++n);
250 >                printf("4 %.6g %.6g %.6g %.6g\n\n",
251 >                                dir[0]*bsdf, dir[1]*bsdf, dir[2]*bsdf,
252 >                                .011*bsdf);
253 >        }
254 > #endif
255 > #if 1                                           /* output continuous surface */
256          puts("void trans tgreen\n0\n0\n7 .7 1 .7 .04 .04 .9 1\n");
257          fflush(stdout);
258          sprintf(buf, "gensurf tgreen bsdf - - - %d %d", GRIDRES-1, GRIDRES-1);
# Line 246 | Line 265 | main(int argc, char *argv[])
265              for (j = 0; j < GRIDRES; j++) {
266                  ovec_from_pos(dir, i, j);
267                  bsdf = eval_rbfrep(dsf_list, dir) / (output_orient*dir[2]);
268 <                bsdf = log(bsdf) - min_log;
268 >                bsdf = log(bsdf + 1e-5) - min_log;
269                  fprintf(pfp, "%.8e %.8e %.8e\n",
270                                  dir[0]*bsdf, dir[1]*bsdf, dir[2]*bsdf);
271              }
272 <        return(pclose(pfp)==0 ? 0 : 1);
272 >        if (pclose(pfp) != 0)
273 >                return(1);
274 > #endif
275 >        return(0);
276   }
277   #endif

Diff Legend

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