simgear/sky/sphere.cxx

   1 // sphere.cxx -- build an ssg sphere object
   2 //
   3 // Pulled straight out of MesaGLU/quadratic.c
   4 //
   5 // Original gluSphere code is Copyright (C) 1999-2000  Brian Paul and
   6 // licensed under the GPL
   7 //
   8 // This program is free software; you can redistribute it and/or
   9 // modify it under the terms of the GNU General Public License as
  10 // published by the Free Software Foundation; either version 2 of the
  11 // License, or (at your option) any later version.
  12 //
  13 // This program is distributed in the hope that it will be useful, but
  14 // WITHOUT ANY WARRANTY; without even the implied warranty of
  15 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  16 // General Public License for more details.
  17 //
  18 // You should have received a copy of the GNU General Public License
  19 // along with this program; if not, write to the Free Software
  20 // Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  21 //
  22 // $Id$
  23
  24
  25 #include <plib/ssg.h>
  26
  27
  28 // return a sphere object as an ssgBranch
  29 ssgBranch *ssgMakeSphere( double radius, int slices, int stacks ) {
  30     float rho, drho, theta, dtheta;
  31     float x, y, z;
  32     float s, t, ds, dt;
  33     int i, j, imin, imax;
  34     float nsign = 1.0;
  35     ssgBranch *sphere = new ssgBranch;
  36     sgVec2 vec2;
  37     sgVec3 vec3;
  38
  39     drho = M_PI / (float) stacks;
  40     dtheta = 2.0 * M_PI / (float) slices;
  41
  42     /* texturing: s goes from 0.0/0.25/0.5/0.75/1.0 at +y/+x/-y/-x/+y
  43        axis t goes from -1.0/+1.0 at z = -radius/+radius (linear along
  44        longitudes) cannot use triangle fan on texturing (s coord. at
  45        top/bottom tip varies) */
  46
  47     ds = 1.0 / slices;
  48     dt = 1.0 / stacks;
  49     t = 1.0;  /* because loop now runs from 0 */
  50     imin = 0;
  51     imax = stacks;
  52
  53     /* build slices as quad strips */
  54     for ( i = imin; i < imax; i++ ) {
  55         ssgVertexArray   *vl = new ssgVertexArray();
  56         ssgNormalArray   *nl = new ssgNormalArray();
  57         ssgTexCoordArray *tl = new ssgTexCoordArray();
  58
  59         rho = i * drho;
  60         s = 0.0;
  61         for ( j = 0; j <= slices; j++ ) {
  62             theta = (j == slices) ? 0.0 : j * dtheta;
  63             x = -sin(theta) * sin(rho);
  64             y = cos(theta) * sin(rho);
  65             z = nsign * cos(rho);
  66
  67             // glNormal3f( x*nsign, y*nsign, z*nsign );
  68             sgSetVec3( vec3, x*nsign, y*nsign, z*nsign );
  69             nl->add( vec3 );
  70
  71             // glTexCoord2f(s,t);
  72             sgSetVec2( vec2, 2, t );
  73             tl->add( vec2 );
  74
  75             // glVertex3f( x*radius, y*radius, z*radius );
  76             sgSetVec3( vec3, x*radius, y*radius, z*radius );
  77             vl->add( vec3 );
  78
  79             x = -sin(theta) * sin(rho+drho);
  80             y = cos(theta) * sin(rho+drho);
  81             z = nsign * cos(rho+drho);
  82
  83             // glNormal3f( x*nsign, y*nsign, z*nsign );
  84             sgSetVec3( vec3, x*nsign, y*nsign, z*nsign );
  85             nl->add( vec3 );
  86
  87             // glTexCoord2f(s,t-dt);
  88             sgSetVec2( vec2, s, t-dt );
  89             tl->add( vec2 );
  90             s += ds;
  91
  92             // glVertex3f( x*radius, y*radius, z*radius );
  93             sgSetVec3( vec3, x*radius, y*radius, z*radius );
  94             vl->add( vec3 );
  95         }
  96
  97         ssgLeaf *slice =
  98             new ssgVtxTable ( GL_TRIANGLE_STRIP, vl, nl, tl, NULL );
  99
 100         sphere->addKid( slice );
 101
 102         t -= dt;
 103     }
 104
 105     return sphere;
 106 }