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.39 by greg, Thu Feb 25 02:13:15 2021 UTC vs.
Revision 2.40 by greg, Mon Apr 5 19:56:18 2021 UTC

# Line 32 | Line 32 | PGINPUT                *inpfile;       /* input files sorted by incidence
32  
33   int             rev_orient = 0; /* shall we reverse surface orientation? */
34  
35 + double          lim_graze = 0;  /* limit scattering near grazing (deg above) */
36 +
37   /* Compare incident angles */
38   static int
39   cmp_indir(const void *p1, const void *p2)
# Line 178 | Line 180 | add_pabopto_inp(const int i)
180                                  fclose(fp);
181                                  return(1);
182                          }
183 <                if (rev_orient) {       /* reverse Z-axis to face outside */
183 >                                        /* check if scatter angle is too low */
184 >                if (fabs(theta_out - 90.) < lim_graze-FTINY)
185 >                        continue;
186 >                if (rev_orient) {       /* reverse Z-axis to face outside? */
187                          theta_out = 180. - theta_out;
188                          phi_out = 360. - phi_out;
189                  }
# Line 224 | Line 229 | int
229   main(int argc, char *argv[])
230   {
231          extern int      nprocs;
232 +        int             auto_grazing = 0;
233          const char      *symmetry = "0";
234          int             ninpfiles, totinc;
235          int             a, i;
# Line 240 | Line 246 | main(int argc, char *argv[])
246                  case 's':
247                          symmetry = argv[++a];
248                          break;
249 +                case 'g':
250 +                        if (toupper(argv[a+1][0]) == 'A')
251 +                                auto_grazing = 1;
252 +                        else
253 +                                lim_graze = atof(argv[a+1]);
254 +                        ++a;
255 +                        break;
256                  default:
257                          goto userr;
258                  }
# Line 251 | Line 264 | main(int argc, char *argv[])
264          inpfile = (PGINPUT *)malloc(sizeof(PGINPUT)*totinc);
265          if (inpfile == NULL)
266                  return(1);
267 <        for (i = 0; i < ninpfiles; i++)
267 >        if (auto_grazing)
268 >                lim_graze = 90.;
269 >        for (i = 0; i < ninpfiles; i++) {
270                  if (!init_pabopto_inp(i, argv[a+i]))
271                          return(1);
272 <
272 >                if (auto_grazing && fabs(inpfile[i].theta - 90.) < lim_graze)
273 >                        lim_graze = fabs(inpfile[i].theta - 90.);
274 >        }
275          for (i = ninpfiles; i < totinc; i++) {  /* copy for "up" symmetry */
276                  inpfile[i] = inpfile[i-ninpfiles];
277                  inpfile[i].phi += 180.;         /* invert duplicate data */
# Line 308 | Line 325 | main(int argc, char *argv[])
325          save_bsdf_rep(stdout);                  /* complete header + data */
326          return(0);
327   userr:
328 <        fprintf(stderr, "Usage: %s [-t][-n nproc][-s symmetry] meas1.dat meas2.dat .. > bsdf.sir\n",
328 >        fprintf(stderr, "Usage: %s [-t][-n nproc][-s symmetry][-g angle|'A'] meas1.dat meas2.dat .. > bsdf.sir\n",
329                                          progname);
330          return(1);
331   }

Diff Legend

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