--- ray/src/cv/mgf2rad.c	2003/07/27 22:12:01	2.27
+++ ray/src/cv/mgf2rad.c	2024/01/17 00:43:45	2.35
@@ -1,5 +1,5 @@
 #ifndef lint
-static const char	RCSid[] = "$Id: mgf2rad.c,v 2.27 2003/07/27 22:12:01 schorsch Exp $";
+static const char	RCSid[] = "$Id: mgf2rad.c,v 2.35 2024/01/17 00:43:45 greg Exp $";
 #endif
 /*
  * Convert MGF (Materials and Geometry Format) to Radiance
@@ -9,9 +9,12 @@ static const char	RCSid[] = "$Id: mgf2rad.c,v 2.27 200
 #include <stdlib.h>
 #include <math.h>
 #include <string.h>
-#include "mgflib/parser.h"
+
+#include "platform.h"
+#include "mgf_parser.h"
 #include "color.h"
 #include "tmesh.h"
+#include "lookup.h"
 
 #define putv(v)		printf("%18.12g %18.12g %18.12g\n",(v)[0],(v)[1],(v)[2])
 
@@ -23,13 +26,31 @@ double  emult = 1.;			/* emitter multiplier */
 
 FILE	*matfp;				/* material output file */
 
-int	r_comment(), r_cone(), r_cyl(), r_face(), r_ies(), r_ring(), r_sph();
-char	*material(), *object(), *addarg();
+int	dospectra = 0;			/* output spectral colors? */
 
 
-main(argc, argv)		/* convert files to stdout */
-int	argc;
-char	*argv[];
+extern int r_comment(int ac, char **av);
+extern int r_color(int ac, char **av);
+extern int r_cone(int ac, char **av);
+extern int r_cyl(int ac, char **av);
+extern int r_sph(int ac, char **av);
+extern int r_ring(int ac, char **av);
+extern int r_face(int ac, char **av);
+extern int r_ies(int ac, char **av);
+extern void putsided(char *mname);
+extern char * material(void);
+extern char * object(void);
+extern char * addarg(char *op, char *arg);
+extern void do_tri(char *mat, C_VERTEX *cv1, C_VERTEX *cv2, C_VERTEX *cv3, int iv);
+extern void cvtcolor(COLOR radrgb, C_COLOR *ciec, double intensity);
+extern char * specolor(COLOR radrgb, C_COLOR *ciec, double intensity);
+
+
+int
+main(
+	int	argc,
+	char	*argv[]
+)
 {
 	int	i;
 
@@ -41,9 +62,9 @@ char	*argv[];
 	mg_ehand[MG_E_COLOR] = c_hcolor;	/* they get color */
 	mg_ehand[MG_E_CONE] = r_cone;		/* we do cones */
 	mg_ehand[MG_E_CMIX] = c_hcolor;		/* they mix colors */
-	mg_ehand[MG_E_CSPEC] = c_hcolor;	/* they get spectra */
 	mg_ehand[MG_E_CXY] = c_hcolor;		/* they get chromaticities */
-	mg_ehand[MG_E_CCT] = c_hcolor;		/* they get color temp's */
+	mg_ehand[MG_E_CSPEC] = r_color;		/* we get spectra */
+	mg_ehand[MG_E_CCT] = r_color;		/* we get color temp's */
 	mg_ehand[MG_E_CYL] = r_cyl;		/* we do cylinders */
 	mg_ehand[MG_E_ED] = c_hmaterial;	/* they get emission */
 	mg_ehand[MG_E_FACE] = r_face;		/* we do faces */
@@ -88,6 +109,9 @@ char	*argv[];
 			}
 			printf(" %s", argv[i]);
 			break;
+		case 's':			/* spectral color output? */
+			dospectra = !dospectra;
+			break;
 		default:
 			goto userr;
 		}
@@ -113,16 +137,17 @@ char	*argv[];
 		}
 	exit(0);
 userr:
