glu.c

Go to the documentation of this file.

#include <stdlib.h>
#include <math.h>
#include <GL/gl.h>
#include "glu.h"
 
 
void drawTorus(float rc, int numc, float rt, int numt)
{
    int i, j, k;
    double s, t;
    double x, y, z;
    double pi, twopi;
 
    pi = 3.14159265358979323846;
    twopi = 2 * pi;
 
    for (i = 0; i < numc; i++) {
    glBegin(GL_QUAD_STRIP);
        for (j = 0; j <= numt; j++) {
        for (k = 1; k >= 0; k--) {
        s = (i + k) % numc + 0.5;
        t = j % numt;
 
        x = cos(t*twopi/numt) * cos(s*twopi/numc);
        y = sin(t*twopi/numt) * cos(s*twopi/numc);
        z = sin(s*twopi/numc);
        glNormal3f(x, y, z);
 
        x = (rt + rc * cos(s*twopi/numc)) * cos(t*twopi/numt);
        y = (rt + rc * cos(s*twopi/numc)) * sin(t*twopi/numt);
        z = rc * sin(s*twopi/numc);
        glVertex3f(x, y, z);
        }
        }
    glEnd();
    }
}
 
static void normal3f( GLfloat x, GLfloat y, GLfloat z )
{
   GLdouble mag;
 
   mag = sqrt( x*x + y*y + z*z );
   if (mag>0.00001F) {
      x /= mag;
      y /= mag;
      z /= mag;
   }
   glNormal3f( x, y, z );
}
 
void gluPerspective( GLdouble fovy, GLdouble aspect,
             GLdouble zNear, GLdouble zFar )
{
   GLdouble xmin, xmax, ymin, ymax;
 
   ymax = zNear * tan( fovy * M_PI / 360.0 );
   ymin = -ymax;
 
   xmin = ymin * aspect;
   xmax = ymax * aspect;
 
   glFrustum( xmin, xmax, ymin, ymax, zNear, zFar );
}
 
GLUquadricObj *gluNewQuadric(void)
{
  return NULL;
}
 
void gluQuadricDrawStyle(GLUquadricObj *obj, int style)
{
}
 
void gluCylinder( GLUquadricObj *qobj,
                  GLdouble baseRadius, GLdouble topRadius, GLdouble height,
                  GLint slices, GLint stacks )
{
   GLdouble da, r, dr, dz;
   GLfloat z, nz, nsign;
   GLint i, j;
   GLfloat du = 1.0 / slices;
   GLfloat dv = 1.0 / stacks;
   GLfloat tcx = 0.0, tcy = 0.0;
   
   nsign = 1.0;
 
   da = 2.0*M_PI / slices;
   dr = (topRadius-baseRadius) / stacks;
   dz = height / stacks;
   nz = (baseRadius-topRadius) / height;  /* Z component of normal vectors */
 
   for (i=0;i<slices;i++) {
     GLfloat x1 = -sin(i*da);
     GLfloat y1 = cos(i*da);
     GLfloat x2 = -sin((i+1)*da);
     GLfloat y2 = cos((i+1)*da);
     z = 0.0;
     r = baseRadius;
     tcy = 0.0;
     glBegin( GL_QUAD_STRIP );
     for (j=0;j<=stacks;j++) {
        if (nsign==1.0) {
           normal3f( x1*nsign, y1*nsign, nz*nsign );
           glTexCoord2f(tcx, tcy);
           glVertex3f( x1*r, y1*r, z );
           normal3f( x2*nsign, y2*nsign, nz*nsign );
           glTexCoord2f(tcx+du, tcy);
           glVertex3f( x2*r, y2*r, z );
        }
        else {
           normal3f( x2*nsign, y2*nsign, nz*nsign );
           glTexCoord2f(tcx, tcy);
           glVertex3f( x2*r, y2*r, z );
           normal3f( x1*nsign, y1*nsign, nz*nsign );
           glTexCoord2f(tcx+du, tcy);
           glVertex3f( x1*r, y1*r, z );
        }
        z += dz;
        r += dr;
        tcy += dv;
     }
     glEnd();
     tcx += du;
      }
}
 
/* Disk (adapted from Mesa) */
 
