--- ray/src/gen/mkillum.h	1991/07/23 16:13:27	1.2
+++ ray/src/gen/mkillum.h	2007/12/08 01:43:09	2.14
@@ -1,38 +1,115 @@
-/* Copyright (c) 1991 Regents of the University of California */
-
-/* SCCSid "$SunId$ LBL" */
-
+/* RCSid: $Id: mkillum.h,v 2.14 2007/12/08 01:43:09 greg Exp $ */
 /*
  * Common definitions for mkillum
  */
+#ifndef _RAD_MKILLUM_H_
+#define _RAD_MKILLUM_H_
 
-#include  "standard.h"
-
-#include  "object.h"
-
+#include  "ray.h"
 #include  "otypes.h"
+#include  "random.h"
 
+#ifdef __cplusplus
+extern "C" {
+#endif
 				/* illum flags */
-#define  IL_FLAT	0x1		/* flat surface */
-#define  IL_LIGHT	0x2		/* light rather than illum */
-#define  IL_COLDST	0x4		/* colored distribution */
-#define  IL_COLAVG	0x8		/* compute average color */
-#define  IL_DATCLB	0x10		/* OK to clobber data file */
+#define  IL_LIGHT	0x1		/* light rather than illum */
+#define  IL_COLDST	0x2		/* use color distribution */
+#define  IL_COLAVG	0x4		/* use average color */
+#define  IL_DATCLB	0x8		/* OK to clobber data file */
 
+				/* up directions */
+typedef enum {
+	UDzneg=-3,
+	UDyneg=-2,
+	UDxneg=-1,
+	UDunknown=0,
+	UDxpos=1,
+	UDypos=2,
+	UDzpos=3
+} UpDir;
+				/* BSDF coordinate calculation routines */
+				/* vectors always point away from surface */
+typedef int	b_vecf(FVECT, int, void *);
+typedef int	b_ndxf(FVECT, void *);
+typedef double	b_radf(int, void *);
+
+/* Bidirectional Scattering Distribution Function */
+struct BSDF_data {
+	int	ninc;			/* number of incoming directions */
+	int	nout;			/* number of outgoing directions */
+	void	*ib_priv;		/* input basis private data */
+	b_vecf	*ib_vec;		/* get input vector from index */
+	b_ndxf	*ib_ndx;		/* get input index from vector */
+	b_radf	*ib_ohm;		/* get input radius for index */
+	void	*ob_priv;		/* output basis private data */
+	b_vecf	*ob_vec;		/* get output vector from index */
+	b_ndxf	*ob_ndx;		/* get output index from vector */
+	b_radf	*ob_ohm;		/* get output radius for index */
+	float	*bsdf;			/* scattering distribution data */
+};				/* bidirectional scattering distrib. func. */
+
 struct illum_args {
 	int	flags;			/* flags from list above */
+	UpDir	udir;			/* up direction */
+	double	thick;			/* object thickness */
 	char	matname[MAXSTR];	/* illum material name */
 	char	datafile[MAXSTR];	/* distribution data file name */
 	int	dfnum;			/* data file number */
 	char	altmat[MAXSTR];		/* alternate material name */
+	int	sampdens;		/* point sample density */
 	int	nsamps;			/* # of samples in each direction */
-	int	nalt, nazi;		/* # of altitude and azimuth angles */
-	FVECT	orient;			/* coordinate system orientation */
+	struct BSDF_data
+		*sd;			/* scattering data (if set) */
+	float	minbrt;			/* minimum average brightness */
+	COLOR	col;			/* computed average color */
 };				/* illum options */
 
-struct rtproc {
-	int	pd[3];			/* rtrace pipe descriptors */
-	float	*buf;			/* rtrace i/o buffer */
-	int	bsiz;			/* maximum rays for rtrace buffer */
-	int	nrays;			/* current length of rtrace buffer */
-};				/* rtrace process */
+#define getBSDF_incvec(v,b,i)	(*(b)->ib_vec)(v,i,(b)->ib_priv)
+#define getBSDF_incndx(b,v)	(*(b)->ib_ndx)(v,(b)->ib_priv)
+#define getBSDF_incohm(b,i)	(*(b)->ib_ohm)(i,(b)->ib_priv)
+#define getBSDF_outvec(v,b,o)	(*(b)->ob_vec)(v,o,(b)->ob_priv)
+#define getBSDF_outndx(b,v)	(*(b)->ob_ndx)(v,(b)->ob_priv)
+#define getBSDF_outohm(b,o)	(*(b)->ob_ohm)(o,(b)->ob_priv)
+#define BSDF_visible(b,i,o)	(b)->bsdf[(o)*(b)->ninc + (i)]
+
+extern struct BSDF_data *load_BSDF(char *fname);
+extern void free_BSDF(struct BSDF_data *b);
+extern int r_BSDF_incvec(FVECT v, struct BSDF_data *b, int i,
+				double rv, MAT4 xm);
+extern int r_BSDF_outvec(FVECT v, struct BSDF_data *b, int o,
+				double rv, MAT4 xm);
+extern int getBSDF_xfm(MAT4 xm, FVECT nrm, UpDir ud);
+extern void redistribute(struct BSDF_data *b, int nalt, int nazi,
+				FVECT u, FVECT v, FVECT w, MAT4 xm);
+
+extern void printobj(char *mod, OBJREC *obj);
+extern int average(struct illum_args *il, COLORV *da, int n);
+extern void flatout(struct illum_args *il, COLORV *da, int n, int m,
+	FVECT u, FVECT v, FVECT w);
+extern void illumout(struct illum_args *il, OBJREC *ob);
+extern void roundout(struct illum_args *il, COLORV *da, int n, int m);
+
+extern void newdist(int siz);
+extern int process_ray(RAY *r, int rv);
+extern void raysamp(int ndx, FVECT org, FVECT dir);
+extern void rayclean(void);
+
+extern void flatdir(FVECT  dv, double  alt, double  azi);
+
+extern int my_default(OBJREC *, struct illum_args *, char *);
+extern int my_face(OBJREC *, struct illum_args *, char *);
+extern int my_sphere(OBJREC *, struct illum_args *, char *);
+extern int my_ring(OBJREC *, struct illum_args *, char *);
+
+extern COLORV *	distarr;		/* distribution array */
+extern int	distsiz;
+
+extern char	*progname;
+
+#ifdef __cplusplus
+}
+#endif
+
+#endif /* _RAD_MKILLUM_H_ */
+