--- ray/src/cv/bsdfrep.h	2012/10/20 07:02:00	2.2
+++ ray/src/cv/bsdfrep.h	2014/03/08 06:37:59	2.14
@@ -1,4 +1,4 @@
-/* RCSid $Id: bsdfrep.h,v 2.2 2012/10/20 07:02:00 greg Exp $ */
+/* RCSid $Id: bsdfrep.h,v 2.14 2014/03/08 06:37:59 greg Exp $ */
 /*
  * Definitions for BSDF representation used to interpolate measured data.
  *
@@ -10,7 +10,7 @@
 #define DEBUG		1
 
 #ifndef GRIDRES
-#define GRIDRES		200		/* grid resolution per side */
+#define GRIDRES		(1<<8)		/* grid resolution per side */
 #endif
 					/* convert to/from coded radians */
 #define ANG2R(r)	(int)((r)*((1<<16)/M_PI))
@@ -18,8 +18,7 @@
 
 typedef struct {
 	float		vsum;		/* DSF sum */
-	unsigned short	nval;		/* number of values in sum */
-	unsigned short	crad;		/* radius (coded angle) */
+	unsigned int	nval;		/* number of values in sum */
 } GRIDVAL;			/* grid value */
 
 typedef struct {
@@ -57,6 +56,12 @@ typedef struct s_rbfnode {
 #define INP_QUAD3	4		/* 180-270 degree quadrant */
 #define INP_QUAD4	8		/* 270-360 degree quadrant */
 
+				/* name and manufacturer if known */
+extern char		bsdf_name[];
+extern char		bsdf_manuf[];
+				/* active grid resolution */
+extern int		grid_res;
+				/* coverage/symmetry using INP_QUAD? flags */
 extern int		inp_coverage;
 
 				/* all incident angles in-plane so far? */
@@ -66,17 +71,26 @@ 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 unsigned long	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;
 
 				/* RBF-linking matrices (edges) */
 extern MIGRATION	*mig_list;
 
-				/* migration edges drawn in raster fashion */
-extern MIGRATION	*mig_grid[GRIDRES][GRIDRES];
-
-#define mtx_nrows(m)	((m)->rbfv[0]->nrbf)
-#define mtx_ncols(m)	((m)->rbfv[1]->nrbf)
+#define mtx_nrows(m)	(m)->rbfv[0]->nrbf
+#define mtx_ncols(m)	(m)->rbfv[1]->nrbf
 #define mtx_coef(m,i,j)	(m)->mtx[(i)*mtx_ncols(m) + (j)]
 #define is_src(rbf,m)	((rbf) == (m)->rbfv[0])
 #define is_dest(rbf,m)	((rbf) == (m)->rbfv[1])
@@ -91,9 +105,9 @@ extern MIGRATION	*mig_grid[GRIDRES][GRIDRES];
 extern char		*progname;
 
 				/* get theta value in degrees [0,180) range */
-#define get_theta180(v)	(180./M_PI)*acos((v)[2])
+#define get_theta180(v)	((180./M_PI)*Acos((v)[2]))
 				/* get phi value in degrees, [0,360) range */
-#define	get_phi360(v)	((180./M_PI)*atan2((v)[1],(v)[0]) + 180.)
+#define	get_phi360(v)	((180./M_PI)*atan2((v)[1],(v)[0]) + 360.*((v)[1]<0))
 
 				/* our loaded grid for this incident angle */
 extern double		theta_in_deg, phi_in_deg;
@@ -114,6 +128,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);
 
@@ -143,6 +160,9 @@ extern int		is_rev_tri(const FVECT v1,
 /* Find vertices completing triangles on either side of the given edge */
 extern int		get_triangles(RBFNODE *rbfv[2], const MIGRATION *mig);
 
+/* Clear our BSDF representation and free memory */
+extern void		clear_bsdf_rep(void);
+
 /* Write our BSDF mesh interpolant out to the given binary stream */
 extern void		save_bsdf_rep(FILE *ofp);
 
@@ -162,11 +182,12 @@ extern RBFNODE *	make_rbfrep(void);
 /* Build our triangle mesh from recorded RBFs */
 extern void		build_mesh(void);
 
-/* Draw edge list into mig_grid array */
-extern void		draw_edges(void);
-
 /* Find edge(s) for interpolating the given vector, applying symmetry */
 extern int		get_interp(MIGRATION *miga[3], FVECT invec);
 
+/* Advect and allocate new RBF along edge (internal call) */
+extern RBFNODE *	e_advect_rbf(const MIGRATION *mig,
+					const FVECT invec, int lobe_lim);
+
 /* 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);