Le source : TP-OndeCirculaire.cpp
/* Auteur: Nicolas JANEY */
/* nicolas.janey@univ-fcomte.fr */
/* Fevrier 2003 */
/* Animation d'un maillage de facettes */
#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include <GL/glut.h>
#include <GL/gl.h>
#include <GL/glu.h>
#include "ModuleMenus.h"
#include "ModuleReshape.h"
#include "ModuleManipulateur.h"
#include "ModuleCouleurs.h"
#ifndef M_PI
#define M_PI 3.14159
#endif
static int mode = 2;
static int disci = 150;
static int discj = 36;
static float d = 0.0F;
static float larg = 20.0F;
static float factz = 5.0F;
static float fact = 5.0F;
struct coordonnee {
float x;
float y;
float z; };
struct coordonnee **tb;
struct coordonnee **no;
struct coordonnee **allocationTableau(int n,int m) {
struct coordonnee **t =(struct coordonnee **) calloc(n,sizeof(struct coordonnee *));
for ( int i = 0 ; i < n ; i++ )
t[i] =(struct coordonnee *) calloc(m,sizeof(struct coordonnee));
return(t);
}
void liberationTableau(struct coordonnee **t,int n,int m) {
for ( int i = 0 ; i < n ; i++ )
free(t[i]);
free(t);
}
void calculTableau(struct coordonnee **t,int n,int m) {
for ( int j = 0 ; j < m ; j++ ) {
float a =(float) j*2.0F*M_PI/m;
float sn = sin(a);
float cs = cos(a);
for ( int i = 0 ; i < n ; i++ ) {
float r =(float) larg/2.0F*i/(n-1);
t[i][j].x = r*cs;
t[i][j].y = r*sn;
double dist = fact*r;
t[i][j].z =(float) factz*sin(dist+d)/(1+dist); } }
}
void calculNormales(struct coordonnee **no,int n,int m) {
for ( int j = 0 ; j < m ; j++ ) {
float a =(float) j*2.0F*M_PI/m;
float csa = cos(a);
float sna = sin(a);
for ( int i = 0 ; i < n ; i++ ) {
double dist =(float) larg/2.0F*i/(n-1);
double cs = cos(d+fact*dist);
double sn = sin(d+fact*dist);
float dx = factz*fact*(sn/(1+fact*dist)-cs)/(1+fact*dist);
float dz = 1.0F;
no[i][j].x = dx*csa;
no[i][j].y = dx*sna;
no[i][j].z = dz; } }
}
void afficheTableauPoints(struct coordonnee **t,struct coordonnee **no,int n,int m) {
glBegin(GL_POINTS);
for ( int i = 0 ; i < n ; i++ )
for ( int j = 0 ; j < m ; j++ ) {
glNormal3f(no[i][j].x,no[i][j].y,no[i][j].z);
glVertex3f(t[i][j].x,t[i][j].y,t[i][j].z); }
glEnd();
}
void afficheTableauLignes(struct coordonnee **t,struct coordonnee **no,int n,int m) {
for ( int i = 0 ; i < n ; i++ ) {
glBegin(GL_LINE_STRIP);
for ( int j = 0 ; j <= m ; j++ ) {
int jj = j%m;
glNormal3f(no[i][jj].x,no[i][jj].y,no[i][jj].z);
glVertex3f(t[i][jj].x,t[i][jj].y,t[i][jj].z); }
glEnd(); }
for ( int j = 0 ; j <= m ; j++ ) {
int jj = j%m;
glBegin(GL_LINE_STRIP);
for ( int i = 0 ; i < n ; i++ ) {
glNormal3f(no[i][jj].x,no[i][jj].y,no[i][jj].z);
glVertex3f(t[i][jj].x,t[i][jj].y,t[i][jj].z); }
glEnd(); }
/* glBegin(GL_LINES);
for ( i = 0 ; i < n ; i++ ) {
for ( j = 0 ; j < m ; j++ ) {
glVertex3f(t[i][j].x,t[i][j].y,t[i][j].z);
glVertex3f(t[i][j].x+no[i][j].x/2,t[i][j].y+no[i][j].y/2,t[i][j].z+no[i][j].z/2);} }
glEnd();*/
}
void afficheTableauFaces(struct coordonnee **t,struct coordonnee **no,int n,int m) {
for ( int i = 0 ; i < n-1 ; i++ ) {
glBegin(GL_QUAD_STRIP);
for ( int j = 0 ; j <= m ; j++ ) {
int jj = j%m;
glNormal3f(no[i][jj].x,no[i][jj].y,no[i][jj].z);
glVertex3f(t[i][jj].x,t[i][jj].y,t[i][jj].z);
glNormal3f(no[i+1][jj].x,no[i+1][jj].y,no[i+1][jj].z);
glVertex3f(t[i+1][jj].x,t[i+1][jj].y,t[i+1][jj].z); }
glEnd(); }
}
void myinit(void) {
GLfloat shinines[] = { 50.0F };
glClearColor(0.5F,0.5F,0.8F,1.0F) ;
glMaterialfv(GL_FRONT,GL_DIFFUSE,couleurBlanc());
glMaterialfv(GL_FRONT,GL_SPECULAR,couleurGrisFonce());
glMaterialfv(GL_FRONT,GL_SHININESS,shinines);
glEnable(GL_LIGHTING);
glEnable(GL_LIGHT0);
glEnable(GL_LIGHT1);
glEnable(GL_LIGHT2);
float dir0[4] = { 1.0F,1.0F,1.0F,0.0F };
glLightfv(GL_LIGHT0,GL_POSITION,dir0) ;
glLightfv(GL_LIGHT0,GL_DIFFUSE,couleurVert()) ;
float dir1[4] = { -1.0F,1.0F,1.0F,0.0F };
glLightfv(GL_LIGHT1,GL_POSITION,dir1) ;
glLightfv(GL_LIGHT1,GL_DIFFUSE,couleurBleu()) ;
float dir2[4] = { 0.0F,-1.0F,1.0F,0.0F };
glLightfv(GL_LIGHT2,GL_POSITION,dir2) ;
glLightfv(GL_LIGHT2,GL_DIFFUSE,couleurRouge()) ;
glEnable(GL_AUTO_NORMAL);
glEnable(GL_NORMALIZE);
glDepthFunc(GL_LESS);
glEnable(GL_DEPTH_TEST);
}
void display(void) {
glClear(GL_COLOR_BUFFER_BIT|GL_DEPTH_BUFFER_BIT);
glPushMatrix();
manipulateurSouris();
manipulateurClavier();
switch ( mode ) {
case 0 : afficheTableauPoints(tb,no,disci,discj);
break;
case 1 : afficheTableauLignes(tb,no,disci,discj);
break;
case 2 : afficheTableauFaces(tb,no,disci,discj);
break; }
glPopMatrix();
glFlush();
glutSwapBuffers();
}
void key(unsigned char key,int x,int y) {
if ( keyManipulateur(key,x,y) )
glutPostRedisplay();
else
switch ( key ) {
case 'a' : liberationTableau(tb,disci,discj);
liberationTableau(no,disci,discj);
disci++ ;
tb = allocationTableau(disci,discj);
no = allocationTableau(disci,discj);
calculTableau(tb,disci,discj);
calculNormales(no,disci,discj);
glutPostRedisplay() ;
break ;
case 'A' : liberationTableau(tb,disci,discj);
liberationTableau(no,disci,discj);
disci-- ;
if ( disci < 3 )
disci = 3;
tb = allocationTableau(disci,discj);
no = allocationTableau(disci,discj);
calculTableau(tb,disci,discj);
calculNormales(no,disci,discj);
glutPostRedisplay() ;
break ;
case 'b' : liberationTableau(tb,disci,discj);
liberationTableau(no,disci,discj);
discj++ ;
tb = allocationTableau(disci,discj);
no = allocationTableau(disci,discj);
calculTableau(tb,disci,discj);
calculNormales(no,disci,discj);
glutPostRedisplay() ;
break ;
case 'B' : liberationTableau(tb,disci,discj);
liberationTableau(no,disci,discj);
discj-- ;
if ( discj < 3 )
discj = 3;
tb = allocationTableau(disci,discj);
no = allocationTableau(disci,discj);
calculTableau(tb,disci,discj);
calculNormales(no,disci,discj);
glutPostRedisplay() ;
break ;
case 43 : fact *= 1.1F;
glutPostRedisplay() ;
break ;
case 45 : fact /= 1.1F;
glutPostRedisplay() ;
break ;
case 'c' : factz *= 1.1F;
glutPostRedisplay() ;
break ;
case 'C' : factz /= 1.1F;
glutPostRedisplay() ;
break ;
case ' ' : mode = (mode+1)%3 ;
glutPostRedisplay() ;
break ; }
}
void idle(void) {
d -=0.05;
calculTableau(tb,disci,discj);
calculNormales(no,disci,discj);
glutPostRedisplay();
}
int main(int argc,char **argv) {
tb = allocationTableau(disci,discj);
no = allocationTableau(disci,discj);
calculTableau(tb,disci,discj);
calculTableau(no,disci,discj);
glutInit(&argc,argv);
glutInitDisplayMode(GLUT_RGBA|GLUT_DEPTH|GLUT_DOUBLE);
glutInitWindowSize(450,300);
glutInitWindowPosition(50,50);
glutCreateWindow("Animation d'un maillage de facettes");
myinit();
creationMenuBasique();
setParametresOrthoBasique(-10.0,10.0,-10.0,10.0,-500.0,500.0);
setManipulateurDistance(1.0F);
glutReshapeFunc(reshapeOrthoBasique);
glutKeyboardFunc(key);
glutIdleFunc(idle);
glutSpecialFunc(specialBasique);
glutMotionFunc(motionBasique);
glutMouseFunc(sourisBasique);
glutDisplayFunc(display);
glutMainLoop();
return(0);
}