--- ray/src/util/rsensor.c	2022/03/15 18:05:03	2.21
+++ ray/src/util/rsensor.c	2025/04/23 02:35:27	2.24
@@ -1,5 +1,5 @@
 #ifndef lint
-static const char RCSid[] = "$Id: rsensor.c,v 2.21 2022/03/15 18:05:03 greg Exp $";
+static const char RCSid[] = "$Id: rsensor.c,v 2.24 2025/04/23 02:35:27 greg Exp $";
 #endif
 
 /*
@@ -11,6 +11,7 @@ static const char RCSid[] = "$Id: rsensor.c,v 2.21 202
 #include "ray.h"
 #include "platform.h"
 #include "source.h"
+#include "func.h"
 #include "view.h"
 #include "random.h"
 
@@ -101,6 +102,8 @@ main(
 	directrelay = 3;
 	ambounce = 1;
 	maxdepth = -10;
+				/* initialize calcomp routines */
+	initfunc();
 				/* get options from command line */
 	for (i = 1; i < argc; i++) {
 		while ((rval = expandarg(&argc, &argv, i)) > 0)
@@ -521,7 +524,7 @@ comp_sensor(
 				ndsamps > 0 ? 1 : 0;
 	char	*err;
 	int	nt, np;
-	COLOR	vsum;
+	SCOLOR	vsum;
 	RAY	rr;
 	double	sf;
 	int	i, j;
@@ -535,7 +538,7 @@ comp_sensor(
 						/* assign probability table */
 	init_ptable(sfile);
 						/* stratified MC sampling */
-	setcolor(vsum, .0f, .0f, .0f);
+	scolorblack(vsum);
 	nt = (int)(sqrt((double)nsamps*ntheta/nphi) + .5);
 	np = nsamps/nt;
 	sf = gscale/nsamps;
@@ -551,9 +554,9 @@ comp_sensor(
 			}
 			rr.rmax = .0;
 			rayorigin(&rr, PRIMARY|SPECULAR, NULL, NULL);
-			scalecolor(rr.rcoef, sf);
+			scalescolor(rr.rcoef, sf);
 			if (ray_pqueue(&rr) == 1)
-				addcolor(vsum, rr.rcol);
+				saddscolor(vsum, rr.rcol);
 		}
 						/* remaining rays pure MC */
 	for (i = nsamps - nt*np; i-- > 0; ) {
@@ -567,14 +570,14 @@ comp_sensor(
 		}
 		rr.rmax = .0;
 		rayorigin(&rr, PRIMARY|SPECULAR, NULL, NULL);
-		scalecolor(rr.rcoef, sf);
+		scalescolor(rr.rcoef, sf);
 		if (ray_pqueue(&rr) == 1)
-			addcolor(vsum, rr.rcol);
+			saddscolor(vsum, rr.rcol);
 	}
 	if (!ray_pnprocs)			/* just printing rays */
 		return;
 						/* scale partial result */
-	scalecolor(vsum, sf);
+	scalescolor(vsum, sf);
 						/* add direct component */
 	for (i = ndirs; i-- > 0; ) {
 		SRCINDEX	si;
@@ -584,19 +587,18 @@ comp_sensor(
 			if (sf <= FTINY)
 				continue;
 			sf *= si.dom/ndirs;
-			scalecolor(rr.rcoef, sf);
+			scalescolor(rr.rcoef, sf);
 			if (ray_pqueue(&rr) == 1) {
-				multcolor(rr.rcol, rr.rcoef);
-				addcolor(vsum, rr.rcol);
+				smultscolor(rr.rcol, rr.rcoef);
+				saddscolor(vsum, rr.rcol);
 			}
 		}
 	}
-						/* finish our calculation */
-	while (ray_presult(&rr, 0) > 0) {
-		multcolor(rr.rcol, rr.rcoef);
-		addcolor(vsum, rr.rcol);
+	while (ray_presult(&rr, 0) > 0) {	/* finish our calculation */
+		smultscolor(rr.rcol, rr.rcoef);
+		saddscolor(vsum, rr.rcol);
 	}
-						/* print our result */
-	printf("%.4e %.4e %.4e\n", colval(vsum,RED),
-				colval(vsum,GRN), colval(vsum,BLU));
+	for (i = 0; i < NCSAMP; i++)		/* print our result */
+		printf(" %.4e", vsum[i]);
+	fputc('\n', stdout);
 }