src/hd/sm.h

/* RCSid: $Id: sm.h,v 3.12 2003/06/20 00:25:49 greg Exp $ */
/*
 *  sm.h
 */

#ifndef _RAD_SM_H_
#define _RAD_SM_H_

#include "rhd_sample.h"
#include "sm_qtree.h"
#include "sm_stree.h"

#ifdef __cplusplus
extern "C" {
#endif

#define NEWSETS


#ifndef TRUE
#define TRUE 1
#define FALSE 0
#endif


#define ON_V 1
#define ON_P 2
#define ON_E 3
#define IN_T 4

#ifdef SMLFLT
#define VERT_EPS 2e-3 /* min edge length in radians */
#define COS_VERT_EPS 0.999998 /* cos min edge length in radians */
#define EDGE_EPS 2e-5 /* min distance until considered "on-edge"*/
#define COLINEAR_EPS 1e-10 /* min sine of between edges angle :amount off PI */
#else
#define VERT_EPS 5e-4 /* min edge length in radians */
#define COS_VERT_EPS 0.999999875 /* cos min edge length in radians */
#define EDGE_EPS 2e-7 /* min distance until considered "on-edge"*/
#define COLINEAR_EPS 1e-10 /* min sine of between edges angle :amount off PI*/
#endif

#define EV_EPS EDGE_EPS      /* Minimum edge-vertex distance */


#define S_REPLACE_SCALE (5.*5.)   /* if (distance to new point squared) is
         > (triangle edge length squared*S_REPLACE_SCALE):for all edges/
         triangle vertices: new point is puncture point: dont add*/
#define S_REPLACE_TRI 2e-8              /* .052 radians to the sixth power */


#define SM_DEFAULT 0
#define SM_BASE_POINTS 162
#define SM_BASE_TRIS 320
#define SM_BASE_VERTS  SM_BASE_POINTS

#define SM_VIEW_FRAC   0.1

#define  SM_ALL_LEVELS -1


typedef int VERT;  /* One triangle that vertex belongs to- the rest
                      are derived by traversing neighbors */

typedef struct _EDGE {
    S_ID verts[2];
    int tris[2];
} EDGE;

#define E_NTH_VERT(e,i)  ((e>0)?Edges[(e)].verts[(i)]:Edges[-(e)].verts[(1-i)])
#define SET_E_NTH_VERT(e,i,v) if(e>0) Edges[(e)].verts[(i)]=(v); \
                               else Edges[-(e)].verts[(1-i)]=(v)
#define E_NTH_TRI(e,i)  ((e>0)?Edges[(e)].tris[(i)]:Edges[-(e)].tris[(1-i)])
#define SET_E_NTH_TRI(e,i,v) if(e>0) Edges[(e)].tris[(i)]=(v); \
                              else Edges[-(e)].tris[(1-i)]=(v)
#define eClear_edges() (Ecnt = 0,free(Edges),Edges=NULL)

#define FOR_ALL_EDGES(i) for((i)=1; (i) <= Ecnt; i++)
#define FOR_ALL_EDGES_FROM(e,i) for((i)=e+1; (i) <= Ecnt; i++)


typedef struct _TRI {
  S_ID verts[3];         /* Ids into sample and vertex array for each vertex*/
  int nbrs[3]; /* Ids for neighboring triangles: -1 if invalid */
}TRI;


#define T_NTH_NBR(t,i) ((t)->nbrs[(i)])
#define T_CLEAR_NBRS(t) (T_NTH_NBR(t,0)=T_NTH_NBR(t,1)=T_NTH_NBR(t,2)=-1)
#define T_NTH_NBR_PTR(t,n) \
         (T_NTH_NBR(n,0)==(t)?0:T_NTH_NBR(n,1)==(t)?1:T_NTH_NBR(n,2)==(t)?2:-1)
#define T_NTH_V(t,i)  ((t)->verts[(i)])
#define T_WHICH_V(t,i)     \
         (T_NTH_V(t,0)==(i)?0:T_NTH_V(t,1)==(i)?1:T_NTH_V(t,2)==(i)?2:-1)
#define T_NEXT_FREE(t) ((t)->nbrs[0])
#define T_VALID_FLAG(t) ((t)->nbrs[1])
#define T_IS_VALID(t)  (T_VALID_FLAG(t)!=-1)
#define T_FLAGS 4

typedef struct _SM {
    FVECT view_center;    /* Canonical view center defining unit sphere */
    SAMP *samples;        /* Sample point information */
    STREE locator;        /* spherical quadtree for point/triangle location */
    int max_tris;         /* Maximum number of triangles */
    int num_tri;          /* Current number of triangles */
    int free_tris;        /* pointer to free_list */
    int max_verts;        /* Maximum number of vertices in the mesh */
    TRI *tris;            /* Pointer to list of triangle structs */
    int32 *flags[T_FLAGS]; /* Bit 0 set if active(in current frustum) */
                          /* Bit 1 set if not rendered since created */
                          /* Bit 2 set if base triangle */
}SM;

#define T_ACTIVE_FLAG   0
#define T_NEW_FLAG      1
#define T_BASE_FLAG     2
#define T_BG_FLAG       3

#define SM_VIEW_CENTER(m)             ((m)->view_center)
#define SM_SAMP(m)                    ((m)->samples)
#define SM_LOCATOR(m)                 (&((m)->locator))
#define SM_MAX_TRIS(m)                ((m)->max_tris)
#define SM_NUM_TRI(m)                 ((m)->num_tri)
#define SM_FREE_TRIS(m)               ((m)->free_tris)
#define SM_MAX_VERTS(m)               ((m)->max_verts)
#define SM_TRIS(m)                    ((m)->tris)
#define SM_NTH_FLAGS(m,n)             ((m)->flags[(n)])
#define SM_FLAGS(m)                   ((m)->flags)


#define SM_IS_NTH_T_FLAG(sm,n,f)  IS_FLAG(SM_NTH_FLAGS(sm,f),n)
#define SM_SET_NTH_T_FLAG(sm,n,f) SET_FLAG(SM_NTH_FLAGS(sm,f),n)
#define SM_CLR_NTH_T_FLAG(sm,n,f) CLR_FLAG(SM_NTH_FLAGS(sm,f),n)

#define SM_IS_NTH_T_ACTIVE(sm,n)     SM_IS_NTH_T_FLAG(sm,n,T_ACTIVE_FLAG)
#define SM_IS_NTH_T_BASE(sm,n)       SM_IS_NTH_T_FLAG(sm,n,T_BASE_FLAG)
#define SM_IS_NTH_T_NEW(sm,n)        SM_IS_NTH_T_FLAG(sm,n,T_NEW_FLAG)
#define SM_IS_NTH_T_BG(sm,n)        SM_IS_NTH_T_FLAG(sm,n,T_BG_FLAG)

#define SM_SET_NTH_T_ACTIVE(sm,n)    SM_SET_NTH_T_FLAG(sm,n,T_ACTIVE_FLAG)
#define SM_SET_NTH_T_BASE(sm,n)      SM_SET_NTH_T_FLAG(sm,n,T_BASE_FLAG)
#define SM_SET_NTH_T_NEW(sm,n)       SM_SET_NTH_T_FLAG(sm,n,T_NEW_FLAG)
#define SM_SET_NTH_T_BG(sm,n)        SM_SET_NTH_T_FLAG(sm,n,T_BG_FLAG)

#define SM_CLR_NTH_T_ACTIVE(sm,n)   SM_CLR_NTH_T_FLAG(sm,n,T_ACTIVE_FLAG)
#define SM_CLR_NTH_T_BASE(sm,n)     SM_CLR_NTH_T_FLAG(sm,n,T_BASE_FLAG)
#define SM_CLR_NTH_T_NEW(sm,n)      SM_CLR_NTH_T_FLAG(sm,n,T_NEW_FLAG)
#define SM_CLR_NTH_T_BG(sm,n)       SM_CLR_NTH_T_FLAG(sm,n,T_BG_FLAG)

#define SM_NTH_TRI(m,n)               (&(SM_TRIS(m)[(n)]))
#define SM_NTH_VERT(m,n)              S_NTH_INFO(SM_SAMP(m),n)

#define SM_T_ID_VALID(s,t_id) T_IS_VALID(SM_NTH_TRI(s,t_id))
 

#define SM_MAX_SAMP(m)                S_MAX_SAMP(SM_SAMP(m))
#define SM_MAX_POINTS(m)              S_MAX_POINTS(SM_SAMP(m))
#define SM_SAMP_BASE(m)               S_BASE(SM_SAMP(m))
#define SM_NTH_WV(m,i)                S_NTH_W_PT(SM_SAMP(m),i)
#define SM_NTH_W_DIR(m,i)             S_NTH_W_DIR(SM_SAMP(m),i)
#define SM_DIR_ID(m,i)            (!SM_BASE_ID(m,i) && SM_NTH_W_DIR(m,i)==-1)
#define SM_NTH_RGB(m,i)               S_NTH_RGB(SM_SAMP(m),i)
#define SM_RGB(m)                    S_RGB(SM_SAMP(m))
#define SM_WP(m)                    S_W_PT(SM_SAMP(m))
#define SM_BRT(m)                     S_BRT(SM_SAMP(m))
#define SM_NTH_BRT(m,i)               S_NTH_BRT(SM_SAMP(m),i)
#define SM_CHR(m)                     S_CHR(SM_SAMP(m))
#define SM_NTH_CHR(m,i)               S_NTH_CHR(SM_SAMP(m),i)
#define SM_NUM_SAMP(m)               S_NUM_SAMP(SM_SAMP(m))
#define SM_TONE_MAP(m)               S_TONE_MAP(SM_SAMP(m))

#define SM_ALLOWED_VIEW_CHANGE(m) (SM_NUM_SAMP(m)/smDist_sum*SM_VIEW_FRAC)

#define SM_FOR_ALL_ADJACENT_TRIS(sm,id,t) for(t=smTri_next_ccw_nbr(sm,t,id); \
t!=SM_NTH_TRI(sm,SM_NTH_VERT(sm,id)); t=smTri_next_ccw_nbr(sm,t,id))

