--- ray/src/cv/bsdfrep.h	2012/10/23 21:09:29	2.4
+++ ray/src/cv/bsdfrep.h	2013/09/26 17:05:00	2.8
@@ -1,4 +1,4 @@
-/* RCSid $Id: bsdfrep.h,v 2.4 2012/10/23 21:09:29 greg Exp $ */
+/* RCSid $Id: bsdfrep.h,v 2.8 2013/09/26 17:05:00 greg Exp $ */
 /*
  * Definitions for BSDF representation used to interpolate measured data.
  *
@@ -57,6 +57,8 @@ typedef struct s_rbfnode {
 #define INP_QUAD3	4		/* 180-270 degree quadrant */
 #define INP_QUAD4	8		/* 270-360 degree quadrant */
 
+				/* active grid resolution */
+extern int		grid_res;
 				/* coverage/symmetry using INP_QUAD? flags */
 extern int		inp_coverage;
 
@@ -67,6 +69,18 @@ extern int		single_plane_incident;
 extern int		input_orient;
 extern int		output_orient;
 
+				/* log BSDF histogram */
+#define	HISTLEN		256
+#define BSDF2BIG	(1./M_PI)
+#define	BSDF2SML	1e-8
+#define	HISTLNR		17.2759509		/* log(BSDF2BIG/BSDF2SML) */
+extern int		bsdf_hist[HISTLEN];
+#define histndx(v)	(int)(log((v)*(1./BSDF2SML))*(HISTLEN/HISTLNR))
+#define histval(i)	(exp(((i)+.5)*(HISTLNR/HISTLEN))*BSDF2SML)
+
+				/* BSDF value for boundary regions */
+extern double		bsdf_min;
+
 				/* processed incident DSF measurements */
 extern RBFNODE		*dsf_list;
 
@@ -112,6 +126,9 @@ extern void		rev_symmetry(FVECT vec, int sym);
 /* Reverse symmetry for an RBF distribution */
 extern void		rev_rbf_symmetry(RBFNODE *rbf, int sym);
 
+/* Rotate RBF to correspond to given incident vector */
+extern void		rotate_rbf(RBFNODE *rbf, const FVECT invec);
+
 /* Compute volume associated with Gaussian lobe */
 extern double		rbf_volume(const RBFVAL *rbfp);
 
@@ -167,4 +184,4 @@ extern void		build_mesh(void);
 extern int		get_interp(MIGRATION *miga[3], FVECT invec);
 
 /* Partially advect between recorded incident angles and allocate new RBF */
-extern RBFNODE *	advect_rbf(const FVECT invec);
+extern RBFNODE *	advect_rbf(const FVECT invec, int lobe_lim);