foam-extend4.1-coherent-io/applications/utilities/surface/surfaceCoarsen/bunnylod/bunnygut.C

/*
 *   Polygon Reduction Demo by Stan Melax (c) 1998
 *  Permission to use any of this code wherever you want is granted..
 *  Although, please do acknowledge authorship if appropriate.
 *
 *  This module initializes the bunny model data and calls
 *  the polygon reduction routine.  At each frame the RenderModel()
 *  routine is called to draw the model.  This module also
 *  animates the parameters (such as number of vertices to
 *  use) to show the model at various levels of detail.
 */

#include <windows.h>
#include <stdio.h>  
#include <math.h>
#include <stdlib.h>
#include <assert.h>
#include <string.h>
#include <GL/gl.h>
#pragma warning(disable : 4244)

#include "vector.h"
#include "font.h"
#include "progmesh.h"
#include "rabdata.h"

extern float DeltaT;  // change in time since last frame
int   render_num;   // number of vertices to draw with
float lodbase=0.5f; // the fraction of vertices used to morph toward
float morph=1.0f;   // where to render between 2 levels of detail
List<Vector> vert;       // global list of vertices
List<tridata> tri;       // global list of triangles
List<int> collapse_map;  // to which neighbor each vertex collapses
int renderpolycount=0;   // polygons rendered in the current frame
Vector model_position;         // position of bunny
Quaternion model_orientation;  // orientation of bunny

// Note that the use of the Map() function and the collapse_map
// list isn't part of the polygon reduction algorithm.
// We just set up this system here in this module
// so that we could retrieve the model at any desired vertex count.
// Therefore if this part of the program confuses you, then
// dont worry about it.  It might help to look over the progmesh.cpp
// module first.

//       Map()
//
// When the model is rendered using a maximum of mx vertices
// then it is vertices 0 through mx-1 that are used.  
// We are able to do this because the vertex list 
// gets sorted according to the collapse order.
// The Map() routine takes a vertex number 'a' and the
// maximum number of vertices 'mx' and returns the 
// appropriate vertex in the range 0 to mx-1.
// When 'a' is greater than 'mx' the Map() routine
// follows the chain of edge collapses until a vertex
// within the limit is reached.
//   An example to make this clear: assume there is
//   a triangle with vertices 1, 3 and 12.  But when
//   rendering the model we limit ourselves to 10 vertices.
//   In that case we find out how vertex 12 was removed
//   by the polygon reduction algorithm.  i.e. which
//   edge was collapsed.  Lets say that vertex 12 was collapsed
//   to vertex number 7.  This number would have been stored
//   in the collapse_map array (i.e. collapse_map[12]==7).
//   Since vertex 7 is in range (less than max of 10) we
//   will want to render the triangle 1,3,7.  
//   Pretend now that we want to limit ourselves to 5 vertices.
//   and vertex 7 was collapsed to vertex 3 
//   (i.e. collapse_map[7]==3).  Then triangle 1,3,12 would now be
//   triangle 1,3,3.  i.e. this polygon was removed by the
//   progressive mesh polygon reduction algorithm by the time
//   it had gotten down to 5 vertices.
//   No need to draw a one dimensional polygon. :-)
int Map(int a,int mx) {
	if(mx<=0) return 0;
	while(a>=mx) {  
		a=collapse_map[a];
	}
	return a;
}

void DrawModelTriangles() {
	assert(collapse_map.num);
	renderpolycount=0;
	int i=0;
	for(i=0;i<tri.num;i++) {
		int p0= Map(tri[i].v[0],render_num);
		int p1= Map(tri[i].v[1],render_num);
		int p2= Map(tri[i].v[2],render_num);
		// note:  serious optimization opportunity here,
		//  by sorting the triangles the following "continue" 
		//  could have been made into a "break" statement.
		if(p0==p1 || p1==p2 || p2==p0) continue;
		renderpolycount++;
		// if we are not currenly morphing between 2 levels of detail
		// (i.e. if morph=1.0) then q0,q1, and q2 are not necessary.
		int q0= Map(p0,(int)(render_num*lodbase));
		int q1= Map(p1,(int)(render_num*lodbase));
		int q2= Map(p2,(int)(render_num*lodbase));
		Vector v0,v1,v2; 
		v0 = vert[p0]*morph + vert[q0]*(1-morph);
		v1 = vert[p1]*morph + vert[q1]*(1-morph);
		v2 = vert[p2]*morph + vert[q2]*(1-morph);
		glBegin(GL_POLYGON);
			// the purpose of the demo is to show polygons
			// therefore just use 1 face normal (flat shading)
			Vector nrml = (v1-v0) * (v2-v1);  // cross product
			if(0<magnitude(nrml)) {
				glNormal3fv(normalize(nrml));
			}
			glVertex3fv(v0);
			glVertex3fv(v1);
			glVertex3fv(v2);
		glEnd();
	}
}


void PermuteVertices(List<int> &permutation) {
	// rearrange the vertex list 
	List<Vector> temp_list;
	int i;
	assert(permutation.num==vert.num);
	for(i=0;i<vert.num;i++) {
		temp_list.Add(vert[i]);
	}
	for(i=0;i<vert.num;i++) {
		vert[permutation[i]]=temp_list[i];
	}
	// update the changes in the entries in the triangle list
	for(i=0;i<tri.num;i++) {
		for(int j=0;j<3;j++) {
			tri[i].v[j] = permutation[tri[i].v[j]];
		}
	}
}