#define SM_INVALID_SAMP_ID(sm,id)  (((id) < 0) || ((id) >= SM_MAX_SAMP(sm)))
#define SM_INVALID_POINT_ID(sm,id)  (((id) < 0) || ((id) >= SM_MAX_POINTS(sm)))
#define SM_T_NTH_WV(sm,t,i)  (SM_NTH_WV(sm,T_NTH_V(t,i)))

#define SM_BASE_ID(s,i)  \
 ((i) >= S_MAX_SAMP(SM_SAMP(s)) && (i) < S_MAX_BASE_PT(SM_SAMP(s)))

#define SM_BG_SAMPLE(sm,i)   (SM_NTH_W_DIR(sm,i)==-1)

#define SM_BG_TRI(sm,i)  (SM_BG_SAMPLE(sm,T_NTH_V(SM_NTH_TRI(sm,i),0)) && \
                          SM_BG_SAMPLE(sm,T_NTH_V(SM_NTH_TRI(sm,i),1)) && \
                          SM_BG_SAMPLE(sm,T_NTH_V(SM_NTH_TRI(sm,i),2)))
#define SM_MIXED_TRI(sm,i)  (SM_BG_SAMPLE(sm,T_NTH_V(SM_NTH_TRI(sm,i),0)) || \
                          SM_BG_SAMPLE(sm,T_NTH_V(SM_NTH_TRI(sm,i),1)) || \
                          SM_BG_SAMPLE(sm,T_NTH_V(SM_NTH_TRI(sm,i),2)))

