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

Comparing ray/src/rt/srcsupp.c (file contents):
Revision 1.4 by greg, Fri Jun 21 16:44:01 1991 UTC vs.
Revision 1.5 by greg, Mon Jun 24 16:00:30 1991 UTC

# Line 168 | Line 168 | FVECT  ocent;
168   OBJREC  *op;
169   {
170          double  maxrad2;
171 <        double  d2;
171 >        double  d;
172          register int  i, j;
173          register FACE  *f;
174          
175          f = getface(op);
176 +        if (f->area == 0.)
177 +                return(0.);
178          for (i = 0; i < 3; i++) {
179                  ocent[i] = 0.;
180                  for (j = 0; j < f->nv; j++)
181                          ocent[i] += VERTEX(f,j)[i];
182                  ocent[i] /= (double)f->nv;
183          }
184 <        if (f->area == 0.)
185 <                return(0.);
184 >        d = DOT(ocent,f->norm);
185 >        for (i = 0; i < 3; i++)
186 >                ocent[i] += (f->offset - d)*f->norm[i];
187          maxrad2 = 0.;
188          for (j = 0; j < f->nv; j++) {
189 <                d2 = dist2(VERTEX(f,j), ocent);
190 <                if (d2 > maxrad2)
191 <                        maxrad2 = d2;
189 >                d = dist2(VERTEX(f,j), ocent);
190 >                if (d > maxrad2)
191 >                        maxrad2 = d;
192          }
193          return(maxrad2);
194   }

Diff Legend

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