--- ray/src/hd/rhdisp3.c	1997/11/26 20:12:19	3.7
+++ ray/src/hd/rhdisp3.c	1998/01/06 15:08:50	3.9
@@ -84,19 +84,20 @@ getback:
  * The ray directions that define the pyramid in visit_cells() needn't
  * be normalized, but they must be given in clockwise order as seen
  * from the pyramid's apex (origin).
+ * If no cell centers fall within the domain, the closest cell is visited.
  */
 int
 visit_cells(orig, pyrd, hp, vf, dp)	/* visit cells within a pyramid */
 FVECT	orig, pyrd[4];		/* pyramid ray directions in clockwise order */
-HOLO	*hp;
+register HOLO	*hp;
 int	(*vf)();
 char	*dp;
 {
-	int	n = 0;
+	int	ncalls = 0, n = 0;
 	int	inflags = 0;
 	FVECT	gp, pn[4], lo, ld;
 	double	po[4], lbeg, lend, d, t;
-	GCOORD	gc;
+	GCOORD	gc, gc2[2];
 	register int	i;
 					/* figure out whose side we're on */
 	hdgrid(gp, hp, orig);
@@ -114,17 +115,17 @@ char	*dp;
 		if (!(inflags & 1<<gc.w))	/* origin on wrong side */
 			continue;
 					/* scanline algorithm */
-		for (gc.i[1] = hp->grid[((gc.w>>1)+2)%3]; gc.i[1]--; ) {
+		for (gc.i[1] = hp->grid[hdwg1[gc.w]]; gc.i[1]--; ) {
 						/* compute scanline */
 			gp[gc.w>>1] = gc.w&1 ? hp->grid[gc.w>>1] : 0;
-			gp[((gc.w>>1)+1)%3] = 0;
-			gp[((gc.w>>1)+2)%3] = gc.i[1] + 0.5;
+			gp[hdwg0[gc.w]] = 0;
+			gp[hdwg1[gc.w]] = gc.i[1] + 0.5;
 			hdworld(lo, hp, gp);
-			gp[((gc.w>>1)+1)%3] = 1;
+			gp[hdwg0[gc.w]] = 1;
 			hdworld(ld, hp, gp);
 			ld[0] -= lo[0]; ld[1] -= lo[1]; ld[2] -= lo[2];
 						/* find scanline limits */
-			lbeg = 0; lend = hp->grid[((gc.w>>1)+1)%3];
+			lbeg = 0; lend = hp->grid[hdwg0[gc.w]];
 			for (i = 0; i < 4; i++) {
 				t = DOT(pn[i], lo) - po[i];
 				d = -DOT(pn[i], ld);
@@ -142,11 +143,26 @@ char	*dp;
 			if (lbeg >= lend)
 				continue;
 			i = lend + .5;		/* visit cells on this scan */
-			for (gc.i[0] = lbeg + .5; gc.i[0] < i; gc.i[0]++)
+			for (gc.i[0] = lbeg + .5; gc.i[0] < i; gc.i[0]++) {
 				n += (*vf)(&gc, dp);
+				ncalls++;
+			}
 		}
 	}
-	return(n);
+	if (ncalls)			/* got one at least */
+		return(n);
+					/* else find closest cell */
+	VSUM(ld, pyrd[0], pyrd[1], 1.);
+	VSUM(ld, ld, pyrd[2], 1.);
+	VSUM(ld, ld, pyrd[3], 1.);
+#if 0
+	if (normalize(ld) == 0.0)	/* technically not necessary */
+		return(0);
+#endif
+	d = hdinter(gc2, NULL, &t, hp, orig, ld);
+	if (d >= FHUGE || t <= 0.)
+		return(0);
+	return((*vf)(gc2+1, dp));	/* visit it */
 }
 
 
@@ -283,8 +299,8 @@ int	(*f)();
 					/* do each wall on each section */
 	for (hd = 0; hdlist[hd] != NULL; hd++)
 		for (w = 0; w < 6; w++) {
-			g0 = ((w>>1)+1)%3;
-			g1 = ((w>>1)+2)%3;
+			g0 = hdwg0[w];
+			g1 = hdwg1[w];
 			d = 1.0/hdlist[hd]->grid[g0];
 			mov[0] = d * hdlist[hd]->xv[g0][0];
 			mov[1] = d * hdlist[hd]->xv[g0][1];