--- ray/src/rt/srcsupp.c	1991/06/25 13:11:55	1.6
+++ ray/src/rt/srcsupp.c	1991/09/13 13:44:10	1.13
@@ -27,7 +27,7 @@ SRCFUNC  sfun[NUMOTYPE];		/* source dispatch table */
 
 initstypes()			/* initialize source dispatch table */
 {
-	extern VSMATERIAL  mirror_vs;
+	extern VSMATERIAL  mirror_vs, direct1_vs, direct2_vs;
 	extern int  fsetsrc(), ssetsrc(), sphsetsrc(), rsetsrc();
 	extern double  fgetplaneq(), rgetplaneq();
 	extern double  fgetmaxdisk(), rgetmaxdisk();
@@ -37,6 +37,8 @@ initstypes()			/* initialize source dispatch table */
 	static SOBJECT  rsobj = {rsetsrc, rgetplaneq, rgetmaxdisk};
 
 	sfun[MAT_MIRROR].mf = &mirror_vs;
+	sfun[MAT_DIRECT1].mf = &direct1_vs;
+	sfun[MAT_DIRECT2].mf = &direct2_vs;
 	sfun[OBJ_FACE].of = &fsobj;
 	sfun[OBJ_SOURCE].of = &ssobj;
 	sfun[OBJ_SPHERE].of = &sphsobj;
@@ -47,17 +49,19 @@ initstypes()			/* initialize source dispatch table */
 int
 newsource()			/* allocate new source in our array */
 {
+#define SRCINC	4
 	if (nsources == 0)
-		source = (SRCREC *)malloc(sizeof(SRCREC));
-	else
+		source = (SRCREC *)malloc(SRCINC*sizeof(SRCREC));
+	else if (nsources%SRCINC == 0)
 		source = (SRCREC *)realloc((char *)source,
-				(unsigned)(nsources+1)*sizeof(SRCREC));
+				(unsigned)(nsources+SRCINC)*sizeof(SRCREC));
 	if (source == NULL)
 		return(-1);
 	source[nsources].sflags = 0;
 	source[nsources].nhits = 1;
 	source[nsources].ntests = 2;	/* initial hit probability = 1/2 */
 	return(nsources++);
+#undef SRCINC
 }
 
 
@@ -105,7 +109,7 @@ register OBJREC  *so;
 	theta = PI/180.0/2.0 * so->oargs.farg[3];
 	if (theta <= FTINY)
 		objerror(so, USER, "zero size");
-	src->ss = theta >= PI/4 ? 1.0 : tan(theta);
+	src->ss = theta >= PI/4.0 ? 1.0 : tan(theta);
 	src->ss2 = 2.0*PI * (1.0 - cos(theta));
 }
 
@@ -288,10 +292,10 @@ FVECT  nrm;		/* source surface normal */
 
 	d = DOT(sp->aim, nrm);
 	if (d > FTINY)			/* center in front? */
-		return(0);
+		return(1);
 					/* else check horizon */
 	d1 = 1. - sp->siz/(2.*PI);
-	return(1.-FTINY-d*d > d1*d1);
+	return(1.-FTINY-d*d < d1*d1);
 }
 
 
@@ -390,7 +394,7 @@ register RAY  *r;
 		first = 0; last = nsources-1;
 	}
 	for (i = first; i <= last; i++)
-		if (source[i].sflags & SDISTANT)
+		if ((source[i].sflags & (SDISTANT|SVIRTUAL)) == SDISTANT)
 			/*
 			 * Check to see if ray is within
 			 * solid angle of source.
@@ -413,17 +417,24 @@ register RAY  *r;
 }
 
 
-#define  wrongsource(m, r)	(m->otype!=MAT_ILLUM && \
-				r->rsrc>=0 && \
-				source[r->rsrc].so!=r->ro)
+#define  wrongsource(m, r)	(r->rsrc>=0 && \
+				source[r->rsrc].so!=r->ro && \
+				(m->otype!=MAT_ILLUM || \
+			objptr(source[r->rsrc].so->omod)->otype==MAT_ILLUM))
 
+#define  distglow(m, r)		(m->otype==MAT_GLOW && \
+				r->rot > m->oargs.farg[3])
+
 #define  badambient(m, r)	((r->crtype&(AMBIENT|SHADOW))==AMBIENT && \
-				!(m->otype==MAT_GLOW&&r->rot>m->oargs.farg[3]))
+				!distglow(m, r))
 
 #define  passillum(m, r)	(m->otype==MAT_ILLUM && \
 				!(r->rsrc>=0&&source[r->rsrc].so==r->ro))
 
+#define  srcignore(m, r)	(directinvis && !(r->crtype&SHADOW) && \
+				!distglow(m, r))
 
+
 m_light(m, r)			/* ray hit a light source */
 register OBJREC  *m;
 register RAY  *r;
@@ -433,24 +444,25 @@ register RAY  *r;
 		return;
 						/* check for passed illum */
 	if (passillum(m, r)) {
-
 		if (m->oargs.nsargs < 1 || !strcmp(m->oargs.sarg[0], VOIDID))
 			raytrans(r);
 		else
 			rayshade(r, modifier(m->oargs.sarg[0]));
-
-						/* otherwise treat as source */
-	} else {
+		return;
+	}
+					/* otherwise treat as source */
 						/* check for behind */
-		if (r->rod < 0.0)
-			return;
+	if (r->rod < 0.0)
+		return;
+						/* check for invisibility */
+	if (srcignore(m, r))
+		return;
 						/* get distribution pattern */
-		raytexture(r, m->omod);
+	raytexture(r, m->omod);
 						/* get source color */
-		setcolor(r->rcol, m->oargs.farg[0],
-				  m->oargs.farg[1],
-				  m->oargs.farg[2]);
+	setcolor(r->rcol, m->oargs.farg[0],
+			  m->oargs.farg[1],
+			  m->oargs.farg[2]);
 						/* modify value */
-		multcolor(r->rcol, r->pcol);
-	}
+	multcolor(r->rcol, r->pcol);
 }