WebSVN – gelsvn – Diff – /branches/cpp11-devel/apps/Raytracer/raytrace.cpp


#include <iostream>
#include <vector>
#include <algorithm>
#include <cmath>
 
#include <GLGraphics/gel_glut.h>
#include <GLGraphics/GLViewController.h>
 
#include <CGLA/Vec2f.h>
#include <CGLA/Vec3i.h>
#include <CGLA/Vec3f.h>
#include <CGLA/Vec3d.h>
#include <CGLA/Mat4x4f.h>
 
#include <HMesh/Manifold.h>
 
#include <Geometry/TriMesh.h>
#include <Geometry/obj_load.h>
#include <Geometry/Ray.h>
#include <Geometry/BSPTree.h>
#include <Geometry/build_bbtree.h>
#include <Geometry/AABox.h>
#include <Geometry/OBox.h>
 
#include <Util/Timer.h>
 
#include "Camera.h"
 
#define USE_BSP
 
using namespace std;
using namespace CGLA;
using namespace Geometry;
using namespace HMesh;
using namespace GLGraphics;
 
/*
 * TODO: 
 * - BBOX remove HMesh dependency - that is crazy.
 * - BBox visit child nodes in order of how far away the intersection point on
 *        the bbox is. Closest nodes visited first. Cull nodes farther than
 *        an actual intersection point.
 * - BBox Smooth interpolation of triangle normals. Straightforward.
 */
					
 
namespace
{
  const int MAX_OBJECTS = 4;   // Maximum number of triangles in a BSP tree node
  const int MAX_LEVEL = 20;    // Maximum number of BSP tree subdivisions
 
  const unsigned int TEX_SIZE = 512;
 
  bool raytrace = false;
  bool done = false;
  bool shadow = false;
 
  unsigned int winx = TEX_SIZE;     // Screen width
  unsigned int winy = TEX_SIZE;     // Screen height
 
  unsigned int PIXEL_SUBDIVS = 1;
 
  int mouse_state = GLUT_UP;
  int mouse_button = 0;
  int spin_timer = 20;
 
  GLViewController *vctrl;
 
  TriMesh mesh;
  vector<const TriMesh*> mesh_vector(1, &mesh);
  vector<Mat4x4f> transforms(1, identity_Mat4x4f());
 
  double light_pow = 1.0;
  Vec3f light_dir = normalize(Vec3f(1.0, 1.0, 1.0));
  Vec3d background(0.8, 0.9, 1.0);
 
  BSPTree tree;
  Camera* cam;
 
  Vec3f image[TEX_SIZE][TEX_SIZE];
  unsigned int image_tex;
 
  // Function to generate a random number between 0 and 1.
  // gel_rand() returns a random number ranging from 0 to GEL_RAND_MAX.
  inline double my_random()
  {
    return gel_rand()/static_cast<double>(GEL_RAND_MAX);
  }
		
  AABBTree bb_tree;		
}
 
void spin(int x);
 
//////////////////////////////////////////////////////////////
//      I N I T I A L I Z A T I O N               
//////////////////////////////////////////////////////////////
 
void init_texture(unsigned int& tex)
{
  glGenTextures(1, reinterpret_cast<GLuint*>(&tex));
 
  glBindTexture(GL_TEXTURE_2D, tex);
 
  glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);
  glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);
 
  // load the texture image
  glTexImage2D(GL_TEXTURE_2D, 0, GL_RGB, 
	             TEX_SIZE, TEX_SIZE,
	             0, GL_RGB, GL_FLOAT, image[0][0].get());
}
 
void initRaytracer()
{
  Vec3f c;
  float r;
  mesh.get_bsphere(c, r);
  r *= 1.5;
 
  // Initialize track ball
  vctrl = new GLViewController(winx, winy, c, r);
 
  // Initialize corresponding camera
  cam = new Camera(c - Vec3f(r), c, Vec3f(0.0f, 1.0f, 0.0f), 1.0f);
 
#ifdef USE_BSP
   cout << "Constructing BSP tree..." << endl;
   tree.init(mesh_vector, transforms, MAX_OBJECTS, MAX_LEVEL);
   tree.build();
#else
	// AABB TREE
	Manifold m;
	vector<int> faces(mesh.geometry.no_faces(), 3);
	cout << "Creating manifold" << endl;
	build_manifold(m, 
				         mesh.geometry.no_vertices(), 
				         &mesh.geometry.vertex(0), 
				         faces.size(), &faces[0], 
				         reinterpret_cast<const int*>(&mesh.geometry.face(0)));
	cout << "Building tree" << endl;
	build_AABBTree(m, bb_tree);
#endif
}
 
void initGL()
{
  glShadeModel(GL_SMOOTH); 
  glDisable(GL_CULL_FACE);
  glFrontFace(GL_CCW);
 
  glClearColor(1.0, 1.0, 1.0, 1.0);
  glColor3f(0.0, 0.0, 0.0);
}
 
 
//////////////////////////////////////////////////////////////
//      S H A D E   F U N C T I O N S
//////////////////////////////////////////////////////////////
 
double shadow_shade(Ray& r)
{
  r.compute_position();
  Ray shadow(r.hit_pos, light_dir);
  double s = tree.intersect(shadow) ? 0.0 : 1.0;
 
  r.compute_normal();
  return s*light_pow*dot(r.hit_normal, light_dir);
}
 
double lambertian_shade(Ray& r)
{
#ifdef USE_BSP
	r.compute_normal();
#endif
  return light_pow*dot(r.hit_normal, light_dir);
}
 
double (*shade_ray[2])(Ray&) = { lambertian_shade,
				                         shadow_shade      };
 
//////////////////////////////////////////////////////////////
//      D R A W   F U N C T I O N S
//////////////////////////////////////////////////////////////
 
/*
void enable_textures(TriMesh& tm)
{
  for(unsigned int i=0;i<tm.materials.size(); ++i)
  {
    Material& mat = tm.materials[i];
    if(mat.tex_name != "")
    {
      string name = mat.tex_path + mat.tex_name;
      
      GLuint tex_id;
      if(load_image_into_texture(name, tex_id))
	mat.tex_id = tex_id;
    }
  }
}
*/
 
void set_perspective_proj()
{
  glMatrixMode(GL_PROJECTION);	 
  glLoadIdentity();            
 
  gluPerspective(64.0, 1.0, 0.1, 1000.0);
 
  glMatrixMode(GL_MODELVIEW);
}
 
