libgender/view.C

#include <cmath>

#include "opengl.h"

#include "view.h"
#include "oct.h"

using namespace gl;

vector<GLuint> occ_query_objects;

static GLuint begin_occ_query ()
{
  GLuint id;

  if (occ_query_objects.size ())
    {
      id = *(occ_query_objects.end () - 1);
      occ_query_objects.pop_back ();
    }
  else
    glGenQueriesARB (1, &id);

  glBeginQueryARB (GL_SAMPLES_PASSED, id);
  return id;
}

inline void end_occ_query ()
{
  glEndQueryARB (GL_SAMPLES_PASSED);
}

static GLuint occ_query_result (GLuint id)
{
  GLuint count;

  glGetQueryObjectuivARB (id, GL_QUERY_RESULT, &count);
  occ_query_objects.push_back (id);

  return count;
}

view::view ()
: gamma(1.0), nz_near (1.F), nz_far (2.F)
{
}

view::~view ()
{
}

void view::begin_occ_query (recv_occ_query &recv, void *id)
{
  occ_queries.push_back (oq_data (&recv, ::begin_occ_query (), id));
}

void view::end_occ_query ()
{
  ::end_occ_query ();
}

bool view::may_draw (const entity *e)
{
  if (drawn.find (e) != drawn.end ())
    return false;

  drawn.insert (e);
  return true;
}

visibility_base *visible::get_visibility (view &ctx)
{
  view::visibility_map::iterator i = ctx.vismap.find (this);
  visibility_base *vs;

  if (i == ctx.vismap.end ())
    {
      vs = new_visibility ();
      ctx.vismap.insert (view::visibility_map::value_type (this, vs));
    }
  else
    {
      vs = i->second;

      if (vs->generation != ctx.generation)
        {
          if (vs->generation + 1 != ctx.generation)
            clear_visibility (vs);

          vs->generation = ctx.generation;
        }
    }

  return vs;

}

void view::reset_projection ()
{
  renormalize (orig, p);

  glViewport (0, 0, w, h);

  glMatrixMode (GL_PROJECTION);
  glLoadIdentity ();

  GLdouble aspect = (GLdouble)w/h;
  GLdouble ftan = tan (fov * (.5 * M_PI / 180.));
  GLdouble ymax = z_near * ftan;
  GLdouble xmax = ymax * aspect;

  glFrustum (-xmax, xmax, -ymax, ymax, z_near, z_far);

  perspfact = z_near / ymax * 0.5F * h;

  d = normalize (d);//D
  u = normalize (u);//D

  vec3 rz = -d;
  vec3 rx = cross (u, rz);
  vec3 ry = cross (rz, rx);

  matrix &m = projection;
  m(0,0) = rx.x; m(0,1) = rx.y; m(0,2) = rx.z; m(0,3) = 0;
  m(1,0) = ry.x; m(1,1) = ry.y; m(1,2) = ry.z; m(1,3) = 0;
  m(2,0) = rz.x; m(2,1) = rz.y; m(2,2) = rz.z; m(2,3) = 0;
  m(3,0) =    0; m(3,1) =    0; m(3,2) =    0; m(3,3) = 1;

  diagfact = abs (rz.x) + abs (rz.y) + abs (rz.z);
  //printf ("diagfact = %f\n", diagfact);
  diagfact = sqrtf (3.);//D WHY???

  glMultMatrixf (m);
  glTranslatef (-p.x, -p.y, -p.z);

  glGetFloatv (GL_PROJECTION_MATRIX, m);

  frustum.l = plane ( m(3,0) + m(0,0), m(3,1) + m(0,1), m(3,2) + m(0,2), m(3,3) + m(0,3) );
  frustum.r = plane ( m(3,0) - m(0,0), m(3,1) - m(0,1), m(3,2) - m(0,2), m(3,3) - m(0,3) );
  frustum.b = plane ( m(3,0) + m(1,0), m(3,1) + m(1,1), m(3,2) + m(1,2), m(3,3) + m(1,3) );
  frustum.t = plane ( m(3,0) - m(1,0), m(3,1) - m(1,1), m(3,2) - m(1,2), m(3,3) - m(1,3) );
  frustum.n = plane ( m(3,0) + m(2,0), m(3,1) + m(2,1), m(3,2) + m(2,2), m(3,3) + m(2,3) );
  frustum.f = plane ( m(3,0) - m(2,0), m(3,1) - m(2,1), m(3,2) - m(2,2), m(3,3) - m(2,3) );

  {
    GLdouble frustlen = z_far - z_near;
    GLdouble fheight  = frustlen * ftan;
    GLdouble fwidth   = fheight * w / h;
    point half (0, 0, z_near + frustlen * .5);
    point corner (fwidth, fheight, frustlen);
    
    frustum.s = sphere (p + d * (.5 * frustlen), length (corner - half));
  }

  {
    GLdouble depth = h / ftan;

    frustum.c = cone (p, d, atan (sqrt (GLdouble (w * w + h * h)) * ftan / h));
  }

  glMatrixMode (GL_MODELVIEW);
  glLoadIdentity ();
}

