ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/hd/sm.h
Revision: 3.8
Committed: Sun Jan 10 10:27:47 1999 UTC (25 years, 3 months ago) by gwlarson
Content type: text/plain
Branch: MAIN
Changes since 3.7: +5 -29 lines
Log Message:
sm.c, sm_del.c sm_ogl.c sm.h
Fixed failure to tone-map on start up
added code for temporary buffer allocation
added bg flag, and fast flag checking routines for drawing
provided graceful backout if deletion triangulation fails

File Contents

# Content
1 /* Copyright (c) 1998 Silicon Graphics, Inc. */
2
3 /* SCCSid "$SunId$ SGI" */
4
5 /*
6 * sm.h
7 */
8
9 #ifndef _SM_H_
10 #define _SM_H_
11
12 #include "rhd_sample.h"
13 #define NEWSETS
14
15
16 #ifndef TRUE
17 #define TRUE 1
18 #define FALSE 0
19 #endif
20
21
22 #define ON_V 1
23 #define ON_P 2
24 #define ON_E 3
25 #define IN_T 4
26
27 #define S_REPLACE_EPS 0.04 /* if (distance on sphere between sample
28 and a base point) < S_REPLACE_EPS,
29 replace base:
30 */
31 #define S_REPLACE_SCALE (5.*5.) /* if (distance to new point squared) is
32 > (triangle edge length squared*
33 S_REPLACE_SCALE):for all edges/triangle
34 vertices: new point is puncture
35 point: dont add
36 */
37 #define S_REPLACE_TRI 2e-8 /* .052 radians to the sixth power */
38
39 #define SQRT3_2 0.8660254
40
41 #define SM_DEFAULT 0
42 #define SM_EXTRA_POINTS 162
43 #define SM_BASE_TRIS 320
44 #define SM_EXTRA_VERTS SM_EXTRA_POINTS
45
46 #define SM_INC_PERCENT 0.60 /* If number of new triangles added
47 since last full redraw is >
48 (SM_INC_PERCENT * total triangles)
49 do full redraw instead of incremental
50 */
51
52 #define SM_VIEW_FRAC 0.1
53
54
55 typedef int VERT; /* One triangle that vertex belongs to- the rest
56 are derived by traversing neighbors */
57
58 typedef struct _EDGE {
59 int verts[2];
60 int tris[2];
61 } EDGE;
62
63 #define E_NTH_VERT(e,i) ((e>0)?Edges[(e)].verts[(i)]:Edges[-(e)].verts[(1-i)])
64 #define SET_E_NTH_VERT(e,i,v) if(e>0) Edges[(e)].verts[(i)]=(v); \
65 else Edges[-(e)].verts[(1-i)]=(v)
66 #define E_NTH_TRI(e,i) ((e>0)?Edges[(e)].tris[(i)]:Edges[-(e)].tris[(1-i)])
67 #define SET_E_NTH_TRI(e,i,v) if(e>0) Edges[(e)].tris[(i)]=(v); \
68 else Edges[-(e)].tris[(1-i)]=(v)
69 #define eClear_edges() (Ecnt = 0)
70
71 #define FOR_ALL_EDGES(i) for((i)=1; (i) <= Ecnt; i++)
72 #define FOR_ALL_EDGES_FROM(e,i) for((i)=++e; (i) <= Ecnt; i++)
73
74
75 typedef struct _TRI {
76 int verts[3]; /* Ids into sample and vertex array for each vertex*/
77 int nbrs[3]; /* Ids for neighboring triangles: -1 if invalid */
78 }TRI;
79
80
81 #define T_NTH_NBR(t,i) ((t)->nbrs[(i)])
82 #define T_CLEAR_NBRS(t) (T_NTH_NBR(t,0)=T_NTH_NBR(t,1)=T_NTH_NBR(t,2)=-1)
83 #define T_NTH_NBR_PTR(t,n) \
84 (T_NTH_NBR(n,0)==(t)?0:T_NTH_NBR(n,1)==(t)?1:T_NTH_NBR(n,2)==(t)?2:-1)
85 #define T_NTH_V(t,i) ((t)->verts[(i)])
86 #define T_WHICH_V(t,i) \
87 (T_NTH_V(t,0)==(i)?0:T_NTH_V(t,1)==(i)?1:T_NTH_V(t,2)==(i)?2:-1)
88 #define T_NEXT_FREE(t) ((t)->nbrs[0])
89 #define T_NEXT_AVAILABLE(t) ((t)->nbrs[0])
90 #define T_VALID_FLAG(t) ((t)->nbrs[1])
91 #define T_IS_VALID(t) (T_VALID_FLAG(t)!=-1)
92 #define T_FLAGS 4
93
94 typedef struct _SM {
95 FVECT view_center; /* Canonical view center defining unit sphere */
96 SAMP *samples; /* Sample point information */
97 STREE locator; /* spherical quadtree for point/triangle location */
98 int max_tris; /* Maximum number of triangles */
99 int num_tri; /* Current number of triangles */
100 int sample_tris; /* Current number of non-base triangles*/
101 int free_tris; /* pointer to free_list */
102 int available_tris; /* pointer to available_list */
103 int max_verts; /* Maximum number of vertices in the mesh */
104 TRI *tris; /* Pointer to list of triangle structs */
105 VERT *verts; /* List of vertices */
106 int4 *flags[T_FLAGS]; /* Bit 0 set if active(in current frustum) */
107 /* Bit 1 set if not rendered since created */
108 /* Bit 2 set if base triangle */
109 }SM;
110
111 #define T_ACTIVE_FLAG 0
112 #define T_NEW_FLAG 1
113 #define T_BASE_FLAG 2
114 #define T_BG_FLAG 3
115
116 #define SM_VIEW_CENTER(m) ((m)->view_center)
117 #define SM_SAMP(m) ((m)->samples)
118 #define SM_LOCATOR(m) (&((m)->locator))
119 #define SM_MAX_TRIS(m) ((m)->max_tris)
120 #define SM_NUM_TRI(m) ((m)->num_tri)
121 #define SM_SAMPLE_TRIS(m) ((m)->sample_tris)
122 #define SM_FREE_TRIS(m) ((m)->free_tris)
123 #define SM_AVAILABLE_TRIS(m) ((m)->available_tris)
124 #define SM_MAX_VERTS(m) ((m)->max_verts)
125 #define SM_TRIS(m) ((m)->tris)
126 #define SM_VERTS(m) ((m)->verts)
127 #define SM_NTH_FLAGS(m,n) ((m)->flags[(n)])
128 #define SM_FLAGS(m) ((m)->flags)
129
130
131 #define SM_IS_NTH_T_FLAG(sm,n,f) IS_FLAG(SM_NTH_FLAGS(sm,f),n)
132 #define SM_SET_NTH_T_FLAG(sm,n,f) SET_FLAG(SM_NTH_FLAGS(sm,f),n)
133 #define SM_CLR_NTH_T_FLAG(sm,n,f) CLR_FLAG(SM_NTH_FLAGS(sm,f),n)
134
135 #define SM_IS_NTH_T_ACTIVE(sm,n) SM_IS_NTH_T_FLAG(sm,n,T_ACTIVE_FLAG)
136 #define SM_IS_NTH_T_BASE(sm,n) SM_IS_NTH_T_FLAG(sm,n,T_BASE_FLAG)
137 #define SM_IS_NTH_T_NEW(sm,n) SM_IS_NTH_T_FLAG(sm,n,T_NEW_FLAG)
138 #define SM_IS_NTH_T_BG(sm,n) SM_IS_NTH_T_FLAG(sm,n,T_BG_FLAG)
139
140 #define SM_SET_NTH_T_ACTIVE(sm,n) SM_SET_NTH_T_FLAG(sm,n,T_ACTIVE_FLAG)
141 #define SM_SET_NTH_T_BASE(sm,n) SM_SET_NTH_T_FLAG(sm,n,T_BASE_FLAG)
142 #define SM_SET_NTH_T_NEW(sm,n) SM_SET_NTH_T_FLAG(sm,n,T_NEW_FLAG)
143 #define SM_SET_NTH_T_BG(sm,n) SM_SET_NTH_T_FLAG(sm,n,T_BG_FLAG)
144
145 #define SM_CLR_NTH_T_ACTIVE(sm,n) SM_CLR_NTH_T_FLAG(sm,n,T_ACTIVE_FLAG)
146 #define SM_CLR_NTH_T_BASE(sm,n) SM_CLR_NTH_T_FLAG(sm,n,T_BASE_FLAG)
147 #define SM_CLR_NTH_T_NEW(sm,n) SM_CLR_NTH_T_FLAG(sm,n,T_NEW_FLAG)
148 #define SM_CLR_NTH_T_BG(sm,n) SM_CLR_NTH_T_FLAG(sm,n,T_BG_FLAG)
149
150 #define SM_NTH_TRI(m,n) (&(SM_TRIS(m)[(n)]))
151 #define SM_NTH_VERT(m,n) (SM_VERTS(m)[(n)])
152
153 #define SM_T_ID_VALID(s,t_id) T_IS_VALID(SM_NTH_TRI(s,t_id))
154
155
156 #define SM_MAX_SAMP(m) S_MAX_SAMP(SM_SAMP(m))
157 #define SM_MAX_POINTS(m) S_MAX_POINTS(SM_SAMP(m))
158 #define SM_SAMP_BASE(m) S_BASE(SM_SAMP(m))
159 #define SM_NTH_WV(m,i) S_NTH_W_PT(SM_SAMP(m),i)
160 #define SM_NTH_W_DIR(m,i) S_NTH_W_DIR(SM_SAMP(m),i)
161 #define SM_DIR_ID(m,i) (!SM_BASE_ID(m,i) && SM_NTH_W_DIR(m,i)==-1)
162 #define SM_NTH_RGB(m,i) S_NTH_RGB(SM_SAMP(m),i)
163 #define SM_RGB(m) S_RGB(SM_SAMP(m))
164 #define SM_WP(m) S_W_PT(SM_SAMP(m))
165 #define SM_BRT(m) S_BRT(SM_SAMP(m))
166 #define SM_NTH_BRT(m,i) S_NTH_BRT(SM_SAMP(m),i)
167 #define SM_CHR(m) S_CHR(SM_SAMP(m))
168 #define SM_NTH_CHR(m,i) S_NTH_CHR(SM_SAMP(m),i)
169 #define SM_NUM_SAMP(m) S_NUM_SAMP(SM_SAMP(m))
170 #define SM_TONE_MAP(m) S_TONE_MAP(SM_SAMP(m))
171
172 #define SM_ALLOWED_VIEW_CHANGE(m) (SM_NUM_SAMP(m)/smDist_sum*SM_VIEW_FRAC)
173
174 #define SM_FOR_ALL_ADJACENT_TRIS(sm,id,t) for(t=smTri_next_ccw_nbr(sm,t,id); \
175 t!=SM_NTH_TRI(sm,SM_NTH_VERT(sm,id)); t=smTri_next_ccw_nbr(sm,t,id))
176
177 #define SM_INVALID_SAMP_ID(sm,id) (((id) < 0) || ((id) >= SM_MAX_SAMP(sm)))
178 #define SM_INVALID_POINT_ID(sm,id) (((id) < 0) || ((id) >= SM_MAX_POINTS(sm)))
179 #define SM_T_NTH_WV(sm,t,i) (SM_NTH_WV(sm,T_NTH_V(t,i)))
180
181 #define SM_BASE_ID(s,i) \
182 ((i) >= S_MAX_SAMP(SM_SAMP(s)) && (i) < S_MAX_BASE_PT(SM_SAMP(s)))
183
184 #define SM_BG_SAMPLE(sm,i) (SM_NTH_W_DIR(sm,i)==-1)
185
186 #define SM_BG_TRI(sm,i) (SM_BG_SAMPLE(sm,T_NTH_V(SM_NTH_TRI(sm,i),0)) && \
187 SM_BG_SAMPLE(sm,T_NTH_V(SM_NTH_TRI(sm,i),1)) && \
188 SM_BG_SAMPLE(sm,T_NTH_V(SM_NTH_TRI(sm,i),2)))
189 #define SM_MIXED_TRI(sm,i) (SM_BG_SAMPLE(sm,T_NTH_V(SM_NTH_TRI(sm,i),0)) || \
190 SM_BG_SAMPLE(sm,T_NTH_V(SM_NTH_TRI(sm,i),1)) || \
191 SM_BG_SAMPLE(sm,T_NTH_V(SM_NTH_TRI(sm,i),2)))
192
193 #define SM_FOR_ALL_SAMPLES(sm,i) for((i)=0;i < SM_NUM_SAMP(sm);(i)++)
194
195 #define smInit_locator(sm,c) (stInit(SM_LOCATOR(sm)), \
196 ST_SET_CENTER(SM_LOCATOR(sm),c))
197 #define smClear_locator(sm) stClear(SM_LOCATOR(sm))
198 #define smAlloc_locator(sm) stAlloc(SM_LOCATOR(sm))
199 #define smFree_locator(sm) stClear(SM_LOCATOR(sm))
200 #define smPointLocateCell(sm,pt) stPoint_locate(SM_LOCATOR(sm),pt)
201 #define smUnalloc_samp(sm,id) sUnalloc_samp(SM_SAMP(sm),id)
202 #define smFree_samples(sm) sFree(SM_SAMP(sm))
203 #define smClear_samples(sm) sClear(SM_SAMP(sm))
204 #define smInit_samples(sm) sInit(SM_SAMP(sm))
205
206
207 #define smClear_vert(sm,id) (SM_NTH_VERT(sm,id) = INVALID)
208
209 typedef struct _RT_ARGS_{
210 FVECT orig,dir;
211 int t_id;
212 OBJECT *os;
213 }RT_ARGS;
214
215
216 typedef struct _ADD_ARGS {
217 int t_id;
218 OBJECT *del_set;
219 }ADD_ARGS;
220
221 extern SM *smMesh;
222 extern int smNew_tri_cnt;
223 extern double smDist_sum;
224
225
226 #ifdef TEST_DRIVER
227 extern VIEW View;
228 extern VIEW Current_View;
229 extern int Pick_tri,Picking,Pick_samp;
230 extern FVECT Pick_point[500],Pick_origin,Pick_dir;
231 extern FVECT Pick_v0[500],Pick_v1[500],Pick_v2[500];
232 extern int Pick_q[500];
233 extern FVECT P0,P1,P2;
234 extern int Pick_cnt;
235 extern FVECT FrustumNear[4],FrustumFar[4];
236 #endif
237
238
239 /*
240 * int
241 * smInit(n) : Initialize/clear data structures for n entries
242 * int n;
243 *
244 * Initialize sampL and other data structures for at least n samples.
245 * If n is 0, then free data structures. Return number actually allocated.
246 *
247 *
248 * int
249 * smNewSamp(c, p, v) : register new sample point and return index
250 * COLR c; : pixel color (RGBE)
251 * FVECT p; : world intersection point
252 * FVECT v; : ray direction vector
253 *
254 * Add new sample point to data structures, removing old values as necessary.
255 * New sample representation will be output in next call to smUpdate().
256 *
257 *
258 * int
259 * smFindSamp(orig, dir): intersect ray with 3D rep. and find closest sample
260 * FVECT orig, dir;
261 *
262 * Find the closest sample to the given ray. Return -1 on failure.
263 *
264 *
265 * smClean() : display has been wiped clean
266 *
267 * Called after display has been effectively cleared, meaning that all
268 * geometry must be resent down the pipeline in the next call to smUpdate().
269 *
270 *
271 * smUpdate(vp, qua) : update OpenGL output geometry for view vp
272 * VIEW *vp; : desired view
273 * int qua; : quality level (percentage on linear time scale)
274 *
275 * Draw new geometric representation using OpenGL calls. Assume that the
276 * view has already been set up and the correct frame buffer has been
277 * selected for drawing. The quality level is on a linear scale, where 100%
278 * is full (final) quality. It is not necessary to redraw geometry that has
279 * been output since the last call to smClean(). (The last view drawn will
280 * be vp==&odev.v each time.)
281 */
282 #endif
283
284
285
286
287
288
289
290