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

Comparing ray/src/common/face.c (file contents):
Revision 2.9 by greg, Fri Mar 14 21:27:45 2003 UTC vs.
Revision 2.13 by greg, Fri Sep 16 15:09:21 2016 UTC

# Line 32 | Line 32 | static const char RCSid[] = "$Id$";
32  
33  
34   FACE *
35 < getface(o)                      /* get arguments for a face */
36 < OBJREC  *o;
35 > getface(                                /* get arguments for a face */
36 >        OBJREC  *o
37 > )
38   {
39          double  d1;
40 <        int  badvert;
40 >        int  smalloff, badvert;
41          FVECT  v1, v2, v3;
42 <        register FACE  *f;
43 <        register int  i;
42 >        FACE  *f;
43 >        int  i;
44  
45          if ((f = (FACE *)o->os) != NULL)
46                  return(f);                      /* already done */
# Line 84 | Line 85 | OBJREC  *o;
85                                                  /* compute offset */
86          badvert = 0;
87          f->offset = DOT(f->norm, VERTEX(f,0));
88 +        smalloff = fabs(f->offset) <= VERTEPS;
89          for (i = 1; i < f->nv; i++) {
90                  d1 = DOT(f->norm, VERTEX(f,i));
91 <                badvert += fabs(1.0 - d1*i/f->offset) > VERTEPS;
91 >                if (smalloff)
92 >                        badvert += fabs(d1 - f->offset/i) > VERTEPS;
93 >                else
94 >                        badvert += fabs(1.0 - d1*i/f->offset) > VERTEPS;
95                  f->offset += d1;
96          }
97          f->offset /= (double)f->nv;
# Line 102 | Line 107 | OBJREC  *o;
107  
108  
109   void
110 < freeface(o)                     /* free memory associated with face */
111 < OBJREC  *o;
110 > freeface(                       /* free memory associated with face */
111 >        OBJREC  *o
112 > )
113   {
114          if (o->os == NULL)
115                  return;
# Line 113 | Line 119 | OBJREC  *o;
119  
120  
121   int
122 < inface(p, f)                    /* determine if point is in face */
123 < FVECT  p;
124 < FACE  *f;
122 > inface(                         /* determine if point is in face */
123 >        FVECT  p,
124 >        FACE  *f
125 > )
126   {
127          int  ncross, n;
128          double  x, y;
129          int  tst;
130 <        register int  xi, yi;
131 <        register FLOAT  *p0, *p1;
130 >        int  xi, yi;
131 >        RREAL  *p0, *p1;
132  
133 <        xi = (f->ax+1)%3;
134 <        yi = (f->ax+2)%3;
133 >        if ((xi = f->ax + 1) >= 3) xi -= 3;
134 >        if ((yi = xi + 1) >= 3) yi -= 3;
135          x = p[xi];
136          y = p[yi];
137          n = f->nv;

Diff Legend

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