void GetRabbitData(){
	// Copy the geometry from the arrays of data in rabdata.cpp into
	// the vert and tri lists which we send to the reduction routine
	int i;
	for(i=0;i<RABBIT_VERTEX_NUM;i++) {
		float *vp=rabbit_vertices[i];
		vert.Add(Vector(vp[0],vp[1],vp[2]));
	}
	for(i=0;i<RABBIT_TRIANGLE_NUM;i++) {
		tridata td;
		td.v[0]=rabbit_triangles[i][0];
		td.v[1]=rabbit_triangles[i][1];
		td.v[2]=rabbit_triangles[i][2];
		tri.Add(td);
	}
	render_num=vert.num;  // by default lets use all the model to render
}


void InitModel() {
	List<int> permutation;
	GetRabbitData();  
	ProgressiveMesh(vert,tri,collapse_map,permutation);
	PermuteVertices(permutation);
	model_position    = Vector(0,0,-3);
	Quaternion yaw(Vector(0,1,0),-3.14f/4);    // 45 degrees
	Quaternion pitch(Vector(1,0,0),3.14f/12);  // 15 degrees 
	model_orientation = pitch*yaw;
}

void StatusDraw() {
	//  Draw a slider type widget looking thing
	// to show portion of vertices being used
	float b = (float)render_num/(float)vert.num;
	float a = b*(lodbase );
	glDisable(GL_LIGHTING);
	glMatrixMode( GL_PROJECTION );
	glPushMatrix();
	glLoadIdentity();
	glOrtho(-0.15,15,-0.1,1.1,-0.1,100);
	glMatrixMode( GL_MODELVIEW );

	glPushMatrix();
	glLoadIdentity();
	glBegin(GL_POLYGON);
	glColor3f(1,0,0);
	glVertex2f(0,0);
	glVertex2f(1,0);
	glVertex2f(1,a);
	glVertex2f(0,a);
	glEnd();
	glBegin(GL_POLYGON);
	glColor3f(1,0,0);
	glVertex2f(0,a);
	glVertex2f(morph,a);
	glVertex2f(morph,b);
	glVertex2f(0,b);
	glEnd();
	glBegin(GL_POLYGON);
	glColor3f(0,0,1);
	glVertex2f(morph,a);
	glVertex2f(1,a);
	glVertex2f(1,b);
	glVertex2f(morph,b);
	glEnd();
	glBegin(GL_POLYGON);
	glColor3f(0,0,1);
	glVertex2f(0,b);
	glVertex2f(1,b);
	glVertex2f(1,1);
	glVertex2f(0,1);
	glEnd();
	glPopMatrix();
	glMatrixMode( GL_PROJECTION );
	glPopMatrix();
	glMatrixMode( GL_MODELVIEW );
}

/*
 *  The following is just a quick hack to animate
 *  the object through various polygon reduced versions.
 */
struct keyframethings {
	float t;   // timestamp
	float n;   // portion of vertices used to start
	float dn;  // rate of change in "n"
	float m;   // morph value
	float dm;  // rate of change in "m"
} keys[]={
	{0 ,1   ,0 ,1, 0},
	{2 ,1   ,-1,1, 0},
	{10,0   ,1 ,1, 0},
	{18,1   ,0 ,1, 0},
	{20,1   ,0 ,1,-1},
	{24,0.5 ,0 ,1, 0},
	{26,0.5 ,0 ,1,-1},
	{30,0.25,0 ,1, 0},
	{32,0.25,0 ,1,-1},
	{36,0.125,0,1, 0},
	{38,0.25,0 ,0, 1},
	{42,0.5 ,0 ,0, 1},
	{46,1   ,0 ,0, 1},
	{50,1   ,0 ,1, 0},
};
void AnimateParameters() {
	static float time=0;  // global time - used for animation
	time+=DeltaT;
	if(time>=50) time=0;  // repeat cycle every so many seconds
	int k=0;
	while(time>keys[k+1].t) {
		k++;
	}
	float interp = (time-keys[k].t)/(keys[k+1].t-keys[k].t);
	render_num = vert.num*(keys[k].n + interp*keys[k].dn);
	morph    = keys[k].m + interp*keys[k].dm;
	morph = (morph>1.0f) ? 1.0f : morph;  // clamp value
	if(render_num>vert.num) render_num=vert.num;
	if(render_num<0       ) render_num=0;
}

void RenderModel() {
	AnimateParameters();
	
	glEnable(GL_LIGHTING);
	glEnable(GL_LIGHT0);
	glColor3f(1,1,1);
	glPushMatrix();
	glTranslatef(model_position.x,model_position.y,model_position.z);
	// Rotate by quaternion: model_orientation
	Vector axis=model_orientation.axis();
	float angle=model_orientation.angle()*180.0f/3.14f;
	glRotatef(angle,axis.x,axis.y,axis.z);
	DrawModelTriangles();
	StatusDraw();
	glPopMatrix();

	char buf[256];
	sprintf(buf,"Polys: %d  Vertices: %d ",renderpolycount,render_num);
	if(morph<1.0) {
		sprintf(buf+strlen(buf),"<-> %d  morph: %4.2f ",
		        (int)(lodbase *render_num),morph); 
	}
	PostString(buf,0,-2,5);
}