--- ray/src/cv/mgf2rad.c	1996/07/24 13:07:48	2.23
+++ ray/src/cv/mgf2rad.c	2024/01/04 01:55:42	2.33
@@ -1,19 +1,20 @@
-/* Copyright (c) 1996 Regents of the University of California */
-
 #ifndef lint
-static char SCCSid[] = "$SunId$ LBL";
+static const char	RCSid[] = "$Id: mgf2rad.c,v 2.33 2024/01/04 01:55:42 greg Exp $";
 #endif
-
 /*
  * Convert MGF (Materials and Geometry Format) to Radiance
  */
 
 #include <stdio.h>
+#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,17 +24,35 @@ double	glowdist = FHUGE;		/* glow test distance */
 
 double  emult = 1.;			/* emitter multiplier */
 
-FILE	*matfp = stdout;		/* material output file */
+FILE	*matfp;				/* material output file */
 
-int	r_comment(), r_cone(), r_cyl(), r_face(), r_ies(), r_ring(), r_sph();
-char	*material(), *object(), *addarg();
 
+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);
 
-main(argc, argv)		/* convert files to stdout */
-int	argc;
-char	*argv[];
+
+int
+main(
+	int	argc,
+	char	*argv[]
+)
 {
 	int	i;
+
+	matfp = stdout;
 				/* print out parser version */
 	printf("## Translated from MGF Version %d.%d\n", MG_VMAJOR, MG_VMINOR);
 				/* initialize dispatch table */
@@ -41,9 +60,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 */
@@ -120,9 +139,10 @@ userr:
 
 
 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 +155,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;
@@ -163,7 +199,7 @@ char	**av;
 		if (r2 == 0.)
 			return(MG_EILL);
 		inv = r2 < 0.;
-	} else if (r2 != 0. && inv ^ r2 < 0.)
+	} else if (r2 != 0. && inv ^ (r2 < 0.))
 		return(MG_EILL);
 	if (inv) {
 		r1 = -r1;
@@ -183,9 +219,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 +258,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 +292,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;
@@ -276,7 +315,7 @@ char	**av;
 	xf_rotvect(norm, cv->n);		/* rotate normal */
 	r1 = xf_scale(atof(av[2]));		/* scale radii */
 	r2 = xf_scale(atof(av[3]));
-	if (r1 < 0. | r2 <= r1)
+	if ((r1 < 0.) | (r2 <= r1))
 		return(MG_EILL);
 	if ((mat = material()) == NULL)		/* get material */
 		return(MG_EBADMAT);
@@ -291,16 +330,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);
@@ -314,16 +355,24 @@ char	**av;
 			if (is0vect(cva[i-1]->n))
 				break;
 		}
-		if (i == ac) {
+		if (i < ac)
+			i = ISFLAT;
+		else
 			i = flat_tri(cva[0]->p, cva[1]->p, cva[2]->p,
 					cva[0]->n, cva[1]->n, cva[2]->n);
-			if (i < 0)
-				return(MG_OK);	/* degenerate (error?) */
+		if (i == DEGEN)
+			return(MG_OK);		/* degenerate (error?) */
+		if (i == RVBENT) {
+			myi = !myi;
+			i = ISBENT;
+		} else if (i == RVFLAT) {
+			myi = !myi;
+			i = ISFLAT;
 		}
-		if (!i) {			/* smoothed triangles */
-			do_tri(mat, cva[0], cva[1], cva[2]);
+		if (i == ISBENT) {		/* smoothed triangles */
+			do_tri(mat, cva[0], cva[1], cva[2], myi);
 			if (ac == 5)
-				do_tri(mat, cva[2], cva[3], cva[0]);
+				do_tri(mat, cva[2], cva[3], cva[0], myi);
 			return(MG_OK);
 		}
 	}
@@ -331,7 +380,7 @@ char	**av;
 	printf("\n%s polygon %sf%d\n", mat, object(), ++nfaces);
 	printf("0\n0\n%d\n", 3*(ac-1));
 	for (i = 1; i < ac; i++) {	/* get, transform, print each vertex */
-		if ((cv = c_getvert(av[invert ? ac-i : i])) == NULL)
+		if ((cv = c_getvert(av[myi ? ac-i : i])) == NULL)
 			return(MG_EUNDEF);
 		xf_xfmpoint(v, cv->p);
 		putv(v);
@@ -341,16 +390,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);
@@ -405,19 +455,24 @@ char	**av;
 }
 
 
-do_tri(mat, cv1, cv2, cv3)		/* put out smoothed triangle */
-char	*mat;
-C_VERTEX	*cv1, *cv2, *cv3;
+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;
-	FLOAT	bcoor[3][3];
+	RREAL	bcoor[3][3];
 	C_VERTEX	*cvt;
 	FVECT	v1, v2, v3;
 	FVECT	n1, n2, n3;
-	register int	i;
+	int	i;
 
-	if (invert) {			/* swap vertex order if inverted */
+	if (iv) {			/* swap vertex order if inverted */
 		cvt = cv1;
 		cv1 = cv3;
 		cv3 = cvt;
@@ -438,7 +493,7 @@ C_VERTEX	*cv1, *cv2, *cv3;
 		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");
@@ -448,13 +503,21 @@ C_VERTEX	*cv1, *cv2, *cv3;
 }
 
 
+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;
@@ -463,15 +526,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));
@@ -480,7 +543,7 @@ material()			/* get (and print) current material */
 	}
 	d = c_cmaterial->rd + c_cmaterial->td +
 			c_cmaterial->rs + c_cmaterial->ts;
-	if (d < 0. | d > 1.)
+	if ((d < 0.) | (d > 1.))
 		return(NULL);
 					/* check for glass/dielectric */
 	if (c_cmaterial->nr > 1.1 &&
@@ -508,73 +571,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;
@@ -584,12 +657,68 @@ 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	i;
+
+	if (!(clr->flags & C_CDSPEC)) {		/* not defined spectrally? */
+		cvtcolor(radrgb, clr, intensity);
+		return("void");
+	}
+	setcolor(radrgb, intensity, intensity, intensity);
+	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 */
+	fprintf(matfp, "\nvoid spectrum %s\n0\n0\n", spname);
+	fprintf(matfp, "%d %d %d", C_CNSS+2, C_CMINWL, C_CMAXWL);
+	mult = (C_CNSS*c_dfcolor.eff)/(clr->ssum*clr->eff);
+	for (i = 0; i < C_CNSS; i++) {
+		if (!((i+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;
@@ -607,11 +736,13 @@ 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++)
+	while ( (*++op = *arg++) )
 		;
 	return(op);
 }