#define SM_FOR_ALL_SAMPLES(sm,i) for((i)=0;i < SM_NUM_SAMP(sm);(i)++)

#define SM_FOR_ALL_VALID_TRIS(m,i) for((i)=0,(i)=smNext_valid_tri(m,i);(i)< \
SM_NUM_TRI(m); (i)++,(i)= smNext_valid_tri(m,i))

#define SM_FOR_ALL_FLAGGED_TRIS(m,i,w,b) for(i=0,i=smNext_tri_flag_set(m,i,w,b);i < SM_NUM_TRI(m);i++,i=smNext_tri_flag_set(m,i,w,b))
#define SM_FOR_ALL_ACTIVE_TRIS(m,i) SM_FOR_ALL_FLAGGED_TRIS(m,i,T_ACTIVE_FLAG,0)
#define smInit_locator(sm)    stInit(SM_LOCATOR(sm))
#define smAlloc_locator(sm)   stAlloc(SM_LOCATOR(sm))
#define smFree_locator(sm)   stFree(SM_LOCATOR(sm))
#define smUnalloc_samp(sm,id)  sUnalloc_samp(SM_SAMP(sm),id)
#define smPoint_locate_cell(sm,pt)  stPoint_locate(SM_LOCATOR(sm),pt)
#define smFree_samples(sm)     sFree(SM_SAMP(sm))
#define smInit_samples(sm)     sInit(SM_SAMP(sm))

