--- ray/src/rt/source.c	1993/04/06 17:55:51	2.11
+++ ray/src/rt/source.c	1993/12/21 15:46:19	2.13
@@ -299,8 +299,10 @@ char  *p;			/* data for f */
 						/* modify threshold */
 	ourthresh = shadthresh / r->rweight;
 						/* test for shadows */
-	nhits = 0;
-	for (sn = 0; sn < ncnts; sn++) {
+	for (nhits = 0, hwt = 0.0, sn = 0; sn < ncnts;
+			hwt += (double)source[scp->sno].nhits /
+				(double)source[scp->sno].ntests,
+			sn++) {
 						/* check threshold */
 		if ((sn+nshadcheck>=ncnts ? cntord[sn].brt :
 				cntord[sn].brt-cntord[sn+nshadcheck].brt)
@@ -329,13 +331,13 @@ char  *p;			/* data for f */
 		nhits++;
 		source[scp->sno].nhits++;
 	}
-					/* surface hit rate */
-	if (sn > 0)
-		hwt = (double)nhits / (double)sn;
+					/* source hit rate */
+	if (hwt > FTINY)
+		hwt = (double)nhits / hwt;
 	else
 		hwt = 0.5;
 #ifdef DEBUG
-	sprintf(errmsg, "%d tested, %d untested, %f hit rate\n",
+	sprintf(errmsg, "%d tested, %d untested, %f conditional hit rate\n",
 			sn, ncnts-sn, hwt);
 	eputs(errmsg);
 #endif
@@ -344,6 +346,8 @@ char  *p;			/* data for f */
 		scp = srccnt + cntord[sn].sndx;
 		prob = hwt * (double)source[scp->sno].nhits /
 				(double)source[scp->sno].ntests;
+		if (prob > 1.0)
+			prob = 1.0;
 		scalecolor(scp->val, prob);
 		addcolor(r->rcol, scp->val);
 	}
@@ -371,6 +375,7 @@ static int weaksrcmod(obj) int obj;	/* efficiency boos
 return(o->otype==MAT_ILLUM|o->otype==MAT_GLOW);}
 
 #define  illumblock(m, r)	(!(source[r->rsrc].sflags&SVIRTUAL) && \
+				r->rod > 0.0 && \
 				weaksrcmod(source[r->rsrc].so->omod))
 
 /* wrongsource *
@@ -407,13 +412,6 @@ return(o->otype==MAT_ILLUM|o->otype==MAT_GLOW);}
 				!(r->crtype&SHADOW || r->rod < 0.0 || \
 					distglow(m, r)))
 
-/* overcount *
- *
- * All overcounting possibilities are contained here.
- */
-
-#define  overcount(m, r)	(badcomponent(m,r) || wrongsource(m,r))
-
 /* passillum *
  *
  * An illum passes to another material type when we didn't hit it
@@ -439,7 +437,9 @@ register OBJREC  *m;
 register RAY  *r;
 {
 						/* check for over-counting */
-	if (overcount(m, r))
+	if (badcomponent(m, r))
+		return;
+	if (wrongsource(m,r))
 		return;
 						/* check for passed illum */
 	if (passillum(m, r)) {