void gluDisk( GLUquadricObj *qobj,
              GLdouble innerRadius, GLdouble outerRadius,
              GLint slices, GLint loops )
{
   GLdouble a, da;
   GLfloat dr;
   GLfloat r1, r2, dtc;
   GLint s, l;
   GLfloat sa,ca;
 
   /* Normal vectors */
     glNormal3f( 0.0, 0.0, +1.0 );
 
   da = 2.0*M_PI / slices;
   dr = (outerRadius-innerRadius) / (GLfloat) loops;
 
   /* texture of a gluDisk is a cut out of the texture unit square */
   /* x, y in [-outerRadius, +outerRadius]; s, t in [0, 1] (linear mapping) */
   dtc = 2.0f * outerRadius;
 
   r1 = innerRadius;
   for (l=0;l<loops;l++) {
        r2 = r1 + dr;
           glBegin( GL_QUAD_STRIP );
           for (s=0;s<=slices;s++) {
          if (s==slices) a = 0.0;
          else  a = s * da;
          sa = sin(a); ca = cos(a);
                  glTexCoord2f(0.5+sa*r2/dtc,0.5+ca*r2/dtc);
                  glVertex2f( r2*sa, r2*ca );
                  glTexCoord2f(0.5+sa*r1/dtc,0.5+ca*r1/dtc);
                  glVertex2f( r1*sa, r1*ca );
           }
           glEnd();
        r1 = r2;
     }
 
}
 
/*
 * Sphère (adapted from Mesa)
 */
 
void gluSphere(GLUquadricObj *qobj,
               float radius,int slices,int stacks)
{
   float rho, drho, theta, dtheta;
   float x, y, z;
   float s, t, ds, dt;
   int i, j, imin, imax;
   int normals;
   float nsign;
  
   normals=1;
   nsign=1;
 
   drho = M_PI / (float) stacks;
   dtheta = 2.0 * M_PI / (float) slices;
 
  /* draw +Z end as a triangle fan */
  glBegin( GL_TRIANGLE_FAN );
  glNormal3f( 0.0, 0.0, 1.0 );
  glTexCoord2f(0.5,0.0);
  glVertex3f( 0.0, 0.0, nsign * radius );
  for (j=0;j<=slices;j++) {
     theta = (j==slices) ? 0.0 : j * dtheta;
     x = -sin(theta) * sin(drho);
     y = cos(theta) * sin(drho);
     z = nsign * cos(drho);
     if (normals)  glNormal3f( x*nsign, y*nsign, z*nsign );
     glVertex3f( x*radius, y*radius, z*radius );
      }
   glEnd();
 
 
      ds = 1.0 / slices;
      dt = 1.0 / stacks;
      t = 1.0;  /* because loop now runs from 0 */
      if (1) {
        imin = 0;
        imax = stacks;
      }
      else {
        imin = 1;
        imax = stacks-1;
      }
 
      /* draw intermediate stacks as quad strips */
      for (i=imin;i<imax;i++) {
     rho = i * drho;
     glBegin( GL_QUAD_STRIP );
         s = 0.0;
     for (j=0;j<=slices;j++) {
        theta = (j==slices) ? 0.0 : j * dtheta;
        x = -sin(theta) * sin(rho);
        y = cos(theta) * sin(rho);
        z = nsign * cos(rho);
        if (normals)  glNormal3f( x*nsign, y*nsign, z*nsign );
        glTexCoord2f(s,1-t);
        glVertex3f( x*radius, y*radius, z*radius );
        x = -sin(theta) * sin(rho+drho);
        y = cos(theta) * sin(rho+drho);
        z = nsign * cos(rho+drho);
        if (normals)  glNormal3f( x*nsign, y*nsign, z*nsign );
        glTexCoord2f(s,1-(t-dt));
            s += ds;
        glVertex3f( x*radius, y*radius, z*radius );
     }
     glEnd();
     t -= dt;
      }
 
/* draw -Z end as a triangle fan */
    glBegin( GL_TRIANGLE_FAN );
    glNormal3f( 0.0, 0.0, -1.0 );
      glTexCoord2f(0.5,1.0);
      glVertex3f( 0.0, 0.0, -radius*nsign );
      rho = M_PI - drho;
      s = 1.0;
      t = dt;
      for (j=slices;j>=0;j--) {
     theta = (j==slices) ? 0.0 : j * dtheta;
     x = -sin(theta) * sin(rho);
     y = cos(theta) * sin(rho);
     z = nsign * cos(rho);
     if (normals)  glNormal3f( x*nsign, y*nsign, z*nsign );
     glTexCoord2f(s,1-t);
         s -= ds;
     glVertex3f( x*radius, y*radius, z*radius );
      }
      glEnd();
}