--- ray/src/rt/p_data.c	1991/01/02 18:22:40	1.4
+++ ray/src/rt/p_data.c	2023/11/15 18:02:53	2.10
@@ -1,22 +1,20 @@
-/* Copyright (c) 1991 Regents of the University of California */
-
 #ifndef lint
-static char SCCSid[] = "$SunId$ LBL";
+static const char	RCSid[] = "$Id: p_data.c,v 2.10 2023/11/15 18:02:53 greg Exp $";
 #endif
-
 /*
  *  p_data.c - routine for stored patterns.
- *
- *     6/4/86
  */
 
-#include  "ray.h"
+#include "copyright.h"
 
+#include  "ray.h"
 #include  "data.h"
+#include  "func.h"
+#include  "rtotypes.h"
 
 /*
- *	A stored pattern can either be brightness or
- *  color data.  Brightness data is specified as:
+ *	A stored pattern can either be brightness,
+ *  color, or spectral data.  Brightness data is specified as:
  *
  *	modifier brightdata name
  *	4+ func dfname vfname v0 v1 .. xf
@@ -37,159 +35,226 @@ static char SCCSid[] = "$SunId$ LBL";
  *	0
  *	n A1 A2 ..
  *
+ *  A simple spectrum is specified as:
+ *
+ *	modifier spectrum name
+ *	0
+ *	0
+ *	5+ nmA nmB s1 s2 s3 ..
+ *
+ *  A constant spectrum from a data file is given as:
+ *
+ *	modifier specfile name
+ *	1 dfname
+ *	0
+ *	0
+ *
  *  Vfname is the name of the file where the variable definitions
  *  can be found.  The list of real arguments can be accessed by
  *  definitions in the file.  The dfnames are the data file
  *  names.  The dimensions of the data files and the number
  *  of variables must match.  The funcs take a single argument
  *  for brightdata, and three for colordata and colorpict to produce
- *  interpolated values from the file.  The xf is a transformation
+ *  interpolated values from the file.  The xf is a transform spec
  *  to get from the original coordinates to the current coordinates.
  */
 
 
