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.13 by greg, Fri Sep 16 15:09:21 2016 UTC vs.
Revision 2.18 by greg, Mon Sep 9 02:48:03 2024 UTC

# Line 57 | Line 57 | getface(                               /* get arguments for a face */
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 98 | Line 98 | getface(                               /* get arguments for a face */
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 117 | Line 117 | freeface(                      /* free memory associated with face */
117          o->os = NULL;
118   }
119  
120 –
120   int
121   inface(                         /* determine if point is in face */
122          FVECT  p,
# Line 140 | Line 139 | inface(                                /* determine if point is in face */
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
> Changed lines