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

Comparing ray/src/util/rsensor.c (file contents):
Revision 2.10 by greg, Thu May 19 16:09:38 2011 UTC vs.
Revision 2.16 by greg, Thu Apr 23 23:54:08 2015 UTC

# Line 15 | Line 15 | static const char RCSid[] = "$Id$";
15  
16   #define DEGREE          (PI/180.)
17  
18 < #define MAXNT           180     /* maximum number of theta divisions */
18 > #define MAXNT           181     /* maximum number of theta divisions */
19   #define MAXNP           360     /* maximum number of phi divisions */
20  
21   extern char     *progname;      /* global argv[0] */
# Line 209 | Line 209 | load_sensor(
209   {
210          int     warnedneg;
211          char    linebuf[8192];
212 +        int     last_pos_val = 0;
213          int     nelem = 1000;
214          float   *sarr = (float *)malloc(sizeof(float)*nelem);
215          FILE    *fp;
# Line 269 | Line 270 | load_sensor(
270                                          error(WARNING, errmsg);
271                                  }
272                                  sarr[i] = .0;
273 <                        }
273 >                        } else if (i > ntp[0]*(ntp[1]+1) && sarr[i] > FTINY)
274 >                                last_pos_val = i;
275                          ++i;
276                  }
277                  if (i == ntp[0]*(ntp[1]+1))
# Line 288 | Line 290 | load_sensor(
290                          error(USER, errmsg);
291                  }
292          }
293 <        nelem = i;
293 >                                                /* truncate zero region */
294 >        ntp[0] = (last_pos_val + ntp[1])/(ntp[1]+1) - 1;
295 >        nelem = (ntp[0]+1)*(ntp[1]+1);
296          fclose(fp);
297          errmsg[0] = '\0';                       /* sanity checks */
298 +        if (!last_pos_val)
299 +                sprintf(errmsg, "no positive sensor values in file '%s'", sfile);
300          if (ntp[0] <= 0)
301                  sprintf(errmsg, "no data in sensor file '%s'", sfile);
302          else if (fabs(sarr[ntp[1]+1]) > FTINY)
# Line 305 | Line 311 | load_sensor(
311                                  sfile);
312          else if (sarr[ntp[0]*(ntp[1]+1)] <= FTINY)
313                  sprintf(errmsg,
314 <                        "maximum theta must be positive in sensor file '%s'",
315 <                                sfile);
314 >                        "maximum theta (%f) must be positive in sensor file '%s'",
315 >                                sarr[ntp[0]*(ntp[1]+1)], sfile);
316          if (errmsg[0])
317                  error(USER, errmsg);
318          return((float *)realloc((void *)sarr, sizeof(float)*nelem));
# Line 398 | Line 404 | init_ptable(
404                  tvals[i] = 1. - ( (1.-frac)*cos(thdiv[t]) +
405                                                  frac*cos(thdiv[t+1]) );
406                                  /* offset b/c sensor values are centered */
407 <                if (!t || (t < sntp[0]-1) & (frac >= 0.5))
407 >                if ((t < sntp[0]-1) & (!t | (frac >= 0.5))) {
408                          frac -= 0.5;
409 <                else {
409 >                } else {
410                          frac += 0.5;
411                          --t;
412                  }
# Line 424 | Line 430 | init_ptable(
430                  }
431                  pvals[i*(nphi+1) + nphi] = phdiv[sntp[1]];
432          }
433 +                                                /* duplicate final row */
434 +        memcpy(pvals+ntheta*(nphi+1), pvals+(ntheta-1)*(nphi+1),
435 +                                sizeof(*pvals)*(nphi+1));
436          tvals[0] = .0f;
437          tvals[ntheta] = (float)tsize;
438   }
# Line 470 | Line 479 | sens_val(
479          int     t, p;
480          
481          dv[2] = DOT(dvec, ourview.vdir);
482 <        theta = (float)((1./DEGREE) * acos(dv[2]));
482 >        theta = acos(dv[2]);
483          if (theta >= maxtheta)
484                  return(.0f);
485          dv[0] = DOT(dvec, ourview.hvec);
486          dv[1] = DOT(dvec, ourview.vvec);
487 <        phi = (float)((1./DEGREE) * atan2(-dv[0], dv[1]));
488 <        while (phi < .0f) phi += 360.f;
487 >        phi = atan2(-dv[0], dv[1]);
488 >        while (phi < .0f) phi += (float)(2.*PI);
489          t = (int)(theta/maxtheta * sntp[0]);
490 <        p = (int)(phi*(1./360.) * sntp[1]);
490 >        p = (int)(phi*(1./(2.*PI)) * sntp[1]);
491                          /* hack for non-uniform sensor grid */
492 +        theta *= (float)(1./DEGREE);
493 +        phi *= (float)(1./DEGREE);
494          while (t+1 < sntp[0] && theta >= s_theta(t+1))
495                  ++t;
496          while (t-1 >= 0 && theta <= s_theta(t-1))
# Line 542 | Line 553 | comp_sensor(
553                                  continue;
554                          }
555                          rr.rmax = .0;
556 <                        rayorigin(&rr, PRIMARY, NULL, NULL);
556 >                        rayorigin(&rr, PRIMARY|SPECULAR, NULL, NULL);
557                          scalecolor(rr.rcoef, sf);
558                          if (ray_pqueue(&rr) == 1)
559                                  addcolor(vsum, rr.rcol);
# Line 558 | Line 569 | comp_sensor(
569                          continue;
570                  }
571                  rr.rmax = .0;
572 <                rayorigin(&rr, PRIMARY, NULL, NULL);
572 >                rayorigin(&rr, PRIMARY|SPECULAR, NULL, NULL);
573                  scalecolor(rr.rcoef, sf);
574                  if (ray_pqueue(&rr) == 1)
575                          addcolor(vsum, rr.rcol);

Diff Legend

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