--- ray/src/common/face.c	2003/02/25 02:47:21	2.7
+++ ray/src/common/face.c	2003/03/14 21:27:45	2.9
@@ -1,5 +1,5 @@
 #ifndef lint
-static const char	RCSid[] = "$Id: face.c,v 2.7 2003/02/25 02:47:21 greg Exp $";
+static const char RCSid[] = "$Id: face.c,v 2.9 2003/03/14 21:27:45 greg Exp $";
 #endif
 /*
  *  face.c - routines dealing with polygonal faces.
@@ -25,9 +25,9 @@ static const char	RCSid[] = "$Id: face.c,v 2.7 2003/02
  */
 
 #ifdef  SMLFLT
-#define  VERTEPS	1e-2		/* allowed vertex error */
+#define  VERTEPS	1e-3		/* allowed vertex error */
 #else
-#define  VERTEPS	1e-4		/* allowed vertex error */
+#define  VERTEPS	1e-5		/* allowed vertex error */
 #endif
 
 
@@ -86,11 +86,11 @@ OBJREC  *o;
 	f->offset = DOT(f->norm, VERTEX(f,0));
 	for (i = 1; i < f->nv; i++) {
 		d1 = DOT(f->norm, VERTEX(f,i));
-		badvert += fabs(d1 - f->offset/i) > VERTEPS;
+		badvert += fabs(1.0 - d1*i/f->offset) > VERTEPS;
 		f->offset += d1;
 	}
 	f->offset /= (double)f->nv;
-	if (badvert)
+	if (f->nv > 3 && badvert)
 		objerror(o, WARNING, "non-planar vertex");
 						/* find axis */
 	f->ax = fabs(f->norm[0]) > fabs(f->norm[1]) ? 0 : 1;
@@ -119,6 +119,7 @@ FACE  *f;
 {
 	int  ncross, n;
 	double  x, y;
+	int  tst;
 	register int  xi, yi;
 	register FLOAT  *p0, *p1;
 
@@ -132,13 +133,15 @@ FACE  *f;
 	ncross = 0;
 					/* positive x axis cross test */
 	while (n--) {
-		if ((p0[yi] > y) ^ (p1[yi] > y))
-			if (p0[xi] > x && p1[xi] > x)
+		if ((p0[yi] > y) ^ (p1[yi] > y)) {
+			tst = (p0[xi] > x) + (p1[xi] > x);
+			if (tst == 2)
 				ncross++;
-			else if (p0[xi] > x || p1[xi] > x)
+			else if (tst)
 				ncross += (p1[yi] > p0[yi]) ^
 						((p0[yi]-y)*(p1[xi]-x) >
 						(p0[xi]-x)*(p1[yi]-y));
+		}
 		p0 = p1;
 		p1 += 3;
 	}