#define SM_S_NTH_QT(sm,s_id)       S_NTH_INFO1(SM_SAMP(sm),s_id)
#define smClear_vert(sm,id)    (SM_NTH_VERT(sm,id) = INVALID)

#define freebuf(b)  tempbuf(-1)

typedef struct _RT_ARGS_{
  FVECT orig,dir;
  int t_id;
  S_ID *os;
}RT_ARGS;

typedef struct _S_ARGS_{
  S_ID s_id;
  S_ID n_id;
}S_ARGS;


typedef struct _ADD_ARGS {
  int t_id;
  S_ID *del_set;
}ADD_ARGS;

extern SM *smMesh;
extern double smDist_sum;


#ifdef TEST_DRIVER
extern VIEW View;
extern VIEW Current_View;
extern int Pick_tri,Picking,Pick_samp;
extern FVECT Pick_point[500],Pick_origin,Pick_dir;
extern FVECT Pick_v0[500],Pick_v1[500],Pick_v2[500];
extern int Pick_q[500];
extern FVECT P0,P1,P2;
extern int Pick_cnt;
extern FVECT FrustumNear[4],FrustumFar[4];
#endif


/*
 * int
 * smInit(n)            : Initialize/clear data structures for n entries
 * int  n;
 *
 * Initialize sampL and other data structures for at least n samples.
 * If n is 0, then free data structures.  Return number actually allocated.
 *
 *
 * int
 * smNewSamp(c, p, v)   : register new sample point and return index
 * COLR c;              : pixel color (RGBE)
 * FVECT        p;      : world intersection point
 * FVECT        v;      : ray direction vector
 *
 * Add new sample point to data structures, removing old values as necessary.
 * New sample representation will be output in next call to smUpdate().
 *
 *
 * int
 * smFindSamp(orig, dir): intersect ray with 3D rep. and find closest sample
 * FVECT        orig, dir;
 *
 * Find the closest sample to the given ray.  Return -1 on failure.
 *
 */

#ifdef __cplusplus
}
#endif
#endif /* define _RAD_SM_H_ */