-	fprintf(stderr, "Usage: %s [-g dist][-e mult][-m matf] [file.mgf] ..\n",
+	fprintf(stderr, "Usage: %s [-s][-g dist][-e mult][-m matf] [file.mgf] ..\n",
 			argv[0]);
 	exit(1);
 }
 
 
 int
-r_comment(ac, av)		/* repeat a comment verbatim */
-register int	ac;
-register char	**av;
+r_comment(		/* repeat a comment verbatim */
+	int	ac,
+	char	**av
+)
 {
 	putchar('#');		/* use Radiance comment character */
 	while (--ac) {			/* pass through verbatim */
@@ -135,10 +160,26 @@ register char	**av;
 
 
 int
-r_cone(ac, av)			/* put out a cone */
-int	ac;
-char	**av;
+r_color(		/* call color handler & remember name */
+	int	ac,
+	char	**av
+)
 {
+	int	rval = c_hcolor(ac, av);
+
+	if (rval == MG_OK)
+		c_ccolor->client_data = c_ccname;
+
+	return(rval);
+}
+
+
+int
+r_cone(			/* put out a cone */
+	int	ac,
+	char	**av
+)
+{
 	static int	ncones;
 	char	*mat;
 	double	r1, r2;
@@ -183,9 +224,10 @@ char	**av;
 
 
 int
-r_cyl(ac, av)			/* put out a cylinder */
-int	ac;
-char	**av;
+r_cyl(			/* put out a cylinder */
+	int	ac,
+	char	**av
+)
 {
 	static int	ncyls;
 	char	*mat;
@@ -221,9 +263,10 @@ char	**av;
 
 
 int
-r_sph(ac, av)			/* put out a sphere */
-int	ac;
-char	**av;
+r_sph(			/* put out a sphere */
+	int	ac,
+	char	**av
+)
 {
 	static int	nsphs;
 	char	*mat;
@@ -254,9 +297,10 @@ char	**av;
 
 
 int
-r_ring(ac, av)			/* put out a ring */
-int	ac;
-char	**av;
+r_ring(			/* put out a ring */
+	int	ac,
+	char	**av
+)
 {
 	static int	nrings;
 	char	*mat;
@@ -291,17 +335,18 @@ char	**av;
 
 
 int
-r_face(ac, av)			/* convert a face */
-int	ac;
-char	**av;
+r_face(			/* convert a face */
+	int	ac,
+	char	**av
+)
 {
 	static int	nfaces;
 	int		myi = invert;
 	char	*mat;
-	register int	i;
-	register C_VERTEX	*cv;
+	int	i;
+	C_VERTEX	*cv;
 	FVECT	v;
-	int	rv;
+
 					/* check argument count and type */
 	if (ac < 4)
 		return(MG_EARGC);
@@ -350,16 +395,17 @@ char	**av;
 
 
 int
-r_ies(ac, av)				/* convert an IES luminaire file */
-int	ac;
-char	**av;
+r_ies(				/* convert an IES luminaire file */
+	int	ac,
+	char	**av
+)
 {
 	int	xa0 = 2;
 	char	combuf[128];
 	char	fname[48];
 	char	*oname;
-	register char	*op;
-	register int	i;
+	char	*op;
+	int	i;
 					/* check argument count */
 	if (ac < 2)
 		return(MG_EARGC);
@@ -414,10 +460,14 @@ char	**av;
 }
 
 
-do_tri(mat, cv1, cv2, cv3, iv)		/* put out smoothed triangle */
-char	*mat;
-C_VERTEX	*cv1, *cv2, *cv3;
-int	iv;
+void
+do_tri(		/* put out smoothed triangle */
+	char	*mat,
+	C_VERTEX	*cv1,
+	C_VERTEX	*cv2,
+	C_VERTEX	*cv3,
+	int	iv
+)
 {
 	static int	ntris;
 	BARYCCM	bvecs;
@@ -425,7 +475,7 @@ int	iv;
 	C_VERTEX	*cvt;
 	FVECT	v1, v2, v3;
 	FVECT	n1, n2, n3;
-	register int	i;
+	int	i;
 
 	if (iv) {			/* swap vertex order if inverted */
 		cvt = cv1;
@@ -448,7 +498,7 @@ int	iv;
 		bcoor[i][1] = n2[i];
 		bcoor[i][2] = n3[i];
 	}
-	put_baryc(&bvecs, bcoor, 3);
+	fput_baryc(&bvecs, bcoor, 3, stdout);
 						/* put out triangle */
 	printf("\nT-nor polygon %st%d\n", object(), ++ntris);
 	printf("0\n0\n9\n");
@@ -458,13 +508,21 @@ int	iv;
 }
 
 
+void
+putsided(char *mname)		/* print out mixfunc for sided material */
+{
+	fprintf(matfp, "\nvoid mixfunc %s\n", mname);
+	fprintf(matfp, "4 %s void if(Rdot,1,0) .\n0\n0\n", mname);
+}
+
+
 char *
-material()			/* get (and print) current material */
+material(void)			/* get (and print) current material */
 {
 	char	*mname = "mat";
+	char	*pname;
 	COLOR	radrgb, c2;
 	double	d;
-	register int	i;
 
 	if (c_cmname != NULL)
 		mname = c_cmname;
@@ -473,15 +531,15 @@ material()			/* get (and print) current material */
 				/* else update output */
 	c_cmaterial->clock = 0;
 	if (c_cmaterial->ed > .1) {	/* emitter */
-		cvtcolor(radrgb, &c_cmaterial->ed_c,
+		pname = specolor(radrgb, &c_cmaterial->ed_c,
 				emult*c_cmaterial->ed/(PI*WHTEFFICACY));
 		if (glowdist < FHUGE) {		/* do a glow */
-			fprintf(matfp, "\nvoid glow %s\n0\n0\n", mname);
+			fprintf(matfp, "\n%s glow %s\n0\n0\n", pname, mname);
 			fprintf(matfp, "4 %f %f %f %f\n", colval(radrgb,RED),
 					colval(radrgb,GRN),
 					colval(radrgb,BLU), glowdist);
 		} else {
-			fprintf(matfp, "\nvoid light %s\n0\n0\n", mname);
+			fprintf(matfp, "\n%s light %s\n0\n0\n", pname, mname);
 			fprintf(matfp, "3 %f %f %f\n", colval(radrgb,RED),
 					colval(radrgb,GRN),
 					colval(radrgb,BLU));
@@ -518,73 +576,83 @@ material()			/* get (and print) current material */
 				colval(radrgb,GRN), colval(radrgb,BLU),
 				c_cmaterial->nr);
 		return(mname);
-		}
+	}
 					/* check for trans */
 	if (c_cmaterial->td > .01 || c_cmaterial->ts > .01) {
-		double	ts, a5, a6;
-
-		if (c_cmaterial->sided) {
-			ts = sqrt(c_cmaterial->ts);	/* approximate */
-			a5 = .5;
-		} else {
-			ts = c_cmaterial->ts;
-			a5 = 1.;
-		}
+		double	a5, a6;
 						/* average colors */
-		d = c_cmaterial->rd + c_cmaterial->td + ts;
+		d = c_cmaterial->rd + c_cmaterial->td + c_cmaterial->ts;
 		cvtcolor(radrgb, &c_cmaterial->rd_c, c_cmaterial->rd/d);
 		cvtcolor(c2, &c_cmaterial->td_c, c_cmaterial->td/d);
 		addcolor(radrgb, c2);
-		cvtcolor(c2, &c_cmaterial->ts_c, ts/d);
+		cvtcolor(c2, &c_cmaterial->ts_c, c_cmaterial->ts/d);
 		addcolor(radrgb, c2);
-		if (c_cmaterial->rs + ts > .0001)
+		if (c_cmaterial->rs + c_cmaterial->ts > .0001)
 			a5 = (c_cmaterial->rs*c_cmaterial->rs_a +
-					ts*a5*c_cmaterial->ts_a) /
-					(c_cmaterial->rs + ts);
-		a6 = (c_cmaterial->td + ts) /
-				(c_cmaterial->rd + c_cmaterial->td + ts);
+					c_cmaterial->ts*c_cmaterial->ts_a) /
+					(c_cmaterial->rs + c_cmaterial->ts);
+		a6 = (c_cmaterial->td + c_cmaterial->ts) /
+				(c_cmaterial->rd + c_cmaterial->td + c_cmaterial->ts);
 		if (a6 < .999)
 			d = c_cmaterial->rd/(1. - c_cmaterial->rs)/(1. - a6);
 		else
-			d = c_cmaterial->td + ts;
+			d = c_cmaterial->td + c_cmaterial->ts;
 		scalecolor(radrgb, d);
 		fprintf(matfp, "\nvoid trans %s\n0\n0\n", mname);
 		fprintf(matfp, "7 %f %f %f\n", colval(radrgb,RED),
 				colval(radrgb,GRN), colval(radrgb,BLU));
 		fprintf(matfp, "\t%f %f %f %f\n", c_cmaterial->rs, a5, a6,
-				ts/(ts + c_cmaterial->td));
+				c_cmaterial->ts/(c_cmaterial->ts + c_cmaterial->td));
+		if (c_cmaterial->sided)
+			putsided(mname);
 		return(mname);
 	}
 					/* check for plastic */
-	if (c_cmaterial->rs < .1) {
-		cvtcolor(radrgb, &c_cmaterial->rd_c,
+	if (c_cmaterial->rs < .1 && (c_cmaterial->rs < .1*c_cmaterial->rd ||
+					c_isgrey(&c_cmaterial->rs_c))) {
+		pname = specolor(radrgb, &c_cmaterial->rd_c,
 					c_cmaterial->rd/(1.-c_cmaterial->rs));
-		fprintf(matfp, "\nvoid plastic %s\n0\n0\n", mname);
+		fprintf(matfp, "\n%s plastic %s\n0\n0\n", pname, mname);
 		fprintf(matfp, "5 %f %f %f %f %f\n", colval(radrgb,RED),
 				colval(radrgb,GRN), colval(radrgb,BLU),
 				c_cmaterial->rs, c_cmaterial->rs_a);
+		if (c_cmaterial->sided)
+			putsided(mname);
 		return(mname);
 	}
 					/* else it's metal */
-						/* average colors */
-	cvtcolor(radrgb, &c_cmaterial->rd_c, c_cmaterial->rd);
-	cvtcolor(c2, &c_cmaterial->rs_c, c_cmaterial->rs);
-	addcolor(radrgb, c2);
-	fprintf(matfp, "\nvoid metal %s\n0\n0\n", mname);
+						/* compute color */
+	if (c_equiv(&c_cmaterial->rd_c, &c_cmaterial->rs_c)) {
+		pname = specolor(radrgb, &c_cmaterial->rs_c, c_cmaterial->rs+c_cmaterial->rd);
+	} else if (c_cmaterial->rd <= .05f) {
+		pname = specolor(radrgb, &c_cmaterial->rs_c, c_cmaterial->rs);
+		cvtcolor(c2, &c_cmaterial->rd_c, c_cmaterial->rd);
+		addcolor(radrgb, c2);
+	} else {
+		pname = "void";
+		cvtcolor(radrgb, &c_cmaterial->rd_c, c_cmaterial->rd);
+		cvtcolor(c2, &c_cmaterial->rs_c, c_cmaterial->rs);
+		addcolor(radrgb, c2);
+	}
+	fprintf(matfp, "\n%s metal %s\n0\n0\n", pname, mname);
 	fprintf(matfp, "5 %f %f %f %f %f\n", colval(radrgb,RED),
 			colval(radrgb,GRN), colval(radrgb,BLU),
 			c_cmaterial->rs/(c_cmaterial->rd + c_cmaterial->rs),
 			c_cmaterial->rs_a);
+	if (c_cmaterial->sided)
+		putsided(mname);
 	return(mname);
 }
 
 
-cvtcolor(radrgb, ciec, intensity)	/* convert a CIE XYZ color to RGB */
-COLOR	radrgb;
-register C_COLOR	*ciec;
-double	intensity;
+void
+cvtcolor(	/* convert a CIE XYZ color to RGB */
+	COLOR	radrgb,
+	C_COLOR	*ciec,
+	double	intensity
+)
 {
-	static COLOR	ciexyz;
+	COLOR	ciexyz;
 
 	c_ccvt(ciec, C_CSXY);		/* get xy representation */
 	ciexyz[1] = intensity;
@@ -594,12 +662,76 @@ double	intensity;
 }
 
 
+static int	/* new spectrum definition? */
+newspecdef(C_COLOR *spc)
+{
+	static LUTAB	spc_tab = LU_SINIT(NULL,free);
+	LUENT	*lp = lu_find(&spc_tab, (const char *)spc->client_data);
+
+	if (lp == NULL)			/* should never happen */
+		return(1);
+	if (lp->data == NULL) {		/* new entry */
+		lp->key = (char *)spc->client_data;
+		lp->data = (char *)malloc(sizeof(C_COLOR));
+	} else if (c_equiv(spc, (C_COLOR *)lp->data))
+		return(0);		/* unchanged */
+
+	if (lp->data != NULL)		/* else remember if we can */
+		*(C_COLOR *)lp->data = *spc;
+	return(1);			/* good as new */
+}
+
+
 char *
-object()			/* return current object name */
+specolor(	/* check if color has spectra and output accordingly */
+	COLOR	radrgb,
+	C_COLOR	*clr,
+	double	intensity
+)
 {
+	static char	spname[128];
+	double	mult;
+	int	cbeg, cend, i;
+
+	if (!dospectra | !(clr->flags & C_CDSPEC)) {
+		cvtcolor(radrgb, clr, intensity);
+		return("void");			/* just use RGB */
+	}
+	setcolor(radrgb, intensity, intensity, intensity);
+	for (cbeg = 0; cbeg < C_CNSS; cbeg++)	/* trim zeros off beginning */
+		if (clr->ssamp[cbeg])
+			break;
+	if (cbeg >= C_CNSS)			/* should never happen! */
+		return("void");
+	if (clr->client_data != NULL) {		/* get name if available */
+		strcpy(spname, (char *)clr->client_data);
+		strcat(spname, "*");		/* make sure it's special */
+		if (!newspecdef(clr))		/* output already? */
+			return(spname);
+	} else
+		strcpy(spname, "spec*");
+	c_ccvt(clr, C_CSEFF);			/* else output spectrum prim */
+	for (cend = 0; !clr->ssamp[C_CNSS-1-cend]; cend++)
+		;				/* trim zeros off end */
+	fprintf(matfp, "\nvoid spectrum %s\n0\n0\n", spname);
+	fprintf(matfp, "%d %d %d", C_CNSS+2-cbeg-cend,
+		C_CMINWL+cbeg*C_CWLI, C_CMAXWL-cend*C_CWLI);
+	mult = (C_CNSS*c_dfcolor.eff)/(clr->ssum*clr->eff);
+	for (i = cbeg; i < C_CNSS-cend; i++) {
+		if (!((i-cbeg+1)%6)) fputc('\n', matfp);
+		fprintf(matfp, "\t%.5f", clr->ssamp[i]*mult);
+	}
+	fputc('\n', matfp);
+	return(spname);
+}
+
+
+char *
+object(void)			/* return current object name */
+{
 	static char	objbuf[64];
-	register int	i;
-	register char	*cp;
+	int	i;
+	char	*cp;
 	int	len;
 						/* tracked by obj_handler */
 	i = obj_nnames - sizeof(objbuf)/16;
@@ -617,8 +749,10 @@ object()			/* return current object name */
 
 
 char *
-addarg(op, arg)				/* add argument and advance pointer */
-register char	*op, *arg;
+addarg(				/* add argument and advance pointer */
+	char *op,
+	char *arg
+)
 {
 	*op = ' ';
 	while ( (*++op = *arg++) )