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

Comparing ray/src/common/face.c (file contents):
Revision 2.12 by greg, Thu Mar 2 16:51:54 2006 UTC vs.
Revision 2.18 by greg, Mon Sep 9 02:48:03 2024 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  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 56 | Line 57 | OBJREC  *o;
57          f->va = o->oargs.farg;
58          f->nv = o->oargs.nfargs / 3;
59                                                  /* check for last==first */
60 <        if (dist2(VERTEX(f,0),VERTEX(f,f->nv-1)) <= FTINY*FTINY)
60 >        if (f->nv > 3 && dist2(VERTEX(f,0),VERTEX(f,f->nv-1)) <= FTINY*FTINY)
61                  f->nv--;
62                                                  /* compute area and normal */
63          f->norm[0] = f->norm[1] = f->norm[2] = 0.0;
# Line 97 | Line 98 | OBJREC  *o;
98          if (f->nv > 3 && badvert)
99                  objerror(o, WARNING, "non-planar vertex");
100                                                  /* find axis */
101 <        f->ax = fabs(f->norm[0]) > fabs(f->norm[1]) ? 0 : 1;
101 >        f->ax = (fabs(f->norm[1]) > fabs(f->norm[0]));
102          if (fabs(f->norm[2]) > fabs(f->norm[f->ax]))
103                  f->ax = 2;
104  
# Line 106 | 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 115 | Line 117 | OBJREC  *o;
117          o->os = NULL;
118   }
119  
118 –
120   int
121 < inface(p, f)                    /* determine if point is in face */
122 < FVECT  p;
123 < FACE  *f;
121 > inface(                         /* determine if point is in face */
122 >        FVECT  p,
123 >        FACE  *f
124 > )
125   {
126          int  ncross, n;
127          double  x, y;
128          int  tst;
129 <        register int  xi, yi;
130 <        register RREAL  *p0, *p1;
129 >        int  xi, yi;
130 >        RREAL  *p0, *p1;
131  
132          if ((xi = f->ax + 1) >= 3) xi -= 3;
133          if ((yi = xi + 1) >= 3) yi -= 3;
# Line 137 | Line 139 | FACE  *f;
139          ncross = 0;
140                                          /* positive x axis cross test */
141          while (n--) {
142 +                if (FABSEQ(p0[yi], y) && FABSEQ(p1[yi], y) &&
143 +                                (p0[xi] > x) ^ (p1[xi] > x))
144 +                        return(1);                      /* edge case #1 */
145                  if ((p0[yi] > y) ^ (p1[yi] > y)) {
146                          tst = (p0[xi] > x) + (p1[xi] > x);
147                          if (tst == 2)
148                                  ncross++;
149 <                        else if (tst)
150 <                                ncross += (p1[yi] > p0[yi]) ^
151 <                                                ((p0[yi]-y)*(p1[xi]-x) >
152 <                                                (p0[xi]-x)*(p1[yi]-y));
149 >                        else if (tst) {
150 >                                double  prodA = (p0[yi]-y)*(p1[xi]-x);
151 >                                double  prodB = (p0[xi]-x)*(p1[yi]-y);
152 >                                if (FABSEQ(prodA, prodB))
153 >                                        return(1);      /* edge case #2 */
154 >                                ncross += (p1[yi] > p0[yi]) ^ (prodA > prodB);
155 >                        } else if (FABSEQ(p0[xi], x) && FABSEQ(p1[xi], x))
156 >                                return(1);              /* edge case #3 */
157                  }
158                  p0 = p1;
159                  p1 += 3;

Diff Legend

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