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

Comparing ray/src/cv/pabopto2xyz.c (file contents):
Revision 2.1 by greg, Wed Jul 27 18:05:53 2016 UTC vs.
Revision 2.5 by greg, Wed Jun 12 00:09:18 2019 UTC

# Line 43 | Line 43 | extern void            SDdisk2square(double sq[2], double diskx,
43  
44   char                    *progname;      /* global argv[0] */
45  
46 + int                     incident_side = 0;      /* signum(intheta) */
47 + int                     exiting_side = 0;       /* signum(outtheta) */
48 +
49   int                     nprocs = 1;     /* number of processes to run */
50  
51   int                     nchild = 0;     /* number of children (-1 in child) */
# Line 240 | Line 243 | init_pabopto_inp(PGINPUT *pgi, const char *fname)
243                  fputs(": unknown incident angle\n", stderr);
244                  return(0);
245          }
246 +        if (!incident_side) {
247 +                incident_side = (pgi->theta < 90.) ? 1 : -1;
248 +        } else if ((incident_side > 0) ^ (pgi->theta < 90.)) {
249 +                fputs(fname, stderr);
250 +                fputs(": incident on opposite side of surface\n", stderr);
251 +                return(0);
252 +        }
253                                  /* convert angle to grid position */
254          pos_from_ang(pgi->igp, pgi->theta, pgi->phi);
255          return(1);
# Line 280 | Line 290 | load_interp(PGINTERP *pgint, const int igp[2], const P
290                                  if ((pgint->ip2 == NULL) | (pgint->va == NULL))
291                                          goto memerr;
292                          }
293 +                        if (!exiting_side) {
294 +                                exiting_side = (theta < 90.) ? 1 : -1;
295 +                        } else if ((exiting_side > 0) ^ (theta < 90.)) {
296 +                                fputs(pginp->fname, stderr);
297 +                                fputs(": exiting measurement on wrong side\n", stderr);
298 +                                return(0);
299 +                        }
300                          sq_from_ang(sq, theta, phi);
301                          pgint->ip2->spt[nread][0] = sq[0]*ANGRES;
302                          pgint->ip2->spt[nread][1] = sq[1]*ANGRES;
# Line 305 | Line 322 | memerr:
322   static double
323   interp2val(const PGINTERP *pgint, double px, double py)
324   {
325 < #define NSMP    12
325 > #define NSMP    36
326          float   wt[NSMP];
327          int     si[NSMP];
328          int     n = interp2_topsamp(wt, si, NSMP, pgint->ip2, px, py);
# Line 371 | Line 388 | interp_xyz(const PGINPUT *inp0, int nf, const PGINTERP
388                                  inp0->fname, outfname);
389   #endif
390                  while (fscanf(ifp, "%lf %lf %lf", &theta, &phi, &val[0]) == 3) {
391 +                        if (!exiting_side) {
392 +                                exiting_side = (theta < 90.) ? 1 : -1;
393 +                        } else if ((exiting_side > 0) ^ (theta < 90.)) {
394 +                                fputs(inp0->fname, stderr);
395 +                                fputs(": exiting measurement on wrong side\n", stderr);
396 +                                return(0);
397 +                        }
398                          sq_from_ang(sq, theta, phi);
399                          val[1] = interp2val(int1, sq[0]*ANGRES, sq[1]*ANGRES);
400                          val[2] = interp2val(int2, sq[0]*ANGRES, sq[1]*ANGRES);
# Line 441 | Line 465 | advance_incidence(PGINPUT *slist[3], int ndx[3])
465   int
466   main(int argc, char *argv[])
467   {
444 –        char    *flist[MAX_INPUTS];
468          PGINPUT *slist[3];
469          int     i, j;
470          int     ndx[3];
# Line 467 | Line 490 | main(int argc, char *argv[])
490                  default:
491                          goto userr;
492                  }
493 <        if (i > argc-3)
493 >        if (i != argc-3)
494                  goto userr;
495          for (j = 0; j < 3; j++) {               /* prep input channels */
496 +                char    *flist[MAX_INPUTS];
497                  int     k, n;
498                  n = wordfile(flist, MAX_INPUTS, argv[i+j]);
499                  if (n <= 0) {

Diff Legend

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