/* RCSid: $Id: glare.h,v 2.7 2008/05/01 16:42:06 greg Exp $ */
/*
 * Common data structures for glare source finding routines
 */
#ifndef _RAD_GLARE_H_
#define _RAD_GLARE_H_

#include "standard.h"
#include "view.h"
#include "color.h"
#include "setscan.h"

#ifdef __cplusplus
extern "C" {
#endif

#define GLAREBR		7.0		/* glare source is this * avg. lum. */

#define SAMPDENS	75		/* default samples per unit in image */
#define TSAMPSTEP	10		/* sample step to compute threshold */

#define SEPS		1		/* sources this close ==> contig. */

#define SAMIN		.005		/* minimum solid angle for source */
#define MAXBUDDY	(4.*sqrt(SAMIN/PI))	/* max separation for pairing */

#define TOOSMALL(s)	((s)->brt*(s)->dom < threshold*SAMIN)

#define SABIG		.025		/* solid angle of splittable source */
#define LCORR		.12		/* linearity of splittable source */

extern VIEW	ourview;		/* our view */
extern VIEW	pictview;		/* picture view */
extern VIEW	leftview, rightview;	/* leftmost and rightmost views */

extern int	verbose;		/* verbose reporting */
extern char	*progname;		/* global argv[0] */

extern double	threshold;		/* threshold value for glare sources */

extern int	sampdens;		/* sample density */
extern ANGLE	glarang[];		/* glare calculation angles */
extern int	nglarangs;
extern double	maxtheta;		/* maximum glare angle (in radians) */
extern int	hsize;			/* horizontal size */

#define nglardirs	(2*nglarangs+1)
#define vsize		(sampdens-1)
#define hscale(v)	sqrt((double)(sampdens*sampdens - (v)*(v)))
#define hlim(v)		(int)(maxtheta*hscale(v))
#define h_theta(h,v)	(-(h)/hscale(v))

extern struct illum {
	float	theta;		/* glare direction */
	float	lcos, lsin;	/* cosine and sine to left view */
	float	rcos, rsin;	/* cosine and sine to right view */
	double	sum;		/* sum of indirect luminances */
	double	n;		/* number of values in sum */
} *indirect;		/* array of indirect illuminances */

struct srcspan {
	short	v;		/* vertical position */
	short	l, r;		/* left and right horizontal limits */
	float	brsum;		/* sum of brightnesses for this span */
	struct srcspan	*next;	/* next source span in list */
};

extern struct source {
	FVECT	dir;		/* source direction */
	double	dom;		/* solid angle of source */
	double	brt;		/* average source brightness */
	struct srcspan	*first;	/* first span for this source */
	struct source	*next;	/* next source in list */
} *donelist;			/* finished sources */


extern long	npixinvw;	/* number of samples in view */
extern long	npixmiss;	/* number of samples missing */

	/* defined in findglare.c */
extern void memerr(char	*s);
extern int compdir(FVECT vd, int x, int y);
extern double pixsize(int x, int y);
	/* defined in glaresrc.c */
extern void comp_thresh(void);
extern void analyze(void);
extern void absorb_specks(void);
	/* defined in glareval.c */
extern void open_pict(char *fn);
extern void fork_rtrace(char *av[]);
extern void close_pict(void);
extern void done_rtrace(void);
extern void getviewspan(int vv, float *vb);
extern double getviewpix(int vh, int vv);

#ifdef __cplusplus
}
#endif
#endif /* _RAD_GLARE_H_ */