-p_bdata(m, r)			/* interpolate brightness data */
-register OBJREC  *m;
-RAY  *r;
+int
+p_bdata(			/* interpolate brightness data */
+	OBJREC  *m,
+	RAY  *r
+)
 {
-	extern double  varvalue(), funvalue(), datavalue();
-	extern int  errno;
-	int  nv;
 	double  bval;
-	double  pt[MAXDIM];
+	double  pt[MAXDDIM];
 	DATARRAY  *dp;
-	register char  **sa;
+	MFUNC  *mf;
+	int  i;
 
-	setfunc(m, r);
-
-	sa = m->oargs.sarg;
-
 	if (m->oargs.nsargs < 4)
 		objerror(m, USER, "bad # arguments");
-	if (!vardefined(sa[3]))
-		loadfunc(sa[2]);
+	dp = getdata(m->oargs.sarg[1]);
+	i = (1 << dp->nd) - 1;
+	mf = getfunc(m, 2, i<<3, 0);
+	setfunc(m, r);
 	errno = 0;
-	for (nv = 0; nv+3 < m->oargs.nsargs &&
-			sa[nv+3][0] != '-'; nv++) {
-		if (nv >= MAXDIM)
-			goto dimerr;
-		pt[nv] = varvalue(sa[nv+3]);
+	for (i = dp->nd; i-- > 0; ) {
+		pt[i] = evalue(mf->ep[i]);
+		if ((errno == EDOM) | (errno == ERANGE))
+			goto computerr;
 	}
-	if (errno)
-		goto computerr;
-	dp = getdata(sa[1]);
-	if (dp->nd != nv)
-		goto dimerr;
 	bval = datavalue(dp, pt);
 	errno = 0;
-	bval = funvalue(sa[0], 1, &bval);
-	if (errno)
+	bval = funvalue(m->oargs.sarg[0], 1, &bval);
+	if ((errno == EDOM) | (errno == ERANGE))
 		goto computerr;
-	scalecolor(r->pcol, bval);
-	return;
-
-dimerr:
-	objerror(m, USER, "dimension error");
-
+	scalescolor(r->pcol, bval);
+	return(0);
 computerr:
 	objerror(m, WARNING, "compute error");
-	return;
+	return(0);
 }
 
 
-p_cdata(m, r)			/* interpolate color data */
-register OBJREC  *m;
-RAY  *r;
+int
+p_cdata(			/* interpolate color data */
+	OBJREC  *m,
+	RAY  *r
+)
 {
-	extern double  varvalue(), funvalue(), datavalue();
-	extern int  errno;
-	int  i, nv;
 	double  col[3];
 	COLOR  cval;
-	double  pt[MAXDIM];
+	double  pt[MAXDDIM];
+	int  nv;
 	DATARRAY  *dp;
-	register char  **sa;
+	MFUNC  *mf;
+	int  i;
 
-	setfunc(m, r);
-
-	sa = m->oargs.sarg;
-
 	if (m->oargs.nsargs < 8)
 		objerror(m, USER, "bad # arguments");
-	if (!vardefined(sa[7]))
-		loadfunc(sa[6]);
-	for (nv = 0; nv+7 < m->oargs.nsargs &&
-			sa[nv+7][0] != '-'; nv++) {
-		if (nv >= MAXDIM)
-			goto dimerr;
-		errno = 0;
-		pt[nv] = varvalue(sa[nv+7]);
-		if (errno)
+	dp = getdata(m->oargs.sarg[3]);
+	i = (1 << (nv = dp->nd)) - 1;
+	mf = getfunc(m, 6, i<<7, 0);
+	setfunc(m, r);
+	errno = 0;
+	for (i = 0; i < nv; i++) {
+		pt[i] = evalue(mf->ep[i]);
+		if ((errno == EDOM) | (errno == ERANGE))
 			goto computerr;
 	}
-	for (i = 0; i < 3; i++) {
-		dp = getdata(sa[i+3]);
+	col[0] = datavalue(dp, pt);
+	for (i = 1; i < 3; i++) {
+		dp = getdata(m->oargs.sarg[i+3]);
 		if (dp->nd != nv)
-			goto dimerr;
+			objerror(m, USER, "dimension error");
 		col[i] = datavalue(dp, pt);
 	}
 	errno = 0;
-	setcolor(cval,	funvalue(sa[0], 3, col),
-			funvalue(sa[1], 3, col),
-			funvalue(sa[2], 3, col));
-	if (errno)
+	for (i = 0; i < 3; i++)
+		if (fundefined(m->oargs.sarg[i]) < 3)
+			colval(cval,i) = funvalue(m->oargs.sarg[i], 1, col+i);
+		else
+			colval(cval,i) = funvalue(m->oargs.sarg[i], 3, col);
+	if ((errno == EDOM) | (errno == ERANGE))
 		goto computerr;
-	multcolor(r->pcol, cval);
-	return;
-
-dimerr:
-	objerror(m, USER, "dimension error");
-
+	smultcolor(r->pcol, cval);
+	return(0);
 computerr:
 	objerror(m, WARNING, "compute error");
-	return;
+	return(0);
 }
 
 
-p_pdata(m, r)			/* interpolate picture data */
-register OBJREC  *m;
-RAY  *r;
+int
+p_pdata(			/* interpolate picture data */
+	OBJREC  *m,
+	RAY  *r
+)
 {
-	extern double  varvalue(), funvalue(), datavalue();
-	extern int  errno;
-	int  i;
 	double  col[3];
 	COLOR  cval;
 	double  pt[2];
 	DATARRAY  *dp;
-	register char  **sa;
+	MFUNC  *mf;
+	int  i;
 
-	setfunc(m, r);
-
-	sa = m->oargs.sarg;
-
 	if (m->oargs.nsargs < 7)
 		objerror(m, USER, "bad # arguments");
-	if (!vardefined(sa[5]))
-		loadfunc(sa[4]);
+	mf = getfunc(m, 4, 0x3<<5, 0);
+	setfunc(m, r);
 	errno = 0;
-	pt[1] = varvalue(sa[5]);	/* y major ordering */
-	pt[0] = varvalue(sa[6]);
-	if (errno)
+	pt[1] = evalue(mf->ep[0]);	/* y major ordering */
+	pt[0] = evalue(mf->ep[1]);
+	if ((errno == EDOM) | (errno == ERANGE))
 		goto computerr;
-	dp = getpict(sa[3]);
+	dp = getpict(m->oargs.sarg[3]);
 	for (i = 0; i < 3; i++)
 		col[i] = datavalue(dp+i, pt);
 	errno = 0;
-	setcolor(cval,	funvalue(sa[0], 3, col),
-			funvalue(sa[1], 3, col),
-			funvalue(sa[2], 3, col));
-	if (errno)
+	for (i = 0; i < 3; i++)
+		if (fundefined(m->oargs.sarg[i]) < 3)
+			colval(cval,i) = funvalue(m->oargs.sarg[i], 1, col+i);
+		else
+			colval(cval,i) = funvalue(m->oargs.sarg[i], 3, col);
+	if ((errno == EDOM) | (errno == ERANGE))
 		goto computerr;
-	multcolor(r->pcol, cval);
-	return;
+	smultcolor(r->pcol, cval);
+	return(0);
 
 computerr:
 	objerror(m, WARNING, "compute error");
-	return;
+	return(0);
+}
+
+
+int
+p_spectrum(			/* simple constant spectrum */
+	OBJREC  *m,
+	RAY  *r
+)
+{
+	COLORV	*scval;
+
+	if ((scval = (COLORV *)m->os) == NULL) {
+		COLORV	*sinp;
+		double	hstep;
+		int	i;
+		if (m->oargs.nfargs < 5)
+			objerror(m, USER, "bad # arguments");
+		sinp = (COLORV *)malloc(sizeof(COLORV)*(m->oargs.nfargs-2));
+		scval = (COLORV *)malloc(sizeof(COLORV)*NCSAMP);
+		if ((sinp == NULL) | (scval == NULL))
+			objerror(m, SYSTEM, "out of memory");
+		for (i = m->oargs.nfargs-2; i--; )
+			sinp[i] = (COLORV)m->oargs.farg[i+2];
+		hstep = 0.5 * (m->oargs.farg[1] - m->oargs.farg[0]) /
+				(m->oargs.nfargs-3.0);
+		convertscolor(scval, NCSAMP, WLPART[0], WLPART[3],
+				sinp, m->oargs.nfargs-2,
+				m->oargs.farg[0]-hstep, m->oargs.farg[1]+hstep);
+		free(sinp);
+		m->os = (void *)scval;
+	}
+	smultscolor(r->pcol, scval);
+	return(0);
+}
+
+
+int
+p_specfile(			/* constant spectrum from 1-D data file */
+	OBJREC  *m,
+	RAY  *r
+)
+{
+	COLORV	*scval;
+
+	if ((scval = (COLORV *)m->os) == NULL) {
+		DATARRAY	*dp;
+		COLORV		*sinp;
+		double		step;
+		int		i;
+		if (m->oargs.nsargs != 1)
+			objerror(m, USER, "bad # arguments");
+		dp = getdata(m->oargs.sarg[0]);
+		if (dp->nd != 1)
+			objerror(m, USER, "data file must be 1-dimensional");
+
+		sinp = (COLORV *)malloc(sizeof(COLORV)*dp->dim[0].ne);
+		scval = (COLORV *)malloc(sizeof(COLORV)*NCSAMP);
+		if ((sinp == NULL) | (scval == NULL))
+			objerror(m, SYSTEM, "out of memory");
+		step = dp->dim[0].siz / (dp->dim[0].ne - 1.0);
+		for (i = dp->dim[0].ne; i-- > 0; ) {
+			double	wl = dp->dim[0].org + i*step;
+			sinp[i] = (COLORV)datavalue(dp, &wl);
+		}
+		convertscolor(scval, NCSAMP, WLPART[0], WLPART[3],
+				sinp, dp->dim[0].ne, dp->dim[0].org-.5*step,
+				dp->dim[0].org+dp->dim[0].siz+.5*step);
+		free(sinp);
+		m->os = (void *)scval;
+	}
+	smultscolor(r->pcol, scval);
+	return(0);
 }