Revision:	3.13
Committed:	Mon Jul 14 22:24:00 2003 UTC (20 years, 9 months ago) by schorsch
Content type:	text/plain
Branch:	MAIN
CVS Tags:	rad3R6P1, rad3R6
Changes since 3.12:	+13 -10 lines
Log Message:	Instrumented headers against multiple inclusion and for use from C++. Moved includes in headers out of "C" scope.
#	User	Rev	Content
1	schorsch	3.13	/* RCSid: $Id: sm.h,v 3.12 2003/06/20 00:25:49 greg Exp $ */
2	gwlarson	3.1	/*
3			* sm.h
4			*/
5
6	schorsch	3.13	#ifndef _RAD_SM_H_
7			#define _RAD_SM_H_
8
9	gwlarson	3.1	#include "rhd_sample.h"
10	gwlarson	3.10	#include "sm_qtree.h"
11			#include "sm_stree.h"
12	schorsch	3.13
13			#ifdef __cplusplus
14			extern "C" {
15			#endif
16
17	gwlarson	3.6	#define NEWSETS
18	gwlarson	3.1
19	gwlarson	3.2
20	gwlarson	3.1	#ifndef TRUE
21			#define TRUE 1
22			#define FALSE 0
23			#endif
24
25	gwlarson	3.6
26			#define ON_V 1
27			#define ON_P 2
28			#define ON_E 3
29			#define IN_T 4
30
31	greg	3.11	#ifdef SMLFLT
32			#define VERT_EPS 2e-3 /* min edge length in radians */
33			#define COS_VERT_EPS 0.999998 /* cos min edge length in radians */
34			#define EDGE_EPS 2e-5 /* min distance until considered "on-edge"*/
35			#define COLINEAR_EPS 1e-10 /* min sine of between edges angle :amount off PI */
36			#else
37	gwlarson	3.10	#define VERT_EPS 5e-4 /* min edge length in radians */
38			#define COS_VERT_EPS 0.999999875 /* cos min edge length in radians */
39			#define EDGE_EPS 2e-7 /* min distance until considered "on-edge"*/
40	greg	3.11	#define COLINEAR_EPS 1e-10 /* min sine of between edges angle :amount off PI*/
41	gwlarson	3.10	#endif
42
43			#define EV_EPS EDGE_EPS /* Minimum edge-vertex distance */
44
45
46	gwlarson	3.5	#define S_REPLACE_SCALE (5.5.) / if (distance to new point squared) is
47	gwlarson	3.10	> (triangle edge length squared*S_REPLACE_SCALE):for all edges/
48			triangle vertices: new point is puncture point: dont add*/
49			#define S_REPLACE_TRI 2e-8 /* .052 radians to the sixth power */
50	gwlarson	3.4
51	gwlarson	3.1
52	gwlarson	3.10
53	gwlarson	3.4	#define SM_DEFAULT 0
54	gwlarson	3.10	#define SM_BASE_POINTS 162
55	gwlarson	3.6	#define SM_BASE_TRIS 320
56	gwlarson	3.10	#define SM_BASE_VERTS SM_BASE_POINTS
57	gwlarson	3.1
58			#define SM_VIEW_FRAC 0.1
59
60	gwlarson	3.9	#define SM_ALL_LEVELS -1
61
62	gwlarson	3.1
63			typedef int VERT; /* One triangle that vertex belongs to- the rest
64			are derived by traversing neighbors */
65
66			typedef struct _EDGE {
67	greg	3.11	S_ID verts[2];
68	gwlarson	3.1	int tris[2];
69			} EDGE;
70	gwlarson	3.4
71	gwlarson	3.1	#define E_NTH_VERT(e,i) ((e>0)?Edges[(e)].verts[(i)]:Edges[-(e)].verts[(1-i)])
72			#define SET_E_NTH_VERT(e,i,v) if(e>0) Edges[(e)].verts[(i)]=(v); \
73			else Edges[-(e)].verts[(1-i)]=(v)
74			#define E_NTH_TRI(e,i) ((e>0)?Edges[(e)].tris[(i)]:Edges[-(e)].tris[(1-i)])
75			#define SET_E_NTH_TRI(e,i,v) if(e>0) Edges[(e)].tris[(i)]=(v); \
76			else Edges[-(e)].tris[(1-i)]=(v)
77	gwlarson	3.10	#define eClear_edges() (Ecnt = 0,free(Edges),Edges=NULL)
78	gwlarson	3.1
79			#define FOR_ALL_EDGES(i) for((i)=1; (i) <= Ecnt; i++)
80	gwlarson	3.10	#define FOR_ALL_EDGES_FROM(e,i) for((i)=e+1; (i) <= Ecnt; i++)
81	gwlarson	3.1
82
83			typedef struct _TRI {
84	greg	3.11	S_ID verts[3]; /* Ids into sample and vertex array for each vertex*/
85	gwlarson	3.1	int nbrs[3]; /* Ids for neighboring triangles: -1 if invalid */
86			}TRI;
87
88
89			#define T_NTH_NBR(t,i) ((t)->nbrs[(i)])
90			#define T_CLEAR_NBRS(t) (T_NTH_NBR(t,0)=T_NTH_NBR(t,1)=T_NTH_NBR(t,2)=-1)
91			#define T_NTH_NBR_PTR(t,n) \
92			(T_NTH_NBR(n,0)==(t)?0:T_NTH_NBR(n,1)==(t)?1:T_NTH_NBR(n,2)==(t)?2:-1)
93			#define T_NTH_V(t,i) ((t)->verts[(i)])
94			#define T_WHICH_V(t,i) \
95			(T_NTH_V(t,0)==(i)?0:T_NTH_V(t,1)==(i)?1:T_NTH_V(t,2)==(i)?2:-1)
96	gwlarson	3.3	#define T_NEXT_FREE(t) ((t)->nbrs[0])
97			#define T_VALID_FLAG(t) ((t)->nbrs[1])
98	gwlarson	3.1	#define T_IS_VALID(t) (T_VALID_FLAG(t)!=-1)
99	gwlarson	3.8	#define T_FLAGS 4
100	gwlarson	3.1
101			typedef struct _SM {
102			FVECT view_center; /* Canonical view center defining unit sphere */
103	gwlarson	3.4	SAMP samples; / Sample point information */
104			STREE locator; /* spherical quadtree for point/triangle location */
105	gwlarson	3.1	int max_tris; /* Maximum number of triangles */
106	gwlarson	3.4	int num_tri; /* Current number of triangles */
107	gwlarson	3.1	int free_tris; /* pointer to free_list */
108			int max_verts; /* Maximum number of vertices in the mesh */
109			TRI tris; / Pointer to list of triangle structs */
110	greg	3.12	int32 flags[T_FLAGS]; / Bit 0 set if active(in current frustum) */
111	gwlarson	3.4	/* Bit 1 set if not rendered since created */
112			/* Bit 2 set if base triangle */
113	gwlarson	3.1	}SM;
114
115			#define T_ACTIVE_FLAG 0
116			#define T_NEW_FLAG 1
117			#define T_BASE_FLAG 2
118	gwlarson	3.8	#define T_BG_FLAG 3
119	gwlarson	3.1
120			#define SM_VIEW_CENTER(m) ((m)->view_center)
121			#define SM_SAMP(m) ((m)->samples)
122			#define SM_LOCATOR(m) (&((m)->locator))
123			#define SM_MAX_TRIS(m) ((m)->max_tris)
124	gwlarson	3.4	#define SM_NUM_TRI(m) ((m)->num_tri)
125	gwlarson	3.1	#define SM_FREE_TRIS(m) ((m)->free_tris)
126			#define SM_MAX_VERTS(m) ((m)->max_verts)
127			#define SM_TRIS(m) ((m)->tris)
128			#define SM_NTH_FLAGS(m,n) ((m)->flags[(n)])
129	gwlarson	3.4	#define SM_FLAGS(m) ((m)->flags)
130	gwlarson	3.1
131
132	gwlarson	3.4	#define SM_IS_NTH_T_FLAG(sm,n,f) IS_FLAG(SM_NTH_FLAGS(sm,f),n)
133			#define SM_SET_NTH_T_FLAG(sm,n,f) SET_FLAG(SM_NTH_FLAGS(sm,f),n)
134			#define SM_CLR_NTH_T_FLAG(sm,n,f) CLR_FLAG(SM_NTH_FLAGS(sm,f),n)
135
136	gwlarson	3.1	#define SM_IS_NTH_T_ACTIVE(sm,n) SM_IS_NTH_T_FLAG(sm,n,T_ACTIVE_FLAG)
137			#define SM_IS_NTH_T_BASE(sm,n) SM_IS_NTH_T_FLAG(sm,n,T_BASE_FLAG)
138			#define SM_IS_NTH_T_NEW(sm,n) SM_IS_NTH_T_FLAG(sm,n,T_NEW_FLAG)
139	gwlarson	3.8	#define SM_IS_NTH_T_BG(sm,n) SM_IS_NTH_T_FLAG(sm,n,T_BG_FLAG)
140	gwlarson	3.1
141			#define SM_SET_NTH_T_ACTIVE(sm,n) SM_SET_NTH_T_FLAG(sm,n,T_ACTIVE_FLAG)
142			#define SM_SET_NTH_T_BASE(sm,n) SM_SET_NTH_T_FLAG(sm,n,T_BASE_FLAG)
143	gwlarson	3.4	#define SM_SET_NTH_T_NEW(sm,n) SM_SET_NTH_T_FLAG(sm,n,T_NEW_FLAG)
144	gwlarson	3.8	#define SM_SET_NTH_T_BG(sm,n) SM_SET_NTH_T_FLAG(sm,n,T_BG_FLAG)
145	gwlarson	3.1
146	gwlarson	3.4	#define SM_CLR_NTH_T_ACTIVE(sm,n) SM_CLR_NTH_T_FLAG(sm,n,T_ACTIVE_FLAG)
147			#define SM_CLR_NTH_T_BASE(sm,n) SM_CLR_NTH_T_FLAG(sm,n,T_BASE_FLAG)
148			#define SM_CLR_NTH_T_NEW(sm,n) SM_CLR_NTH_T_FLAG(sm,n,T_NEW_FLAG)
149	gwlarson	3.8	#define SM_CLR_NTH_T_BG(sm,n) SM_CLR_NTH_T_FLAG(sm,n,T_BG_FLAG)
150	gwlarson	3.1
151			#define SM_NTH_TRI(m,n) (&(SM_TRIS(m)[(n)]))
152	gwlarson	3.10	#define SM_NTH_VERT(m,n) S_NTH_INFO(SM_SAMP(m),n)
153	gwlarson	3.6
154			#define SM_T_ID_VALID(s,t_id) T_IS_VALID(SM_NTH_TRI(s,t_id))
155	gwlarson	3.1
156	gwlarson	3.7
157	gwlarson	3.4	#define SM_MAX_SAMP(m) S_MAX_SAMP(SM_SAMP(m))
158			#define SM_MAX_POINTS(m) S_MAX_POINTS(SM_SAMP(m))
159			#define SM_SAMP_BASE(m) S_BASE(SM_SAMP(m))
160			#define SM_NTH_WV(m,i) S_NTH_W_PT(SM_SAMP(m),i)
161			#define SM_NTH_W_DIR(m,i) S_NTH_W_DIR(SM_SAMP(m),i)
162	gwlarson	3.7	#define SM_DIR_ID(m,i) (!SM_BASE_ID(m,i) && SM_NTH_W_DIR(m,i)==-1)
163	gwlarson	3.4	#define SM_NTH_RGB(m,i) S_NTH_RGB(SM_SAMP(m),i)
164			#define SM_RGB(m) S_RGB(SM_SAMP(m))
165	gwlarson	3.7	#define SM_WP(m) S_W_PT(SM_SAMP(m))
166	gwlarson	3.4	#define SM_BRT(m) S_BRT(SM_SAMP(m))
167			#define SM_NTH_BRT(m,i) S_NTH_BRT(SM_SAMP(m),i)
168			#define SM_CHR(m) S_CHR(SM_SAMP(m))
169			#define SM_NTH_CHR(m,i) S_NTH_CHR(SM_SAMP(m),i)
170			#define SM_NUM_SAMP(m) S_NUM_SAMP(SM_SAMP(m))
171			#define SM_TONE_MAP(m) S_TONE_MAP(SM_SAMP(m))
172	gwlarson	3.1
173			#define SM_ALLOWED_VIEW_CHANGE(m) (SM_NUM_SAMP(m)/smDist_sum*SM_VIEW_FRAC)
174
175			#define SM_FOR_ALL_ADJACENT_TRIS(sm,id,t) for(t=smTri_next_ccw_nbr(sm,t,id); \
176			t!=SM_NTH_TRI(sm,SM_NTH_VERT(sm,id)); t=smTri_next_ccw_nbr(sm,t,id))
177
178			#define SM_INVALID_SAMP_ID(sm,id) (((id) < 0) \|\| ((id) >= SM_MAX_SAMP(sm)))
179	gwlarson	3.4	#define SM_INVALID_POINT_ID(sm,id) (((id) < 0) \|\| ((id) >= SM_MAX_POINTS(sm)))
180	gwlarson	3.1	#define SM_T_NTH_WV(sm,t,i) (SM_NTH_WV(sm,T_NTH_V(t,i)))
181
182			#define SM_BASE_ID(s,i) \
183	gwlarson	3.4	((i) >= S_MAX_SAMP(SM_SAMP(s)) && (i) < S_MAX_BASE_PT(SM_SAMP(s)))
184	gwlarson	3.1
185			#define SM_BG_SAMPLE(sm,i) (SM_NTH_W_DIR(sm,i)==-1)
186
187	gwlarson	3.2	#define SM_BG_TRI(sm,i) (SM_BG_SAMPLE(sm,T_NTH_V(SM_NTH_TRI(sm,i),0)) && \
188			SM_BG_SAMPLE(sm,T_NTH_V(SM_NTH_TRI(sm,i),1)) && \
189			SM_BG_SAMPLE(sm,T_NTH_V(SM_NTH_TRI(sm,i),2)))
190			#define SM_MIXED_TRI(sm,i) (SM_BG_SAMPLE(sm,T_NTH_V(SM_NTH_TRI(sm,i),0)) \|\| \
191	gwlarson	3.1	SM_BG_SAMPLE(sm,T_NTH_V(SM_NTH_TRI(sm,i),1)) \|\| \
192			SM_BG_SAMPLE(sm,T_NTH_V(SM_NTH_TRI(sm,i),2)))
193
194			#define SM_FOR_ALL_SAMPLES(sm,i) for((i)=0;i < SM_NUM_SAMP(sm);(i)++)
195
196	gwlarson	3.10	#define SM_FOR_ALL_VALID_TRIS(m,i) for((i)=0,(i)=smNext_valid_tri(m,i);(i)< \
197			SM_NUM_TRI(m); (i)++,(i)= smNext_valid_tri(m,i))
198
199			#define SM_FOR_ALL_FLAGGED_TRIS(m,i,w,b) for(i=0,i=smNext_tri_flag_set(m,i,w,b);i < SM_NUM_TRI(m);i++,i=smNext_tri_flag_set(m,i,w,b))
200			#define SM_FOR_ALL_ACTIVE_TRIS(m,i) SM_FOR_ALL_FLAGGED_TRIS(m,i,T_ACTIVE_FLAG,0)
201			#define smInit_locator(sm) stInit(SM_LOCATOR(sm))
202	gwlarson	3.4	#define smAlloc_locator(sm) stAlloc(SM_LOCATOR(sm))
203	gwlarson	3.10	#define smFree_locator(sm) stFree(SM_LOCATOR(sm))
204	gwlarson	3.4	#define smUnalloc_samp(sm,id) sUnalloc_samp(SM_SAMP(sm),id)
205	gwlarson	3.10	#define smPoint_locate_cell(sm,pt) stPoint_locate(SM_LOCATOR(sm),pt)
206	gwlarson	3.4	#define smFree_samples(sm) sFree(SM_SAMP(sm))
207			#define smInit_samples(sm) sInit(SM_SAMP(sm))
208
209	gwlarson	3.10	#define SM_S_NTH_QT(sm,s_id) S_NTH_INFO1(SM_SAMP(sm),s_id)
210	gwlarson	3.4	#define smClear_vert(sm,id) (SM_NTH_VERT(sm,id) = INVALID)
211
212	gwlarson	3.10	#define freebuf(b) tempbuf(-1)
213
214	gwlarson	3.4	typedef struct _RT_ARGS_{
215			FVECT orig,dir;
216			int t_id;
217	greg	3.11	S_ID *os;
218	gwlarson	3.4	}RT_ARGS;
219	gwlarson	3.1
220	gwlarson	3.10	typedef struct _S_ARGS_{
221	greg	3.11	S_ID s_id;
222			S_ID n_id;
223	gwlarson	3.10	}S_ARGS;
224	gwlarson	3.4
225	gwlarson	3.10
226	gwlarson	3.4	typedef struct _ADD_ARGS {
227			int t_id;
228	greg	3.11	S_ID *del_set;
229	gwlarson	3.4	}ADD_ARGS;
230
231	gwlarson	3.1	extern SM *smMesh;
232			extern double smDist_sum;
233
234	gwlarson	3.4
235	gwlarson	3.1	#ifdef TEST_DRIVER
236			extern VIEW View;
237			extern VIEW Current_View;
238	gwlarson	3.2	extern int Pick_tri,Picking,Pick_samp;
239	gwlarson	3.1	extern FVECT Pick_point[500],Pick_origin,Pick_dir;
240			extern FVECT Pick_v0[500],Pick_v1[500],Pick_v2[500];
241	gwlarson	3.5	extern int Pick_q[500];
242	gwlarson	3.1	extern FVECT P0,P1,P2;
243			extern int Pick_cnt;
244			extern FVECT FrustumNear[4],FrustumFar[4];
245			#endif
246
247	gwlarson	3.4
248			/*
249			* int
250			* smInit(n) : Initialize/clear data structures for n entries
251			* int n;
252			*
253			* Initialize sampL and other data structures for at least n samples.
254			* If n is 0, then free data structures. Return number actually allocated.
255			*
256			*
257			* int
258			* smNewSamp(c, p, v) : register new sample point and return index
259			* COLR c; : pixel color (RGBE)
260			* FVECT p; : world intersection point
261			* FVECT v; : ray direction vector
262			*
263			* Add new sample point to data structures, removing old values as necessary.
264			* New sample representation will be output in next call to smUpdate().
265			*
266			*
267			* int
268			* smFindSamp(orig, dir): intersect ray with 3D rep. and find closest sample
269			* FVECT orig, dir;
270			*
271			* Find the closest sample to the given ray. Return -1 on failure.
272			*
273			*/
274	schorsch	3.13
275			#ifdef __cplusplus
276			}
277	gwlarson	3.1	#endif
278	schorsch	3.13	#endif /* define _RAD_SM_H_ */
279	gwlarson	3.1