void set_ortho_proj()
{
  glMatrixMode(GL_PROJECTION);	 
  glLoadIdentity();             
    
  glOrtho(0.0, 1.0, 0.0, 1.0, -1.0, 1.0);
 
  glMatrixMode(GL_MODELVIEW);  
}
 
void draw_texture(unsigned int tex)
{
  static GLfloat verts[] = { 0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 1.0, 1.0 };
 
  glColor4f(1.0, 1.0, 1.0, 1.0);
 
  glBindTexture(GL_TEXTURE_2D, tex);
  glEnable(GL_TEXTURE_2D);
 
  glEnableClientState(GL_VERTEX_ARRAY);
  glEnableClientState(GL_TEXTURE_COORD_ARRAY);
  
  glVertexPointer(2, GL_FLOAT, 0, verts);
  glTexCoordPointer(2, GL_FLOAT, 0, verts);
 
  glDrawArrays(GL_TRIANGLE_STRIP, 0, 4);
 
  glDisableClientState(GL_VERTEX_ARRAY);
  glDisableClientState(GL_TEXTURE_COORD_ARRAY);
 
  glDisable(GL_TEXTURE_2D);
}
 
void drawOBJ()
{
  static bool washere = false;
  static unsigned int disp_list;
 
  if(!washere)
  {
    disp_list = glGenLists(1);
    glNewList(disp_list, GL_COMPILE);
 
    glBegin(GL_TRIANGLES);
    for(int i = 0; i < mesh.geometry.no_faces(); ++i)
    {
      Vec3i n_face = mesh.normals.face(i);
      Vec3i g_face = mesh.geometry.face(i);
      for(int j=0;j<3;j++)
      {
        double shade = 0.5;
 
        if(n_face != Geometry::NULL_FACE)
        {
          Vec3f norm = normalize(mesh.normals.vertex(n_face[j]));
          glNormal3fv(norm.get());
          shade = light_pow*dot(norm, light_dir);
        }
 
        glColor3d(shade, shade, shade);
        Vec3f vert = mesh.geometry.vertex(g_face[j]);	  
        glVertex3fv(vert.get());
      }
    }
    glEnd();
 
    glEndList();
    washere = true;
  }
  glCallList(disp_list);
}
 
 
//////////////////////////////////////////////////////////////
//      G L U T   C A L L B A C K   F U N C T I O N S                 
//////////////////////////////////////////////////////////////
 
void display()
{
  static bool first = true;
 
  if(first)
  {
    first = false;
    glutTimerFunc(spin_timer, spin, 0);
  }
 
  Vec3f eye, focus, up;
	vctrl->get_view_param(eye, focus, up);
	cam->set(eye, focus, up, cam->get_focal_dist());
	
  if(raytrace)
  {
    raytrace = false;
    
    cout << "Raytracing";
    float win_to_vp = 1.0f/static_cast<float>(TEX_SIZE);
    float lowerleft = 0.5 + win_to_vp*0.5;
    float step = win_to_vp/static_cast<float>(PIXEL_SUBDIVS);
    
    vector<Vec2f> jitter(PIXEL_SUBDIVS*PIXEL_SUBDIVS); 
    for(unsigned int i = 0; i < PIXEL_SUBDIVS; ++i)
      for(unsigned int j = 0; j < PIXEL_SUBDIVS; ++j)
      {
	      jitter[i*PIXEL_SUBDIVS + j][0] = (my_random() + (j%PIXEL_SUBDIVS))*step; 
	      jitter[i*PIXEL_SUBDIVS + j][1] = (my_random() + (i%PIXEL_SUBDIVS))*step; 
      }
 
    Util::Timer tim;
    tim.start();
    for(unsigned int i = 0; i < TEX_SIZE; ++i)
    {
	    for(unsigned int j = 0; j < TEX_SIZE; ++j)
	    {
		    Vec3d sum(0.0f);
		    Vec2f vp_pos(j*win_to_vp - lowerleft, i*win_to_vp - lowerleft);
    	
		    for(unsigned int ky = 0; ky < PIXEL_SUBDIVS; ++ky)
		      for(unsigned int kx = 0; kx < PIXEL_SUBDIVS; ++kx)
		      {
			      Ray r = cam->get_ray(vp_pos + jitter[ky*PIXEL_SUBDIVS + kx]);
 
#ifdef USE_BSP
			      if(tree.intersect(r))
			        sum += Vec3d(shade_ray[shadow](r));
			      else
			        sum += background;
#else
			      float t = FLT_MAX;
			      bb_tree.intersect(r);
			      if(r.has_hit)
			        sum += Vec3d(shade_ray[0](r));
			      else 
			        sum += background;
#endif		
		      }
	      image[i][j] = Vec3f(sum/static_cast<double>(PIXEL_SUBDIVS*PIXEL_SUBDIVS));
      }
      if(((i + 1) % 50) == 0) cerr << ".";
    }
	  cout << " - " << tim.get_secs() << " secs " << endl;
    cout << endl;
 
    init_texture(image_tex);
 
    done = true;
  }
 
  if(done)
  {
    set_ortho_proj();
    glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
    glLoadIdentity();
 
    draw_texture(image_tex);
  }
  else
  {
    glEnable(GL_DEPTH_TEST);
 
    cam->glSetPerspective(winx, winy);
 
    glClearColor(background[0], background[1], background[2], 1.0);
    glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
    glLoadIdentity();
 
	  cam->glSetCamera();
    
    glColor3f(0.5, 0.5, 0.5);
    drawOBJ();
 
    glDisable(GL_DEPTH_TEST);
  }
 
  glutSwapBuffers();  
}
 
void reshape(int w, int h)
{
  winx = w; winy = h;
 
  vctrl->reshape(winx, winy);
 
  glViewport(0, 0, winx, winy);
}
 
void keyboard(unsigned char key, int x, int y)
{
  switch(key)
  {
  case '+':
    ++PIXEL_SUBDIVS;
    cout << "Rays per pixel: " << PIXEL_SUBDIVS*PIXEL_SUBDIVS << endl;
    break;
  case '-':
    if(PIXEL_SUBDIVS > 1)
      --PIXEL_SUBDIVS;
    cout << "Rays per pixel: " << PIXEL_SUBDIVS*PIXEL_SUBDIVS << endl;
    break;
  case 'r':
    if(done) done = false;
    else raytrace = true;
    break;
  case 's':
    shadow = !shadow;
    cout << "Shadow " << (shadow ? "on." : "off.") << endl;
    break;
  case 27:
    delete vctrl;
    delete cam;
    exit(0);
  }
}
 
