Remove demos that have moved to git+ssh://git.freedesktop.org/git/mesa/demos.
[mesa.git] / progs / demos / isosurf.c
diff --git a/progs/demos/isosurf.c b/progs/demos/isosurf.c
deleted file mode 100644 (file)
index a5b21ff..0000000
+++ /dev/null
@@ -1,1146 +0,0 @@
-
-/*
- * Display an isosurface of 3-D wind speed volume.
- *
- * Command line options:
- *    -info      print GL implementation information
- *
- * Brian Paul  This file in public domain.
- */
-
-
-/* Keys:
- * =====
- *
- *   - Arrow keys to rotate
- *   - 's' toggles smooth shading
- *   - 'l' toggles lighting
- *   - 'f' toggles fog
- *   - 'I' and 'i' zoom in and out
- *   - 'c' toggles a user clip plane
- *   - 'm' toggles colorful materials in GL_TRIANGLES modes.
- *   - '+' and '-' move the user clip plane
- *
- * Other options are available via the popup menu.
- */
-
-#include <stdio.h>
-#include <string.h>
-#include <stdlib.h>
-#include <math.h>
-#ifdef _WIN32
-#include <windows.h>
-#undef CLIP_MASK
-#endif
-#include <GL/glew.h>
-#include "GL/glut.h"
-
-#include "readtex.h"
-#define TEXTURE_FILE "../images/reflect.rgb"
-
-#define LIT            0x00000001
-#define UNLIT          0x00000002
-#define REFLECT                0x00000004
-#define POINT_FILTER   0x00000008
-#define LINEAR_FILTER  0x00000010
-#define GLVERTEX       0x00000020
-#define DRAW_ELTS      0x00000040 
-#define DRAW_ARRAYS    0x00000080 
-#define ARRAY_ELT      0x00000100
-#define LOCKED         0x00000200
-#define UNLOCKED       0x00000400 
-#define IMMEDIATE      0x00000800
-#define DISPLAYLIST    0x00001000
-#define SHADE_SMOOTH   0x00002000
-#define SHADE_FLAT     0x00004000
-#define TRIANGLES      0x00008000
-#define STRIPS         0x00010000
-#define POINTS         0x00020000
-#define USER_CLIP      0x00040000
-#define NO_USER_CLIP   0x00080000
-#define MATERIALS      0x00100000
-#define NO_MATERIALS   0x00200000
-#define FOG            0x00400000
-#define NO_FOG         0x00800000
-#define QUIT           0x01000000
-#define GLINFO         0x02000000
-#define STIPPLE                0x04000000
-#define NO_STIPPLE     0x08000000
-#define POLYGON_FILL   0x10000000
-#define POLYGON_LINE   0x20000000
-#define POLYGON_POINT  0x40000000
-
-#define LIGHT_MASK             (LIT|UNLIT|REFLECT)
-#define FILTER_MASK            (POINT_FILTER|LINEAR_FILTER)
-#define RENDER_STYLE_MASK      (GLVERTEX|DRAW_ARRAYS|DRAW_ELTS|ARRAY_ELT)
-#define DLIST_MASK             (IMMEDIATE|DISPLAYLIST)
-#define LOCK_MASK              (LOCKED|UNLOCKED)
-#define MATERIAL_MASK          (MATERIALS|NO_MATERIALS)
-#define PRIMITIVE_MASK         (TRIANGLES|STRIPS|POINTS)
-#define CLIP_MASK              (USER_CLIP|NO_USER_CLIP)
-#define SHADE_MASK             (SHADE_SMOOTH|SHADE_FLAT)
-#define FOG_MASK               (FOG|NO_FOG)
-#define STIPPLE_MASK           (STIPPLE|NO_STIPPLE)
-#define POLYGON_MASK           (POLYGON_FILL|POLYGON_LINE|POLYGON_POINT)
-
-#define MAXVERTS 10000
-static GLint maxverts = MAXVERTS;
-static float data[MAXVERTS][6];
-static float compressed_data[MAXVERTS][6];
-static float expanded_data[MAXVERTS*3][6];
-static GLuint indices[MAXVERTS];
-static GLuint tri_indices[MAXVERTS*3];
-static GLuint strip_indices[MAXVERTS];
-static GLfloat col[100][4];
-static GLint numverts, num_tri_verts, numuniq;
-
-static GLfloat xrot;
-static GLfloat yrot;
-static GLfloat dist;
-static GLint state, allowed = ~0;
-static GLboolean doubleBuffer = GL_TRUE;
-static GLdouble plane[4];
-static GLuint surf1, dlist_state;
-
-static GLboolean PrintInfo = GL_FALSE;
-
-
-static GLubyte halftone[] = {
-   0xAA, 0xAA, 0xAA, 0xAA, 0x55, 0x55, 0x55, 0x55, 0xAA, 0xAA, 0xAA, 0xAA,
-   0x55, 0x55, 0x55, 0x55, 0xAA, 0xAA, 0xAA, 0xAA, 0x55, 0x55, 0x55, 0x55,
-   0xAA, 0xAA, 0xAA, 0xAA, 0x55, 0x55, 0x55, 0x55, 0xAA, 0xAA, 0xAA, 0xAA,
-   0x55, 0x55, 0x55, 0x55, 0xAA, 0xAA, 0xAA, 0xAA, 0x55, 0x55, 0x55, 0x55,
-   0xAA, 0xAA, 0xAA, 0xAA, 0x55, 0x55, 0x55, 0x55, 0xAA, 0xAA, 0xAA, 0xAA,
-   0x55, 0x55, 0x55, 0x55, 0xAA, 0xAA, 0xAA, 0xAA, 0x55, 0x55, 0x55, 0x55,
-   0xAA, 0xAA, 0xAA, 0xAA, 0x55, 0x55, 0x55, 0x55, 0xAA, 0xAA, 0xAA, 0xAA,
-   0x55, 0x55, 0x55, 0x55, 0xAA, 0xAA, 0xAA, 0xAA, 0x55, 0x55, 0x55, 0x55,
-   0xAA, 0xAA, 0xAA, 0xAA, 0x55, 0x55, 0x55, 0x55, 0xAA, 0xAA, 0xAA, 0xAA,
-   0x55, 0x55, 0x55, 0x55, 0xAA, 0xAA, 0xAA, 0xAA, 0x55, 0x55, 0x55, 0x55,
-   0xAA, 0xAA, 0xAA, 0xAA, 0x55, 0x55, 0x55, 0x55};
-
-
-static void read_surface( char *filename )
-{
-   FILE *f;
-
-   f = fopen(filename,"r");
-   if (!f) {
-      printf("couldn't read %s\n", filename);
-      exit(1);
-   }
-
-   numverts = 0;
-   while (!feof(f) && numverts<maxverts) {
-      int result;
-      result = fscanf( f, "%f %f %f  %f %f %f",
-                      &data[numverts][0], &data[numverts][1], &data[numverts][2],
-                      &data[numverts][3], &data[numverts][4], &data[numverts][5] );
-      (void) result;
-      numverts++;
-   }
-   numverts--;
-
-   printf("%d vertices, %d triangles\n", numverts, numverts-2);
-   fclose(f);
-}
-
-
-
-static void print_flags( const char *msg, GLuint flags ) 
-{
-   fprintf(stderr, 
-          "%s (0x%x): %s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s\n",
-          msg, flags,
-          (flags & GLVERTEX) ? "glVertex, " : "",
-          (flags & DRAW_ARRAYS) ? "glDrawArrays, " : "",
-          (flags & DRAW_ELTS) ? "glDrawElements, " : "",
-          (flags & ARRAY_ELT) ? "glArrayElement, " : "",
-          (flags & LOCKED) ? "locked arrays, " : "",
-          (flags & TRIANGLES) ? "GL_TRIANGLES, " : "",
-          (flags & STRIPS) ? "GL_TRIANGLE_STRIP, " : "",
-          (flags & POINTS) ? "GL_POINTS, " : "",
-          (flags & DISPLAYLIST) ? "as a displaylist, " : "",
-          (flags & LIT) ? "lit, " : "",
-          (flags & UNLIT) ? "unlit, " : "",
-          (flags & REFLECT) ? "reflect, " : "",
-          (flags & SHADE_FLAT) ? "flat-shaded, " : "",
-          (flags & USER_CLIP) ? "user_clip, " : "",
-          (flags & MATERIALS) ? "materials, " : "",
-          (flags & FOG) ? "fog, " : "",
-          (flags & STIPPLE) ? "stipple, " : "",
-          (flags & POLYGON_LINE) ? "polygon mode line, " : "",
-          (flags & POLYGON_POINT) ? "polygon mode point, " : "");
-}
-
-
-
-struct data_idx {
-   float *data;
-   int idx;
-   int uniq_idx;
-};
-
-
-#define COMPARE_FUNC( AXIS )                            \
-static int compare_axis_##AXIS( const void *a, const void *b ) \
-{                                                      \
-   float t = ( (*(struct data_idx *)a).data[AXIS] -    \
-              (*(struct data_idx *)b).data[AXIS] );    \
-                                                       \
-   if (t < 0) return -1;                               \
-   if (t > 0) return 1;                                        \
-   return 0;                                           \
-}
-
-COMPARE_FUNC(0)
-COMPARE_FUNC(1)
-COMPARE_FUNC(2)
-COMPARE_FUNC(3)
-COMPARE_FUNC(4)
-COMPARE_FUNC(5)
-COMPARE_FUNC(6)
-
-int (*(compare[7]))( const void *a, const void *b ) =
-{
-   compare_axis_0,
-   compare_axis_1,
-   compare_axis_2,
-   compare_axis_3,
-   compare_axis_4,
-   compare_axis_5,
-   compare_axis_6,
-};
-
-
-#define VEC_ELT(f, s, i)  (float *)(((char *)f) + s * i)
-
-static int sort_axis( int axis,
-                     int vec_size,
-                     int vec_stride,
-                     struct data_idx *indices,
-                     int start,
-                     int finish,
-                     float *out,
-                     int uniq,
-                     const float fudge )
-{
-   int i;
-
-   if (finish-start > 2)
-   {
-      qsort( indices+start, finish-start, sizeof(*indices), compare[axis] );
-   }
-   else if (indices[start].data[axis] > indices[start+1].data[axis])
-   {
-      struct data_idx tmp = indices[start];
-      indices[start] = indices[start+1];
-      indices[start+1] = tmp;
-   }
-
-   if (axis == vec_size-1) {
-      for (i = start ; i < finish ; ) {
-        float max = indices[i].data[axis] + fudge;
-        float *dest = VEC_ELT(out, vec_stride, uniq);
-        int j;
-
-        for (j = 0 ; j < vec_size ; j++)
-           dest[j] = indices[i].data[j];
-
-        for ( ; i < finish && max >= indices[i].data[axis]; i++)
-           indices[i].uniq_idx = uniq;
-
-        uniq++;
-      }
-   } else {
-      for (i = start ; i < finish ; ) {
-        int j = i + 1;
-        float max = indices[i].data[axis] + fudge;
-        while (j < finish && max >= indices[j].data[axis]) j++;
-        if (j == i+1) {
-           float *dest = VEC_ELT(out, vec_stride, uniq);
-           int k;
-
-           indices[i].uniq_idx = uniq;
-
-           for (k = 0 ; k < vec_size ; k++)
-              dest[k] = indices[i].data[k];
-
-           uniq++;
-        } else {
-           uniq = sort_axis( axis+1, vec_size, vec_stride,
-                             indices, i, j, out, uniq, fudge );
-        }
-        i = j;
-      }
-   }
-
-   return uniq;
-}
-
-
-static void extract_indices1( const struct data_idx *in, unsigned int *out,
-                             int n )
-{
-   int i;
-   for ( i = 0 ; i < n ; i++ ) {
-      out[in[i].idx] = in[i].uniq_idx;
-   }
-}
-
-
-static void compactify_arrays(void)
-{
-   int i;
-   struct data_idx *ind;
-
-   ind = (struct data_idx *) malloc( sizeof(struct data_idx) * numverts );
-
-   for (i = 0 ; i < numverts ; i++) {
-      ind[i].idx = i;
-      ind[i].data = data[i];
-   }
-
-   numuniq = sort_axis(0,
-                      sizeof(compressed_data[0])/sizeof(float),
-                      sizeof(compressed_data[0]),
-                      ind,
-                      0,
-                      numverts,
-                      (float *)compressed_data,
-                      0,
-                      1e-6);
-
-   printf("Nr unique vertex/normal pairs: %d\n", numuniq);
-
-   extract_indices1( ind, indices, numverts );
-   free( ind );
-}
-
-static void expand_arrays(void) 
-{
-   int i;
-   int parity = 0;
-   for (i = 2 ; i < numverts ; i++, parity ^= 1) {
-      int v0 = i-2+parity;
-      int v1 = i-1-parity;
-      int v2 = i;
-      memcpy( expanded_data[(i-2)*3+0], data[v0], sizeof(data[0]) );
-      memcpy( expanded_data[(i-2)*3+1], data[v1], sizeof(data[0]) );
-      memcpy( expanded_data[(i-2)*3+2], data[v2], sizeof(data[0]) );
-   }
-}
-
-static float myrand( float max )
-{
-   return max*rand()/(RAND_MAX+1.0);
-}
-
-
-static void make_tri_indices( void )
-{
-   unsigned int *v = tri_indices;
-   unsigned int parity = 0;
-   int i, j;
-
-   for (j=2;j<numverts;j++,parity^=1) {
-      if (parity) {
-        *v++ = indices[j-1];
-        *v++ = indices[j-2];
-        *v++ = indices[j];
-      } else {
-        *v++ = indices[j-2];
-        *v++ = indices[j-1];
-        *v++ = indices[j];
-      }
-   }
-
-   num_tri_verts = v - tri_indices;
-   printf("num_tri_verts: %d\n", num_tri_verts);
-
-   for (i = j = 0 ; i < num_tri_verts ; i += 600, j++) {
-      col[j][3] = 1;
-      col[j][2] = myrand(1);
-      col[j][1] = myrand(1);
-      col[j][0] = myrand(1);
-   }
-
-   for (i = 0; i < numverts ; i++)
-      strip_indices[i] = i;
-}
-
-#define MIN(x,y) (x < y) ? x : y
-
-static void draw_surface( unsigned int with_state )
-{
-   GLint i, j;
-   
-   if (with_state & DISPLAYLIST) {
-      if ((with_state & (RENDER_STYLE_MASK|PRIMITIVE_MASK|MATERIAL_MASK)) != 
-         dlist_state) {
-        /* 
-         */
-        fprintf(stderr, "rebuilding displaylist\n");
-
-        if (dlist_state)
-           glDeleteLists( surf1, 1 );
-
-        dlist_state = with_state & (RENDER_STYLE_MASK|PRIMITIVE_MASK|
-                                    MATERIAL_MASK);
-        surf1 = glGenLists(1);
-        glNewList(surf1, GL_COMPILE);
-        draw_surface( dlist_state );
-        glEndList();
-      }
-
-      glCallList( surf1 );
-      return;
-   }
-
-   switch (with_state & (RENDER_STYLE_MASK|PRIMITIVE_MASK)) {
-#ifdef GL_EXT_vertex_array
-
-   case (DRAW_ELTS|TRIANGLES):
-      if (with_state & MATERIALS) {
-        for (j = i = 0 ; i < num_tri_verts ; i += 600, j++) {
-           GLuint nr = MIN(num_tri_verts-i, 600);
-           glMaterialfv(GL_FRONT_AND_BACK, GL_SPECULAR, col[j]);
-           glMaterialfv(GL_FRONT_AND_BACK, GL_DIFFUSE, col[j]);
-           glDrawElements( GL_TRIANGLES, nr, GL_UNSIGNED_INT, tri_indices+i );
-        }
-      } else {
-        glDrawElements( GL_TRIANGLES, num_tri_verts, GL_UNSIGNED_INT,
-                        tri_indices );
-      }
-      break;
-
-   case (DRAW_ARRAYS|TRIANGLES):
-      glDrawArraysEXT( GL_TRIANGLES, 0, (numverts-2)*3 );
-      break;
-
-   case (ARRAY_ELT|TRIANGLES):
-      if (with_state & MATERIALS) {
-        for (j = i = 0 ; i < num_tri_verts ; i += 600, j++) {
-           GLuint nr = MIN(num_tri_verts-i, 600);
-           GLuint k;
-           glMaterialfv(GL_FRONT_AND_BACK, GL_SPECULAR, col[j]);
-           glMaterialfv(GL_FRONT_AND_BACK, GL_DIFFUSE, col[j]);
-           glBegin( GL_TRIANGLES );
-           for (k = 0 ; k < nr ; k++)
-              glArrayElement( tri_indices[i+k] );
-           glEnd();
-        }
-      } else {
-        glBegin( GL_TRIANGLES );
-        for (i = 0 ; i < num_tri_verts ; i++)
-           glArrayElement( tri_indices[i] );
-
-        glEnd();
-      }
-      break;
-
-
-      /* Uses the original arrays (including duplicate elements):
-       */
-   case (DRAW_ARRAYS|STRIPS):
-      glDrawArraysEXT( GL_TRIANGLE_STRIP, 0, numverts );
-      break;
-   case (DRAW_ELTS|STRIPS):
-      glDrawElements( GL_TRIANGLE_STRIP, numverts,
-                     GL_UNSIGNED_INT, strip_indices );
-      break;
-
-      /* Uses the original arrays (including duplicate elements):
-       */
-   case (ARRAY_ELT|STRIPS):
-      glBegin( GL_TRIANGLE_STRIP );
-      for (i = 0 ; i < numverts ; i++)
-        glArrayElement( i );
-      glEnd();
-      break;
-
-   case (DRAW_ARRAYS|POINTS):
-      glDrawArraysEXT( GL_POINTS, 0, numuniq );
-      break;
-   case (DRAW_ELTS|POINTS):
-      /* can use numuniq with strip_indices as strip_indices[i] == i.
-       */
-      glDrawElements( GL_POINTS, numuniq, 
-                     GL_UNSIGNED_INT, strip_indices );
-      break;
-   case (ARRAY_ELT|POINTS):
-      /* just emit each unique element once:
-       */
-      glBegin( GL_POINTS );
-      for (i = 0 ; i < numuniq ; i++)
-        glArrayElement( i );
-      glEnd();
-      break;
-#endif
-
-   case (GLVERTEX|TRIANGLES):
-      if (with_state & MATERIALS) {
-        for (j = i = 0 ; i < num_tri_verts ; i += 600, j++) {
-           GLuint nr = MIN(num_tri_verts-i, 600);
-           GLuint k;
-           glMaterialfv(GL_FRONT_AND_BACK, GL_SPECULAR, col[j]);
-           glMaterialfv(GL_FRONT_AND_BACK, GL_DIFFUSE, col[j]);
-           glBegin( GL_TRIANGLES );
-           for (k = 0 ; k < nr ; k++) {
-              glNormal3fv( &compressed_data[tri_indices[i+k]][3] );
-              glVertex3fv( &compressed_data[tri_indices[i+k]][0] );
-           }
-           glEnd();
-        }
-      } else {
-        glBegin( GL_TRIANGLES );
-        for (i = 0 ; i < num_tri_verts ; i++) {
-           glNormal3fv( &compressed_data[tri_indices[i]][3] );
-           glVertex3fv( &compressed_data[tri_indices[i]][0] );
-        }
-        glEnd();
-      }
-      break;
-
-   case (GLVERTEX|POINTS):
-      /* Renders all points, but not in strip order...  Shouldn't be a
-       * problem, but people may be confused as to why points are so
-       * much faster in this demo...  And why cva doesn't help them...
-       */
-      glBegin( GL_POINTS );
-      for ( i = 0 ; i < numuniq ; i++ ) {
-         glNormal3fv( &compressed_data[i][3] );
-         glVertex3fv( &compressed_data[i][0] );
-      }
-      glEnd();
-      break;
-
-   case (GLVERTEX|STRIPS):
-      if (with_state & MATERIALS) {
-         glBegin( GL_TRIANGLE_STRIP );
-         for (i=0;i<numverts;i++) {
-            if (i % 600 == 0 && i != 0) {
-               unsigned j = i / 600;
-               glMaterialfv(GL_FRONT_AND_BACK, GL_SPECULAR, col[j]);
-               glMaterialfv(GL_FRONT_AND_BACK, GL_DIFFUSE, col[j]);
-            }
-            glNormal3fv( &data[i][3] );
-            glVertex3fv( &data[i][0] );
-         }
-         glEnd();
-      }
-      else {
-         glBegin( GL_TRIANGLE_STRIP );
-         for (i=0;i<numverts;i++) {
-            glNormal3fv( &data[i][3] );
-            glVertex3fv( &data[i][0] );
-         }
-         glEnd();
-      }
-      break;
-
-   default:
-      fprintf(stderr, "unimplemented mode %x...\n", 
-             (with_state & (RENDER_STYLE_MASK|PRIMITIVE_MASK)));
-      break;
-   }
-}
-
-
-
-static void Display(void)
-{
-    glClear( GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT );
-    draw_surface( state );
-    glFlush();
-    if (doubleBuffer) glutSwapBuffers();
-}
-
-
-/* KW: only do this when necessary, so CVA can re-use results.
- */
-static void set_matrix( void )
-{
-   glMatrixMode(GL_MODELVIEW);
-   glLoadIdentity();
-   glTranslatef( 0.0, 0.0, dist );
-   glRotatef( yrot, 0.0, 1.0, 0.0 );
-   glRotatef( xrot, 1.0, 0.0, 0.0 );
-}
-
-static void Benchmark( float xdiff, float ydiff )
-{
-   int startTime, endTime;
-   int draws;
-   double seconds, fps, triPerSecond;
-
-   printf("Benchmarking...\n");
-
-   draws = 0;
-   startTime = glutGet(GLUT_ELAPSED_TIME);
-   xrot = 0.0;
-   do {
-      xrot += xdiff;
-      yrot += ydiff;
-      set_matrix();
-      Display();
-      draws++;
-      endTime = glutGet(GLUT_ELAPSED_TIME);
-   } while (endTime - startTime < 5000);   /* 5 seconds */
-
-   /* Results */
-   seconds = (double) (endTime - startTime) / 1000.0;
-   triPerSecond = (numverts - 2) * draws / seconds;
-   fps = draws / seconds;
-   printf("Result:  triangles/sec: %g  fps: %g\n", triPerSecond, fps);
-}
-
-
-static void InitMaterials(void)
-{
-    static float ambient[] = {0.1, 0.1, 0.1, 1.0};
-    static float diffuse[] = {0.5, 1.0, 1.0, 1.0};
-    static float position0[] = {0.0, 0.0, 20.0, 0.0};
-    static float position1[] = {0.0, 0.0, -20.0, 0.0};
-    static float front_mat_shininess[] = {60.0};
-    static float front_mat_specular[] = {0.2, 0.2, 0.2, 1.0};
-    static float front_mat_diffuse[] = {0.5, 0.28, 0.38, 1.0};
-    /*
-    static float back_mat_shininess[] = {60.0};
-    static float back_mat_specular[] = {0.5, 0.5, 0.2, 1.0};
-    static float back_mat_diffuse[] = {1.0, 1.0, 0.2, 1.0};
-    */
-    static float lmodel_ambient[] = {1.0, 1.0, 1.0, 1.0};
-    static float lmodel_twoside[] = {GL_FALSE};
-
-    glLightfv(GL_LIGHT0, GL_AMBIENT, ambient);
-    glLightfv(GL_LIGHT0, GL_DIFFUSE, diffuse);
-    glLightfv(GL_LIGHT0, GL_POSITION, position0);
-    glEnable(GL_LIGHT0);
-
-    glLightfv(GL_LIGHT1, GL_AMBIENT, ambient);
-    glLightfv(GL_LIGHT1, GL_DIFFUSE, diffuse);
-    glLightfv(GL_LIGHT1, GL_POSITION, position1);
-    glEnable(GL_LIGHT1);
-
-    glLightModelfv(GL_LIGHT_MODEL_AMBIENT, lmodel_ambient);
-    glLightModelfv(GL_LIGHT_MODEL_TWO_SIDE, lmodel_twoside);
-
-    glMaterialfv(GL_FRONT_AND_BACK, GL_SHININESS, front_mat_shininess);
-    glMaterialfv(GL_FRONT_AND_BACK, GL_SPECULAR, front_mat_specular);
-    glMaterialfv(GL_FRONT_AND_BACK, GL_DIFFUSE, front_mat_diffuse);
-
-    glPolygonStipple (halftone);
-}
-
-
-
-#define UPDATE(o,n,mask) (o&=~mask, o|=n&mask)
-#define CHANGED(o,n,mask) ((n&mask) && (n&mask) != (o&mask) )
-
-static void ModeMenu(int m)
-{
-   m &= allowed;
-
-   if (!m) return;
-
-   if (m==QUIT)
-      exit(0);
-
-   if (m==GLINFO) {
-      printf("GL_VERSION: %s\n", (char *) glGetString(GL_VERSION));
-      printf("GL_EXTENSIONS: %s\n", (char *) glGetString(GL_EXTENSIONS));
-      printf("GL_RENDERER: %s\n", (char *) glGetString(GL_RENDERER));
-      return;
-   }
-
-   if (CHANGED(state, m, FILTER_MASK)) {
-      UPDATE(state, m, FILTER_MASK);
-      if (m & LINEAR_FILTER) {
-        glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_LINEAR);
-        glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_LINEAR);
-      } else {
-        glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);
-        glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);
-      }
-   }
-
-   if (CHANGED(state, m, LIGHT_MASK)) {
-      UPDATE(state, m, LIGHT_MASK);
-      if (m & LIT) {
-        glEnable(GL_LIGHTING);
-        glDisable(GL_TEXTURE_GEN_S);
-        glDisable(GL_TEXTURE_GEN_T);
-        glDisable(GL_TEXTURE_2D);
-      }
-      else if (m & UNLIT) {
-        glDisable(GL_LIGHTING);
-        glDisable(GL_TEXTURE_GEN_S);
-        glDisable(GL_TEXTURE_GEN_T);
-        glDisable(GL_TEXTURE_2D);
-      }
-      else if (m & REFLECT) {
-        glDisable(GL_LIGHTING);
-        glEnable(GL_TEXTURE_GEN_S);
-        glEnable(GL_TEXTURE_GEN_T);
-        glEnable(GL_TEXTURE_2D);
-      }
-   }
-
-   if (CHANGED(state, m, SHADE_MASK)) {
-      UPDATE(state, m, SHADE_MASK);
-      if (m & SHADE_SMOOTH)
-        glShadeModel(GL_SMOOTH);
-      else
-        glShadeModel(GL_FLAT);
-   }
-
-
-   if (CHANGED(state, m, CLIP_MASK)) {
-      UPDATE(state, m, CLIP_MASK);
-      if (m & USER_CLIP) {
-        glEnable(GL_CLIP_PLANE0);
-      } else {
-        glDisable(GL_CLIP_PLANE0);
-      }
-   }
-
-   if (CHANGED(state, m, FOG_MASK)) {
-      UPDATE(state, m, FOG_MASK);
-      if (m & FOG) {
-        glEnable(GL_FOG);
-      }
-      else {
-        glDisable(GL_FOG);
-      }
-   }
-
-   if (CHANGED(state, m, STIPPLE_MASK)) {
-      UPDATE(state, m, STIPPLE_MASK);
-      if (m & STIPPLE) {
-        glEnable(GL_POLYGON_STIPPLE);
-      }
-      else {
-        glDisable(GL_POLYGON_STIPPLE);
-      }
-   }
-
-   if (CHANGED(state, m, POLYGON_MASK)) {
-      UPDATE(state, m, POLYGON_MASK);
-      if (m & POLYGON_FILL) {
-        glPolygonMode(GL_FRONT_AND_BACK, GL_FILL);
-      }
-      else if (m & POLYGON_LINE) {
-        glPolygonMode(GL_FRONT_AND_BACK, GL_LINE);
-      }
-      else {
-        glPolygonMode(GL_FRONT_AND_BACK, GL_POINT);
-      }
-   }
-
-#ifdef GL_EXT_vertex_array
-   if (CHANGED(state, m, (LOCK_MASK|RENDER_STYLE_MASK|PRIMITIVE_MASK)))
-   {
-      if (m & (PRIMITIVE_MASK)) {
-        UPDATE(state, m, (PRIMITIVE_MASK));
-      }
-
-      if (m & (RENDER_STYLE_MASK)) {
-        UPDATE(state, m, (RENDER_STYLE_MASK));
-      }
-
-      if (m & LOCK_MASK) {
-        UPDATE(state, m, (LOCK_MASK));
-      }
-
-
-      print_flags("primitive", state & PRIMITIVE_MASK);
-      print_flags("render style", state & RENDER_STYLE_MASK);
-
-      if ((state & PRIMITIVE_MASK) != STRIPS &&
-         ((state & RENDER_STYLE_MASK) == DRAW_ELTS ||
-          (state & RENDER_STYLE_MASK) == ARRAY_ELT || 
-          (state & PRIMITIVE_MASK) == POINTS))
-      {
-        fprintf(stderr, "enabling small arrays\n");
-        /* Rendering any primitive with draw-element/array-element
-         *  --> Can't do strips here as ordering has been lost in
-         *  compaction process...
-         */
-        glVertexPointerEXT( 3, GL_FLOAT, sizeof(data[0]), numuniq,
-                            compressed_data );
-        glNormalPointerEXT( GL_FLOAT, sizeof(data[0]), numuniq,
-                            &compressed_data[0][3]);
-#ifdef GL_EXT_compiled_vertex_array
-        if (allowed & LOCKED) {
-           if (state & LOCKED) {
-              glLockArraysEXT( 0, numuniq );
-           } else {
-              glUnlockArraysEXT();
-           }
-        }
-#endif
-      }
-      else if ((state & PRIMITIVE_MASK) == TRIANGLES &&
-              (state & RENDER_STYLE_MASK) == DRAW_ARRAYS) {
-        fprintf(stderr, "enabling big arrays\n");
-        /* Only get here for TRIANGLES and drawarrays
-         */
-        glVertexPointerEXT( 3, GL_FLOAT, sizeof(data[0]), (numverts-2) * 3,
-                            expanded_data );
-        glNormalPointerEXT( GL_FLOAT, sizeof(data[0]), (numverts-2) * 3,
-                            &expanded_data[0][3]);
-
-#ifdef GL_EXT_compiled_vertex_array
-        if (allowed & LOCKED) {
-           if (state & LOCKED) {
-              glLockArraysEXT( 0, (numverts-2)*3 );
-           } else {
-              glUnlockArraysEXT();
-           }
-        }
-#endif
-      }
-      else {
-        fprintf(stderr, "enabling normal arrays\n");
-        glVertexPointerEXT( 3, GL_FLOAT, sizeof(data[0]), numverts, data );
-        glNormalPointerEXT( GL_FLOAT, sizeof(data[0]), numverts, &data[0][3]);
-#ifdef GL_EXT_compiled_vertex_array
-        if (allowed & LOCKED) {
-           if (state & LOCKED) {
-              glLockArraysEXT( 0, numverts );
-           } else {
-              glUnlockArraysEXT();
-           }
-        }
-#endif
-      }
-
-   }
-#endif
-
-
-   if (m & DLIST_MASK) {
-      UPDATE(state, m, DLIST_MASK);
-   }
-
-   if (m & MATERIAL_MASK) {
-      UPDATE(state, m, MATERIAL_MASK);
-   }
-
-   print_flags("new flags", state);
-
-   glutPostRedisplay();
-}
-
-
-
-static void Init(int argc, char *argv[])
-{
-   GLfloat fogColor[4] = {0.5,1.0,0.5,1.0};
-
-   xrot = 0;
-   yrot = 0;
-   dist = -6;
-   plane[0] = 1.0;
-   plane[1] = 0.0;
-   plane[2] = -1.0;
-   plane[3] = 0.0;
-
-   glClearColor(0.0, 0.0, 1.0, 0.0);
-   glEnable( GL_DEPTH_TEST );
-   glEnableClientState( GL_VERTEX_ARRAY );
-   glEnableClientState( GL_NORMAL_ARRAY );
-
-   glMatrixMode(GL_PROJECTION);
-   glLoadIdentity();
-   glFrustum( -1.0, 1.0, -1.0, 1.0, 5, 25 );
-
-   glMatrixMode(GL_MODELVIEW);
-   glLoadIdentity();
-   glClipPlane(GL_CLIP_PLANE0, plane);
-
-   InitMaterials();
-
-   set_matrix();
-
-   glTexEnvi(GL_TEXTURE_ENV, GL_TEXTURE_ENV_MODE, GL_REPLACE);
-   glHint(GL_PERSPECTIVE_CORRECTION_HINT, GL_FASTEST);
-
-   glTexGeni(GL_S, GL_TEXTURE_GEN_MODE, GL_SPHERE_MAP);
-   glTexGeni(GL_T, GL_TEXTURE_GEN_MODE, GL_SPHERE_MAP);
-
-
-   /* Green fog is easy to see */
-   glFogi(GL_FOG_MODE,GL_EXP2);
-   glFogfv(GL_FOG_COLOR,fogColor);
-   glFogf(GL_FOG_DENSITY,0.15);
-   glHint(GL_FOG_HINT,GL_DONT_CARE);
-
-   {
-      static int firsttime = 1;
-      if (firsttime) {
-        firsttime = 0;
-        compactify_arrays();
-        expand_arrays();
-        make_tri_indices();
-
-        if (!LoadRGBMipmaps(TEXTURE_FILE, GL_RGB)) {
-           printf("Error: couldn't load texture image\n");
-           exit(1);
-        }
-      }
-   }
-
-   ModeMenu(SHADE_SMOOTH|
-           LIT|
-           POINT_FILTER|
-           NO_USER_CLIP|
-           NO_MATERIALS|
-           NO_FOG|
-           NO_STIPPLE|
-           IMMEDIATE|
-           STRIPS|
-           UNLOCKED|
-           GLVERTEX);
-
-   if (PrintInfo) {
-      printf("GL_RENDERER   = %s\n", (char *) glGetString(GL_RENDERER));
-      printf("GL_VERSION    = %s\n", (char *) glGetString(GL_VERSION));
-      printf("GL_VENDOR     = %s\n", (char *) glGetString(GL_VENDOR));
-      printf("GL_EXTENSIONS = %s\n", (char *) glGetString(GL_EXTENSIONS));
-   }
-}
-
-
-
-static void Reshape(int width, int height)
-{
-    glViewport(0, 0, (GLint)width, (GLint)height);
-}
-
-
-
-static void Key( unsigned char key, int x, int y )
-{
-   (void) x;
-   (void) y;
-   switch (key) {
-   case 27:
-      exit(0);
-   case 'f':
-      ModeMenu((state ^ FOG_MASK) & FOG_MASK);
-      break;
-   case 's':
-      ModeMenu((state ^ SHADE_MASK) & SHADE_MASK);
-      break;
-   case 't':
-      ModeMenu((state ^ STIPPLE_MASK) & STIPPLE_MASK);
-      break;
-   case 'l':
-      ModeMenu((state ^ LIGHT_MASK) & (LIT|UNLIT));
-      break;
-   case 'm':
-      ModeMenu((state ^ MATERIAL_MASK) & MATERIAL_MASK);
-      break;
-   case 'c':
-      ModeMenu((state ^ CLIP_MASK) & CLIP_MASK);
-      break;
-   case 'v':
-      ModeMenu((LOCKED|IMMEDIATE|DRAW_ELTS|TRIANGLES) & allowed);
-      break;
-   case 'V':
-      ModeMenu(UNLOCKED|IMMEDIATE|GLVERTEX|STRIPS);
-      break;
-   case 'b':
-      Benchmark(5.0, 0);
-      break;
-   case 'B':
-      Benchmark(0, 5.0);
-      break;
-   case 'i':
-      dist += .25;
-      set_matrix();
-      glutPostRedisplay();
-      break;
-   case 'I':
-      dist -= .25;
-      set_matrix();
-      glutPostRedisplay();
-      break;
-   case '-':
-   case '_':
-      plane[3] += 2.0;
-      glMatrixMode(GL_MODELVIEW);
-      glLoadIdentity();
-      glClipPlane(GL_CLIP_PLANE0, plane);
-      set_matrix();
-      glutPostRedisplay();
-      break;
-   case '+':
-   case '=':
-      plane[3] -= 2.0;
-      glMatrixMode(GL_MODELVIEW);
-      glLoadIdentity();
-      glClipPlane(GL_CLIP_PLANE0, plane);
-      set_matrix();
-      glutPostRedisplay();
-      break;
-   case ' ':
-      Init(0,0);
-      break;
-   }
-}
-
-
-static void SpecialKey( int key, int x, int y )
-{
-   (void) x;
-   (void) y;
-   switch (key) {
-   case GLUT_KEY_LEFT:
-      yrot -= 15.0;
-      break;
-   case GLUT_KEY_RIGHT:
-      yrot += 15.0;
-      break;
-   case GLUT_KEY_UP:
-      xrot += 15.0;
-      break;
-   case GLUT_KEY_DOWN:
-      xrot -= 15.0;
-      break;
-   default:
-      return;
-   }
-   set_matrix();
-   glutPostRedisplay();
-}
-
-
-
-static GLint Args(int argc, char **argv)
-{
-   GLint i;
-   GLint mode = 0;
-
-   for (i = 1; i < argc; i++) {
-      if (strcmp(argv[i], "-sb") == 0) {
-         doubleBuffer = GL_FALSE;
-      }
-      else if (strcmp(argv[i], "-db") == 0) {
-         doubleBuffer = GL_TRUE;
-      }
-      else if (strcmp(argv[i], "-info") == 0) {
-         PrintInfo = GL_TRUE;
-      }
-      else if (strcmp(argv[i], "-10") == 0) {
-         maxverts = 10;
-      }
-      else if (strcmp(argv[i], "-100") == 0) {
-        maxverts = 100;
-      }
-      else if (strcmp(argv[i], "-1000") == 0) {
-        maxverts = 1000;
-      }
-      else {
-         printf("%s (Bad option).\n", argv[i]);
-        return QUIT;
-      }
-   }
-
-   return mode;
-}
-
-int main(int argc, char **argv)
-{
-   GLenum type;
-
-   GLuint arg_mode = Args(argc, argv);
-
-   if (arg_mode & QUIT)
-      exit(0);
-
-   read_surface( "isosurf.dat" );
-
-   glutInitWindowSize(400, 400);
-   glutInit( &argc, argv);
-
-   type = GLUT_DEPTH;
-   type |= GLUT_RGB;
-   type |= (doubleBuffer) ? GLUT_DOUBLE : GLUT_SINGLE;
-   glutInitDisplayMode(type);
-
-   if (glutCreateWindow("Isosurface") <= 0) {
-      exit(0);
-   }
-
-   glewInit();
-
-   /* Make sure server supports the vertex array extension */
-   if (!GLEW_EXT_vertex_array)
-   {
-      printf("Vertex arrays not supported by this renderer\n");
-      allowed &= ~(LOCKED|DRAW_ARRAYS|DRAW_ELTS|ARRAY_ELT);
-   }
-   else if (!GLEW_EXT_compiled_vertex_array)
-   {
-      printf("Compiled vertex arrays not supported by this renderer\n");
-      allowed &= ~LOCKED;
-   }
-
-   Init(argc, argv);
-   ModeMenu(arg_mode);
-
-   glutCreateMenu(ModeMenu);
-   glutAddMenuEntry("GL info",               GLINFO);
-   glutAddMenuEntry("", 0);
-   glutAddMenuEntry("Lit",                   LIT);
-   glutAddMenuEntry("Unlit",                 UNLIT);
-   glutAddMenuEntry("Reflect",               REFLECT);
-   glutAddMenuEntry("", 0);
-   glutAddMenuEntry("Smooth",                SHADE_SMOOTH);
-   glutAddMenuEntry("Flat",                  SHADE_FLAT);
-   glutAddMenuEntry("", 0);
-   glutAddMenuEntry("Fog",                   FOG);
-   glutAddMenuEntry("No Fog",                NO_FOG);
-   glutAddMenuEntry("", 0);
-   glutAddMenuEntry("Stipple",               STIPPLE);
-   glutAddMenuEntry("No Stipple",            NO_STIPPLE);
-   glutAddMenuEntry("", 0);
-   glutAddMenuEntry("Polygon Mode Fill",     POLYGON_FILL);
-   glutAddMenuEntry("Polygon Mode Line",     POLYGON_LINE);
-   glutAddMenuEntry("Polygon Mode Points",   POLYGON_POINT);
-   glutAddMenuEntry("", 0);
-   glutAddMenuEntry("Point Filtered",        POINT_FILTER);
-   glutAddMenuEntry("Linear Filtered",       LINEAR_FILTER);
-   glutAddMenuEntry("", 0);
-   glutAddMenuEntry("GL_TRIANGLES",          TRIANGLES);
-   glutAddMenuEntry("GL_TRIANGLE_STRIPS",    STRIPS);
-   glutAddMenuEntry("GL_POINTS",             POINTS);
-   glutAddMenuEntry("", 0);
-   glutAddMenuEntry("Displaylist",           DISPLAYLIST);
-   glutAddMenuEntry("Immediate",             IMMEDIATE);
-   glutAddMenuEntry("", 0);
-   if (allowed & LOCKED) {
-      glutAddMenuEntry("Locked Arrays (CVA)", LOCKED);
-      glutAddMenuEntry("Unlocked Arrays",     UNLOCKED);
-      glutAddMenuEntry("", 0);
-   }
-   glutAddMenuEntry("glVertex",               GLVERTEX);
-   if (allowed & DRAW_ARRAYS) {
-      glutAddMenuEntry("glDrawElements",      DRAW_ELTS);
-      glutAddMenuEntry("glDrawArrays",       DRAW_ARRAYS);
-      glutAddMenuEntry("glArrayElement",      ARRAY_ELT);
-   }
-   glutAddMenuEntry("", 0);
-   glutAddMenuEntry("Quit",                   QUIT);
-   glutAttachMenu(GLUT_RIGHT_BUTTON);
-
-   glutReshapeFunc(Reshape);
-   glutKeyboardFunc(Key);
-   glutSpecialFunc(SpecialKey);
-   glutDisplayFunc(Display);
-
-   glutMainLoop();
-   return 0;
-}