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

Comparing ray/src/rt/ambcomp.c (file contents):
Revision 1.14 by greg, Mon Aug 26 10:11:39 1991 UTC vs.
Revision 2.1 by greg, Tue Nov 12 17:10:29 1991 UTC

# Line 206 | Line 206 | FVECT  pg, dg;
206          b = 1.0/ndivs;
207          scalecolor(acol, b);
208          if (arad <= FTINY)
209 –                arad = FHUGE;
210 –        else
211 –                arad = (ndivs+ns)/arad;
212 –        if (arad > maxarad)
209                  arad = maxarad;
210 <        else if (arad < minarad)
210 >        else {
211 >                arad = (ndivs+ns)/arad;
212 >                if (arad > maxarad)
213 >                        arad = maxarad;
214 >        }
215 >        if (pg != NULL) {               /* reduce radius if gradient large */
216 >                d = DOT(pg,pg);
217 >                if (d*arad*arad > 1.0)
218 >                        arad = 1.0/sqrt(d);
219 >        }
220 >        if (arad < minarad) {
221                  arad = minarad;
222 <        arad /= sqrt(wt);
223 <        if (pg != NULL) {               /* clip pos. gradient if too large */
218 <                d = 4.0*DOT(pg,pg)*arad*arad;
219 <                if (d > 1.0) {
220 <                        d = 1.0/sqrt(d);
222 >                if (pg != NULL && d*arad*arad > 1.0) {  /* cap gradient */
223 >                        d = 1.0/arad/sqrt(d);
224                          for (i = 0; i < 3; i++)
225                                  pg[i] *= d;
226                  }
227          }
228 <        return(arad);
228 >        return(arad/sqrt(wt));
229   oopsy:
230          if (div != NULL)
231                  free((char *)div);
# Line 333 | Line 336 | AMBHEMI  *hp;
336                          if (i > 0) {
337                                  d = dp[-hp->np].r;
338                                  if (dp[0].r > d) d = dp[0].r;
339 <                                d *= 1.0 - sqrt((double)i/hp->nt);
339 >                                d *= 1.0 - (double)i/hp->nt;    /* cos(t)^2 */
340                                  mag0 += d*(b - bright(dp[-hp->np].v));
341                          }
342                          if (j > 0) {
# Line 381 | Line 384 | AMBHEMI  *hp;
384                                  error(CONSISTENCY,
385                                          "division order in dirgradient");
386   #endif
387 <                        mag += sqrt((i+.5)/hp->nt)*bright(dp->v);
387 >                        mag += sqrt((i+.5)/hp->nt)*bright(dp->v);  /* sin(t) */
388                          dp += hp->np;
389                  }
390                  phi = 2.0*PI * (j+.5)/hp->np + PI/2.0;

Diff Legend

– Removed lines
+ Added lines
< Changed lines (old)
> Changed lines (new)