--- ray/src/common/face.c	1989/02/02 10:34:31	1.1
+++ ray/src/common/face.c	1995/08/24 20:54:54	2.4
@@ -1,4 +1,4 @@
-/* Copyright (c) 1986 Regents of the University of California */
+/* Copyright (c) 1995 Regents of the University of California */
 
 #ifndef lint
 static char SCCSid[] = "$SunId$ LBL";
@@ -27,14 +27,17 @@ static char SCCSid[] = "$SunId$ LBL";
  *  The last vertex is automatically connected to the first.
  */
 
+#ifdef  SMLFLT
+#define  VERTEPS	1e-2		/* allowed vertex error */
+#else
 #define  VERTEPS	1e-4		/* allowed vertex error */
+#endif
 
 
 FACE *
 getface(o)			/* get arguments for a face */
 OBJREC  *o;
 {
-	double  fabs();
 	double  d1;
 	int  badvert;
 	FVECT  v1, v2, v3;
@@ -51,8 +54,13 @@ OBJREC  *o;
 	if (o->oargs.nfargs < 9 || o->oargs.nfargs % 3)
 		objerror(o, USER, "bad # arguments");
 
+	o->os = (char *)f;			/* save face */
+
 	f->va = o->oargs.farg;
 	f->nv = o->oargs.nfargs / 3;
+						/* check for last==first */
+	if (dist2(VERTEX(f,0),VERTEX(f,f->nv-1)) <= FTINY*FTINY)
+		f->nv--;
 						/* compute area and normal */
 	f->norm[0] = f->norm[1] = f->norm[2] = 0.0;
 	v1[0] = v1[1] = v1[2] = 0.0;
@@ -69,20 +77,20 @@ OBJREC  *o;
 	f->area = normalize(f->norm);
 	if (f->area == 0.0) {
 		objerror(o, WARNING, "zero area");	/* used to be fatal */
-		f->const = 0.0;
+		f->offset = 0.0;
 		f->ax = 0;
 		return(f);
 	}
 	f->area *= 0.5;
-						/* compute constant */
+						/* compute offset */
 	badvert = 0;
-	f->const = DOT(f->norm, VERTEX(f,0));
+	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->const/i) > VERTEPS;
-		f->const += d1;
+		badvert += fabs(d1 - f->offset/i) > VERTEPS;
+		f->offset += d1;
 	}
-	f->const /= f->nv;
+	f->offset /= (double)f->nv;
 	if (badvert)
 		objerror(o, WARNING, "non-planar vertex");
 						/* find axis */
@@ -90,7 +98,6 @@ OBJREC  *o;
 	if (fabs(f->norm[2]) > fabs(f->norm[f->ax]))
 		f->ax = 2;
 
-	(FACE *)o->os = f;			/* save face */
 	return(f);
 }
 
@@ -98,6 +105,8 @@ OBJREC  *o;
 freeface(o)			/* free memory associated with face */
 OBJREC  *o;
 {
+	if (o->os == NULL)
+		return;
 	free(o->os);
 	o->os = NULL;
 }
@@ -110,7 +119,7 @@ FACE  *f;
 	int  ncross, n;
 	double  x, y;
 	register int  xi, yi;
-	register double  *p0, *p1;
+	register FLOAT  *p0, *p1;
 
 	xi = (f->ax+1)%3;
 	yi = (f->ax+2)%3;