void mouse(int btn, int state, int x, int y)
{
  if(state == GLUT_DOWN) 
  {
    if(btn == GLUT_LEFT_BUTTON) 
      vctrl->grab_ball(ROTATE_ACTION, Vec2i(x, y));
    else if(btn == GLUT_MIDDLE_BUTTON) 
      vctrl->grab_ball(ZOOM_ACTION, Vec2i(x, y));
    else if(btn == GLUT_RIGHT_BUTTON) 
      vctrl->grab_ball(PAN_ACTION, Vec2i(x, y));
  }
  else if(state == GLUT_UP)
			vctrl->release_ball();
 
  mouse_state = state;
  mouse_button = btn;
 
  glutPostRedisplay();
}
 
void move(int x, int y)
{
    vctrl->roll_ball(Vec2i(x, y));
 
  glutPostRedisplay();
}
 
void spin(int x)
{
  vctrl->try_spin();
  glutTimerFunc(spin_timer, spin, 0);  
  glutPostRedisplay();
}
 
 
//////////////////////////////////////////////////////////////
//                        M A I N
//////////////////////////////////////////////////////////////
 
int main(int argc, char** argv)
{
#ifdef __BORLANDC__
  _control87(MCW_EM, MCW_EM);  // Borland C++ will crash OpenGL if this
                               // magic line is not inserted
#endif
 
  glutInit(&argc, argv);
  glutInitDisplayMode(GLUT_RGBA | GLUT_DOUBLE | GLUT_DEPTH);
  glutInitWindowSize(winx, winy);
  glutCreateWindow("Press 'r' to raytrace");
  glutDisplayFunc(display);
  glutReshapeFunc(reshape);
  glutKeyboardFunc(keyboard);
  glutMouseFunc(mouse);
  glutMotionFunc(move);
 
  // LOAD OBJ
  string filename;
  if(argc > 1)
  {
    filename = argv[1];
    cout << "Loading " << filename << endl;
 
    obj_load(filename, mesh);
    //if(!mesh.has_normals())
    //{
      cout << "Computing normals" << endl;
      mesh.compute_normals();
    //}
 
    cout << "No. of triangles: " << mesh.geometry.no_faces() << endl;
  }
  else
  {
    obj_load("../../../data/dolphins.obj", mesh);
 
  	cout << "Computing normals" << endl;
    mesh.compute_normals();
 
    //cout << "Usage: raytrace any_object.obj";
    //exit(0);
  }
 
  initRaytracer();
  initGL();    
 
  glutMainLoop();
 
  return 0;
}
 
 