void view::begin ()
{
  generation++;

  glColorMask (1, 1, 1, 1);
  glDepthMask (1);
  glClear (GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT | GL_STENCIL_BUFFER_BIT);

  vislist.clear ();

  z_near = max (nz_near, 1.F);
  z_far = max (nz_far, z_near * 2.F);
  c_far = nc_far;

  reset_projection ();

  nz_near = 1e37;
  nc_far = nz_far = 1.F;
}

void view::end ()
{
  vislist.clear ();
}

#define DEPTH_OFFSET (1. / (GLdouble)(1L << 16))

void view::render (enum pass p)
{
  pass = p;

  switch (pass)
    {
      case DEPTH:
        //glEnable (GL_POLYGON_OFFSET_FILL);
        //glPolygonOffset (0, 5);
        glDisable (GL_MINMAX);
        glDepthRange (DEPTH_OFFSET, 1.);
        glDepthFunc (GL_LESS);
        glEnable (GL_DEPTH_TEST);
        glColorMask (0, 0, 0, 0);
        glDepthMask (1);
        glDisable (GL_DEPTH_CLAMP_NV);

        world.detect_visibility (*this);

        printf ("fps %f NF %f:%f vis %d CAM (%d,%d,%d)\n", timer.fps, z_near, z_far, drawn.size (), orig.x, orig.y, orig.z);//D
        break;

      case POSTDEPTH:
        glDepthRange (0., 1. - DEPTH_OFFSET);
        glDepthFunc (GL_LESS);
        glColorMask (0, 0, 0, 0);
        glDepthMask (0);
        glEnable (GL_DEPTH_CLAMP_NV);
       
        // check occlusion queries
        for (vector<oq_data>::iterator i = occ_queries.begin (); i != occ_queries.end (); ++i)
          {
            occ_query oq(*this, i->data, ::occ_query_result (i->id));
            i->recv->event (oq);
          }

        occ_queries.clear ();

        break;

      case LIGHTED:
        glEnable (GL_MINMAX);
        glDepthRange (0., 1. - DEPTH_OFFSET);
        glDepthFunc (GL_LESS);
        glColorMask (1, 1, 1, 1);
        glDepthMask (0);
        glDisable (GL_DEPTH_CLAMP_NV);

        break;
    }
  
  for (vector<octant *>::iterator i = vislist.begin (); i != vislist.end (); ++i)
    (*i)->display (*this);

  drawn.clear ();

#if 0
  if (pass == view::DEPTH)
    {
      glEnable (GL_DEPTH_CLAMP_NV);
      glDepthMask (0);
      glDepthFunc (GL_LESS);

      for (vector<octant *>::iterator i = farlist.begin (); i != farlist.end (); ++i)
        {
          begin_occ_query (**i);
          (*i)->draw_bbox (*this);
          end_occ_query ();
       }

      glDisable (GL_DEPTH_CLAMP_NV);
    }
#endif
}


