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 1.5 by greg, Thu Jun 13 09:41:17 1991 UTC vs.
Revision 1.8 by greg, Thu Jun 13 13:58:24 1991 UTC

# Line 86 | Line 86 | RAY  *r;
86          rayvalue(&ar);
87          ndims--;
88          addcolor(dp->v, ar.rcol);
89 <        if (ar.rt < FHUGE)
89 >        if (ar.rt > FTINY && ar.rt < FHUGE)
90                  dp->r += 1.0/ar.rt;
91                                          /* (re)initialize error */
92          if (dp->n++) {
# Line 180 | Line 180 | FVECT  pg, dg;
180                          addcolor(acol, dp->v);
181                  }
182                  b = bright(acol);
183 <                if (b > FTINY)
183 >                if (b > FTINY) {
184                          b = ndivs/b;
185 <                else
186 <                        b = 0.0;
187 <                if (pg != NULL) {
188 <                        posgradient(pg, div, &hemi);
189 <                        for (i = 0; i < 3; i++)
190 <                                pg[i] *= b;
185 >                        if (pg != NULL) {
186 >                                posgradient(pg, div, &hemi);
187 >                                for (i = 0; i < 3; i++)
188 >                                        pg[i] *= b;
189 >                        }
190 >                        if (dg != NULL) {
191 >                                dirgradient(dg, div, &hemi);
192 >                                for (i = 0; i < 3; i++)
193 >                                        dg[i] *= b;
194 >                        }
195 >                } else {
196 >                        if (pg != NULL)
197 >                                for (i = 0; i < 3; i++)
198 >                                        pg[i] = 0.0;
199 >                        if (dg != NULL)
200 >                                for (i = 0; i < 3; i++)
201 >                                        dg[i] = 0.0;
202                  }
192 –                if (dg != NULL) {
193 –                        dirgradient(dg, div, &hemi);
194 –                        for (i = 0; i < 3; i++)
195 –                                dg[i] *= b;
196 –                }
203                  free((char *)div);
204          }
205          b = 1.0/ndivs;
# Line 207 | Line 213 | FVECT  pg, dg;
213          else if (arad < minarad)
214                  arad = minarad;
215          arad /= sqrt(r->rweight);
216 +        if (pg != NULL) {               /* clip pos. gradient if too large */
217 +                d = 4.0*DOT(pg,pg)*arad*arad;
218 +                if (d > 1.0) {
219 +                        d = 1.0/sqrt(d);
220 +                        for (i = 0; i < 3; i++)
221 +                                pg[i] *= d;
222 +                }
223 +        }
224          return(arad);
225   oopsy:
226          if (div != NULL)
# Line 249 | Line 263 | register AMBHEMI  *hp;
263          dp = da;
264          for (i = 0; i < hp->nt; i++)
265                  for (j = 0; j < hp->np; j++) {
266 + #ifdef  DEBUG
267 +                        if (dp->t != i || dp->p != j)
268 +                                error(CONSISTENCY,
269 +                                        "division order in comperrs");
270 + #endif
271                          b = bright(dp[0].v);
272                          if (i > 0) {            /* from above */
273                                  b2 = bright(dp[-hp->np].v) - b;

Diff Legend

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