Rev 630	Rev 638
1	`#include <iostream>`	1	`#include <iostream>`
2	`#include <vector>`	2	`#include <vector>`
3	`#include <algorithm>`	3	`#include <algorithm>`
4	`#include <cmath>`	4	`#include <cmath>`
5		5
6	`#include <GLGraphics/gel_glut.h>`	6	`#include <GLGraphics/gel_glut.h>`
7	`#include <GLGraphics/GLViewController.h>`	7	`#include <GLGraphics/GLViewController.h>`
8		8
9	`#include <CGLA/Vec2f.h>`	9	`#include <CGLA/Vec2f.h>`
10	`#include <CGLA/Vec3i.h>`	10	`#include <CGLA/Vec3i.h>`
11	`#include <CGLA/Vec3f.h>`	11	`#include <CGLA/Vec3f.h>`
12	`#include <CGLA/Vec3d.h>`	12	`#include <CGLA/Vec3d.h>`
13	`#include <CGLA/Mat4x4f.h>`	13	`#include <CGLA/Mat4x4f.h>`
14		14
15	`#include <HMesh/Manifold.h>`	15	`#include <HMesh/Manifold.h>`
16		16
17	`#include <Geometry/TriMesh.h>`	17	`#include <Geometry/TriMesh.h>`
18	`#include <Geometry/obj_load.h>`	18	`#include <Geometry/obj_load.h>`
19	`#include <Geometry/Ray.h>`	19	`#include <Geometry/Ray.h>`
20	`#include <Geometry/BSPTree.h>`	20	`#include <Geometry/BSPTree.h>`
21	`#include <Geometry/build_bbtree.h>`	21	`#include <Geometry/build_bbtree.h>`
22	`#include <Geometry/AABox.h>`	22	`#include <Geometry/AABox.h>`
23	`#include <Geometry/OBox.h>`	23	`#include <Geometry/OBox.h>`
24		24
25	`#include <Util/Timer.h>`	25	`#include <Util/Timer.h>`
26		26
27	`#include "Camera.h"`	27	`#include "Camera.h"`
28		28
29	`#define USE_BSP`	29	`#define USE_BSP`
30		30
31	`using namespace std;`	31	`using namespace std;`
32	`using namespace CGLA;`	32	`using namespace CGLA;`
33	`using namespace Geometry;`	33	`using namespace Geometry;`
34	`using namespace HMesh;`	34	`using namespace HMesh;`
35	`using namespace GLGraphics;`	35	`using namespace GLGraphics;`
36		36
37	`/*`	37	`/*`
38	`* TODO:`	38	`* TODO:`
39	`* - BBOX remove HMesh dependency - that is crazy.`	39	`* - BBOX remove HMesh dependency - that is crazy.`
40	`* - BBox visit child nodes in order of how far away the intersection point on`	40	`* - BBox visit child nodes in order of how far away the intersection point on`
41	`* the bbox is. Closest nodes visited first. Cull nodes farther than`	41	`* the bbox is. Closest nodes visited first. Cull nodes farther than`
42	`* an actual intersection point.`	42	`* an actual intersection point.`
43	`* - BBox Smooth interpolation of triangle normals. Straightforward.`	43	`* - BBox Smooth interpolation of triangle normals. Straightforward.`
44	`*/`	44	`*/`
45		45
46		46
47	`namespace`	47	`namespace`
48	`{`	48	`{`
49	`const int MAX_OBJECTS = 4; // Maximum number of triangles in a BSP tree node`	49	`const int MAX_OBJECTS = 4; // Maximum number of triangles in a BSP tree node`
50	`const int MAX_LEVEL = 20; // Maximum number of BSP tree subdivisions`	50	`const int MAX_LEVEL = 20; // Maximum number of BSP tree subdivisions`
51		51
52	`const unsigned int TEX_SIZE = 512;`	52	`const unsigned int TEX_SIZE = 512;`
53		53
54	`bool raytrace = false;`	54	`bool raytrace = false;`
55	`bool done = false;`	55	`bool done = false;`
56	`bool shadow = false;`	56	`bool shadow = false;`
57		57
58	`unsigned int winx = TEX_SIZE; // Screen width`	58	`unsigned int winx = TEX_SIZE; // Screen width`
59	`unsigned int winy = TEX_SIZE; // Screen height`	59	`unsigned int winy = TEX_SIZE; // Screen height`
60		60
61	`unsigned int PIXEL_SUBDIVS = 1;`	61	`unsigned int PIXEL_SUBDIVS = 1;`
62		62
63	`int mouse_state = GLUT_UP;`	63	`int mouse_state = GLUT_UP;`
64	`int mouse_button = 0;`	64	`int mouse_button = 0;`
65	`int spin_timer = 20;`	65	`int spin_timer = 20;`
66		66
67	`GLViewController *vctrl;`	67	`GLViewController *vctrl;`
68		68
69	`TriMesh mesh;`	69	`TriMesh mesh;`
70	`vector<const TriMesh*> mesh_vector(1, &mesh);`	70	`vector<const TriMesh*> mesh_vector(1, &mesh);`
71	`vector<Mat4x4f> transforms(1, identity_Mat4x4f());`	71	`vector<Mat4x4f> transforms(1, identity_Mat4x4f());`
72		72
73	`double light_pow = 1.0;`	73	`double light_pow = 1.0;`
74	`Vec3f light_dir = normalize(Vec3f(1.0, 1.0, 1.0));`	74	`Vec3f light_dir = normalize(Vec3f(1.0, 1.0, 1.0));`
75	`Vec3d background(0.8, 0.9, 1.0);`	75	`Vec3d background(0.8, 0.9, 1.0);`
76		76
77	`BSPTree tree;`	77	`BSPTree tree;`
78	`Camera* cam;`	78	`Camera* cam;`
79		79
80	`Vec3f image[TEX_SIZE][TEX_SIZE];`	80	`Vec3f image[TEX_SIZE][TEX_SIZE];`
81	`unsigned int image_tex;`	81	`unsigned int image_tex;`
82		82
83	`// Function to generate a random number between 0 and 1.`	83	`// Function to generate a random number between 0 and 1.`
84	`// gel_rand() returns a random number ranging from 0 to GEL_RAND_MAX.`	84	`// gel_rand() returns a random number ranging from 0 to GEL_RAND_MAX.`
85	`inline double my_random()`	85	`inline double my_random()`
86	`{`	86	`{`
87	`return gel_rand()/static_cast<double>(GEL_RAND_MAX);`	87	`return gel_rand()/static_cast<double>(GEL_RAND_MAX);`
88	`}`	88	`}`
89		89
90	`AABBTree bb_tree;`	90	`AABBTree bb_tree;`
91	`}`	91	`}`
92		92
93	`void spin(int x);`	93	`void spin(int x);`
94		94
95	`//////////////////////////////////////////////////////////////`	95	`//////////////////////////////////////////////////////////////`
96	`// I N I T I A L I Z A T I O N`	96	`// I N I T I A L I Z A T I O N`
97	`//////////////////////////////////////////////////////////////`	97	`//////////////////////////////////////////////////////////////`
98		98
99	`void init_texture(unsigned int& tex)`	99	`void init_texture(unsigned int& tex)`
100	`{`	100	`{`
101	`glGenTextures(1, reinterpret_cast<GLuint*>(&tex));`	101	`glGenTextures(1, reinterpret_cast<GLuint*>(&tex));`
102		102
103	`glBindTexture(GL_TEXTURE_2D, tex);`	103	`glBindTexture(GL_TEXTURE_2D, tex);`
104		104
105	`glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);`	105	`glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);`
106	`glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);`	106	`glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);`
107		107
108	`// load the texture image`	108	`// load the texture image`
109	`glTexImage2D(GL_TEXTURE_2D, 0, GL_RGB,`	109	`glTexImage2D(GL_TEXTURE_2D, 0, GL_RGB,`
110	`TEX_SIZE, TEX_SIZE,`	110	`TEX_SIZE, TEX_SIZE,`
111	`0, GL_RGB, GL_FLOAT, image[0][0].get());`	111	`0, GL_RGB, GL_FLOAT, image[0][0].get());`
112	`}`	112	`}`
113		113
114	`void initRaytracer()`	114	`void initRaytracer()`
115	`{`	115	`{`
116	`Vec3f c;`	116	`Vec3f c;`
117	`float r;`	117	`float r;`
118	`mesh.get_bsphere(c, r);`	118	`mesh.get_bsphere(c, r);`
119	`r *= 1.5;`	119	`r *= 1.5;`
120		120
121	`// Initialize track ball`	121	`// Initialize track ball`
122	`vctrl = new GLViewController(winx, winy, c, r);`	122	`vctrl = new GLViewController(winx, winy, c, r);`
123		123
124	`// Initialize corresponding camera`	124	`// Initialize corresponding camera`
125	`cam = new Camera(c - Vec3f(r), c, Vec3f(0.0f, 1.0f, 0.0f), 1.0f);`	125	`cam = new Camera(c - Vec3f(r), c, Vec3f(0.0f, 1.0f, 0.0f), 1.0f);`
126		126
127	`#ifdef USE_BSP`	127	`#ifdef USE_BSP`
128	`cout << "Constructing BSP tree..." << endl;`	128	`cout << "Constructing BSP tree..." << endl;`
129	`tree.init(mesh_vector, transforms, MAX_OBJECTS, MAX_LEVEL);`	129	`tree.init(mesh_vector, transforms, MAX_OBJECTS, MAX_LEVEL);`
130	`tree.build();`	130	`tree.build();`
131	`#else`	131	`#else`
132	`// AABB TREE`	132	`// AABB TREE`
133	`Manifold m;`	133	`Manifold m;`
134	`vector<int> faces(mesh.geometry.no_faces(), 3);`	134	`vector<int> faces(mesh.geometry.no_faces(), 3);`
135	`cout << "Creating manifold" << endl;`	135	`cout << "Creating manifold" << endl;`
136	`build_manifold(m,`	136	`build_manifold(m,`
137	`mesh.geometry.no_vertices(),`	137	`mesh.geometry.no_vertices(),`
138	`&mesh.geometry.vertex(0),`	138	`&mesh.geometry.vertex(0),`
139	`faces.size(), &faces[0],`	139	`faces.size(), &faces[0],`
140	`reinterpret_cast<const int*>(&mesh.geometry.face(0)));`	140	`reinterpret_cast<const int*>(&mesh.geometry.face(0)));`
141	`cout << "Building tree" << endl;`	141	`cout << "Building tree" << endl;`
142	`build_AABBTree(m, bb_tree);`	142	`build_AABBTree(m, bb_tree);`
143	`#endif`	143	`#endif`
144	`}`	144	`}`
145		145
146	`void initGL()`	146	`void initGL()`
147	`{`	147	`{`
148	`glShadeModel(GL_SMOOTH);`	148	`glShadeModel(GL_SMOOTH);`
149	`glDisable(GL_CULL_FACE);`	149	`glDisable(GL_CULL_FACE);`
150	`glFrontFace(GL_CCW);`	150	`glFrontFace(GL_CCW);`
151		151
152	`glClearColor(1.0, 1.0, 1.0, 1.0);`	152	`glClearColor(1.0, 1.0, 1.0, 1.0);`
153	`glColor3f(0.0, 0.0, 0.0);`	153	`glColor3f(0.0, 0.0, 0.0);`
154	`}`	154	`}`
155		155
156		156
157	`//////////////////////////////////////////////////////////////`	157	`//////////////////////////////////////////////////////////////`
158	`// S H A D E F U N C T I O N S`	158	`// S H A D E F U N C T I O N S`
159	`//////////////////////////////////////////////////////////////`	159	`//////////////////////////////////////////////////////////////`
160		160
161	`double shadow_shade(Ray& r)`	161	`double shadow_shade(Ray& r)`
162	`{`	162	`{`
163	`r.compute_position();`	163	`r.compute_position();`
164	`Ray shadow(r.hit_pos, light_dir);`	164	`Ray shadow(r.hit_pos, light_dir);`
165	`double s = tree.intersect(shadow) ? 0.0 : 1.0;`	165	`double s = tree.intersect(shadow) ? 0.0 : 1.0;`
166		166
167	`r.compute_normal();`	167	`r.compute_normal();`
168	`return slight_powdot(r.hit_normal, light_dir);`	168	`return slight_powdot(r.hit_normal, light_dir);`
169	`}`	169	`}`
170		170
171	`double lambertian_shade(Ray& r)`	171	`double lambertian_shade(Ray& r)`
172	`{`	172	`{`
173	`#ifdef USE_BSP`	173	`#ifdef USE_BSP`
174	`r.compute_normal();`	174	`r.compute_normal();`
175	`#endif`	175	`#endif`
176	`return light_pow*dot(r.hit_normal, light_dir);`	176	`return light_pow*dot(r.hit_normal, light_dir);`
177	`}`	177	`}`
178		178
179	`double (*shade_ray[2])(Ray&) = { lambertian_shade,`	179	`double (*shade_ray[2])(Ray&) = { lambertian_shade,`
180	`shadow_shade };`	180	`shadow_shade };`
181		181
182	`//////////////////////////////////////////////////////////////`	182	`//////////////////////////////////////////////////////////////`
183	`// D R A W F U N C T I O N S`	183	`// D R A W F U N C T I O N S`
184	`//////////////////////////////////////////////////////////////`	184	`//////////////////////////////////////////////////////////////`
185		185
186	`/*`	186	`/*`
187	`void enable_textures(TriMesh& tm)`	187	`void enable_textures(TriMesh& tm)`
188	`{`	188	`{`
189	`for(unsigned int i=0;i<tm.materials.size(); ++i)`	189	`for(unsigned int i=0;i<tm.materials.size(); ++i)`
190	`{`	190	`{`
191	`Material& mat = tm.materials[i];`	191	`Material& mat = tm.materials[i];`
192	`if(mat.tex_name != "")`	192	`if(mat.tex_name != "")`
193	`{`	193	`{`
194	`string name = mat.tex_path + mat.tex_name;`	194	`string name = mat.tex_path + mat.tex_name;`
195		195
196	`GLuint tex_id;`	196	`GLuint tex_id;`
197	`if(load_image_into_texture(name, tex_id))`	197	`if(load_image_into_texture(name, tex_id))`
198	`mat.tex_id = tex_id;`	198	`mat.tex_id = tex_id;`
199	`}`	199	`}`
200	`}`	200	`}`
201	`}`	201	`}`
202	`*/`	202	`*/`
203		203
204	`void set_perspective_proj()`	204	`void set_perspective_proj()`
205	`{`	205	`{`
206	`glMatrixMode(GL_PROJECTION);`	206	`glMatrixMode(GL_PROJECTION);`
207	`glLoadIdentity();`	207	`glLoadIdentity();`
208		208
209	`gluPerspective(64.0, 1.0, 0.1, 1000.0);`	209	`gluPerspective(64.0, 1.0, 0.1, 1000.0);`
210		210
211	`glMatrixMode(GL_MODELVIEW);`	211	`glMatrixMode(GL_MODELVIEW);`
212	`}`	212	`}`
213		213
214	`void set_ortho_proj()`	214	`void set_ortho_proj()`
215	`{`	215	`{`
216	`glMatrixMode(GL_PROJECTION);`	216	`glMatrixMode(GL_PROJECTION);`
217	`glLoadIdentity();`	217	`glLoadIdentity();`
218		218
219	`glOrtho(0.0, 1.0, 0.0, 1.0, -1.0, 1.0);`	219	`glOrtho(0.0, 1.0, 0.0, 1.0, -1.0, 1.0);`
220		220
221	`glMatrixMode(GL_MODELVIEW);`	221	`glMatrixMode(GL_MODELVIEW);`
222	`}`	222	`}`
223		223
224	`void draw_texture(unsigned int tex)`	224	`void draw_texture(unsigned int tex)`
225	`{`	225	`{`
226	`static GLfloat verts[] = { 0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 1.0, 1.0 };`	226	`static GLfloat verts[] = { 0.0, 0.0, 1.0, 0.0, 0.0, 1.0, 1.0, 1.0 };`
227		227
228	`glColor4f(1.0, 1.0, 1.0, 1.0);`	228	`glColor4f(1.0, 1.0, 1.0, 1.0);`
229		229
230	`glBindTexture(GL_TEXTURE_2D, tex);`	230	`glBindTexture(GL_TEXTURE_2D, tex);`
231	`glEnable(GL_TEXTURE_2D);`	231	`glEnable(GL_TEXTURE_2D);`
232		232
233	`glEnableClientState(GL_VERTEX_ARRAY);`	233	`glEnableClientState(GL_VERTEX_ARRAY);`
234	`glEnableClientState(GL_TEXTURE_COORD_ARRAY);`	234	`glEnableClientState(GL_TEXTURE_COORD_ARRAY);`
235		235
236	`glVertexPointer(2, GL_FLOAT, 0, verts);`	236	`glVertexPointer(2, GL_FLOAT, 0, verts);`
237	`glTexCoordPointer(2, GL_FLOAT, 0, verts);`	237	`glTexCoordPointer(2, GL_FLOAT, 0, verts);`
238		238
239	`glDrawArrays(GL_TRIANGLE_STRIP, 0, 4);`	239	`glDrawArrays(GL_TRIANGLE_STRIP, 0, 4);`
240		240
241	`glDisableClientState(GL_VERTEX_ARRAY);`	241	`glDisableClientState(GL_VERTEX_ARRAY);`
242	`glDisableClientState(GL_TEXTURE_COORD_ARRAY);`	242	`glDisableClientState(GL_TEXTURE_COORD_ARRAY);`
243		243
244	`glDisable(GL_TEXTURE_2D);`	244	`glDisable(GL_TEXTURE_2D);`
245	`}`	245	`}`
246		246
247	`void drawOBJ()`	247	`void drawOBJ()`
248	`{`	248	`{`
249	`static bool washere = false;`	249	`static bool washere = false;`
250	`static unsigned int disp_list;`	250	`static unsigned int disp_list;`
251		251
252	`if(!washere)`	252	`if(!washere)`
253	`{`	253	`{`
254	`disp_list = glGenLists(1);`	254	`disp_list = glGenLists(1);`
255	`glNewList(disp_list, GL_COMPILE);`	255	`glNewList(disp_list, GL_COMPILE);`
256		256
257	`glBegin(GL_TRIANGLES);`	257	`glBegin(GL_TRIANGLES);`
258	`for(int i = 0; i < mesh.geometry.no_faces(); ++i)`	258	`for(int i = 0; i < mesh.geometry.no_faces(); ++i)`
259	`{`	259	`{`
260	`Vec3i n_face = mesh.normals.face(i);`	260	`Vec3i n_face = mesh.normals.face(i);`
261	`Vec3i g_face = mesh.geometry.face(i);`	261	`Vec3i g_face = mesh.geometry.face(i);`
262	`for(int j=0;j<3;j++)`	262	`for(int j=0;j<3;j++)`
263	`{`	263	`{`
264	`double shade = 0.5;`	264	`double shade = 0.5;`
265		265
266	`if(n_face != Geometry::NULL_FACE)`	266	`if(n_face != Geometry::NULL_FACE)`
267	`{`	267	`{`
268	`Vec3f norm = normalize(mesh.normals.vertex(n_face[j]));`	268	`Vec3f norm = normalize(mesh.normals.vertex(n_face[j]));`
269	`glNormal3fv(norm.get());`	269	`glNormal3fv(norm.get());`
270	`shade = light_pow*dot(norm, light_dir);`	270	`shade = light_pow*dot(norm, light_dir);`
271	`}`	271	`}`
272		272
273	`glColor3d(shade, shade, shade);`	273	`glColor3d(shade, shade, shade);`
274	`Vec3f vert = mesh.geometry.vertex(g_face[j]);`	274	`Vec3f vert = mesh.geometry.vertex(g_face[j]);`
275	`glVertex3fv(vert.get());`	275	`glVertex3fv(vert.get());`
276	`}`	276	`}`
277	`}`	277	`}`
278	`glEnd();`	278	`glEnd();`
279		279
280	`glEndList();`	280	`glEndList();`
281	`washere = true;`	281	`washere = true;`
282	`}`	282	`}`
283	`glCallList(disp_list);`	283	`glCallList(disp_list);`
284	`}`	284	`}`
285		285
286		286
287	`//////////////////////////////////////////////////////////////`	287	`//////////////////////////////////////////////////////////////`
288	`// G L U T C A L L B A C K F U N C T I O N S`	288	`// G L U T C A L L B A C K F U N C T I O N S`
289	`//////////////////////////////////////////////////////////////`	289	`//////////////////////////////////////////////////////////////`
290		290
291	`void display()`	291	`void display()`
292	`{`	292	`{`
293	`static bool first = true;`	293	`static bool first = true;`
294		294
295	`if(first)`	295	`if(first)`
296	`{`	296	`{`
297	`first = false;`	297	`first = false;`
298	`glutTimerFunc(spin_timer, spin, 0);`	298	`glutTimerFunc(spin_timer, spin, 0);`
299	`}`	299	`}`
300		300
301	`Vec3f eye, focus, up;`	301	`Vec3f eye, focus, up;`
302	`vctrl->get_view_param(eye, focus, up);`	302	`vctrl->get_view_param(eye, focus, up);`
303	`cam->set(eye, focus, up, cam->get_focal_dist());`	303	`cam->set(eye, focus, up, cam->get_focal_dist());`
304		304
305	`if(raytrace)`	305	`if(raytrace)`
306	`{`	306	`{`
307	`raytrace = false;`	307	`raytrace = false;`
308		308
309	`cout << "Raytracing";`	309	`cout << "Raytracing";`
310	`float win_to_vp = 1.0f/static_cast<float>(TEX_SIZE);`	310	`float win_to_vp = 1.0f/static_cast<float>(TEX_SIZE);`
311	`float lowerleft = 0.5 + win_to_vp*0.5;`	311	`float lowerleft = 0.5 + win_to_vp*0.5;`
312	`float step = win_to_vp/static_cast<float>(PIXEL_SUBDIVS);`	312	`float step = win_to_vp/static_cast<float>(PIXEL_SUBDIVS);`
313		313
314	`vector<Vec2f> jitter(PIXEL_SUBDIVS*PIXEL_SUBDIVS);`	314	`vector<Vec2f> jitter(PIXEL_SUBDIVS*PIXEL_SUBDIVS);`
315	`for(unsigned int i = 0; i < PIXEL_SUBDIVS; ++i)`	315	`for(unsigned int i = 0; i < PIXEL_SUBDIVS; ++i)`
316	`for(unsigned int j = 0; j < PIXEL_SUBDIVS; ++j)`	316	`for(unsigned int j = 0; j < PIXEL_SUBDIVS; ++j)`
317	`{`	317	`{`
318	`jitter[iPIXEL_SUBDIVS + j][0] = (my_random() + (j%PIXEL_SUBDIVS))step;`	318	`jitter[iPIXEL_SUBDIVS + j][0] = (my_random() + (j%PIXEL_SUBDIVS))step;`
319	`jitter[iPIXEL_SUBDIVS + j][1] = (my_random() + (i%PIXEL_SUBDIVS))step;`	319	`jitter[iPIXEL_SUBDIVS + j][1] = (my_random() + (i%PIXEL_SUBDIVS))step;`
320	`}`	320	`}`
321		321
322	`Util::Timer tim;`	322	`Util::Timer tim;`
323	`tim.start();`	323	`tim.start();`
324	`for(unsigned int i = 0; i < TEX_SIZE; ++i)`	324	`for(unsigned int i = 0; i < TEX_SIZE; ++i)`
325	`{`	325	`{`
326	`for(unsigned int j = 0; j < TEX_SIZE; ++j)`	326	`for(unsigned int j = 0; j < TEX_SIZE; ++j)`
327	`{`	327	`{`
328	`Vec3d sum(0.0f);`	328	`Vec3d sum(0.0f);`
329	`Vec2f vp_pos(jwin_to_vp - lowerleft, iwin_to_vp - lowerleft);`	329	`Vec2f vp_pos(jwin_to_vp - lowerleft, iwin_to_vp - lowerleft);`
330		330
331	`for(unsigned int ky = 0; ky < PIXEL_SUBDIVS; ++ky)`	331	`for(unsigned int ky = 0; ky < PIXEL_SUBDIVS; ++ky)`
332	`for(unsigned int kx = 0; kx < PIXEL_SUBDIVS; ++kx)`	332	`for(unsigned int kx = 0; kx < PIXEL_SUBDIVS; ++kx)`
333	`{`	333	`{`
334	`Ray r = cam->get_ray(vp_pos + jitter[ky*PIXEL_SUBDIVS + kx]);`	334	`Ray r = cam->get_ray(vp_pos + jitter[ky*PIXEL_SUBDIVS + kx]);`
335		335
336	`#ifdef USE_BSP`	336	`#ifdef USE_BSP`
337	`if(tree.intersect(r))`	337	`if(tree.intersect(r))`
338	`sum += Vec3d(shade_ray[shadow](r));`	338	`sum += Vec3d(shade_ray[shadow](r));`
339	`else`	339	`else`
340	`sum += background;`	340	`sum += background;`
341	`#else`	341	`#else`
342	`float t = FLT_MAX;`	342	`float t = FLT_MAX;`
343	`bb_tree.intersect(r);`	343	`bb_tree.intersect(r);`
344	`if(r.has_hit)`	344	`if(r.has_hit)`
345	`sum += Vec3d(shade_ray[0](r));`	345	`sum += Vec3d(shade_ray[0](r));`
346	`else`	346	`else`
347	`sum += background;`	347	`sum += background;`
348	`#endif`	348	`#endif`
349	`}`	349	`}`
350	`image[i][j] = Vec3f(sum/static_cast<double>(PIXEL_SUBDIVS*PIXEL_SUBDIVS));`	350	`image[i][j] = Vec3f(sum/static_cast<double>(PIXEL_SUBDIVS*PIXEL_SUBDIVS));`
351	`}`	351	`}`
352	`if(((i + 1) % 50) == 0) cerr << ".";`	352	`if(((i + 1) % 50) == 0) cerr << ".";`
353	`}`	353	`}`
354	`cout << " - " << tim.get_secs() << " secs " << endl;`	354	`cout << " - " << tim.get_secs() << " secs " << endl;`
355	`cout << endl;`	355	`cout << endl;`
356		356
357	`init_texture(image_tex);`	357	`init_texture(image_tex);`
358		358
359	`done = true;`	359	`done = true;`
360	`}`	360	`}`
361		361
362	`if(done)`	362	`if(done)`
363	`{`	363	`{`
364	`set_ortho_proj();`	364	`set_ortho_proj();`
365	`glClear(GL_COLOR_BUFFER_BIT \| GL_DEPTH_BUFFER_BIT);`	365	`glClear(GL_COLOR_BUFFER_BIT \| GL_DEPTH_BUFFER_BIT);`
366	`glLoadIdentity();`	366	`glLoadIdentity();`
367		367
368	`draw_texture(image_tex);`	368	`draw_texture(image_tex);`
369	`}`	369	`}`
370	`else`	370	`else`
371	`{`	371	`{`
372	`glEnable(GL_DEPTH_TEST);`	372	`glEnable(GL_DEPTH_TEST);`
373		373
374	`cam->glSetPerspective(winx, winy);`	374	`cam->glSetPerspective(winx, winy);`
375		375
376	`glClearColor(background[0], background[1], background[2], 1.0);`	376	`glClearColor(background[0], background[1], background[2], 1.0);`
377	`glClear(GL_COLOR_BUFFER_BIT \| GL_DEPTH_BUFFER_BIT);`	377	`glClear(GL_COLOR_BUFFER_BIT \| GL_DEPTH_BUFFER_BIT);`
378	`glLoadIdentity();`	378	`glLoadIdentity();`
379		379
380	`cam->glSetCamera();`	380	`cam->glSetCamera();`
381		381
382	`glColor3f(0.5, 0.5, 0.5);`	382	`glColor3f(0.5, 0.5, 0.5);`
383	`drawOBJ();`	383	`drawOBJ();`
384		384
385	`glDisable(GL_DEPTH_TEST);`	385	`glDisable(GL_DEPTH_TEST);`
386	`}`	386	`}`
387		387
388	`glutSwapBuffers();`	388	`glutSwapBuffers();`
389	`}`	389	`}`
390		390
391	`void reshape(int w, int h)`	391	`void reshape(int w, int h)`
392	`{`	392	`{`
393	`winx = w; winy = h;`	393	`winx = w; winy = h;`
394		394
395	`vctrl->reshape(winx, winy);`	395	`vctrl->reshape(winx, winy);`
396		396
397	`glViewport(0, 0, winx, winy);`	397	`glViewport(0, 0, winx, winy);`
398	`}`	398	`}`
399		399
400	`void keyboard(unsigned char key, int x, int y)`	400	`void keyboard(unsigned char key, int x, int y)`
401	`{`	401	`{`
402	`switch(key)`	402	`switch(key)`
403	`{`	403	`{`
404	`case '+':`	404	`case '+':`
405	`++PIXEL_SUBDIVS;`	405	`++PIXEL_SUBDIVS;`
406	`cout << "Rays per pixel: " << PIXEL_SUBDIVS*PIXEL_SUBDIVS << endl;`	406	`cout << "Rays per pixel: " << PIXEL_SUBDIVS*PIXEL_SUBDIVS << endl;`
407	`break;`	407	`break;`
408	`case '-':`	408	`case '-':`
409	`if(PIXEL_SUBDIVS > 1)`	409	`if(PIXEL_SUBDIVS > 1)`
410	`--PIXEL_SUBDIVS;`	410	`--PIXEL_SUBDIVS;`
411	`cout << "Rays per pixel: " << PIXEL_SUBDIVS*PIXEL_SUBDIVS << endl;`	411	`cout << "Rays per pixel: " << PIXEL_SUBDIVS*PIXEL_SUBDIVS << endl;`
412	`break;`	412	`break;`
413	`case 'r':`	413	`case 'r':`
414	`if(done) done = false;`	414	`if(done) done = false;`
415	`else raytrace = true;`	415	`else raytrace = true;`
416	`break;`	416	`break;`
417	`case 's':`	417	`case 's':`
418	`shadow = !shadow;`	418	`shadow = !shadow;`
419	`cout << "Shadow " << (shadow ? "on." : "off.") << endl;`	419	`cout << "Shadow " << (shadow ? "on." : "off.") << endl;`
420	`break;`	420	`break;`
421	`case 27:`	421	`case 27:`
422	`delete vctrl;`	422	`delete vctrl;`
423	`delete cam;`	423	`delete cam;`
424	`exit(0);`	424	`exit(0);`
425	`}`	425	`}`
426	`}`	426	`}`
427		427
428	`void mouse(int btn, int state, int x, int y)`	428	`void mouse(int btn, int state, int x, int y)`
429	`{`	429	`{`
430	`if(state == GLUT_DOWN)`	430	`if(state == GLUT_DOWN)`
431	`{`	431	`{`
432	`if(btn == GLUT_LEFT_BUTTON)`	432	`if(btn == GLUT_LEFT_BUTTON)`
433	`vctrl->grab_ball(ROTATE_ACTION, Vec2i(x, y));`	433	`vctrl->grab_ball(ROTATE_ACTION, Vec2i(x, y));`
434	`else if(btn == GLUT_MIDDLE_BUTTON)`	434	`else if(btn == GLUT_MIDDLE_BUTTON)`
435	`vctrl->grab_ball(ZOOM_ACTION, Vec2i(x, y));`	435	`vctrl->grab_ball(ZOOM_ACTION, Vec2i(x, y));`
436	`else if(btn == GLUT_RIGHT_BUTTON)`	436	`else if(btn == GLUT_RIGHT_BUTTON)`
437	`vctrl->grab_ball(PAN_ACTION, Vec2i(x, y));`	437	`vctrl->grab_ball(PAN_ACTION, Vec2i(x, y));`
438	`}`	438	`}`
439	`else if(state == GLUT_UP)`	439	`else if(state == GLUT_UP)`
440	`vctrl->release_ball();`	440	`vctrl->release_ball();`
441		441
442	`mouse_state = state;`	442	`mouse_state = state;`
443	`mouse_button = btn;`	443	`mouse_button = btn;`
444		444
445	`glutPostRedisplay();`	445	`glutPostRedisplay();`
446	`}`	446	`}`
447		447
448	`void move(int x, int y)`	448	`void move(int x, int y)`
449	`{`	449	`{`
450	`vctrl->roll_ball(Vec2i(x, y));`	450	`vctrl->roll_ball(Vec2i(x, y));`
451		451
452	`glutPostRedisplay();`	452	`glutPostRedisplay();`
453	`}`	453	`}`
454		454
455	`void spin(int x)`	455	`void spin(int x)`
456	`{`	456	`{`
457	`vctrl->try_spin();`	457	`vctrl->try_spin();`
458	`glutTimerFunc(spin_timer, spin, 0);`	458	`glutTimerFunc(spin_timer, spin, 0);`
459	`glutPostRedisplay();`	459	`glutPostRedisplay();`
460	`}`	460	`}`
461		461
462		462
463	`//////////////////////////////////////////////////////////////`	463	`//////////////////////////////////////////////////////////////`
464	`// M A I N`	464	`// M A I N`
465	`//////////////////////////////////////////////////////////////`	465	`//////////////////////////////////////////////////////////////`
466		466
467	`int main(int argc, char** argv)`	467	`int main(int argc, char** argv)`
468	`{`	468	`{`
469	`#ifdef __BORLANDC__`	469	`#ifdef __BORLANDC__`
470	`_control87(MCW_EM, MCW_EM); // Borland C++ will crash OpenGL if this`	470	`_control87(MCW_EM, MCW_EM); // Borland C++ will crash OpenGL if this`
471	`// magic line is not inserted`	471	`// magic line is not inserted`
472	`#endif`	472	`#endif`
473		473
474	`glutInit(&argc, argv);`	474	`glutInit(&argc, argv);`
475	`glutInitDisplayMode(GLUT_RGBA \| GLUT_DOUBLE \| GLUT_DEPTH);`	475	`glutInitDisplayMode(GLUT_RGBA \| GLUT_DOUBLE \| GLUT_DEPTH);`
476	`glutInitWindowSize(winx, winy);`	476	`glutInitWindowSize(winx, winy);`
477	`glutCreateWindow("Press 'r' to raytrace");`	477	`glutCreateWindow("Press 'r' to raytrace");`
478	`glutDisplayFunc(display);`	478	`glutDisplayFunc(display);`
479	`glutReshapeFunc(reshape);`	479	`glutReshapeFunc(reshape);`
480	`glutKeyboardFunc(keyboard);`	480	`glutKeyboardFunc(keyboard);`
481	`glutMouseFunc(mouse);`	481	`glutMouseFunc(mouse);`
482	`glutMotionFunc(move);`	482	`glutMotionFunc(move);`
483		483
484	`// LOAD OBJ`	484	`// LOAD OBJ`
485	`string filename;`	485	`string filename;`
486	`if(argc > 1)`	486	`if(argc > 1)`
487	`{`	487	`{`
488	`filename = argv[1];`	488	`filename = argv[1];`
489	`cout << "Loading " << filename << endl;`	489	`cout << "Loading " << filename << endl;`
490		490
491	`obj_load(filename, mesh);`	491	`obj_load(filename, mesh);`
492	`//if(!mesh.has_normals())`	492	`//if(!mesh.has_normals())`
493	`//{`	493	`//{`
494	`cout << "Computing normals" << endl;`	494	`cout << "Computing normals" << endl;`
495	`mesh.compute_normals();`	495	`mesh.compute_normals();`
496	`//}`	496	`//}`
497		497
498	`cout << "No. of triangles: " << mesh.geometry.no_faces() << endl;`	498	`cout << "No. of triangles: " << mesh.geometry.no_faces() << endl;`
499	`}`	499	`}`
500	`else`	500	`else`

Subversion Repositories gelsvn

(root)/branches/cpp11-devel/apps/Raytracer/raytrace.cpp – Rev 630 → 638