glsurfeval.h revision f220fa62
1/*
2 * SGI FREE SOFTWARE LICENSE B (Version 2.0, Sept. 18, 2008)
3 * Copyright (C) 1991-2000 Silicon Graphics, Inc. All Rights Reserved.
4 *
5 * Permission is hereby granted, free of charge, to any person obtaining a
6 * copy of this software and associated documentation files (the "Software"),
7 * to deal in the Software without restriction, including without limitation
8 * the rights to use, copy, modify, merge, publish, distribute, sublicense,
9 * and/or sell copies of the Software, and to permit persons to whom the
10 * Software is furnished to do so, subject to the following conditions:
11 *
12 * The above copyright notice including the dates of first publication and
13 * either this permission notice or a reference to
14 * http://oss.sgi.com/projects/FreeB/
15 * shall be included in all copies or substantial portions of the Software.
16 *
17 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
18 * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
19 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
20 * SILICON GRAPHICS, INC. BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY,
21 * WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF
22 * OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
23 * SOFTWARE.
24 *
25 * Except as contained in this notice, the name of Silicon Graphics, Inc.
26 * shall not be used in advertising or otherwise to promote the sale, use or
27 * other dealings in this Software without prior written authorization from
28 * Silicon Graphics, Inc.
29 */
30
31/*
32 * glsurfeval.h
33 *
34 */
35
36#ifndef __gluglsurfeval_h_
37#define __gluglsurfeval_h_
38
39#include "basicsurfeval.h"
40#include "bezierPatchMesh.h" //in case output triangles
41#include <GL/gl.h>
42#include <GL/glu.h>
43
44class SurfaceMap;
45class OpenGLSurfaceEvaluator;
46class StoredVertex;
47
48#define TYPECOORD	1
49#define TYPEPOINT	2
50
51/* Cache up to 3 vertices from tmeshes */
52#define VERTEX_CACHE_SIZE	3
53
54/*for internal evaluator callback stuff*/
55#ifndef IN_MAX_BEZIER_ORDER
56#define IN_MAX_BEZIER_ORDER 40 /*XXX should be bigger than machine order*/
57#endif
58
59#ifndef IN_MAX_DIMENSION
60#define IN_MAX_DIMENSION 4
61#endif
62
63typedef struct surfEvalMachine{
64  REAL uprime;//cached previusly evaluated uprime.
65  REAL vprime;
66  int k; /*the dimension*/
67  REAL u1;
68  REAL u2;
69  int ustride;
70  int uorder;
71  REAL v1;
72  REAL v2;
73  int vstride;
74  int vorder;
75  REAL ctlPoints[IN_MAX_BEZIER_ORDER*IN_MAX_BEZIER_ORDER*IN_MAX_DIMENSION];
76  REAL ucoeff[IN_MAX_BEZIER_ORDER]; /*cache the polynomial values*/
77  REAL vcoeff[IN_MAX_BEZIER_ORDER];
78  REAL ucoeffDeriv[IN_MAX_BEZIER_ORDER]; /*cache the polynomial derivatives*/
79  REAL vcoeffDeriv[IN_MAX_BEZIER_ORDER];
80} surfEvalMachine;
81
82
83
84class StoredVertex {
85public:
86    		StoredVertex() { type = 0; coord[0] = 0; coord[1] = 0; point[0] = 0; point[1] = 0; }
87		~StoredVertex(void) {}
88    void	saveEvalCoord(REAL x, REAL y)
89		    {coord[0] = x; coord[1] = y; type = TYPECOORD; }
90    void	saveEvalPoint(long x, long y)
91		    {point[0] = x; point[1] = y; type = TYPEPOINT; }
92    void	invoke(OpenGLSurfaceEvaluator *eval);
93
94private:
95    int		type;
96    REAL	coord[2];
97    long	point[2];
98};
99
100class OpenGLSurfaceEvaluator : public BasicSurfaceEvaluator {
101public:
102			OpenGLSurfaceEvaluator();
103    			virtual ~OpenGLSurfaceEvaluator( void );
104    void		polymode( long style );
105    void		range2f( long, REAL *, REAL * );
106    void		domain2f( REAL, REAL, REAL, REAL );
107    void		addMap( SurfaceMap * ) { }
108
109    void		enable( long );
110    void		disable( long );
111    void		bgnmap2f( long );
112    void		map2f( long, REAL, REAL, long, long,
113				     REAL, REAL, long, long, REAL * );
114    void		mapgrid2f( long, REAL, REAL, long, REAL, REAL );
115    void		mapmesh2f( long, long, long, long, long );
116    void		evalcoord2f( long, REAL, REAL );
117    void		evalpoint2i( long, long );
118    void		endmap2f( void );
119
120    void	 	bgnline( void );
121    void	 	endline( void );
122    void	 	bgnclosedline( void );
123    void	 	endclosedline( void );
124    void	 	bgntmesh( void );
125    void	 	swaptmesh( void );
126    void	 	endtmesh( void );
127    void	 	bgnqstrip( void );
128    void	 	endqstrip( void );
129
130    void                bgntfan( void );
131    void                endtfan( void );
132    void                evalUStrip(int n_upper, REAL v_upper, REAL* upper_val,
133                                   int n_lower, REAL v_lower, REAL* lower_val);
134    void                evalVStrip(int n_left, REAL u_left, REAL* left_val,
135                                   int n_right, REAL u_right, REAL* right_val);
136
137    void		coord2f( REAL, REAL );
138    void		point2i( long, long );
139
140    void		newtmeshvert( REAL, REAL );
141    void		newtmeshvert( long, long );
142
143#ifdef _WIN32
144    void 	        putCallBack(GLenum which, void (GLAPIENTRY *fn)() );
145#else
146    void 	        putCallBack(GLenum which, _GLUfuncptr fn );
147#endif
148
149    int                 get_vertices_call_back()
150      {
151	return output_triangles;
152      }
153    void                put_vertices_call_back(int flag)
154      {
155	output_triangles = flag;
156      }
157
158    void                 put_callback_auto_normal(int flag)
159      {
160        callback_auto_normal = flag;
161      }
162
163   int                   get_callback_auto_normal()
164     {
165        return callback_auto_normal;
166      }
167
168   void                  set_callback_userData(void* data)
169     {
170       userData = data;
171     }
172
173    /**************begin for LOD_eval_list***********/
174    void LOD_eval_list(int level);
175
176
177
178
179private:
180    StoredVertex	*vertexCache[VERTEX_CACHE_SIZE];
181    int			tmeshing;
182    int			which;
183    int			vcount;
184
185    GLint              gl_polygon_mode[2];/*to save and restore so that
186					 *no side effect
187					 */
188    bezierPatchMesh        *global_bpm; //for output triangles
189    int                output_triangles; //true 1 or false 0
190
191
192
193    void (GLAPIENTRY *beginCallBackN) (GLenum type);
194    void (GLAPIENTRY *endCallBackN)   (void);
195    void (GLAPIENTRY *vertexCallBackN) (const GLfloat *vert);
196    void (GLAPIENTRY *normalCallBackN) (const GLfloat *normal);
197    void (GLAPIENTRY *colorCallBackN) (const GLfloat *color);
198    void (GLAPIENTRY *texcoordCallBackN) (const GLfloat *texcoord);
199
200    void (GLAPIENTRY *beginCallBackData) (GLenum type, void* data);
201    void (GLAPIENTRY *endCallBackData)   (void* data);
202    void (GLAPIENTRY *vertexCallBackData) (const GLfloat *vert, void* data);
203    void (GLAPIENTRY *normalCallBackData) (const GLfloat *normal, void* data);
204    void (GLAPIENTRY *colorCallBackData) (const GLfloat *color, void* data);
205    void (GLAPIENTRY *texcoordCallBackData) (const GLfloat *texcoord, void* data);
206
207    void               beginCallBack (GLenum type, void* data);
208    void               endCallBack   (void* data);
209    void               vertexCallBack (const GLfloat *vert, void* data);
210    void               normalCallBack (const GLfloat *normal, void* data);
211    void               colorCallBack (const GLfloat *color, void* data);
212    void               texcoordCallBack (const GLfloat *texcoord, void* data);
213
214
215    void* userData; //the opaque pointer for Data callback functions.
216
217   /*LOD evaluation*/
218   void LOD_triangle(REAL A[2], REAL B[2], REAL C[2],
219		     int level);
220   void LOD_eval(int num_vert, REAL* verts, int type, int level);
221
222  int LOD_eval_level; //set by LOD_eval_list()
223
224   /*************begin for internal evaluators*****************/
225
226 /*the following global variables are only defined in this file.
227 *They are used to cache the precomputed Bezier polynomial values.
228 *These calues may be used consecutively in which case we don't have
229 *recompute these values again.
230 */
231 int global_uorder; /*store the uorder in the previous evaluation*/
232 int global_vorder; /*store the vorder in the previous evaluation*/
233 REAL global_uprime;
234 REAL global_vprime;
235 REAL global_vprime_BV;
236 REAL global_uprime_BU;
237 int global_uorder_BV; /*store the uorder in the previous evaluation*/
238 int global_vorder_BV; /*store the vorder in the previous evaluation*/
239 int global_uorder_BU; /*store the uorder in the previous evaluation*/
240 int global_vorder_BU; /*store the vorder in the previous evaluation*/
241
242 REAL global_ucoeff[IN_MAX_BEZIER_ORDER]; /*cache the polynomial values*/
243 REAL global_vcoeff[IN_MAX_BEZIER_ORDER];
244 REAL global_ucoeffDeriv[IN_MAX_BEZIER_ORDER]; /*cache the polynomial derivatives*/
245 REAL global_vcoeffDeriv[IN_MAX_BEZIER_ORDER];
246
247 REAL global_BV[IN_MAX_BEZIER_ORDER][IN_MAX_DIMENSION];
248 REAL global_PBV[IN_MAX_BEZIER_ORDER][IN_MAX_DIMENSION];
249 REAL global_BU[IN_MAX_BEZIER_ORDER][IN_MAX_DIMENSION];
250 REAL global_PBU[IN_MAX_BEZIER_ORDER][IN_MAX_DIMENSION];
251 REAL* global_baseData;
252
253 int    global_ev_k; /*the dimension*/
254 REAL global_ev_u1;
255 REAL global_ev_u2;
256 int    global_ev_ustride;
257 int    global_ev_uorder;
258 REAL global_ev_v1;
259 REAL global_ev_v2;
260 int    global_ev_vstride;
261 int    global_ev_vorder;
262 REAL global_ev_ctlPoints[IN_MAX_BEZIER_ORDER*IN_MAX_BEZIER_ORDER*IN_MAX_DIMENSION];
263
264 REAL  global_grid_u0;
265 REAL  global_grid_u1;
266 int     global_grid_nu;
267 REAL  global_grid_v0;
268 REAL  global_grid_v1;
269 int     global_grid_nv;
270
271/*functions*/
272 void inDoDomain2WithDerivs(int k, REAL u, REAL v,
273				REAL u1, REAL u2, int uorder,
274				REAL v1,  REAL v2, int vorder,
275				REAL *baseData,
276				REAL *retPoint, REAL *retdu, REAL *retdv);
277 void inPreEvaluate(int order, REAL vprime, REAL *coeff);
278 void inPreEvaluateWithDeriv(int order, REAL vprime, REAL *coeff, REAL *coeffDeriv);
279 void inComputeFirstPartials(REAL *p, REAL *pu, REAL *pv);
280 void inComputeNormal2(REAL *pu, REAL *pv, REAL *n);
281 void inDoEvalCoord2(REAL u, REAL v,
282		     REAL *retPoint, REAL *retNormal);
283 void inDoEvalCoord2NOGE(REAL u, REAL v,
284		     REAL *retPoint, REAL *retNormal);
285 void inMap2f(int k,
286	      REAL ulower,
287	      REAL uupper,
288	      int ustride,
289	      int uorder,
290	      REAL vlower,
291	      REAL vupper,
292	      int vstride,
293	      int vorder,
294	      REAL *ctlPoints);
295
296 void inMapGrid2f(int nu, REAL u0, REAL u1,
297		  int nv, REAL v0, REAL v1);
298
299 void inEvalMesh2(int lowU, int lowV, int highU, int highV);
300 void inEvalPoint2(int i, int j);
301 void inEvalCoord2f(REAL u, REAL v);
302
303void inEvalULine(int n_points, REAL v, REAL* u_vals,
304	int stride, REAL ret_points[][3], REAL ret_normals[][3]);
305
306void inEvalVLine(int n_points, REAL u, REAL* v_vals,
307	int stride, REAL ret_points[][3], REAL ret_normals[][3]);
308
309void inEvalUStrip(int n_upper, REAL v_upper, REAL* upper_val,
310                       int n_lower, REAL v_lower, REAL* lower_val
311                       );
312void inEvalVStrip(int n_left, REAL u_left, REAL* left_val, int n_right, REAL u_right, REAL* right_val);
313
314void inPreEvaluateBV(int k, int uorder, int vorder, REAL vprime, REAL *baseData);
315void inPreEvaluateBU(int k, int uorder, int vorder, REAL uprime, REAL *baseData);
316void inPreEvaluateBV_intfac(REAL v )
317  {
318   inPreEvaluateBV(global_ev_k, global_ev_uorder, global_ev_vorder, (v-global_ev_v1)/(global_ev_v2-global_ev_v1), global_ev_ctlPoints);
319  }
320
321void inPreEvaluateBU_intfac(REAL u)
322  {
323    inPreEvaluateBU(global_ev_k, global_ev_uorder, global_ev_vorder, (u-global_ev_u1)/(global_ev_u2-global_ev_u1), global_ev_ctlPoints);
324  }
325
326void inDoDomain2WithDerivsBV(int k, REAL u, REAL v,
327			     REAL u1, REAL u2, int uorder,
328			     REAL v1, REAL v2, int vorder,
329			     REAL *baseData,
330			     REAL *retPoint, REAL* retdu, REAL *retdv);
331
332void inDoDomain2WithDerivsBU(int k, REAL u, REAL v,
333			     REAL u1, REAL u2, int uorder,
334			     REAL v1, REAL v2, int vorder,
335			     REAL *baseData,
336			     REAL *retPoint, REAL* retdu, REAL *retdv);
337
338
339void inDoEvalCoord2NOGE_BV(REAL u, REAL v,
340			   REAL *retPoint, REAL *retNormal);
341
342void inDoEvalCoord2NOGE_BU(REAL u, REAL v,
343			   REAL *retPoint, REAL *retNormal);
344
345void inBPMEval(bezierPatchMesh* bpm);
346void inBPMListEval(bezierPatchMesh* list);
347
348/*-------------begin for surfEvalMachine -------------*/
349surfEvalMachine em_vertex;
350surfEvalMachine em_normal;
351surfEvalMachine em_color;
352surfEvalMachine em_texcoord;
353
354int auto_normal_flag; //whether to output normla or not in callback
355                      //determined by GL_AUTO_NORMAL and callback_auto_normal
356int callback_auto_normal; //GLU_CALLBACK_AUTO_NORMAL_EXT
357int vertex_flag;
358int normal_flag;
359int color_flag;
360int texcoord_flag;
361
362void inMap2fEM(int which, //0:vert,1:norm,2:color,3:tex
363	       int dimension,
364	      REAL ulower,
365	      REAL uupper,
366	      int ustride,
367	      int uorder,
368	      REAL vlower,
369	      REAL vupper,
370	      int vstride,
371	      int vorder,
372	      REAL *ctlPoints);
373
374void inDoDomain2WithDerivsEM(surfEvalMachine *em, REAL u, REAL v,
375				REAL *retPoint, REAL *retdu, REAL *retdv);
376void inDoDomain2EM(surfEvalMachine *em, REAL u, REAL v,
377				REAL *retPoint);
378 void inDoEvalCoord2EM(REAL u, REAL v);
379
380void inBPMEvalEM(bezierPatchMesh* bpm);
381void inBPMListEvalEM(bezierPatchMesh* list);
382
383/*-------------end for surfEvalMachine -------------*/
384
385
386   /*************end for internal evaluators*****************/
387
388};
389
390inline void StoredVertex::invoke(OpenGLSurfaceEvaluator *eval)
391{
392    switch(type) {
393      case TYPECOORD:
394	eval->coord2f(coord[0], coord[1]);
395	break;
396      case TYPEPOINT:
397	eval->point2i(point[0], point[1]);
398	break;
399      default:
400	break;
401    }
402}
403
404#endif /* __gluglsurfeval_h_ */
405