src/glu/sgi/libnurbs/interface/glsurfeval.h

   1 /*
   2 ** License Applicability. Except to the extent portions of this file are
   3 ** made subject to an alternative license as permitted in the SGI Free
   4 ** Software License B, Version 1.1 (the "License"), the contents of this
   5 ** file are subject only to the provisions of the License. You may not use
   6 ** this file except in compliance with the License. You may obtain a copy
   7 ** of the License at Silicon Graphics, Inc., attn: Legal Services, 1600
   8 ** Amphitheatre Parkway, Mountain View, CA 94043-1351, or at:
   9 **
  10 ** http://oss.sgi.com/projects/FreeB
  11 **
  12 ** Note that, as provided in the License, the Software is distributed on an
  13 ** "AS IS" basis, with ALL EXPRESS AND IMPLIED WARRANTIES AND CONDITIONS
  14 ** DISCLAIMED, INCLUDING, WITHOUT LIMITATION, ANY IMPLIED WARRANTIES AND
  15 ** CONDITIONS OF MERCHANTABILITY, SATISFACTORY QUALITY, FITNESS FOR A
  16 ** PARTICULAR PURPOSE, AND NON-INFRINGEMENT.
  17 **
  18 ** Original Code. The Original Code is: OpenGL Sample Implementation,
  19 ** Version 1.2.1, released January 26, 2000, developed by Silicon Graphics,
  20 ** Inc. The Original Code is Copyright (c) 1991-2000 Silicon Graphics, Inc.
  21 ** Copyright in any portions created by third parties is as indicated
  22 ** elsewhere herein. All Rights Reserved.
  23 **
  24 ** Additional Notice Provisions: The application programming interfaces
  25 ** established by SGI in conjunction with the Original Code are The
  26 ** OpenGL(R) Graphics System: A Specification (Version 1.2.1), released
  27 ** April 1, 1999; The OpenGL(R) Graphics System Utility Library (Version
  28 ** 1.3), released November 4, 1998; and OpenGL(R) Graphics with the X
  29 ** Window System(R) (Version 1.3), released October 19, 1998. This software
  30 ** was created using the OpenGL(R) version 1.2.1 Sample Implementation
  31 ** published by SGI, but has not been independently verified as being
  32 ** compliant with the OpenGL(R) version 1.2.1 Specification.
  33 */
  34
  35 /*
  36  * glsurfeval.h
  37  *
  38  * $Date: 2001/07/16 15:46:42 $ $Revision: 1.2 $
  39  * $Header: /home/krh/git/sync/mesa-cvs-repo/Mesa/src/glu/sgi/libnurbs/interface/glsurfeval.h,v 1.2 2001/07/16 15:46:42 brianp Exp $
  40  */
  41
  42 #ifndef __gluglsurfeval_h_
  43 #define __gluglsurfeval_h_
  44
  45 #include "basicsurfeval.h"
  46 #include "bezierPatchMesh.h" //in case output triangles
  47 #include "GL/gl.h"
  48
  49 class SurfaceMap;
  50 class OpenGLSurfaceEvaluator;
  51 class StoredVertex;
  52
  53 #define TYPECOORD       1
  54 #define TYPEPOINT       2
  55
  56 /* Cache up to 3 vertices from tmeshes */
  57 #define VERTEX_CACHE_SIZE       3
  58
  59 /*for internal evaluator callback stuff*/
  60 #ifndef IN_MAX_BEZIER_ORDER
  61 #define IN_MAX_BEZIER_ORDER 40 /*XXX should be bigger than machine order*/
  62 #endif
  63
  64 #ifndef IN_MAX_DIMENSION
  65 #define IN_MAX_DIMENSION 4
  66 #endif
  67
  68 typedef struct surfEvalMachine{
  69   REAL uprime;//cached previusly evaluated uprime.
  70   REAL vprime;
  71   int k; /*the dimension*/
  72   REAL u1;
  73   REAL u2;
  74   int ustride;
  75   int uorder;
  76   REAL v1;
  77   REAL v2;
  78   int vstride;
  79   int vorder;
  80   REAL ctlPoints[IN_MAX_BEZIER_ORDER*IN_MAX_BEZIER_ORDER*IN_MAX_DIMENSION];
  81   REAL ucoeff[IN_MAX_BEZIER_ORDER]; /*cache the polynomial values*/
  82   REAL vcoeff[IN_MAX_BEZIER_ORDER];
  83   REAL ucoeffDeriv[IN_MAX_BEZIER_ORDER]; /*cache the polynomial derivatives*/
  84   REAL vcoeffDeriv[IN_MAX_BEZIER_ORDER];
  85 } surfEvalMachine;
  86
  87
  88
  89 class StoredVertex {
  90 public:
  91                 StoredVertex() { type = 0; }
  92                 ~StoredVertex(void) {}
  93     void        saveEvalCoord(REAL x, REAL y)
  94                     {coord[0] = x; coord[1] = y; type = TYPECOORD; }
  95     void        saveEvalPoint(long x, long y)
  96                     {point[0] = x; point[1] = y; type = TYPEPOINT; }
  97     void        invoke(OpenGLSurfaceEvaluator *eval);
  98
  99 private:
 100     int         type;
 101     REAL        coord[2];
 102     long        point[2];
 103 };
 104
 105 class OpenGLSurfaceEvaluator : public BasicSurfaceEvaluator {
 106 public:
 107                         OpenGLSurfaceEvaluator();
 108                         ~OpenGLSurfaceEvaluator( void );
 109     void                polymode( long style );
 110     void                range2f( long, REAL *, REAL * );
 111     void                domain2f( REAL, REAL, REAL, REAL );
 112     void                addMap( SurfaceMap * ) { }
 113
 114     void                enable( long );
 115     void                disable( long );
 116     void                bgnmap2f( long );
 117     void                map2f( long, REAL, REAL, long, long,
 118                                      REAL, REAL, long, long, REAL * );
 119     void                mapgrid2f( long, REAL, REAL, long, REAL, REAL );
 120     void                mapmesh2f( long, long, long, long, long );
 121     void                evalcoord2f( long, REAL, REAL );
 122     void                evalpoint2i( long, long );
 123     void                endmap2f( void );
 124
 125     void                bgnline( void );
 126     void                endline( void );
 127     void                bgnclosedline( void );
 128     void                endclosedline( void );
 129     void                bgntmesh( void );
 130     void                swaptmesh( void );
 131     void                endtmesh( void );
 132     void                bgnqstrip( void );
 133     void                endqstrip( void );
 134
 135     void                bgntfan( void );
 136     void                endtfan( void );
 137     void                evalUStrip(int n_upper, REAL v_upper, REAL* upper_val,
 138                                    int n_lower, REAL v_lower, REAL* lower_val);
 139     void                evalVStrip(int n_left, REAL u_left, REAL* left_val,
 140                                    int n_right, REAL u_right, REAL* right_val);
 141
 142     void                coord2f( REAL, REAL );
 143     void                point2i( long, long );
 144
 145     void                newtmeshvert( REAL, REAL );
 146     void                newtmeshvert( long, long );
 147
 148     void                putCallBack(GLenum which, _GLUfuncptr fn );
 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
 179 private:
 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
 303 void inEvalULine(int n_points, REAL v, REAL* u_vals,
 304         int stride, REAL ret_points[][3], REAL ret_normals[][3]);
 305
 306 void inEvalVLine(int n_points, REAL u, REAL* v_vals,
 307         int stride, REAL ret_points[][3], REAL ret_normals[][3]);
 308
 309 void inEvalUStrip(int n_upper, REAL v_upper, REAL* upper_val,
 310                        int n_lower, REAL v_lower, REAL* lower_val
 311                        );
 312 void inEvalVStrip(int n_left, REAL u_left, REAL* left_val, int n_right, REAL u_right, REAL* right_val);
 313
 314 void inPreEvaluateBV(int k, int uorder, int vorder, REAL vprime, REAL *baseData);
 315 void inPreEvaluateBU(int k, int uorder, int vorder, REAL uprime, REAL *baseData);
 316 void 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
 321 void 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
 326 void 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
 332 void 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
 339 void inDoEvalCoord2NOGE_BV(REAL u, REAL v,
 340                            REAL *retPoint, REAL *retNormal);
 341
 342 void inDoEvalCoord2NOGE_BU(REAL u, REAL v,
 343                            REAL *retPoint, REAL *retNormal);
 344
 345 void inBPMEval(bezierPatchMesh* bpm);
 346 void inBPMListEval(bezierPatchMesh* list);
 347
 348 /*-------------begin for surfEvalMachine -------------*/
 349 surfEvalMachine em_vertex;
 350 surfEvalMachine em_normal;
 351 surfEvalMachine em_color;
 352 surfEvalMachine em_texcoord;
 353
 354 int auto_normal_flag; //whether to output normla or not in callback
 355                       //determined by GL_AUTO_NORMAL and callback_auto_normal
 356 int callback_auto_normal; //GLU_CALLBACK_AUTO_NORMAL_EXT
 357 int vertex_flag;
 358 int normal_flag;
 359 int color_flag;
 360 int texcoord_flag;
 361
 362 void 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
 374 void inDoDomain2WithDerivsEM(surfEvalMachine *em, REAL u, REAL v,
 375                                 REAL *retPoint, REAL *retdu, REAL *retdv);
 376 void inDoDomain2EM(surfEvalMachine *em, REAL u, REAL v,
 377                                 REAL *retPoint);
 378  void inDoEvalCoord2EM(REAL u, REAL v);
 379
 380 void inBPMEvalEM(bezierPatchMesh* bpm);
 381 void inBPMListEvalEM(bezierPatchMesh* list);
 382
 383 /*-------------end for surfEvalMachine -------------*/
 384
 385
 386    /*************end for internal evaluators*****************/
 387
 388 };
 389
 390 inline 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_ */