Revision:	1.61
Committed:	Tue Oct 19 15:34:58 2004 UTC (19 years, 7 months ago) by root
Content type:	text/plain
Branch:	MAIN
Changes since 1.60:	+5 -8 lines
Log Message:	* empty log message *
#	Content
1	#include <cmath>
2
3	#include "opengl.h"
4
5	#include "view.h"
6	#include "oct.h"
7
8	using namespace gl;
9
10	vector<GLuint> occ_query_objects;
11
12	static GLuint begin_occ_query ()
13	{
14	GLuint id;
15
16	if (occ_query_objects.size ())
17	{
18	id = *(occ_query_objects.end () - 1);
19	occ_query_objects.pop_back ();
20	}
21	else
22	glGenQueriesARB (1, &id);
23
24	glBeginQueryARB (GL_SAMPLES_PASSED, id);
25	return id;
26	}
27
28	inline void end_occ_query ()
29	{
30	glEndQueryARB (GL_SAMPLES_PASSED);
31	}
32
33	static GLuint occ_query_result (GLuint id)
34	{
35	GLuint count;
36
37	glGetQueryObjectuivARB (id, GL_QUERY_RESULT, &count);
38	occ_query_objects.push_back (id);
39
40	return count;
41	}
42
43	view::view ()
44	: gamma(1.0), nz_near (1.F), nz_far (2.F)
45	{
46	}
47
48	view::~view ()
49	{
50	}
51
52	void view::begin_occ_query (recv_occ_query &recv, void *id)
53	{
54	occ_queries.push_back (oq_data (&recv, ::begin_occ_query (), id));
55	}
56
57	void view::end_occ_query ()
58	{
59	::end_occ_query ();
60	}
61
62	bool view::may_draw (const entity *e)
63	{
64	if (drawn.find (e) != drawn.end ())
65	return false;
66
67	drawn.insert (e);
68	return true;
69	}
70
71	visibility_base *visible::get_visibility (view &ctx)
72	{
73	view::visibility_map::iterator i = ctx.vismap.find (this);
74	visibility_base *vs;
75
76	if (i == ctx.vismap.end ())
77	{
78	vs = new_visibility ();
79	ctx.vismap.insert (view::visibility_map::value_type (this, vs));
80	}
81	else
82	{
83	vs = i->second;
84
85	if (vs->generation != ctx.generation)
86	{
87	if (vs->generation + 1 != ctx.generation)
88	clear_visibility (vs);
89
90	vs->generation = ctx.generation;
91	}
92	}
93
94	return vs;
95
96	}
97
98	void view::reset_projection ()
99	{
100	renormalize (orig, p);
101
102	glViewport (0, 0, w, h);
103
104	glMatrixMode (GL_PROJECTION);
105	glLoadIdentity ();
106
107	GLdouble aspect = (GLdouble)w/h;
108	GLdouble ftan = tan (fov * (.5 * M_PI / 180.));
109	GLdouble ymax = z_near * ftan;
110	GLdouble xmax = ymax * aspect;
111
112	glFrustum (-xmax, xmax, -ymax, ymax, z_near, z_far);
113
114	perspfact = z_near / ymax * 0.5F * h;
115
116	d = normalize (d);//D
117	u = normalize (u);//D
118
119	vec3 rz = -d;
120	vec3 rx = cross (u, rz);
121	vec3 ry = cross (rz, rx);
122
123	matrix &m = projection;
124	m(0,0) = rx.x; m(0,1) = rx.y; m(0,2) = rx.z; m(0,3) = 0;
125	m(1,0) = ry.x; m(1,1) = ry.y; m(1,2) = ry.z; m(1,3) = 0;
126	m(2,0) = rz.x; m(2,1) = rz.y; m(2,2) = rz.z; m(2,3) = 0;
127	m(3,0) = 0; m(3,1) = 0; m(3,2) = 0; m(3,3) = 1;
128
129	diagfact = abs (rz.x) + abs (rz.y) + abs (rz.z);
130	//printf ("diagfact = %f\n", diagfact);
131	diagfact = sqrtf (3.);//D WHY???
132
133	glMultMatrixf (m);
134	glTranslatef (-p.x, -p.y, -p.z);
135
136	glGetFloatv (GL_PROJECTION_MATRIX, m);
137
138	frustum.l = plane ( m(3,0) + m(0,0), m(3,1) + m(0,1), m(3,2) + m(0,2), m(3,3) + m(0,3) );
139	frustum.r = plane ( m(3,0) - m(0,0), m(3,1) - m(0,1), m(3,2) - m(0,2), m(3,3) - m(0,3) );
140	frustum.b = plane ( m(3,0) + m(1,0), m(3,1) + m(1,1), m(3,2) + m(1,2), m(3,3) + m(1,3) );
141	frustum.t = plane ( m(3,0) - m(1,0), m(3,1) - m(1,1), m(3,2) - m(1,2), m(3,3) - m(1,3) );
142	frustum.n = plane ( m(3,0) + m(2,0), m(3,1) + m(2,1), m(3,2) + m(2,2), m(3,3) + m(2,3) );
143	frustum.f = plane ( m(3,0) - m(2,0), m(3,1) - m(2,1), m(3,2) - m(2,2), m(3,3) - m(2,3) );
144
145	{
146	GLdouble frustlen = z_far - z_near;
147	GLdouble fheight = frustlen * ftan;
148	GLdouble fwidth = fheight * w / h;
149	point half (0, 0, z_near + frustlen * .5);
150	point corner (fwidth, fheight, frustlen);
151
152	frustum.s = sphere (p + d * (.5 * frustlen), length (corner - half));
153	}
154
155	{
156	GLdouble depth = h / ftan;
157
158	frustum.c = cone (p, d, atan (sqrt (GLdouble (w * w + h * h)) * ftan / h));
159	}
160
161	glMatrixMode (GL_MODELVIEW);
162	glLoadIdentity ();
163	}
164
165	void view::begin ()
166	{
167	generation++;
168
169	glColorMask (1, 1, 1, 1);
170	glDepthMask (1);
171	glClear (GL_COLOR_BUFFER_BIT \| GL_DEPTH_BUFFER_BIT \| GL_STENCIL_BUFFER_BIT);
172
173	vislist.clear ();
174
175	z_near = max (nz_near, 1.F);
176	z_far = max (nz_far, z_near * 2.F);
177	c_far = nc_far;
178
179	reset_projection ();
180
181	nz_near = 1e37;
182	nc_far = nz_far = 1.F;
183	}
184
185	void view::end ()
186	{
187	vislist.clear ();
188	}
189
190	#define DEPTH_OFFSET (1. / (GLdouble)(1L << 16))
191
192	void view::render (enum pass p)
193	{
194	pass = p;
195
196	switch (pass)
197	{
198	case DEPTH:
199	//glEnable (GL_POLYGON_OFFSET_FILL);
200	//glPolygonOffset (0, 5);
201	glDisable (GL_MINMAX);
202	glDepthRange (DEPTH_OFFSET, 1.);
203	glDepthFunc (GL_LESS);
204	glEnable (GL_DEPTH_TEST);
205	glColorMask (0, 0, 0, 0);
206	glDepthMask (1);
207	glDisable (GL_DEPTH_CLAMP_NV);
208
209	world.detect_visibility (*this);
210
211	printf ("fps %f NF %f:%f vis %d CAM (%d,%d,%d)\n", timer.fps, z_near, z_far, drawn.size (), orig.x, orig.y, orig.z);//D
212	break;
213
214	case POSTDEPTH:
215	glDepthRange (0., 1. - DEPTH_OFFSET);
216	glDepthFunc (GL_LESS);
217	glColorMask (0, 0, 0, 0);
218	glDepthMask (0);
219	glEnable (GL_DEPTH_CLAMP_NV);
220
221	// check occlusion queries
222	for (vector<oq_data>::iterator i = occ_queries.begin (); i != occ_queries.end (); ++i)
223	{
224	occ_query oq(*this, i->data, ::occ_query_result (i->id));
225	i->recv->event (oq);
226	}
227
228	occ_queries.clear ();
229
230	break;
231
232	case LIGHTED:
233	glEnable (GL_MINMAX);
234	glDepthRange (0., 1. - DEPTH_OFFSET);
235	glDepthFunc (GL_LESS);
236	glColorMask (1, 1, 1, 1);
237	glDepthMask (0);
238	glDisable (GL_DEPTH_CLAMP_NV);
239
240	break;
241	}
242
243	for (vector<octant *>::iterator i = vislist.begin (); i != vislist.end (); ++i)
244	(i)->display (this);
245
246	drawn.clear ();
247
248	#if 0
249	if (pass == view::DEPTH)
250	{
251	glEnable (GL_DEPTH_CLAMP_NV);
252	glDepthMask (0);
253	glDepthFunc (GL_LESS);
254
255	for (vector<octant *>::iterator i = farlist.begin (); i != farlist.end (); ++i)
256	{
257	begin_occ_query (**i);
258	(i)->draw_bbox (this);
259	end_occ_query ();
260	}
261
262	glDisable (GL_DEPTH_CLAMP_NV);
263	}
264	#endif
265	}
266
267