halomagic.c
上传用户:xk288cn
上传日期:2007-05-28
资源大小:4876k
文件大小:31k
源码类别:

GIS编程

开发平台:

Visual C++

  1. /* Copyright (c) Mark J. Kilgard, 1994, 1997.  */
  2. /* This program is freely distributable without licensing fees 
  3.    and is provided without guarantee or warrantee expressed or 
  4.    implied. This program is -not- in the public domain. */
  5. /* dinoshade.c with an added "magic halo" effect when you hit the
  6.    space bar.  You can use an overlaid or blended halo (blending is
  7.    slower). */
  8. /* Example for PC game developers to show how to *combine* texturing,
  9.    reflections, and projected shadows all in real-time with OpenGL.
  10.    Robust reflections use stenciling.  Robust projected shadows
  11.    use both stenciling and polygon offset.  PC game programmers
  12.    should realize that neither stenciling nor polygon offset are 
  13.    supported by Direct3D, so these real-time rendering algorithms
  14.    are only really viable with OpenGL. 
  15.    
  16.    The program has modes for disabling the stenciling and polygon
  17.    offset uses.  It is worth running this example with these features
  18.    toggled off so you can see the sort of artifacts that result.
  19.    
  20.    Notice that the floor texturing, reflections, and shadowing
  21.    all co-exist properly. */
  22. /* When you run this program:  Left mouse button controls the
  23.    view.  Middle mouse button controls light position (left &
  24.    right rotates light around dino; up & down moves light
  25.    position up and down).  Right mouse button pops up menu. */
  26. /* Check out the comments in the "redraw" routine to see how the
  27.    reflection blending and surface stenciling is done.  You can
  28.    also see in "redraw" how the projected shadows are rendered,
  29.    including the use of stenciling and polygon offset. */
  30. /* This program is derived from glutdino.c */
  31. /* Compile: cc -o halomagic halomagic.c -lglut -lGLU -lGL -lXmu -lXext -lX11 -lm */
  32. #include <stdio.h>
  33. #include <stdlib.h>
  34. #include <string.h>
  35. #include <math.h>       /* for cos(), sin(), and sqrt() */
  36. #include <GL/glut.h>    /* OpenGL Utility Toolkit header */
  37. /* Some <math.h> files do not define M_PI... */
  38. #ifndef M_PI
  39. #define M_PI 3.14159265358979323846
  40. #endif
  41. /* Variable controlling various rendering modes. */
  42. static int stencilReflection = 1, stencilShadow = 1, offsetShadow = 1;
  43. static int renderShadow = 0, renderDinosaur = 1, renderReflection = 0;
  44. static int linearFiltering = 0, useMipmaps = 0, useTexture = 0;
  45. static int reportSpeed = 0;
  46. static int animation = 0;
  47. static GLboolean lightSwitch = GL_TRUE;
  48. static int directionalLight = 1;
  49. static int forceExtension = 0;
  50. static int haloMagic = 0, blendedHalo = 0;
  51. static GLfloat haloScale = 1.0, haloTime = 0.0;
  52. /* Time varying or user-controled variables. */
  53. static float jump = 0.0;
  54. static float lightAngle = 0.0, lightHeight = 20;
  55. GLfloat angle = -150;   /* in degrees */
  56. GLfloat angle2 = 30;   /* in degrees */
  57. int moving, startx, starty;
  58. int lightMoving = 0, lightStartX, lightStartY;
  59. enum {
  60.   MISSING, EXTENSION, ONE_DOT_ONE
  61. };
  62. int polygonOffsetVersion;
  63. static GLdouble bodyWidth = 3.0;
  64. /* *INDENT-OFF* */
  65. static GLfloat body[][2] = { {0, 3}, {1, 1}, {5, 1}, {8, 4}, {10, 4}, {11, 5},
  66.   {11, 11.5}, {13, 12}, {13, 13}, {10, 13.5}, {13, 14}, {13, 15}, {11, 16},
  67.   {8, 16}, {7, 15}, {7, 13}, {8, 12}, {7, 11}, {6, 6}, {4, 3}, {3, 2},
  68.   {1, 2} };
  69. static GLfloat arm[][2] = { {8, 10}, {9, 9}, {10, 9}, {13, 8}, {14, 9}, {16, 9},
  70.   {15, 9.5}, {16, 10}, {15, 10}, {15.5, 11}, {14.5, 10}, {14, 11}, {14, 10},
  71.   {13, 9}, {11, 11}, {9, 11} };
  72. static GLfloat leg[][2] = { {8, 6}, {8, 4}, {9, 3}, {9, 2}, {8, 1}, {8, 0.5}, {9, 0},
  73.   {12, 0}, {10, 1}, {10, 2}, {12, 4}, {11, 6}, {10, 7}, {9, 7} };
  74. static GLfloat eye[][2] = { {8.75, 15}, {9, 14.7}, {9.6, 14.7}, {10.1, 15},
  75.   {9.6, 15.25}, {9, 15.25} };
  76. static GLfloat lightPosition[4];
  77. static GLfloat lightColor[] = {0.8, 1.0, 0.8, 1.0}; /* green-tinted */
  78. static GLfloat skinColor[] = {0.1, 1.0, 0.1, 1.0}, eyeColor[] = {1.0, 0.2, 0.2, 1.0};
  79. /* *INDENT-ON* */
  80. /* Nice floor texture tiling pattern. */
  81. static char *circles[] = {
  82.   "....xxxx........",
  83.   "..xxxxxxxx......",
  84.   ".xxxxxxxxxx.....",
  85.   ".xxx....xxx.....",
  86.   "xxx......xxx....",
  87.   "xxx......xxx....",
  88.   "xxx......xxx....",
  89.   "xxx......xxx....",
  90.   ".xxx....xxx.....",
  91.   ".xxxxxxxxxx.....",
  92.   "..xxxxxxxx......",
  93.   "....xxxx........",
  94.   "................",
  95.   "................",
  96.   "................",
  97.   "................",
  98. };
  99. static void
  100. makeFloorTexture(void)
  101. {
  102.   GLubyte floorTexture[16][16][3];
  103.   GLubyte *loc;
  104.   int s, t;
  105.   /* Setup RGB image for the texture. */
  106.   loc = (GLubyte*) floorTexture;
  107.   for (t = 0; t < 16; t++) {
  108.     for (s = 0; s < 16; s++) {
  109.       if (circles[t][s] == 'x') {
  110. /* Nice blue. */
  111.         loc[0] = 0x1f;
  112.         loc[1] = 0x1f;
  113.         loc[2] = 0x8f;
  114.       } else {
  115. /* Light gray. */
  116.         loc[0] = 0xca;
  117.         loc[1] = 0xca;
  118.         loc[2] = 0xca;
  119.       }
  120.       loc += 3;
  121.     }
  122.   }
  123.   glPixelStorei(GL_UNPACK_ALIGNMENT, 1);
  124.   if (useMipmaps) {
  125.     glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER,
  126.       GL_LINEAR_MIPMAP_LINEAR);
  127.     gluBuild2DMipmaps(GL_TEXTURE_2D, 3, 16, 16,
  128.       GL_RGB, GL_UNSIGNED_BYTE, floorTexture);
  129.   } else {
  130.     if (linearFiltering) {
  131.       glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_LINEAR);
  132.     } else {
  133.       glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);
  134.     }
  135.     glTexImage2D(GL_TEXTURE_2D, 0, 3, 16, 16, 0,
  136.       GL_RGB, GL_UNSIGNED_BYTE, floorTexture);
  137.   }
  138. }
  139. enum {
  140.   X, Y, Z, W
  141. };
  142. enum {
  143.   A, B, C, D
  144. };
  145. /* Create a matrix that will project the desired shadow. */
  146. void
  147. shadowMatrix(GLfloat shadowMat[4][4],
  148.   GLfloat groundplane[4],
  149.   GLfloat lightpos[4])
  150. {
  151.   GLfloat dot;
  152.   /* Find dot product between light position vector and ground plane normal. */
  153.   dot = groundplane[X] * lightpos[X] +
  154.     groundplane[Y] * lightpos[Y] +
  155.     groundplane[Z] * lightpos[Z] +
  156.     groundplane[W] * lightpos[W];
  157.   shadowMat[0][0] = dot - lightpos[X] * groundplane[X];
  158.   shadowMat[1][0] = 0.f - lightpos[X] * groundplane[Y];
  159.   shadowMat[2][0] = 0.f - lightpos[X] * groundplane[Z];
  160.   shadowMat[3][0] = 0.f - lightpos[X] * groundplane[W];
  161.   shadowMat[X][1] = 0.f - lightpos[Y] * groundplane[X];
  162.   shadowMat[1][1] = dot - lightpos[Y] * groundplane[Y];
  163.   shadowMat[2][1] = 0.f - lightpos[Y] * groundplane[Z];
  164.   shadowMat[3][1] = 0.f - lightpos[Y] * groundplane[W];
  165.   shadowMat[X][2] = 0.f - lightpos[Z] * groundplane[X];
  166.   shadowMat[1][2] = 0.f - lightpos[Z] * groundplane[Y];
  167.   shadowMat[2][2] = dot - lightpos[Z] * groundplane[Z];
  168.   shadowMat[3][2] = 0.f - lightpos[Z] * groundplane[W];
  169.   shadowMat[X][3] = 0.f - lightpos[W] * groundplane[X];
  170.   shadowMat[1][3] = 0.f - lightpos[W] * groundplane[Y];
  171.   shadowMat[2][3] = 0.f - lightpos[W] * groundplane[Z];
  172.   shadowMat[3][3] = dot - lightpos[W] * groundplane[W];
  173. }
  174. /* Find the plane equation given 3 points. */
  175. void
  176. findPlane(GLfloat plane[4],
  177.   GLfloat v0[3], GLfloat v1[3], GLfloat v2[3])
  178. {
  179.   GLfloat vec0[3], vec1[3];
  180.   /* Need 2 vectors to find cross product. */
  181.   vec0[X] = v1[X] - v0[X];
  182.   vec0[Y] = v1[Y] - v0[Y];
  183.   vec0[Z] = v1[Z] - v0[Z];
  184.   vec1[X] = v2[X] - v0[X];
  185.   vec1[Y] = v2[Y] - v0[Y];
  186.   vec1[Z] = v2[Z] - v0[Z];
  187.   /* find cross product to get A, B, and C of plane equation */
  188.   plane[A] = vec0[Y] * vec1[Z] - vec0[Z] * vec1[Y];
  189.   plane[B] = -(vec0[X] * vec1[Z] - vec0[Z] * vec1[X]);
  190.   plane[C] = vec0[X] * vec1[Y] - vec0[Y] * vec1[X];
  191.   plane[D] = -(plane[A] * v0[X] + plane[B] * v0[Y] + plane[C] * v0[Z]);
  192. }
  193. void
  194. extrudeSolidFromPolygon(GLfloat data[][2], unsigned int dataSize,
  195.   GLdouble thickness, GLuint side, GLuint edge, GLuint whole)
  196. {
  197.   static GLUtriangulatorObj *tobj = NULL;
  198.   GLdouble vertex[3], dx, dy, len;
  199.   int i;
  200.   int count = dataSize / (int) (2 * sizeof(GLfloat));
  201.   if (tobj == NULL) {
  202.     tobj = gluNewTess();  /* create and initialize a GLU
  203.                              polygon tesselation object */
  204.     gluTessCallback(tobj, GLU_BEGIN, glBegin);
  205.     gluTessCallback(tobj, GLU_VERTEX, glVertex2fv);  /* semi-tricky */
  206.     gluTessCallback(tobj, GLU_END, glEnd);
  207.   }
  208.   glNewList(side, GL_COMPILE);
  209.   glShadeModel(GL_SMOOTH);  /* smooth minimizes seeing
  210.                                tessellation */
  211.   gluBeginPolygon(tobj);
  212.   for (i = 0; i < count; i++) {
  213.     vertex[0] = data[i][0];
  214.     vertex[1] = data[i][1];
  215.     vertex[2] = 0;
  216.     gluTessVertex(tobj, vertex, data[i]);
  217.   }
  218.   gluEndPolygon(tobj);
  219.   glEndList();
  220.   glNewList(edge, GL_COMPILE);
  221.   glShadeModel(GL_FLAT);  /* flat shade keeps angular hands
  222.                              from being "smoothed" */
  223.   glBegin(GL_QUAD_STRIP);
  224.   for (i = 0; i <= count; i++) {
  225.     /* mod function handles closing the edge */
  226.     glVertex3f(data[i % count][0], data[i % count][1], 0.0);
  227.     glVertex3f(data[i % count][0], data[i % count][1], thickness);
  228.     /* Calculate a unit normal by dividing by Euclidean
  229.        distance. We * could be lazy and use
  230.        glEnable(GL_NORMALIZE) so we could pass in * arbitrary
  231.        normals for a very slight performance hit. */
  232.     dx = data[(i + 1) % count][1] - data[i % count][1];
  233.     dy = data[i % count][0] - data[(i + 1) % count][0];
  234.     len = sqrt(dx * dx + dy * dy);
  235.     glNormal3f(dx / len, dy / len, 0.0);
  236.   }
  237.   glEnd();
  238.   glEndList();
  239.   glNewList(whole, GL_COMPILE);
  240.   glFrontFace(GL_CW);
  241.   glCallList(edge);
  242.   glNormal3f(0.0, 0.0, -1.0);  /* constant normal for side */
  243.   glCallList(side);
  244.   glPushMatrix();
  245.   glTranslatef(0.0, 0.0, thickness);
  246.   glFrontFace(GL_CCW);
  247.   glNormal3f(0.0, 0.0, 1.0);  /* opposite normal for other side */
  248.   glCallList(side);
  249.   glPopMatrix();
  250.   glEndList();
  251. }
  252. /* Enumerants for refering to display lists. */
  253. typedef enum {
  254.   RESERVED, BODY_SIDE, BODY_EDGE, BODY_WHOLE, ARM_SIDE, ARM_EDGE, ARM_WHOLE,
  255.   LEG_SIDE, LEG_EDGE, LEG_WHOLE, EYE_SIDE, EYE_EDGE, EYE_WHOLE
  256. } displayLists;
  257. static void
  258. makeDinosaur(void)
  259. {
  260.   extrudeSolidFromPolygon(body, sizeof(body), bodyWidth,
  261.     BODY_SIDE, BODY_EDGE, BODY_WHOLE);
  262.   extrudeSolidFromPolygon(arm, sizeof(arm), bodyWidth / 4,
  263.     ARM_SIDE, ARM_EDGE, ARM_WHOLE);
  264.   extrudeSolidFromPolygon(leg, sizeof(leg), bodyWidth / 2,
  265.     LEG_SIDE, LEG_EDGE, LEG_WHOLE);
  266.   extrudeSolidFromPolygon(eye, sizeof(eye), bodyWidth + 0.2,
  267.     EYE_SIDE, EYE_EDGE, EYE_WHOLE);
  268. }
  269. static void
  270. drawDinosaur(void)
  271. {
  272.   glPushMatrix();
  273.   /* Translate the dinosaur to be at (0,8,0). */
  274.   glTranslatef(-8, -8, -bodyWidth / 2);
  275.   glTranslatef(0.0, jump, 0.0);
  276.   glMaterialfv(GL_FRONT, GL_DIFFUSE, skinColor);
  277.   glCallList(BODY_WHOLE);
  278.   glTranslatef(0.0, 0.0, bodyWidth);
  279.   glCallList(ARM_WHOLE);
  280.   glCallList(LEG_WHOLE);
  281.   glTranslatef(0.0, 0.0, -bodyWidth - bodyWidth / 4);
  282.   glCallList(ARM_WHOLE);
  283.   glTranslatef(0.0, 0.0, -bodyWidth / 4);
  284.   glCallList(LEG_WHOLE);
  285.   glTranslatef(0.0, 0.0, bodyWidth / 2 - 0.1);
  286.   glMaterialfv(GL_FRONT, GL_DIFFUSE, eyeColor);
  287.   glCallList(EYE_WHOLE);
  288.   glPopMatrix();
  289. }
  290. enum {
  291.   MOD_DINO, MOD_SPHERE, MOD_CUBE, MOD_ICO
  292. };
  293. static int currentModel = MOD_DINO;
  294. static GLfloat blueMaterial[] = {0.0, 0.2, 1.0, 1.0},
  295.   redMaterial[] = {0.6, 0.1, 0.0, 1.0},
  296.   purpleMaterial[] = {0.3, 0.0, 0.3, 1.0},
  297.   greenMaterial[] = {1.0, 0.2, 0.0, 1.0};
  298. static void
  299. drawModel(void)
  300. {
  301.   switch(currentModel) {
  302.   case MOD_DINO:
  303.     drawDinosaur();
  304.     break;
  305.   case MOD_SPHERE:
  306.     glMaterialfv(GL_FRONT, GL_DIFFUSE, blueMaterial);
  307.     glutSolidSphere(6.0, 15, 15);
  308.     break;
  309.   case MOD_CUBE:
  310.     glMaterialfv(GL_FRONT, GL_DIFFUSE, redMaterial);
  311.     glutSolidCube(6.0);
  312.     break;
  313.   case MOD_ICO:
  314.     glMaterialfv(GL_FRONT, GL_DIFFUSE, purpleMaterial);
  315.     glPushMatrix();
  316.       glEnable(GL_NORMALIZE);
  317.       glScalef(7.0, 7.0, 7.0);
  318.       glutSolidIcosahedron();
  319.       glDisable(GL_NORMALIZE);
  320.     glPopMatrix();
  321.     break;
  322.   }
  323. }
  324. static void
  325. drawBox(GLfloat xsize, GLfloat ysize, GLfloat zsize)
  326. {
  327.   static GLfloat n[6][3] =
  328.   {
  329.     {-1.0, 0.0, 0.0},
  330.     {0.0, 1.0, 0.0},
  331.     {1.0, 0.0, 0.0},
  332.     {0.0, -1.0, 0.0},
  333.     {0.0, 0.0, 1.0},
  334.     {0.0, 0.0, -1.0}
  335.   };
  336.   static GLint faces[6][4] =
  337.   {
  338.     {0, 1, 2, 3},
  339.     {3, 2, 6, 7},
  340.     {7, 6, 5, 4},
  341.     {4, 5, 1, 0},
  342.     {5, 6, 2, 1},
  343.     {7, 4, 0, 3}
  344.   };
  345.   GLfloat v[8][3];
  346.   GLint i;
  347.   v[0][0] = v[1][0] = v[2][0] = v[3][0] = -xsize / 2;
  348.   v[4][0] = v[5][0] = v[6][0] = v[7][0] = xsize / 2;
  349.   v[0][1] = v[1][1] = v[4][1] = v[5][1] = -ysize / 2;
  350.   v[2][1] = v[3][1] = v[6][1] = v[7][1] = ysize / 2;
  351.   v[0][2] = v[3][2] = v[4][2] = v[7][2] = -zsize / 2;
  352.   v[1][2] = v[2][2] = v[5][2] = v[6][2] = zsize / 2;
  353.   for (i = 0; i < 6; i++) {
  354.     glBegin(GL_QUADS);
  355.     glNormal3fv(&n[i][0]);
  356.     glVertex3fv(&v[faces[i][0]][0]);
  357.     glVertex3fv(&v[faces[i][1]][0]);
  358.     glVertex3fv(&v[faces[i][2]][0]);
  359.     glVertex3fv(&v[faces[i][3]][0]);
  360.     glEnd();
  361.   }
  362. }
  363. static void
  364. drawPillar(void)
  365. {
  366.   glEnable(GL_NORMALIZE);
  367.   glMaterialfv(GL_FRONT, GL_DIFFUSE, greenMaterial);
  368.   glPushMatrix();
  369.     glTranslatef(8.0, 4.01, 8.0);
  370.     drawBox(2.0, 8.0, 2.0);
  371.     glutSolidCube(2.0);
  372.   glPopMatrix();
  373.   glDisable(GL_NORMALIZE);
  374. }
  375. static GLfloat floorVertices[4][3] = {
  376.   { -20.0, 0.0, 20.0 },
  377.   { 20.0, 0.0, 20.0 },
  378.   { 20.0, 0.0, -20.0 },
  379.   { -20.0, 0.0, -20.0 },
  380. };
  381. /* Draw a floor (possibly textured). */
  382. static void
  383. drawFloor(void)
  384. {
  385.   glDisable(GL_LIGHTING);
  386.   if (useTexture) {
  387.     glEnable(GL_TEXTURE_2D);
  388.   }
  389.   glBegin(GL_QUADS);
  390.     glTexCoord2f(0.0, 0.0);
  391.     glVertex3fv(floorVertices[0]);
  392.     glTexCoord2f(0.0, 16.0);
  393.     glVertex3fv(floorVertices[1]);
  394.     glTexCoord2f(16.0, 16.0);
  395.     glVertex3fv(floorVertices[2]);
  396.     glTexCoord2f(16.0, 0.0);
  397.     glVertex3fv(floorVertices[3]);
  398.   glEnd();
  399.   if (useTexture) {
  400.     glDisable(GL_TEXTURE_2D);
  401.   }
  402.   glEnable(GL_LIGHTING);
  403. }
  404. static GLfloat floorPlane[4];
  405. static GLfloat floorShadow[4][4];
  406. static void
  407. redraw(void)
  408. {
  409.   int start, end;
  410.   if (reportSpeed) {
  411.     start = glutGet(GLUT_ELAPSED_TIME);
  412.   }
  413.   /* Clear; default stencil clears to zero. */
  414.   if ((stencilReflection && renderReflection) || (stencilShadow && renderShadow) || (haloScale > 1.0)) {
  415.     glStencilMask(0xffffffff);
  416.     glClearStencil(0x4);
  417.     glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT | GL_STENCIL_BUFFER_BIT);
  418.   } else {
  419.     /* Avoid clearing stencil when not using it. */
  420.     glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
  421.   }
  422.   /* Reposition the light source. */
  423.   lightPosition[0] = 15*cos(lightAngle);
  424.   lightPosition[1] = lightHeight;
  425.   lightPosition[2] = 15*sin(lightAngle);
  426.   if (directionalLight) {
  427.     lightPosition[3] = 0.0;
  428.   } else {
  429.     lightPosition[3] = 1.0;
  430.   }
  431.   shadowMatrix(floorShadow, floorPlane, lightPosition);
  432.   glPushMatrix();
  433.     /* Perform scene rotations based on user mouse input. */
  434.     glRotatef(angle2, 1.0, 0.0, 0.0);
  435.     glRotatef(angle, 0.0, 1.0, 0.0);
  436.      
  437.     /* Tell GL new light source position. */
  438.     glLightfv(GL_LIGHT0, GL_POSITION, lightPosition);
  439.     if (renderReflection) {
  440.       if (stencilReflection) {
  441.         /* We can eliminate the visual "artifact" of seeing the "flipped"
  442.       model underneath the floor by using stencil.  The idea is
  443.    draw the floor without color or depth update but so that 
  444.    a stencil value of one is where the floor will be.  Later when
  445.    rendering the model reflection, we will only update pixels
  446.    with a stencil value of 1 to make sure the reflection only
  447.    lives on the floor, not below the floor. */
  448.         /* Don't update color or depth. */
  449.         glDisable(GL_DEPTH_TEST);
  450.         glColorMask(GL_FALSE, GL_FALSE, GL_FALSE, GL_FALSE);
  451.         /* Draw 1 into the stencil buffer. */
  452.         glEnable(GL_STENCIL_TEST);
  453.         glStencilOp(GL_KEEP, GL_KEEP, GL_REPLACE);
  454.         glStencilFunc(GL_ALWAYS, 1, 0x1);
  455. glStencilMask(0x1);
  456.         /* Now render floor; floor pixels just get their stencil set to 1. */
  457.         drawFloor();
  458.         /* Re-enable update of color and depth. */ 
  459.         glColorMask(GL_TRUE, GL_TRUE, GL_TRUE, GL_TRUE);
  460.         glEnable(GL_DEPTH_TEST);
  461.         /* Now, only render where stencil is set to 1. */
  462.         glStencilFunc(GL_EQUAL, 1, 0x1);  /* draw if ==1 */
  463.         glStencilOp(GL_KEEP, GL_KEEP, GL_KEEP);
  464.       }
  465.       glPushMatrix();
  466.         /* The critical reflection step: Reflect 3D model through the floor
  467.            (the Y=0 plane) to make a relection. */
  468.         glScalef(1.0, -1.0, 1.0);
  469. /* Reflect the light position. */
  470.         glLightfv(GL_LIGHT0, GL_POSITION, lightPosition);
  471.         /* To avoid our normals getting reversed and hence botched lighting
  472.    on the reflection, turn on normalize.  */
  473.         glEnable(GL_NORMALIZE);
  474.         glCullFace(GL_FRONT);
  475.         /* Draw the reflected model. */
  476. glPushMatrix();
  477.   glTranslatef(0, 8.01, 0);
  478.           drawModel();
  479. glPopMatrix();
  480. drawPillar();
  481.         /* Disable noramlize again and re-enable back face culling. */
  482.         glDisable(GL_NORMALIZE);
  483.         glCullFace(GL_BACK);
  484.       glPopMatrix();
  485.       /* Switch back to the unreflected light position. */
  486.       glLightfv(GL_LIGHT0, GL_POSITION, lightPosition);
  487.       if (stencilReflection) {
  488.         glDisable(GL_STENCIL_TEST);
  489.       }
  490.     }
  491.     /* Back face culling will get used to only draw either the top or the
  492.        bottom floor.  This let's us get a floor with two distinct
  493.        appearances.  The top floor surface is reflective and kind of red.
  494.        The bottom floor surface is not reflective and blue. */
  495.     /* Draw "bottom" of floor in blue. */
  496.     glFrontFace(GL_CW);  /* Switch face orientation. */
  497.     glColor4f(0.1, 0.1, 0.7, 1.0);
  498.     drawFloor();
  499.     glFrontFace(GL_CCW);
  500.     if (renderShadow && stencilShadow) {
  501.      /* Draw the floor with stencil value 2.  This helps us only 
  502. draw the shadow once per floor pixel (and only on the
  503.         floor pixels). */
  504.       glEnable(GL_STENCIL_TEST);
  505.       glStencilFunc(GL_ALWAYS, 0x2, 0x2);
  506.       glStencilMask(0x2);
  507.       glStencilOp(GL_KEEP, GL_KEEP, GL_REPLACE);
  508.     }
  509.     /* Draw "top" of floor.  Use blending to blend in reflection. */
  510.     glEnable(GL_BLEND);
  511.     glBlendFunc(GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
  512.     glColor4f(1.0, 1.0, 1.0, 0.3);
  513.     drawFloor();
  514.     glDisable(GL_BLEND);
  515.     if (renderShadow && stencilShadow) {
  516.       glDisable(GL_STENCIL_TEST);
  517.     }
  518.     if (renderDinosaur) {
  519.       drawPillar();
  520.       if (haloScale > 1.0) {
  521. /* If halo effect is enabled, draw the model with its stencil set to 6
  522.    (arbitary value); later, we'll make sure not to update pixels tagged
  523.    as 6. */
  524.         glEnable(GL_STENCIL_TEST);
  525.         glStencilFunc(GL_ALWAYS, 0x0, 0x4);
  526.         glStencilMask(0x4);
  527.         glStencilOp(GL_KEEP, GL_KEEP, GL_REPLACE);
  528.       } 
  529.       /* Draw "actual" dinosaur (or other model), not its reflection. */
  530.       glPushMatrix();
  531.         glTranslatef(0, 8.01, 0);
  532.         drawModel();
  533.       glPopMatrix();
  534.     }
  535.     /* Begin shadow render. */
  536.     if (renderShadow) {
  537.       /* Render the projected shadow. */
  538.       if (stencilShadow) {
  539.         /* Now, only render where stencil is set above 5 (ie, 6 where
  540.    the top floor is).  Update stencil with 2 where the shadow
  541.    gets drawn so we don't redraw (and accidently reblend) the
  542.    shadow). */
  543. glEnable(GL_STENCIL_TEST);
  544.         glStencilFunc(GL_NOTEQUAL, 0x0, 0x2);
  545.         glStencilMask(0x2);
  546.         glStencilOp(GL_KEEP, GL_KEEP, GL_REPLACE);
  547.       }
  548.       /* To eliminate depth buffer artifacts, we use polygon offset
  549.  to raise the depth of the projected shadow slightly so
  550.  that it does not depth buffer alias with the floor. */
  551.       if (offsetShadow) {
  552. switch (polygonOffsetVersion) {
  553. case EXTENSION:
  554. #ifdef GL_EXT_polygon_offset
  555.   glEnable(GL_POLYGON_OFFSET_EXT);
  556.   break;
  557. #endif
  558. #ifdef GL_VERSION_1_1
  559. case ONE_DOT_ONE:
  560.           glEnable(GL_POLYGON_OFFSET_FILL);
  561.   break;
  562. #endif
  563. case MISSING:
  564.   /* Oh well. */
  565.   break;
  566. }
  567.       }
  568.       /* Render 50% black shadow color on top of whatever the
  569.          floor appareance is. */
  570.       glEnable(GL_BLEND);
  571.       glBlendFunc(GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
  572.       glDisable(GL_LIGHTING);  /* Force the 50% black. */
  573.       glColor4f(0.0, 0.0, 0.0, 0.5);
  574.       glPushMatrix();
  575. /* Project the shadow. */
  576.         glMultMatrixf((GLfloat *) floorShadow);
  577. glPushMatrix();
  578.           glTranslatef(0, 8.01, 0);
  579.           drawModel();
  580. glPopMatrix();
  581. drawPillar();
  582.       glPopMatrix();
  583.       glDisable(GL_BLEND);
  584.       glEnable(GL_LIGHTING);
  585.       if (offsetShadow) {
  586. switch (polygonOffsetVersion) {
  587. #ifdef GL_EXT_polygon_offset
  588. case EXTENSION:
  589.   glDisable(GL_POLYGON_OFFSET_EXT);
  590.   break;
  591. #endif
  592. #ifdef GL_VERSION_1_1
  593. case ONE_DOT_ONE:
  594.           glDisable(GL_POLYGON_OFFSET_FILL);
  595.   break;
  596. #endif
  597. case MISSING:
  598.   /* Oh well. */
  599.   break;
  600. }
  601.       }
  602.       if (stencilShadow) {
  603.         glDisable(GL_STENCIL_TEST);
  604.       }
  605.     } /* End shadow render. */
  606.     /* Begin light source location render. */
  607.     glPushMatrix();
  608.       glDisable(GL_LIGHTING);
  609.       glColor3f(1.0, 1.0, 0.0);
  610.       if (directionalLight) {
  611.         /* Draw an arrowhead. */
  612.         glDisable(GL_CULL_FACE);
  613.         glTranslatef(lightPosition[0], lightPosition[1], lightPosition[2]);
  614.         glRotatef(lightAngle * -180.0 / M_PI, 0, 1, 0);
  615.         glRotatef(atan(lightHeight/12) * 180.0 / M_PI, 0, 0, 1);
  616.         glBegin(GL_TRIANGLE_FAN);
  617.      glVertex3f(0, 0, 0);
  618.   glVertex3f(2, 1, 1);
  619.   glVertex3f(2, -1, 1);
  620.   glVertex3f(2, -1, -1);
  621.   glVertex3f(2, 1, -1);
  622.   glVertex3f(2, 1, 1);
  623.         glEnd();
  624.         /* Draw a white line from light direction. */
  625.         glColor3f(1.0, 1.0, 1.0);
  626.         glBegin(GL_LINES);
  627.   glVertex3f(0.1, 0, 0);
  628.   glVertex3f(5, 0, 0);
  629.         glEnd();
  630.         glEnable(GL_CULL_FACE);
  631.       } else {
  632.         /* Draw a yellow ball at the light source. */
  633.         glTranslatef(lightPosition[0], lightPosition[1], lightPosition[2]);
  634.         glutSolidSphere(1.0, 5, 5);
  635.       }
  636.       glEnable(GL_LIGHTING);
  637.     glPopMatrix();
  638.     /* End light source location render. */
  639.     /* Add a halo effect around the 3D model. */
  640.     if (haloScale > 1.0) {
  641.       glDisable(GL_LIGHTING);
  642.       if (blendedHalo) {
  643. /* If we are doing a nice blended halo, enable blending and
  644.    make sure we only blend a halo pixel once and that we do
  645.    not draw to pixels tagged as 6 (where the model is). */
  646. glEnable(GL_BLEND);
  647. glEnable(GL_STENCIL_TEST);
  648. glColor4f(0.8, 0.8, 0.0, 0.3);  /* 30% sorta yellow. */
  649.         glStencilFunc(GL_EQUAL, 0x4, 0x4);
  650.         glStencilMask(0x4);
  651.         glStencilOp(GL_KEEP, GL_KEEP, GL_INVERT);
  652.       } else {
  653. /* Be cheap; no blending.  Just draw yellow halo but not updating
  654.    pixels where the model is.  We don't update stencil at all. */
  655.         glDisable(GL_BLEND);
  656. glEnable(GL_STENCIL_TEST);
  657. glColor3f(0.5, 0.5, 0.0);  /* Half yellow. */
  658.         glStencilFunc(GL_EQUAL, 0x4, 0x4);
  659.         glStencilMask(0x4);
  660.         glStencilOp(GL_KEEP, GL_KEEP, GL_INVERT);
  661.       }
  662.       glPushMatrix();
  663.         glTranslatef(0, 8.01, 0);
  664.         glScalef(haloScale, haloScale, haloScale);
  665.         drawModel();
  666.       glPopMatrix();
  667.       if (blendedHalo) {
  668.         glDisable(GL_BLEND);
  669.       }
  670.       glDisable(GL_STENCIL_TEST);
  671.       glEnable(GL_LIGHTING);
  672.     }
  673.     /* End halo effect render. */
  674.   glPopMatrix();
  675.   if (reportSpeed) {
  676.     glFinish();
  677.     end = glutGet(GLUT_ELAPSED_TIME);
  678.     printf("Speed %.3g frames/sec (%d ms)n", 1000.0/(end-start), end-start);
  679.   }
  680.   glutSwapBuffers();
  681. }
  682. /* ARGSUSED2 */
  683. static void
  684. mouse(int button, int state, int x, int y)
  685. {
  686.   if (button == GLUT_LEFT_BUTTON) {
  687.     if (state == GLUT_DOWN) {
  688.       moving = 1;
  689.       startx = x;
  690.       starty = y;
  691.     }
  692.     if (state == GLUT_UP) {
  693.       moving = 0;
  694.     }
  695.   }
  696.   if (button == GLUT_MIDDLE_BUTTON) {
  697.     if (state == GLUT_DOWN) {
  698.       lightMoving = 1;
  699.       lightStartX = x;
  700.       lightStartY = y;
  701.     }
  702.     if (state == GLUT_UP) {
  703.       lightMoving = 0;
  704.     }
  705.   }
  706. }
  707. /* ARGSUSED1 */
  708. static void
  709. motion(int x, int y)
  710. {
  711.   if (moving) {
  712.     angle = angle + (x - startx);
  713.     angle2 = angle2 + (y - starty);
  714.     startx = x;
  715.     starty = y;
  716.     glutPostRedisplay();
  717.   }
  718.   if (lightMoving) {
  719.     lightAngle += (x - lightStartX)/40.0;
  720.     lightHeight += (lightStartY - y)/20.0;
  721.     lightStartX = x;
  722.     lightStartY = y;
  723.     glutPostRedisplay();
  724.   }
  725. }
  726. static const float maxHalo[] = { 0.2, 0.35, 0.3, 0.5 };
  727. /* Advance time varying state when idle callback registered. */
  728. static void
  729. idle(void)
  730. {
  731.   static float time = 0.0;
  732.   if (animation) {
  733.     time = glutGet(GLUT_ELAPSED_TIME) / 500.0;
  734.     jump = 4.0 * fabs(sin(time)*0.8);
  735.     if (!lightMoving) {
  736.       lightAngle += 0.03;
  737.     }
  738.   }
  739.   if (haloMagic) {
  740.     haloTime += 0.1;
  741.     haloScale = 1.0 + maxHalo[currentModel] * sin(haloTime);
  742.     if (haloScale <= 1.0) {
  743.       haloMagic = 0;
  744.       if (!animation) {
  745.         glutIdleFunc(NULL);
  746.       }
  747.     }
  748.   }
  749.   glutPostRedisplay();
  750. }
  751. enum {
  752.   M_NONE, M_BLENDED_HALO, M_SHOW_HALO, M_SWITCH_MODEL, M_MOTION, M_LIGHT,
  753.   M_TEXTURE, M_SHADOWS, M_REFLECTION, M_DINOSAUR,
  754.   M_STENCIL_REFLECTION, M_STENCIL_SHADOW, M_OFFSET_SHADOW,
  755.   M_POSITIONAL, M_DIRECTIONAL, M_PERFORMANCE
  756. };
  757. static void
  758. controlLights(int value)
  759. {
  760.   switch (value) {
  761.   case M_NONE:
  762.     return;
  763.   case M_SWITCH_MODEL:
  764.     currentModel = (currentModel + 1) % 4;
  765.     break;
  766.   case M_SHOW_HALO:
  767.     haloScale = 1.0 + maxHalo[currentModel];
  768.     break;
  769.   case M_BLENDED_HALO:
  770.     blendedHalo = 1 - blendedHalo;
  771.     break;
  772.   case M_MOTION:
  773.     animation = 1 - animation;
  774.     if (animation || haloMagic) {
  775.       glutIdleFunc(idle);
  776.     } else {
  777.       glutIdleFunc(NULL);
  778.     }
  779.     break;
  780.   case M_LIGHT:
  781.     lightSwitch = !lightSwitch;
  782.     if (lightSwitch) {
  783.       glEnable(GL_LIGHT0);
  784.     } else {
  785.       glDisable(GL_LIGHT0);
  786.     }
  787.     break;
  788.   case M_TEXTURE:
  789.     useTexture = !useTexture;
  790.     break;
  791.   case M_SHADOWS:
  792.     renderShadow = 1 - renderShadow;
  793.     break;
  794.   case M_REFLECTION:
  795.     renderReflection = 1 - renderReflection;
  796.     break;
  797.   case M_DINOSAUR:
  798.     renderDinosaur = 1 - renderDinosaur;
  799.     break;
  800.   case M_STENCIL_REFLECTION:
  801.     stencilReflection = 1 - stencilReflection;
  802.     break;
  803.   case M_STENCIL_SHADOW:
  804.     stencilShadow = 1 - stencilShadow;
  805.     break;
  806.   case M_OFFSET_SHADOW:
  807.     offsetShadow = 1 - offsetShadow;
  808.     break;
  809.   case M_POSITIONAL:
  810.     directionalLight = 0;
  811.     break;
  812.   case M_DIRECTIONAL:
  813.     directionalLight = 1;
  814.     break;
  815.   case M_PERFORMANCE:
  816.     reportSpeed = 1 - reportSpeed;
  817.     break;
  818.   }
  819.   glutPostRedisplay();
  820. }
  821. /* When not visible, stop animating.  Restart when visible again. */
  822. static void 
  823. visible(int vis)
  824. {
  825.   if (vis == GLUT_VISIBLE) {
  826.     if (animation || haloMagic)
  827.       glutIdleFunc(idle);
  828.   } else {
  829.     if (!animation && !haloMagic)
  830.       glutIdleFunc(NULL);
  831.   }
  832. }
  833. /* Press any key to redraw; good when motion stopped and
  834.    performance reporting on. */
  835. /* ARGSUSED */
  836. static void
  837. key(unsigned char c, int x, int y)
  838. {
  839.   if (c == 27) {
  840.     exit(0);  /* IRIS GLism, Escape quits. */
  841.   }
  842.   if (c == ' ') {
  843.     haloMagic = 1;
  844.     haloTime = 0.0;
  845.     glutIdleFunc(idle);
  846.   }
  847.   glutPostRedisplay();
  848. }
  849. /* Press any key to redraw; good when motion stopped and
  850.    performance reporting on. */
  851. /* ARGSUSED */
  852. static void
  853. special(int k, int x, int y)
  854. {
  855.   glutPostRedisplay();
  856. }
  857. static int
  858. supportsOneDotOne(void)
  859. {
  860.   const char *version;
  861.   int major, minor;
  862.   version = (char *) glGetString(GL_VERSION);
  863.   if (sscanf(version, "%d.%d", &major, &minor) == 2)
  864.     return major >= 1 && minor >= 1;
  865.   return 0;            /* OpenGL version string malformed! */
  866. }
  867. int
  868. main(int argc, char **argv)
  869. {
  870.   int i;
  871.   glutInit(&argc, argv);
  872.   for (i=1; i<argc; i++) {
  873.     if (!strcmp("-linear", argv[i])) {
  874.       linearFiltering = 1;
  875.     } else if (!strcmp("-mipmap", argv[i])) {
  876.       useMipmaps = 1;
  877.     } else if (!strcmp("-ext", argv[i])) {
  878.       forceExtension = 1;
  879.     }
  880.   }
  881.   glutInitDisplayMode(GLUT_RGB | GLUT_DOUBLE | GLUT_DEPTH | GLUT_STENCIL | GLUT_MULTISAMPLE);
  882. #if 1
  883.   /* In GLUT 4.0, you'll be able to do this an be sure to
  884.      get 2 bits of stencil if the machine has it for you. */
  885.   glutInitDisplayString("samples stencil~3 rgb double depth");
  886. #endif
  887.   glutCreateWindow("OpenGL Halo Magic (hit Space)");
  888.   if (glutGet(GLUT_WINDOW_STENCIL_SIZE) < 3) {
  889.     printf("halomagic: Sorry, I need at least 3 bits of stencil.n");
  890.     exit(1);
  891.   }
  892.   /* Register GLUT callbacks. */
  893.   glutDisplayFunc(redraw);
  894.   glutMouseFunc(mouse);
  895.   glutMotionFunc(motion);
  896.   glutVisibilityFunc(visible);
  897.   glutKeyboardFunc(key);
  898.   glutSpecialFunc(special);
  899.   glutCreateMenu(controlLights);
  900.   glutAddMenuEntry("Toggle halo blending", M_BLENDED_HALO);
  901.   glutAddMenuEntry("Show halo", M_SHOW_HALO);
  902.   glutAddMenuEntry("Switch model", M_SWITCH_MODEL);
  903.   glutAddMenuEntry("Toggle motion", M_MOTION);
  904.   glutAddMenuEntry("-----------------------", M_NONE);
  905.   glutAddMenuEntry("Toggle light", M_LIGHT);
  906.   glutAddMenuEntry("Toggle texture", M_TEXTURE);
  907.   glutAddMenuEntry("Toggle shadows", M_SHADOWS);
  908.   glutAddMenuEntry("Toggle reflection", M_REFLECTION);
  909.   glutAddMenuEntry("Toggle object", M_DINOSAUR);
  910.   glutAddMenuEntry("-----------------------", M_NONE);
  911.   glutAddMenuEntry("Toggle reflection stenciling", M_STENCIL_REFLECTION);
  912.   glutAddMenuEntry("Toggle shadow stenciling", M_STENCIL_SHADOW);
  913.   glutAddMenuEntry("Toggle shadow offset", M_OFFSET_SHADOW);
  914.   glutAddMenuEntry("----------------------", M_NONE);
  915.   glutAddMenuEntry("Positional light", M_POSITIONAL);
  916.   glutAddMenuEntry("Directional light", M_DIRECTIONAL);
  917.   glutAddMenuEntry("-----------------------", M_NONE);
  918.   glutAddMenuEntry("Toggle performance", M_PERFORMANCE);
  919.   glutAttachMenu(GLUT_RIGHT_BUTTON);
  920.   makeDinosaur();
  921. #ifdef GL_VERSION_1_1
  922.   if (supportsOneDotOne() && !forceExtension) {
  923.     polygonOffsetVersion = ONE_DOT_ONE;
  924.     glPolygonOffset(-2.0, -1.0);
  925.   } else
  926. #endif
  927.   {
  928. #ifdef GL_EXT_polygon_offset
  929.   /* check for the polygon offset extension */
  930.   if (glutExtensionSupported("GL_EXT_polygon_offset")) {
  931.     polygonOffsetVersion = EXTENSION;
  932.     glPolygonOffsetEXT(-0.1, -0.002);
  933.   } else 
  934. #endif
  935.     {
  936.       polygonOffsetVersion = MISSING;
  937.       printf("ndinoshine: Missing polygon offset.n");
  938.       printf("           Expect shadow depth aliasing artifacts.nn");
  939.     }
  940.   }
  941.   glEnable(GL_CULL_FACE);
  942.   glEnable(GL_DEPTH_TEST);
  943.   glLineWidth(3.0);
  944.   glMatrixMode(GL_PROJECTION);
  945.   gluPerspective( /* field of view in degree */ 40.0,
  946.   /* aspect ratio */ 1.0,
  947.     /* Z near */ 20.0, /* Z far */ 100.0);
  948.   glMatrixMode(GL_MODELVIEW);
  949.   gluLookAt(0.0, 8.0, 60.0,  /* eye is at (0,0,30) */
  950.     0.0, 8.0, 0.0,      /* center is at (0,0,0) */
  951.     0.0, 1.0, 0.);      /* up is in postivie Y direction */
  952.   glLightModeli(GL_LIGHT_MODEL_LOCAL_VIEWER, 1);
  953.   glLightfv(GL_LIGHT0, GL_DIFFUSE, lightColor);
  954.   glLightf(GL_LIGHT0, GL_CONSTANT_ATTENUATION, 0.1);
  955.   glLightf(GL_LIGHT0, GL_LINEAR_ATTENUATION, 0.05);
  956.   glEnable(GL_LIGHT0);
  957.   glEnable(GL_LIGHTING);
  958.   makeFloorTexture();
  959.   /* Setup floor plane for projected shadow calculations. */
  960.   findPlane(floorPlane, floorVertices[1], floorVertices[2], floorVertices[3]);
  961.   glutMainLoop();
  962.   return 0;             /* ANSI C requires main to return int. */
  963. }