2 * freeglut_geometry.c
\r
4 * Freeglut geometry rendering methods.
\r
6 * Copyright (c) 1999-2000 Pawel W. Olszta. All Rights Reserved.
\r
7 * Written by Pawel W. Olszta, <olszta@sourceforge.net>
\r
8 * Creation date: Fri Dec 3 1999
\r
10 * Permission is hereby granted, free of charge, to any person obtaining a
\r
11 * copy of this software and associated documentation files (the "Software"),
\r
12 * to deal in the Software without restriction, including without limitation
\r
13 * the rights to use, copy, modify, merge, publish, distribute, sublicense,
\r
14 * and/or sell copies of the Software, and to permit persons to whom the
\r
15 * Software is furnished to do so, subject to the following conditions:
\r
17 * The above copyright notice and this permission notice shall be included
\r
18 * in all copies or substantial portions of the Software.
\r
20 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
\r
21 * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
\r
22 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
\r
23 * PAWEL W. OLSZTA BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER
\r
24 * IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN
\r
25 * CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
\r
28 #include <GL/freeglut.h>
\r
29 #include "freeglut_internal.h"
\r
32 * TODO BEFORE THE STABLE RELEASE:
\r
34 * Following functions have been contributed by Andreas Umbach.
\r
36 * glutWireCube() -- looks OK
\r
37 * glutSolidCube() -- OK
\r
39 * Those functions have been implemented by John Fay.
\r
41 * glutWireTorus() -- looks OK
\r
42 * glutSolidTorus() -- looks OK
\r
43 * glutWireDodecahedron() -- looks OK
\r
44 * glutSolidDodecahedron() -- looks OK
\r
45 * glutWireOctahedron() -- looks OK
\r
46 * glutSolidOctahedron() -- looks OK
\r
47 * glutWireTetrahedron() -- looks OK
\r
48 * glutSolidTetrahedron() -- looks OK
\r
49 * glutWireIcosahedron() -- looks OK
\r
50 * glutSolidIcosahedron() -- looks OK
\r
52 * The Following functions have been updated by Nigel Stewart, based
\r
53 * on FreeGLUT 2.0.0 implementations:
\r
55 * glutWireSphere() -- looks OK
\r
56 * glutSolidSphere() -- looks OK
\r
57 * glutWireCone() -- looks OK
\r
58 * glutSolidCone() -- looks OK
\r
62 /* -- INTERFACE FUNCTIONS -------------------------------------------------- */
\r
65 * Draws a wireframed cube. Code contributed by Andreas Umbach <marvin@dataway.ch>
\r
67 void FGAPIENTRY glutWireCube( GLdouble dSize )
\r
69 double size = dSize * 0.5;
\r
71 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutWireCube" );
\r
73 # define V(a,b,c) glVertex3d( a size, b size, c size );
\r
74 # define N(a,b,c) glNormal3d( a, b, c );
\r
76 /* PWO: I dared to convert the code to use macros... */
\r
77 glBegin( GL_LINE_LOOP ); N( 1.0, 0.0, 0.0); V(+,-,+); V(+,-,-); V(+,+,-); V(+,+,+); glEnd();
\r
78 glBegin( GL_LINE_LOOP ); N( 0.0, 1.0, 0.0); V(+,+,+); V(+,+,-); V(-,+,-); V(-,+,+); glEnd();
\r
79 glBegin( GL_LINE_LOOP ); N( 0.0, 0.0, 1.0); V(+,+,+); V(-,+,+); V(-,-,+); V(+,-,+); glEnd();
\r
80 glBegin( GL_LINE_LOOP ); N(-1.0, 0.0, 0.0); V(-,-,+); V(-,+,+); V(-,+,-); V(-,-,-); glEnd();
\r
81 glBegin( GL_LINE_LOOP ); N( 0.0,-1.0, 0.0); V(-,-,+); V(-,-,-); V(+,-,-); V(+,-,+); glEnd();
\r
82 glBegin( GL_LINE_LOOP ); N( 0.0, 0.0,-1.0); V(-,-,-); V(-,+,-); V(+,+,-); V(+,-,-); glEnd();
\r
89 * Draws a solid cube. Code contributed by Andreas Umbach <marvin@dataway.ch>
\r
91 void FGAPIENTRY glutSolidCube( GLdouble dSize )
\r
93 double size = dSize * 0.5;
\r
95 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutSolidCube" );
\r
97 # define V(a,b,c) glVertex3d( a size, b size, c size );
\r
98 # define N(a,b,c) glNormal3d( a, b, c );
\r
100 /* PWO: Again, I dared to convert the code to use macros... */
\r
101 glBegin( GL_QUADS );
\r
102 N( 1.0, 0.0, 0.0); V(+,-,+); V(+,-,-); V(+,+,-); V(+,+,+);
\r
103 N( 0.0, 1.0, 0.0); V(+,+,+); V(+,+,-); V(-,+,-); V(-,+,+);
\r
104 N( 0.0, 0.0, 1.0); V(+,+,+); V(-,+,+); V(-,-,+); V(+,-,+);
\r
105 N(-1.0, 0.0, 0.0); V(-,-,+); V(-,+,+); V(-,+,-); V(-,-,-);
\r
106 N( 0.0,-1.0, 0.0); V(-,-,+); V(-,-,-); V(+,-,-); V(+,-,+);
\r
107 N( 0.0, 0.0,-1.0); V(-,-,-); V(-,+,-); V(+,+,-); V(+,-,-);
\r
115 * Compute lookup table of cos and sin values forming a cirle
\r
118 * It is the responsibility of the caller to free these tables
\r
119 * The size of the table is (n+1) to form a connected loop
\r
120 * The last entry is exactly the same as the first
\r
121 * The sign of n can be flipped to get the reverse loop
\r
124 static void fghCircleTable(double **sint,double **cost,const int n)
\r
128 /* Table size, the sign of n flips the circle direction */
\r
130 const int size = abs(n);
\r
132 /* Determine the angle between samples */
\r
134 const double angle = 2*M_PI/(double)( ( n == 0 ) ? 1 : n );
\r
136 /* Allocate memory for n samples, plus duplicate of first entry at the end */
\r
138 *sint = (double *) calloc(sizeof(double), size+1);
\r
139 *cost = (double *) calloc(sizeof(double), size+1);
\r
141 /* Bail out if memory allocation fails, fgError never returns */
\r
143 if (!(*sint) || !(*cost))
\r
147 fgError("Failed to allocate memory in fghCircleTable");
\r
150 /* Compute cos and sin around the circle */
\r
155 for (i=1; i<size; i++)
\r
157 (*sint)[i] = sin(angle*i);
\r
158 (*cost)[i] = cos(angle*i);
\r
161 /* Last sample is duplicate of the first */
\r
163 (*sint)[size] = (*sint)[0];
\r
164 (*cost)[size] = (*cost)[0];
\r
168 * Draws a solid sphere
\r
170 void FGAPIENTRY glutSolidSphere(GLdouble radius, GLint slices, GLint stacks)
\r
174 /* Adjust z and radius as stacks are drawn. */
\r
179 /* Pre-computed circle */
\r
181 double *sint1,*cost1;
\r
182 double *sint2,*cost2;
\r
184 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutSolidSphere" );
\r
186 fghCircleTable(&sint1,&cost1,-slices);
\r
187 fghCircleTable(&sint2,&cost2,stacks*2);
\r
189 /* The top stack is covered with a triangle fan */
\r
192 z1 = cost2[(stacks>0)?1:0];
\r
194 r1 = sint2[(stacks>0)?1:0];
\r
196 glBegin(GL_TRIANGLE_FAN);
\r
199 glVertex3d(0,0,radius);
\r
201 for (j=slices; j>=0; j--)
\r
203 glNormal3d(cost1[j]*r1, sint1[j]*r1, z1 );
\r
204 glVertex3d(cost1[j]*r1*radius, sint1[j]*r1*radius, z1*radius);
\r
209 /* Cover each stack with a quad strip, except the top and bottom stacks */
\r
211 for( i=1; i<stacks-1; i++ )
\r
213 z0 = z1; z1 = cost2[i+1];
\r
214 r0 = r1; r1 = sint2[i+1];
\r
216 glBegin(GL_QUAD_STRIP);
\r
218 for(j=0; j<=slices; j++)
\r
220 glNormal3d(cost1[j]*r1, sint1[j]*r1, z1 );
\r
221 glVertex3d(cost1[j]*r1*radius, sint1[j]*r1*radius, z1*radius);
\r
222 glNormal3d(cost1[j]*r0, sint1[j]*r0, z0 );
\r
223 glVertex3d(cost1[j]*r0*radius, sint1[j]*r0*radius, z0*radius);
\r
229 /* The bottom stack is covered with a triangle fan */
\r
234 glBegin(GL_TRIANGLE_FAN);
\r
236 glNormal3d(0,0,-1);
\r
237 glVertex3d(0,0,-radius);
\r
239 for (j=0; j<=slices; j++)
\r
241 glNormal3d(cost1[j]*r0, sint1[j]*r0, z0 );
\r
242 glVertex3d(cost1[j]*r0*radius, sint1[j]*r0*radius, z0*radius);
\r
247 /* Release sin and cos tables */
\r
256 * Draws a wire sphere
\r
258 void FGAPIENTRY glutWireSphere(GLdouble radius, GLint slices, GLint stacks)
\r
262 /* Adjust z and radius as stacks and slices are drawn. */
\r
267 /* Pre-computed circle */
\r
269 double *sint1,*cost1;
\r
270 double *sint2,*cost2;
\r
272 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutWireSphere" );
\r
274 fghCircleTable(&sint1,&cost1,-slices );
\r
275 fghCircleTable(&sint2,&cost2, stacks*2);
\r
277 /* Draw a line loop for each stack */
\r
279 for (i=1; i<stacks; i++)
\r
284 glBegin(GL_LINE_LOOP);
\r
286 for(j=0; j<=slices; j++)
\r
292 glVertex3d(x*r*radius,y*r*radius,z*radius);
\r
298 /* Draw a line loop for each slice */
\r
300 for (i=0; i<slices; i++)
\r
302 glBegin(GL_LINE_STRIP);
\r
304 for(j=0; j<=stacks; j++)
\r
306 x = cost1[i]*sint2[j];
\r
307 y = sint1[i]*sint2[j];
\r
311 glVertex3d(x*radius,y*radius,z*radius);
\r
317 /* Release sin and cos tables */
\r
326 * Draws a solid cone
\r
328 void FGAPIENTRY glutSolidCone( GLdouble base, GLdouble height, GLint slices, GLint stacks )
\r
332 /* Step in z and radius as stacks are drawn. */
\r
337 const double zStep = height / ( ( stacks > 0 ) ? stacks : 1 );
\r
338 const double rStep = base / ( ( stacks > 0 ) ? stacks : 1 );
\r
340 /* Scaling factors for vertex normals */
\r
342 const double cosn = ( height / sqrt ( height * height + base * base ));
\r
343 const double sinn = ( base / sqrt ( height * height + base * base ));
\r
345 /* Pre-computed circle */
\r
347 double *sint,*cost;
\r
349 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutSolidCone" );
\r
351 fghCircleTable(&sint,&cost,-slices);
\r
353 /* Cover the circular base with a triangle fan... */
\r
361 glBegin(GL_TRIANGLE_FAN);
\r
363 glNormal3d(0.0,0.0,-1.0);
\r
364 glVertex3d(0.0,0.0, z0 );
\r
366 for (j=0; j<=slices; j++)
\r
367 glVertex3d(cost[j]*r0, sint[j]*r0, z0);
\r
371 /* Cover each stack with a quad strip, except the top stack */
\r
373 for( i=0; i<stacks-1; i++ )
\r
375 glBegin(GL_QUAD_STRIP);
\r
377 for(j=0; j<=slices; j++)
\r
379 glNormal3d(cost[j]*cosn, sint[j]*cosn, sinn);
\r
380 glVertex3d(cost[j]*r0, sint[j]*r0, z0 );
\r
381 glVertex3d(cost[j]*r1, sint[j]*r1, z1 );
\r
384 z0 = z1; z1 += zStep;
\r
385 r0 = r1; r1 -= rStep;
\r
390 /* The top stack is covered with individual triangles */
\r
392 glBegin(GL_TRIANGLES);
\r
394 glNormal3d(cost[0]*sinn, sint[0]*sinn, cosn);
\r
396 for (j=0; j<slices; j++)
\r
398 glVertex3d(cost[j+0]*r0, sint[j+0]*r0, z0 );
\r
399 glVertex3d(0, 0, height);
\r
400 glNormal3d(cost[j+1]*sinn, sint[j+1]*sinn, cosn );
\r
401 glVertex3d(cost[j+1]*r0, sint[j+1]*r0, z0 );
\r
406 /* Release sin and cos tables */
\r
413 * Draws a wire cone
\r
415 void FGAPIENTRY glutWireCone( GLdouble base, GLdouble height, GLint slices, GLint stacks)
\r
419 /* Step in z and radius as stacks are drawn. */
\r
424 const double zStep = height / ( ( stacks > 0 ) ? stacks : 1 );
\r
425 const double rStep = base / ( ( stacks > 0 ) ? stacks : 1 );
\r
427 /* Scaling factors for vertex normals */
\r
429 const double cosn = ( height / sqrt ( height * height + base * base ));
\r
430 const double sinn = ( base / sqrt ( height * height + base * base ));
\r
432 /* Pre-computed circle */
\r
434 double *sint,*cost;
\r
436 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutWireCone" );
\r
438 fghCircleTable(&sint,&cost,-slices);
\r
440 /* Draw the stacks... */
\r
442 for (i=0; i<stacks; i++)
\r
444 glBegin(GL_LINE_LOOP);
\r
446 for( j=0; j<slices; j++ )
\r
448 glNormal3d(cost[j]*sinn, sint[j]*sinn, cosn);
\r
449 glVertex3d(cost[j]*r, sint[j]*r, z );
\r
458 /* Draw the slices */
\r
464 for (j=0; j<slices; j++)
\r
466 glNormal3d(cost[j]*sinn, sint[j]*sinn, cosn );
\r
467 glVertex3d(cost[j]*r, sint[j]*r, 0.0 );
\r
468 glVertex3d(0.0, 0.0, height);
\r
473 /* Release sin and cos tables */
\r
481 * Draws a solid cylinder
\r
483 void FGAPIENTRY glutSolidCylinder(GLdouble radius, GLdouble height, GLint slices, GLint stacks)
\r
487 /* Step in z and radius as stacks are drawn. */
\r
490 const double zStep = height / ( ( stacks > 0 ) ? stacks : 1 );
\r
492 /* Pre-computed circle */
\r
494 double *sint,*cost;
\r
496 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutSolidCylinder" );
\r
498 fghCircleTable(&sint,&cost,-slices);
\r
500 /* Cover the base and top */
\r
502 glBegin(GL_TRIANGLE_FAN);
\r
503 glNormal3d(0.0, 0.0, -1.0 );
\r
504 glVertex3d(0.0, 0.0, 0.0 );
\r
505 for (j=0; j<=slices; j++)
\r
506 glVertex3d(cost[j]*radius, sint[j]*radius, 0.0);
\r
509 glBegin(GL_TRIANGLE_FAN);
\r
510 glNormal3d(0.0, 0.0, 1.0 );
\r
511 glVertex3d(0.0, 0.0, height);
\r
512 for (j=slices; j>=0; j--)
\r
513 glVertex3d(cost[j]*radius, sint[j]*radius, height);
\r
516 /* Do the stacks */
\r
521 for (i=1; i<=stacks; i++)
\r
526 glBegin(GL_QUAD_STRIP);
\r
527 for (j=0; j<=slices; j++ )
\r
529 glNormal3d(cost[j], sint[j], 0.0 );
\r
530 glVertex3d(cost[j]*radius, sint[j]*radius, z0 );
\r
531 glVertex3d(cost[j]*radius, sint[j]*radius, z1 );
\r
535 z0 = z1; z1 += zStep;
\r
538 /* Release sin and cos tables */
\r
545 * Draws a wire cylinder
\r
547 void FGAPIENTRY glutWireCylinder(GLdouble radius, GLdouble height, GLint slices, GLint stacks)
\r
551 /* Step in z and radius as stacks are drawn. */
\r
554 const double zStep = height / ( ( stacks > 0 ) ? stacks : 1 );
\r
556 /* Pre-computed circle */
\r
558 double *sint,*cost;
\r
560 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutWireCylinder" );
\r
562 fghCircleTable(&sint,&cost,-slices);
\r
564 /* Draw the stacks... */
\r
566 for (i=0; i<=stacks; i++)
\r
571 glBegin(GL_LINE_LOOP);
\r
573 for( j=0; j<slices; j++ )
\r
575 glNormal3d(cost[j], sint[j], 0.0);
\r
576 glVertex3d(cost[j]*radius, sint[j]*radius, z );
\r
584 /* Draw the slices */
\r
588 for (j=0; j<slices; j++)
\r
590 glNormal3d(cost[j], sint[j], 0.0 );
\r
591 glVertex3d(cost[j]*radius, sint[j]*radius, 0.0 );
\r
592 glVertex3d(cost[j]*radius, sint[j]*radius, height);
\r
597 /* Release sin and cos tables */
\r
604 * Draws a wire torus
\r
606 void FGAPIENTRY glutWireTorus( GLdouble dInnerRadius, GLdouble dOuterRadius, GLint nSides, GLint nRings )
\r
608 double iradius = dInnerRadius, oradius = dOuterRadius, phi, psi, dpsi, dphi;
\r
609 double *vertex, *normal;
\r
611 double spsi, cpsi, sphi, cphi ;
\r
613 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutWireTorus" );
\r
615 if ( nSides < 1 ) nSides = 1;
\r
616 if ( nRings < 1 ) nRings = 1;
\r
618 /* Allocate the vertices array */
\r
619 vertex = (double *)calloc( sizeof(double), 3 * nSides * nRings );
\r
620 normal = (double *)calloc( sizeof(double), 3 * nSides * nRings );
\r
624 dpsi = 2.0 * M_PI / (double)nRings ;
\r
625 dphi = -2.0 * M_PI / (double)nSides ;
\r
628 for( j=0; j<nRings; j++ )
\r
630 cpsi = cos ( psi ) ;
\r
631 spsi = sin ( psi ) ;
\r
634 for( i=0; i<nSides; i++ )
\r
636 int offset = 3 * ( j * nSides + i ) ;
\r
637 cphi = cos ( phi ) ;
\r
638 sphi = sin ( phi ) ;
\r
639 *(vertex + offset + 0) = cpsi * ( oradius + cphi * iradius ) ;
\r
640 *(vertex + offset + 1) = spsi * ( oradius + cphi * iradius ) ;
\r
641 *(vertex + offset + 2) = sphi * iradius ;
\r
642 *(normal + offset + 0) = cpsi * cphi ;
\r
643 *(normal + offset + 1) = spsi * cphi ;
\r
644 *(normal + offset + 2) = sphi ;
\r
651 for( i=0; i<nSides; i++ )
\r
653 glBegin( GL_LINE_LOOP );
\r
655 for( j=0; j<nRings; j++ )
\r
657 int offset = 3 * ( j * nSides + i ) ;
\r
658 glNormal3dv( normal + offset );
\r
659 glVertex3dv( vertex + offset );
\r
665 for( j=0; j<nRings; j++ )
\r
667 glBegin(GL_LINE_LOOP);
\r
669 for( i=0; i<nSides; i++ )
\r
671 int offset = 3 * ( j * nSides + i ) ;
\r
672 glNormal3dv( normal + offset );
\r
673 glVertex3dv( vertex + offset );
\r
685 * Draws a solid torus
\r
687 void FGAPIENTRY glutSolidTorus( GLdouble dInnerRadius, GLdouble dOuterRadius, GLint nSides, GLint nRings )
\r
689 double iradius = dInnerRadius, oradius = dOuterRadius, phi, psi, dpsi, dphi;
\r
690 double *vertex, *normal;
\r
692 double spsi, cpsi, sphi, cphi ;
\r
694 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutSolidTorus" );
\r
696 if ( nSides < 1 ) nSides = 1;
\r
697 if ( nRings < 1 ) nRings = 1;
\r
699 /* Increment the number of sides and rings to allow for one more point than surface */
\r
703 /* Allocate the vertices array */
\r
704 vertex = (double *)calloc( sizeof(double), 3 * nSides * nRings );
\r
705 normal = (double *)calloc( sizeof(double), 3 * nSides * nRings );
\r
709 dpsi = 2.0 * M_PI / (double)(nRings - 1) ;
\r
710 dphi = -2.0 * M_PI / (double)(nSides - 1) ;
\r
713 for( j=0; j<nRings; j++ )
\r
715 cpsi = cos ( psi ) ;
\r
716 spsi = sin ( psi ) ;
\r
719 for( i=0; i<nSides; i++ )
\r
721 int offset = 3 * ( j * nSides + i ) ;
\r
722 cphi = cos ( phi ) ;
\r
723 sphi = sin ( phi ) ;
\r
724 *(vertex + offset + 0) = cpsi * ( oradius + cphi * iradius ) ;
\r
725 *(vertex + offset + 1) = spsi * ( oradius + cphi * iradius ) ;
\r
726 *(vertex + offset + 2) = sphi * iradius ;
\r
727 *(normal + offset + 0) = cpsi * cphi ;
\r
728 *(normal + offset + 1) = spsi * cphi ;
\r
729 *(normal + offset + 2) = sphi ;
\r
736 glBegin( GL_QUADS );
\r
737 for( i=0; i<nSides-1; i++ )
\r
739 for( j=0; j<nRings-1; j++ )
\r
741 int offset = 3 * ( j * nSides + i ) ;
\r
742 glNormal3dv( normal + offset );
\r
743 glVertex3dv( vertex + offset );
\r
744 glNormal3dv( normal + offset + 3 );
\r
745 glVertex3dv( vertex + offset + 3 );
\r
746 glNormal3dv( normal + offset + 3 * nSides + 3 );
\r
747 glVertex3dv( vertex + offset + 3 * nSides + 3 );
\r
748 glNormal3dv( normal + offset + 3 * nSides );
\r
749 glVertex3dv( vertex + offset + 3 * nSides );
\r
763 void FGAPIENTRY glutWireDodecahedron( void )
\r
765 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutWireDodecahedron" );
\r
767 /* Magic Numbers: It is possible to create a dodecahedron by attaching two pentagons to each face of
\r
768 * of a cube. The coordinates of the points are:
\r
769 * (+-x,0, z); (+-1, 1, 1); (0, z, x )
\r
770 * where x = (-1 + sqrt(5))/2, z = (1 + sqrt(5))/2 or
\r
771 * x = 0.61803398875 and z = 1.61803398875.
\r
773 glBegin ( GL_LINE_LOOP ) ;
\r
774 glNormal3d ( 0.0, 0.525731112119, 0.850650808354 ) ; glVertex3d ( 0.0, 1.61803398875, 0.61803398875 ) ; glVertex3d ( -1.0, 1.0, 1.0 ) ; glVertex3d ( -0.61803398875, 0.0, 1.61803398875 ) ; glVertex3d ( 0.61803398875, 0.0, 1.61803398875 ) ; glVertex3d ( 1.0, 1.0, 1.0 ) ;
\r
776 glBegin ( GL_LINE_LOOP ) ;
\r
777 glNormal3d ( 0.0, 0.525731112119, -0.850650808354 ) ; glVertex3d ( 0.0, 1.61803398875, -0.61803398875 ) ; glVertex3d ( 1.0, 1.0, -1.0 ) ; glVertex3d ( 0.61803398875, 0.0, -1.61803398875 ) ; glVertex3d ( -0.61803398875, 0.0, -1.61803398875 ) ; glVertex3d ( -1.0, 1.0, -1.0 ) ;
\r
779 glBegin ( GL_LINE_LOOP ) ;
\r
780 glNormal3d ( 0.0, -0.525731112119, 0.850650808354 ) ; glVertex3d ( 0.0, -1.61803398875, 0.61803398875 ) ; glVertex3d ( 1.0, -1.0, 1.0 ) ; glVertex3d ( 0.61803398875, 0.0, 1.61803398875 ) ; glVertex3d ( -0.61803398875, 0.0, 1.61803398875 ) ; glVertex3d ( -1.0, -1.0, 1.0 ) ;
\r
782 glBegin ( GL_LINE_LOOP ) ;
\r
783 glNormal3d ( 0.0, -0.525731112119, -0.850650808354 ) ; glVertex3d ( 0.0, -1.61803398875, -0.61803398875 ) ; glVertex3d ( -1.0, -1.0, -1.0 ) ; glVertex3d ( -0.61803398875, 0.0, -1.61803398875 ) ; glVertex3d ( 0.61803398875, 0.0, -1.61803398875 ) ; glVertex3d ( 1.0, -1.0, -1.0 ) ;
\r
786 glBegin ( GL_LINE_LOOP ) ;
\r
787 glNormal3d ( 0.850650808354, 0.0, 0.525731112119 ) ; glVertex3d ( 0.61803398875, 0.0, 1.61803398875 ) ; glVertex3d ( 1.0, -1.0, 1.0 ) ; glVertex3d ( 1.61803398875, -0.61803398875, 0.0 ) ; glVertex3d ( 1.61803398875, 0.61803398875, 0.0 ) ; glVertex3d ( 1.0, 1.0, 1.0 ) ;
\r
789 glBegin ( GL_LINE_LOOP ) ;
\r
790 glNormal3d ( -0.850650808354, 0.0, 0.525731112119 ) ; glVertex3d ( -0.61803398875, 0.0, 1.61803398875 ) ; glVertex3d ( -1.0, 1.0, 1.0 ) ; glVertex3d ( -1.61803398875, 0.61803398875, 0.0 ) ; glVertex3d ( -1.61803398875, -0.61803398875, 0.0 ) ; glVertex3d ( -1.0, -1.0, 1.0 ) ;
\r
792 glBegin ( GL_LINE_LOOP ) ;
\r
793 glNormal3d ( 0.850650808354, 0.0, -0.525731112119 ) ; glVertex3d ( 0.61803398875, 0.0, -1.61803398875 ) ; glVertex3d ( 1.0, 1.0, -1.0 ) ; glVertex3d ( 1.61803398875, 0.61803398875, 0.0 ) ; glVertex3d ( 1.61803398875, -0.61803398875, 0.0 ) ; glVertex3d ( 1.0, -1.0, -1.0 ) ;
\r
795 glBegin ( GL_LINE_LOOP ) ;
\r
796 glNormal3d ( -0.850650808354, 0.0, -0.525731112119 ) ; glVertex3d ( -0.61803398875, 0.0, -1.61803398875 ) ; glVertex3d ( -1.0, -1.0, -1.0 ) ; glVertex3d ( -1.61803398875, -0.61803398875, 0.0 ) ; glVertex3d ( -1.61803398875, 0.61803398875, 0.0 ) ; glVertex3d ( -1.0, 1.0, -1.0 ) ;
\r
799 glBegin ( GL_LINE_LOOP ) ;
\r
800 glNormal3d ( 0.525731112119, 0.850650808354, 0.0 ) ; glVertex3d ( 1.61803398875, 0.61803398875, 0.0 ) ; glVertex3d ( 1.0, 1.0, -1.0 ) ; glVertex3d ( 0.0, 1.61803398875, -0.61803398875 ) ; glVertex3d ( 0.0, 1.61803398875, 0.61803398875 ) ; glVertex3d ( 1.0, 1.0, 1.0 ) ;
\r
802 glBegin ( GL_LINE_LOOP ) ;
\r
803 glNormal3d ( 0.525731112119, -0.850650808354, 0.0 ) ; glVertex3d ( 1.61803398875, -0.61803398875, 0.0 ) ; glVertex3d ( 1.0, -1.0, 1.0 ) ; glVertex3d ( 0.0, -1.61803398875, 0.61803398875 ) ; glVertex3d ( 0.0, -1.61803398875, -0.61803398875 ) ; glVertex3d ( 1.0, -1.0, -1.0 ) ;
\r
805 glBegin ( GL_LINE_LOOP ) ;
\r
806 glNormal3d ( -0.525731112119, 0.850650808354, 0.0 ) ; glVertex3d ( -1.61803398875, 0.61803398875, 0.0 ) ; glVertex3d ( -1.0, 1.0, 1.0 ) ; glVertex3d ( 0.0, 1.61803398875, 0.61803398875 ) ; glVertex3d ( 0.0, 1.61803398875, -0.61803398875 ) ; glVertex3d ( -1.0, 1.0, -1.0 ) ;
\r
808 glBegin ( GL_LINE_LOOP ) ;
\r
809 glNormal3d ( -0.525731112119, -0.850650808354, 0.0 ) ; glVertex3d ( -1.61803398875, -0.61803398875, 0.0 ) ; glVertex3d ( -1.0, -1.0, -1.0 ) ; glVertex3d ( 0.0, -1.61803398875, -0.61803398875 ) ; glVertex3d ( 0.0, -1.61803398875, 0.61803398875 ) ; glVertex3d ( -1.0, -1.0, 1.0 ) ;
\r
816 void FGAPIENTRY glutSolidDodecahedron( void )
\r
818 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutSolidDodecahedron" );
\r
820 /* Magic Numbers: It is possible to create a dodecahedron by attaching two pentagons to each face of
\r
821 * of a cube. The coordinates of the points are:
\r
822 * (+-x,0, z); (+-1, 1, 1); (0, z, x )
\r
823 * where x = (-1 + sqrt(5))/2, z = (1 + sqrt(5))/2 or
\r
824 * x = 0.61803398875 and z = 1.61803398875.
\r
826 glBegin ( GL_POLYGON ) ;
\r
827 glNormal3d ( 0.0, 0.525731112119, 0.850650808354 ) ; glVertex3d ( 0.0, 1.61803398875, 0.61803398875 ) ; glVertex3d ( -1.0, 1.0, 1.0 ) ; glVertex3d ( -0.61803398875, 0.0, 1.61803398875 ) ; glVertex3d ( 0.61803398875, 0.0, 1.61803398875 ) ; glVertex3d ( 1.0, 1.0, 1.0 ) ;
\r
829 glBegin ( GL_POLYGON ) ;
\r
830 glNormal3d ( 0.0, 0.525731112119, -0.850650808354 ) ; glVertex3d ( 0.0, 1.61803398875, -0.61803398875 ) ; glVertex3d ( 1.0, 1.0, -1.0 ) ; glVertex3d ( 0.61803398875, 0.0, -1.61803398875 ) ; glVertex3d ( -0.61803398875, 0.0, -1.61803398875 ) ; glVertex3d ( -1.0, 1.0, -1.0 ) ;
\r
832 glBegin ( GL_POLYGON ) ;
\r
833 glNormal3d ( 0.0, -0.525731112119, 0.850650808354 ) ; glVertex3d ( 0.0, -1.61803398875, 0.61803398875 ) ; glVertex3d ( 1.0, -1.0, 1.0 ) ; glVertex3d ( 0.61803398875, 0.0, 1.61803398875 ) ; glVertex3d ( -0.61803398875, 0.0, 1.61803398875 ) ; glVertex3d ( -1.0, -1.0, 1.0 ) ;
\r
835 glBegin ( GL_POLYGON ) ;
\r
836 glNormal3d ( 0.0, -0.525731112119, -0.850650808354 ) ; glVertex3d ( 0.0, -1.61803398875, -0.61803398875 ) ; glVertex3d ( -1.0, -1.0, -1.0 ) ; glVertex3d ( -0.61803398875, 0.0, -1.61803398875 ) ; glVertex3d ( 0.61803398875, 0.0, -1.61803398875 ) ; glVertex3d ( 1.0, -1.0, -1.0 ) ;
\r
839 glBegin ( GL_POLYGON ) ;
\r
840 glNormal3d ( 0.850650808354, 0.0, 0.525731112119 ) ; glVertex3d ( 0.61803398875, 0.0, 1.61803398875 ) ; glVertex3d ( 1.0, -1.0, 1.0 ) ; glVertex3d ( 1.61803398875, -0.61803398875, 0.0 ) ; glVertex3d ( 1.61803398875, 0.61803398875, 0.0 ) ; glVertex3d ( 1.0, 1.0, 1.0 ) ;
\r
842 glBegin ( GL_POLYGON ) ;
\r
843 glNormal3d ( -0.850650808354, 0.0, 0.525731112119 ) ; glVertex3d ( -0.61803398875, 0.0, 1.61803398875 ) ; glVertex3d ( -1.0, 1.0, 1.0 ) ; glVertex3d ( -1.61803398875, 0.61803398875, 0.0 ) ; glVertex3d ( -1.61803398875, -0.61803398875, 0.0 ) ; glVertex3d ( -1.0, -1.0, 1.0 ) ;
\r
845 glBegin ( GL_POLYGON ) ;
\r
846 glNormal3d ( 0.850650808354, 0.0, -0.525731112119 ) ; glVertex3d ( 0.61803398875, 0.0, -1.61803398875 ) ; glVertex3d ( 1.0, 1.0, -1.0 ) ; glVertex3d ( 1.61803398875, 0.61803398875, 0.0 ) ; glVertex3d ( 1.61803398875, -0.61803398875, 0.0 ) ; glVertex3d ( 1.0, -1.0, -1.0 ) ;
\r
848 glBegin ( GL_POLYGON ) ;
\r
849 glNormal3d ( -0.850650808354, 0.0, -0.525731112119 ) ; glVertex3d ( -0.61803398875, 0.0, -1.61803398875 ) ; glVertex3d ( -1.0, -1.0, -1.0 ) ; glVertex3d ( -1.61803398875, -0.61803398875, 0.0 ) ; glVertex3d ( -1.61803398875, 0.61803398875, 0.0 ) ; glVertex3d ( -1.0, 1.0, -1.0 ) ;
\r
852 glBegin ( GL_POLYGON ) ;
\r
853 glNormal3d ( 0.525731112119, 0.850650808354, 0.0 ) ; glVertex3d ( 1.61803398875, 0.61803398875, 0.0 ) ; glVertex3d ( 1.0, 1.0, -1.0 ) ; glVertex3d ( 0.0, 1.61803398875, -0.61803398875 ) ; glVertex3d ( 0.0, 1.61803398875, 0.61803398875 ) ; glVertex3d ( 1.0, 1.0, 1.0 ) ;
\r
855 glBegin ( GL_POLYGON ) ;
\r
856 glNormal3d ( 0.525731112119, -0.850650808354, 0.0 ) ; glVertex3d ( 1.61803398875, -0.61803398875, 0.0 ) ; glVertex3d ( 1.0, -1.0, 1.0 ) ; glVertex3d ( 0.0, -1.61803398875, 0.61803398875 ) ; glVertex3d ( 0.0, -1.61803398875, -0.61803398875 ) ; glVertex3d ( 1.0, -1.0, -1.0 ) ;
\r
858 glBegin ( GL_POLYGON ) ;
\r
859 glNormal3d ( -0.525731112119, 0.850650808354, 0.0 ) ; glVertex3d ( -1.61803398875, 0.61803398875, 0.0 ) ; glVertex3d ( -1.0, 1.0, 1.0 ) ; glVertex3d ( 0.0, 1.61803398875, 0.61803398875 ) ; glVertex3d ( 0.0, 1.61803398875, -0.61803398875 ) ; glVertex3d ( -1.0, 1.0, -1.0 ) ;
\r
861 glBegin ( GL_POLYGON ) ;
\r
862 glNormal3d ( -0.525731112119, -0.850650808354, 0.0 ) ; glVertex3d ( -1.61803398875, -0.61803398875, 0.0 ) ; glVertex3d ( -1.0, -1.0, -1.0 ) ; glVertex3d ( 0.0, -1.61803398875, -0.61803398875 ) ; glVertex3d ( 0.0, -1.61803398875, 0.61803398875 ) ; glVertex3d ( -1.0, -1.0, 1.0 ) ;
\r
869 void FGAPIENTRY glutWireOctahedron( void )
\r
871 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutWireOctahedron" );
\r
873 #define RADIUS 1.0f
\r
874 glBegin( GL_LINE_LOOP );
\r
875 glNormal3d( 0.577350269189, 0.577350269189, 0.577350269189); glVertex3d( RADIUS, 0.0, 0.0 ); glVertex3d( 0.0, RADIUS, 0.0 ); glVertex3d( 0.0, 0.0, RADIUS );
\r
876 glNormal3d( 0.577350269189, 0.577350269189,-0.577350269189); glVertex3d( RADIUS, 0.0, 0.0 ); glVertex3d( 0.0, 0.0,-RADIUS ); glVertex3d( 0.0, RADIUS, 0.0 );
\r
877 glNormal3d( 0.577350269189,-0.577350269189, 0.577350269189); glVertex3d( RADIUS, 0.0, 0.0 ); glVertex3d( 0.0, 0.0, RADIUS ); glVertex3d( 0.0,-RADIUS, 0.0 );
\r
878 glNormal3d( 0.577350269189,-0.577350269189,-0.577350269189); glVertex3d( RADIUS, 0.0, 0.0 ); glVertex3d( 0.0,-RADIUS, 0.0 ); glVertex3d( 0.0, 0.0,-RADIUS );
\r
879 glNormal3d(-0.577350269189, 0.577350269189, 0.577350269189); glVertex3d(-RADIUS, 0.0, 0.0 ); glVertex3d( 0.0, 0.0, RADIUS ); glVertex3d( 0.0, RADIUS, 0.0 );
\r
880 glNormal3d(-0.577350269189, 0.577350269189,-0.577350269189); glVertex3d(-RADIUS, 0.0, 0.0 ); glVertex3d( 0.0, RADIUS, 0.0 ); glVertex3d( 0.0, 0.0,-RADIUS );
\r
881 glNormal3d(-0.577350269189,-0.577350269189, 0.577350269189); glVertex3d(-RADIUS, 0.0, 0.0 ); glVertex3d( 0.0,-RADIUS, 0.0 ); glVertex3d( 0.0, 0.0, RADIUS );
\r
882 glNormal3d(-0.577350269189,-0.577350269189,-0.577350269189); glVertex3d(-RADIUS, 0.0, 0.0 ); glVertex3d( 0.0, 0.0,-RADIUS ); glVertex3d( 0.0,-RADIUS, 0.0 );
\r
890 void FGAPIENTRY glutSolidOctahedron( void )
\r
892 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutSolidOctahedron" );
\r
894 #define RADIUS 1.0f
\r
895 glBegin( GL_TRIANGLES );
\r
896 glNormal3d( 0.577350269189, 0.577350269189, 0.577350269189); glVertex3d( RADIUS, 0.0, 0.0 ); glVertex3d( 0.0, RADIUS, 0.0 ); glVertex3d( 0.0, 0.0, RADIUS );
\r
897 glNormal3d( 0.577350269189, 0.577350269189,-0.577350269189); glVertex3d( RADIUS, 0.0, 0.0 ); glVertex3d( 0.0, 0.0,-RADIUS ); glVertex3d( 0.0, RADIUS, 0.0 );
\r
898 glNormal3d( 0.577350269189,-0.577350269189, 0.577350269189); glVertex3d( RADIUS, 0.0, 0.0 ); glVertex3d( 0.0, 0.0, RADIUS ); glVertex3d( 0.0,-RADIUS, 0.0 );
\r
899 glNormal3d( 0.577350269189,-0.577350269189,-0.577350269189); glVertex3d( RADIUS, 0.0, 0.0 ); glVertex3d( 0.0,-RADIUS, 0.0 ); glVertex3d( 0.0, 0.0,-RADIUS );
\r
900 glNormal3d(-0.577350269189, 0.577350269189, 0.577350269189); glVertex3d(-RADIUS, 0.0, 0.0 ); glVertex3d( 0.0, 0.0, RADIUS ); glVertex3d( 0.0, RADIUS, 0.0 );
\r
901 glNormal3d(-0.577350269189, 0.577350269189,-0.577350269189); glVertex3d(-RADIUS, 0.0, 0.0 ); glVertex3d( 0.0, RADIUS, 0.0 ); glVertex3d( 0.0, 0.0,-RADIUS );
\r
902 glNormal3d(-0.577350269189,-0.577350269189, 0.577350269189); glVertex3d(-RADIUS, 0.0, 0.0 ); glVertex3d( 0.0,-RADIUS, 0.0 ); glVertex3d( 0.0, 0.0, RADIUS );
\r
903 glNormal3d(-0.577350269189,-0.577350269189,-0.577350269189); glVertex3d(-RADIUS, 0.0, 0.0 ); glVertex3d( 0.0, 0.0,-RADIUS ); glVertex3d( 0.0,-RADIUS, 0.0 );
\r
908 /* Magic Numbers: r0 = ( 1, 0, 0 )
\r
909 * r1 = ( -1/3, 2 sqrt(2) / 3, 0 )
\r
910 * r2 = ( -1/3, -sqrt(2) / 3, sqrt(6) / 3 )
\r
911 * r3 = ( -1/3, -sqrt(2) / 3, -sqrt(6) / 3 )
\r
912 * |r0| = |r1| = |r2| = |r3| = 1
\r
913 * Distance between any two points is 2 sqrt(6) / 3
\r
915 * Normals: The unit normals are simply the negative of the coordinates of the point not on the surface.
\r
918 #define NUM_TETR_FACES 4
\r
920 static GLdouble tet_r[4][3] = { { 1.0, 0.0, 0.0 },
\r
921 { -0.333333333333, 0.942809041582, 0.0 },
\r
922 { -0.333333333333, -0.471404520791, 0.816496580928 },
\r
923 { -0.333333333333, -0.471404520791, -0.816496580928 } } ;
\r
925 static GLint tet_i[4][3] = /* Vertex indices */
\r
927 { 1, 3, 2 }, { 0, 2, 3 }, { 0, 3, 1 }, { 0, 1, 2 }
\r
933 void FGAPIENTRY glutWireTetrahedron( void )
\r
935 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutWireTetrahedron" );
\r
937 glBegin( GL_LINE_LOOP ) ;
\r
938 glNormal3d ( -tet_r[0][0], -tet_r[0][1], -tet_r[0][2] ) ; glVertex3dv ( tet_r[1] ) ; glVertex3dv ( tet_r[3] ) ; glVertex3dv ( tet_r[2] ) ;
\r
939 glNormal3d ( -tet_r[1][0], -tet_r[1][1], -tet_r[1][2] ) ; glVertex3dv ( tet_r[0] ) ; glVertex3dv ( tet_r[2] ) ; glVertex3dv ( tet_r[3] ) ;
\r
940 glNormal3d ( -tet_r[2][0], -tet_r[2][1], -tet_r[2][2] ) ; glVertex3dv ( tet_r[0] ) ; glVertex3dv ( tet_r[3] ) ; glVertex3dv ( tet_r[1] ) ;
\r
941 glNormal3d ( -tet_r[3][0], -tet_r[3][1], -tet_r[3][2] ) ; glVertex3dv ( tet_r[0] ) ; glVertex3dv ( tet_r[1] ) ; glVertex3dv ( tet_r[2] ) ;
\r
948 void FGAPIENTRY glutSolidTetrahedron( void )
\r
950 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutSolidTetrahedron" );
\r
952 glBegin( GL_TRIANGLES ) ;
\r
953 glNormal3d ( -tet_r[0][0], -tet_r[0][1], -tet_r[0][2] ) ; glVertex3dv ( tet_r[1] ) ; glVertex3dv ( tet_r[3] ) ; glVertex3dv ( tet_r[2] ) ;
\r
954 glNormal3d ( -tet_r[1][0], -tet_r[1][1], -tet_r[1][2] ) ; glVertex3dv ( tet_r[0] ) ; glVertex3dv ( tet_r[2] ) ; glVertex3dv ( tet_r[3] ) ;
\r
955 glNormal3d ( -tet_r[2][0], -tet_r[2][1], -tet_r[2][2] ) ; glVertex3dv ( tet_r[0] ) ; glVertex3dv ( tet_r[3] ) ; glVertex3dv ( tet_r[1] ) ;
\r
956 glNormal3d ( -tet_r[3][0], -tet_r[3][1], -tet_r[3][2] ) ; glVertex3dv ( tet_r[0] ) ; glVertex3dv ( tet_r[1] ) ; glVertex3dv ( tet_r[2] ) ;
\r
963 static double icos_r[12][3] = {
\r
965 { 0.447213595500, 0.894427191000, 0.0 },
\r
966 { 0.447213595500, 0.276393202252, 0.850650808354 },
\r
967 { 0.447213595500, -0.723606797748, 0.525731112119 },
\r
968 { 0.447213595500, -0.723606797748, -0.525731112119 },
\r
969 { 0.447213595500, 0.276393202252, -0.850650808354 },
\r
970 { -0.447213595500, -0.894427191000, 0.0 },
\r
971 { -0.447213595500, -0.276393202252, 0.850650808354 },
\r
972 { -0.447213595500, 0.723606797748, 0.525731112119 },
\r
973 { -0.447213595500, 0.723606797748, -0.525731112119 },
\r
974 { -0.447213595500, -0.276393202252, -0.850650808354 },
\r
978 static int icos_v [20][3] = {
\r
1001 void FGAPIENTRY glutWireIcosahedron( void )
\r
1005 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutWireIcosahedron" );
\r
1007 for ( i = 0; i < 20; i++ )
\r
1009 double normal[3] ;
\r
1010 normal[0] = ( icos_r[icos_v[i][1]][1] - icos_r[icos_v[i][0]][1] ) * ( icos_r[icos_v[i][2]][2] - icos_r[icos_v[i][0]][2] ) - ( icos_r[icos_v[i][1]][2] - icos_r[icos_v[i][0]][2] ) * ( icos_r[icos_v[i][2]][1] - icos_r[icos_v[i][0]][1] ) ;
\r
1011 normal[1] = ( icos_r[icos_v[i][1]][2] - icos_r[icos_v[i][0]][2] ) * ( icos_r[icos_v[i][2]][0] - icos_r[icos_v[i][0]][0] ) - ( icos_r[icos_v[i][1]][0] - icos_r[icos_v[i][0]][0] ) * ( icos_r[icos_v[i][2]][2] - icos_r[icos_v[i][0]][2] ) ;
\r
1012 normal[2] = ( icos_r[icos_v[i][1]][0] - icos_r[icos_v[i][0]][0] ) * ( icos_r[icos_v[i][2]][1] - icos_r[icos_v[i][0]][1] ) - ( icos_r[icos_v[i][1]][1] - icos_r[icos_v[i][0]][1] ) * ( icos_r[icos_v[i][2]][0] - icos_r[icos_v[i][0]][0] ) ;
\r
1013 glBegin ( GL_LINE_LOOP ) ;
\r
1014 glNormal3dv ( normal ) ;
\r
1015 glVertex3dv ( icos_r[icos_v[i][0]] ) ;
\r
1016 glVertex3dv ( icos_r[icos_v[i][1]] ) ;
\r
1017 glVertex3dv ( icos_r[icos_v[i][2]] ) ;
\r
1025 void FGAPIENTRY glutSolidIcosahedron( void )
\r
1029 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutSolidIcosahedron" );
\r
1031 glBegin ( GL_TRIANGLES ) ;
\r
1032 for ( i = 0; i < 20; i++ )
\r
1034 double normal[3] ;
\r
1035 normal[0] = ( icos_r[icos_v[i][1]][1] - icos_r[icos_v[i][0]][1] ) * ( icos_r[icos_v[i][2]][2] - icos_r[icos_v[i][0]][2] ) - ( icos_r[icos_v[i][1]][2] - icos_r[icos_v[i][0]][2] ) * ( icos_r[icos_v[i][2]][1] - icos_r[icos_v[i][0]][1] ) ;
\r
1036 normal[1] = ( icos_r[icos_v[i][1]][2] - icos_r[icos_v[i][0]][2] ) * ( icos_r[icos_v[i][2]][0] - icos_r[icos_v[i][0]][0] ) - ( icos_r[icos_v[i][1]][0] - icos_r[icos_v[i][0]][0] ) * ( icos_r[icos_v[i][2]][2] - icos_r[icos_v[i][0]][2] ) ;
\r
1037 normal[2] = ( icos_r[icos_v[i][1]][0] - icos_r[icos_v[i][0]][0] ) * ( icos_r[icos_v[i][2]][1] - icos_r[icos_v[i][0]][1] ) - ( icos_r[icos_v[i][1]][1] - icos_r[icos_v[i][0]][1] ) * ( icos_r[icos_v[i][2]][0] - icos_r[icos_v[i][0]][0] ) ;
\r
1038 glNormal3dv ( normal ) ;
\r
1039 glVertex3dv ( icos_r[icos_v[i][0]] ) ;
\r
1040 glVertex3dv ( icos_r[icos_v[i][1]] ) ;
\r
1041 glVertex3dv ( icos_r[icos_v[i][2]] ) ;
\r
1050 static double rdod_r[14][3] = {
\r
1051 { 0.0, 0.0, 1.0 },
\r
1052 { 0.707106781187, 0.000000000000, 0.5 },
\r
1053 { 0.000000000000, 0.707106781187, 0.5 },
\r
1054 { -0.707106781187, 0.000000000000, 0.5 },
\r
1055 { 0.000000000000, -0.707106781187, 0.5 },
\r
1056 { 0.707106781187, 0.707106781187, 0.0 },
\r
1057 { -0.707106781187, 0.707106781187, 0.0 },
\r
1058 { -0.707106781187, -0.707106781187, 0.0 },
\r
1059 { 0.707106781187, -0.707106781187, 0.0 },
\r
1060 { 0.707106781187, 0.000000000000, -0.5 },
\r
1061 { 0.000000000000, 0.707106781187, -0.5 },
\r
1062 { -0.707106781187, 0.000000000000, -0.5 },
\r
1063 { 0.000000000000, -0.707106781187, -0.5 },
\r
1064 { 0.0, 0.0, -1.0 }
\r
1067 static int rdod_v [12][4] = {
\r
1077 { 6, 10, 13, 11 },
\r
1078 { 7, 11, 13, 12 },
\r
1082 static double rdod_n[12][3] = {
\r
1083 { 0.353553390594, 0.353553390594, 0.5 },
\r
1084 { -0.353553390594, 0.353553390594, 0.5 },
\r
1085 { -0.353553390594, -0.353553390594, 0.5 },
\r
1086 { 0.353553390594, -0.353553390594, 0.5 },
\r
1087 { 0.000000000000, 1.000000000000, 0.0 },
\r
1088 { -1.000000000000, 0.000000000000, 0.0 },
\r
1089 { 0.000000000000, -1.000000000000, 0.0 },
\r
1090 { 1.000000000000, 0.000000000000, 0.0 },
\r
1091 { 0.353553390594, 0.353553390594, -0.5 },
\r
1092 { -0.353553390594, 0.353553390594, -0.5 },
\r
1093 { -0.353553390594, -0.353553390594, -0.5 },
\r
1094 { 0.353553390594, -0.353553390594, -0.5 }
\r
1097 void FGAPIENTRY glutWireRhombicDodecahedron( void )
\r
1101 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutWireRhombicDodecahedron" );
\r
1103 for ( i = 0; i < 12; i++ )
\r
1105 glBegin ( GL_LINE_LOOP ) ;
\r
1106 glNormal3dv ( rdod_n[i] ) ;
\r
1107 glVertex3dv ( rdod_r[rdod_v[i][0]] ) ;
\r
1108 glVertex3dv ( rdod_r[rdod_v[i][1]] ) ;
\r
1109 glVertex3dv ( rdod_r[rdod_v[i][2]] ) ;
\r
1110 glVertex3dv ( rdod_r[rdod_v[i][3]] ) ;
\r
1118 void FGAPIENTRY glutSolidRhombicDodecahedron( void )
\r
1122 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutSolidRhombicDodecahedron" );
\r
1124 glBegin ( GL_QUADS ) ;
\r
1125 for ( i = 0; i < 12; i++ )
\r
1127 glNormal3dv ( rdod_n[i] ) ;
\r
1128 glVertex3dv ( rdod_r[rdod_v[i][0]] ) ;
\r
1129 glVertex3dv ( rdod_r[rdod_v[i][1]] ) ;
\r
1130 glVertex3dv ( rdod_r[rdod_v[i][2]] ) ;
\r
1131 glVertex3dv ( rdod_r[rdod_v[i][3]] ) ;
\r
1137 void FGAPIENTRY glutWireSierpinskiSponge ( int num_levels, GLdouble offset[3], GLdouble scale )
\r
1141 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutWireSierpinskiSponge" );
\r
1143 if ( num_levels == 0 )
\r
1146 for ( i = 0 ; i < NUM_TETR_FACES ; i++ )
\r
1148 glBegin ( GL_LINE_LOOP ) ;
\r
1149 glNormal3d ( -tet_r[i][0], -tet_r[i][1], -tet_r[i][2] ) ;
\r
1150 for ( j = 0; j < 3; j++ )
\r
1152 double x = offset[0] + scale * tet_r[tet_i[i][j]][0] ;
\r
1153 double y = offset[1] + scale * tet_r[tet_i[i][j]][1] ;
\r
1154 double z = offset[2] + scale * tet_r[tet_i[i][j]][2] ;
\r
1155 glVertex3d ( x, y, z ) ;
\r
1161 else if ( num_levels > 0 )
\r
1163 GLdouble local_offset[3] ; /* Use a local variable to avoid buildup of roundoff errors */
\r
1166 for ( i = 0 ; i < NUM_TETR_FACES ; i++ )
\r
1168 local_offset[0] = offset[0] + scale * tet_r[i][0] ;
\r
1169 local_offset[1] = offset[1] + scale * tet_r[i][1] ;
\r
1170 local_offset[2] = offset[2] + scale * tet_r[i][2] ;
\r
1171 glutWireSierpinskiSponge ( num_levels, local_offset, scale ) ;
\r
1176 void FGAPIENTRY glutSolidSierpinskiSponge ( int num_levels, GLdouble offset[3], GLdouble scale )
\r
1180 FREEGLUT_EXIT_IF_NOT_INITIALISED ( "glutSolidSierpinskiSponge" );
\r
1182 if ( num_levels == 0 )
\r
1184 glBegin ( GL_TRIANGLES ) ;
\r
1186 for ( i = 0 ; i < NUM_TETR_FACES ; i++ )
\r
1188 glNormal3d ( -tet_r[i][0], -tet_r[i][1], -tet_r[i][2] ) ;
\r
1189 for ( j = 0; j < 3; j++ )
\r
1191 double x = offset[0] + scale * tet_r[tet_i[i][j]][0] ;
\r
1192 double y = offset[1] + scale * tet_r[tet_i[i][j]][1] ;
\r
1193 double z = offset[2] + scale * tet_r[tet_i[i][j]][2] ;
\r
1194 glVertex3d ( x, y, z ) ;
\r
1200 else if ( num_levels > 0 )
\r
1202 GLdouble local_offset[3] ; /* Use a local variable to avoid buildup of roundoff errors */
\r
1205 for ( i = 0 ; i < NUM_TETR_FACES ; i++ )
\r
1207 local_offset[0] = offset[0] + scale * tet_r[i][0] ;
\r
1208 local_offset[1] = offset[1] + scale * tet_r[i][1] ;
\r
1209 local_offset[2] = offset[2] + scale * tet_r[i][2] ;
\r
1210 glutSolidSierpinskiSponge ( num_levels, local_offset, scale ) ;
\r
1215 /*** END OF FILE ***/
\r