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

Comparing ray/src/rt/ambcomp.c (file contents):
Revision 2.61 by greg, Sun May 18 18:59:55 2014 UTC vs.
Revision 2.68 by greg, Thu Oct 23 18:19:14 2014 UTC

# Line 21 | Line 21 | static const char      RCSid[] = "$Id$";
21   #include  "ambient.h"
22   #include  "random.h"
23  
24 < #ifdef NEWAMB
24 > #ifndef OLDAMB
25  
26   extern void             SDsquare2disk(double ds[2], double seedx, double seedy);
27  
# Line 68 | Line 68 | ambsample(                             /* initial ambient division sample */
68                  setcolor(ar.rcoef, AVGREFL, AVGREFL, AVGREFL);
69          else
70                  copycolor(ar.rcoef, hp->acoef);
71 <        if (rayorigin(&ar, AMBIENT, hp->rp, ar.rcoef) < 0) {
72 <                if (!n) memset(ap, 0, sizeof(AMBSAMP));
71 >        if (rayorigin(&ar, AMBIENT, hp->rp, ar.rcoef) < 0)
72                  return(0);
74 –        }
73          if (ambacc > FTINY) {
74                  multcolor(ar.rcoef, hp->acoef);
75                  scalecolor(ar.rcoef, 1./AVGREFL);
# Line 80 | Line 78 | ambsample(                             /* initial ambient division sample */
78          hlist[1] = j;
79          hlist[2] = i;
80          multisamp(spt, 2, urand(ilhash(hlist,3)+n));
81 <        if (!n) {                       /* avoid border samples for n==0 */
81 >                                        /* avoid coincident samples */
82 >        if (!n && (0 < i) & (i < hp->ns-1) &&
83 >                        (0 < j) & (j < hp->ns-1)) {
84                  if ((spt[0] < 0.1) | (spt[0] >= 0.9))
85                          spt[0] = 0.1 + 0.8*frandom();
86                  if ((spt[1] < 0.1) | (spt[1] >= 0.9))
# Line 99 | Line 99 | ambsample(                             /* initial ambient division sample */
99          if (ar.rt <= FTINY)
100                  return(0);              /* should never happen */
101          multcolor(ar.rcol, ar.rcoef);   /* apply coefficient */
102 <        if (!n || ar.rt*ap->d < 1.0)    /* new/closer distance? */
102 >        if (ar.rt*ap->d < 1.0)          /* new/closer distance? */
103                  ap->d = 1.0/ar.rt;
104          if (!n) {                       /* record first vertex & value */
105                  if (ar.rt > 10.0*thescene.cusize)
# Line 180 | Line 180 | ambsupersamp(AMBHEMI *hp, int cnt)
180          float   *earr = getambdiffs(hp);
181          double  e2rem = 0;
182          AMBSAMP *ap;
183 –        RAY     ar;
183          float   *ep;
184          int     i, j, n, nss;
185  
# Line 195 | Line 194 | ambsupersamp(AMBHEMI *hp, int cnt)
194                  if (e2rem <= FTINY)
195                          goto done;      /* nothing left to do */
196                  nss = *ep/e2rem*cnt + frandom();
197 <                for (n = 1; n <= nss; n++)
198 <                        cnt -= ambsample(hp, i, j, n);
197 >                for (n = 1; n <= nss && ambsample(hp,i,j,n); n++)
198 >                        --cnt;
199                  e2rem -= *ep++;         /* update remainder */
200          }
201   done:
# Line 229 | Line 228 | samp_hemi(                             /* sample indirect hemisphere */
228          hp->rp = r;
229          hp->ns = n;
230          hp->acol[RED] = hp->acol[GRN] = hp->acol[BLU] = 0.0;
231 +        memset(hp->sa, 0, sizeof(AMBSAMP)*n*n);
232          hp->sampOK = 0;
233                                          /* assign coefficient */
234          copycolor(hp->acoef, rcol);
# Line 609 | Line 609 | static uint32
609   ambcorral(AMBHEMI *hp, FVECT uv[2], const double r0, const double r1)
610   {
611          const double    max_d = 1.0/(minarad*ambacc + 0.001);
612 <        const double    ang_res = 0.5*PI/(hp->ns-1);
613 <        const double    ang_step = ang_res/((int)(16/PI*ang_res) + (1+FTINY));
612 >        const double    ang_res = 0.5*PI/hp->ns;
613 >        const double    ang_step = ang_res/((int)(16/PI*ang_res) + 1.01);
614          double          avg_d = 0;
615          uint32          flgs = 0;
616          FVECT           vec;
617 <        double          d, u, v;
617 >        double          u, v;
618          double          ang, a1;
619          int             i, j;
620                                          /* don't bother for a few samples */
# Line 636 | Line 636 | ambcorral(AMBHEMI *hp, FVECT uv[2], const double r0, c
636                  if ((ap->d <= FTINY) | (ap->d >= max_d))
637                          continue;       /* too far or too near */
638                  VSUB(vec, ap->p, hp->rp->rop);
639 <                d = DOT(vec, hp->rp->ron);
640 <                d = 1.0/sqrt(DOT(vec,vec) - d*d);
641 <                u = DOT(vec, uv[0]) * d;
642 <                v = DOT(vec, uv[1]) * d;
643 <                if ((r0*r0*u*u + r1*r1*v*v) * ap->d*ap->d <= 1.0)
639 >                u = DOT(vec, uv[0]);
640 >                v = DOT(vec, uv[1]);
641 >                if ((r0*r0*u*u + r1*r1*v*v) * ap->d*ap->d <= u*u + v*v)
642                          continue;       /* occluder outside ellipse */
643                  ang = atan2a(v, u);     /* else set direction flags */
644 <                for (a1 = ang-.5*ang_res; a1 <= ang+.5*ang_res; a1 += ang_step)
644 >                for (a1 = ang-ang_res; a1 <= ang+ang_res; a1 += ang_step)
645                          flgs |= 1L<<(int)(16/PI*(a1 + 2.*PI*(a1 < 0)));
646              }
647                                          /* add low-angle incident (< 20deg) */
# Line 696 | Line 694 | doambient(                             /* compute ambient component */
694                  return(0);
695  
696          if ((ra == NULL) & (pg == NULL) & (dg == NULL) ||
697 <                        (hp->sampOK < 0) | (hp->ns < 4)) {
697 >                        (hp->sampOK < 0) | (hp->ns < 6)) {
698                  free(hp);               /* Hessian not requested/possible */
699                  return(-1);             /* value-only return value */
700          }

Diff Legend

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