--- ray/src/rt/rtrace.c	2005/09/19 11:30:10	2.49
+++ ray/src/rt/rtrace.c	2006/04/05 06:22:57	2.53
@@ -1,5 +1,5 @@
 #ifndef lint
-static const char	RCSid[] = "$Id: rtrace.c,v 2.49 2005/09/19 11:30:10 schorsch Exp $";
+static const char	RCSid[] = "$Id: rtrace.c,v 2.53 2006/04/05 06:22:57 greg Exp $";
 #endif
 /*
  *  rtrace.c - program and variables for individual ray tracing.
@@ -54,7 +54,7 @@ void  (*trace)() = NULL;		/* trace call */
 char  *tralist[128];			/* list of modifers to trace (or no) */
 int  traincl = -1;			/* include == 1, exclude == 0 */
 #ifndef  MAXTSET
-#define	 MAXTSET	1024		/* maximum number in trace set */
+#define	 MAXTSET	8192		/* maximum number in trace set */
 #endif
 OBJECT	traset[MAXTSET+1]={0};		/* trace include/exclude set */
 
@@ -101,7 +101,7 @@ typedef void putf_t(double v);
 static putf_t puta, putd, putf;
 
 typedef void oputf_t(RAY *r);
-static oputf_t  oputo, oputd, oputv, oputl, oputL, oputc, oputp,
+static oputf_t  oputo, oputd, oputv, oputV, oputl, oputL, oputc, oputp,
 		oputn, oputN, oputs, oputw, oputW, oputm, oputM, oputtilde;
 
 static void setoutput(char *vs);
@@ -198,7 +198,7 @@ rtrace(				/* trace rays from file */
 				nextflush = hresolu;
 			}
 		} else {
-			samplendx = rand_samp ? random() : samplendx+1;
+			samplendx++;
 							/* compute and print */
 			if (imm_irrad)
 				irrad(orig, direc);
@@ -215,7 +215,8 @@ rtrace(				/* trace rays from file */
 		if (vcount && !--vcount)		/* check for end */
 			break;
 	}
-	fflush(stdout);
+	if (fflush(stdout) < 0)
+		error(SYSTEM, "write error");
 	if (vcount)
 		error(USER, "unexpected EOF on input");
 	if (fname != NULL)
@@ -264,6 +265,12 @@ setoutput(				/* set up output tables */
 			*table++ = oputv;
 			castonly = 0;
 			break;
+		case 'V':				/* contribution */
+			*table++ = oputV;
+			if (ambounce > 0 && (ambacc > FTINY || ambssamp > 0))
+				error(WARNING,
+					"-otV accuracy depends on -aa 0 -as 0");
+			break;
 		case 'l':				/* effective distance */
 			*table++ = oputl;
 			castonly = 0;
@@ -531,6 +538,21 @@ oputv(				/* print value */
 
 
 static void
+oputV(				/* print value contribution */
+	RAY *r
+)
+{
+	double	contr[3];
+
+	raycontrib(contr, r, PRIMARY);
+	multcolor(contr, r->rcol);
+	(*putreal)(contr[RED]);
+	(*putreal)(contr[GRN]);
+	(*putreal)(contr[BLU]);
+}
+
+
+static void
 oputl(				/* print effective distance */
 	RAY  *r
 )
@@ -635,7 +657,7 @@ oputw(				/* print weight */
 
 
 static void
-oputW(				/* print contribution */
+oputW(				/* print coefficient */
 	RAY  *r
 )
 {