libgender/util.h

#ifndef UTIL_H
#define UTIL_H

#include <cmath>
#include <cfloat>
#include <cstdlib>
#include <vector>
#include <string>

#include "opengl.h"

#include <SDL/SDL_image.h>

using namespace std;

extern CGcontext cgc;
extern CGprogram vsh;
extern CGprogram fsh;
extern CGparameter mv, mvp, lightpos;
extern CGprofile vsh_profile, fsh_profile;

typedef long long soffs; // 32 bit
typedef unsigned long long uoffs;
#define OFFS_BITS 63
#define SOFFS_MIN (soffs)-(1LL << (OFFS_BITS - 2))
#define SOFFS_MAX (soffs)+(1LL << (OFFS_BITS - 2))
#define MAXEXTENT (1ULL << (OFFS_BITS - 1))

#define ABS(n) ((n) < 0 ? -(n) : (n))

struct sector
{
  soffs x, y, z;

  sector (soffs x, soffs y, soffs z) : x(x), y(y), z(z) { };
  sector (soffs xyz = 0) : x(xyz), y(xyz), z(xyz) { };

  void offset (int subindex, uoffs extent)
  {
    if (subindex & 1) x += extent;
    if (subindex & 2) y += extent;
    if (subindex & 4) z += extent;
  }
};

inline sector operator +(const sector &a, const sector &b)
{
  return sector (a.x + b.x, a.y + b.y, a.z + b.z);
}

inline sector operator -(const sector &a, const sector &b)
{
  return sector (a.x - b.x, a.y - b.y, a.z - b.z);
}

inline sector operator /(const sector &a, soffs b)
{
  return sector (a.x / b, a.y / b, a.z / b);
}

inline sector operator >>(const sector &a, unsigned int b)
{
  return sector (a.x >> b, a.y >> b, a.z >> b);
}

inline bool operator <=(const sector &a, const sector &b)
{
  return a.x <= b.x && a.y <= b.y && a.z <= b.z;
}

inline soffs max (const sector &a)
{
  return max (a.x, max (a.y, a.z));
}

inline sector translate (const sector &p, const sector &src, const sector &dst)
{
  return p + (dst - src);
}

inline sector abs (const sector &s)
{
  return sector (ABS (s.x), ABS (s.y), ABS (s.z));
}

struct vec3
{
  GLfloat x, y, z;
  vec3 () { };
  vec3 (GLfloat s) : x(s), y(s), z(s) { };
  vec3 (GLfloat x, GLfloat y, GLfloat z) : x(x), y(y), z(z) { };
  vec3 (const sector &s) : x(s.x), y(s.y), z(s.z) { };

  const vec3 operator -() const
  { return vec3 (-x, -y, -z); }
};

const vec3 normalize (const vec3 &v);
const vec3 cross (const vec3 &a, const vec3 &b);

inline const vec3 operator *(const vec3 &a, GLfloat s)
{
  return vec3 (a.x * s, a.y * s, a.z * s);
}

inline const vec3 operator +(const vec3 &a, const vec3 &b)
{
  return vec3 (a.x + b.x, a.y + b.y, a.z + b.z);
}

inline const vec3 operator -(const vec3 &a, const vec3 &b)
{
  return vec3 (a.x - b.x, a.y - b.y, a.z - b.z);
}

inline GLfloat dot (const vec3 &a, const vec3 &b)
{
  return a.x * b.x + a.y * b.y + a.z * b.z;
}

// squared length
inline const GLfloat length2 (const vec3 &v)
{
  return dot (v, v);
}

inline const GLfloat length (const vec3 &v)
{
  return sqrtf (length2 (v));
}

typedef vec3 point;

// a generic plane
struct plane
{
  vec3 n;
  GLfloat d;

  plane () { };
  plane (GLfloat a, GLfloat b, GLfloat c, GLfloat d);
};

inline GLfloat distance (const plane &a, const point &p)
{
  return dot (a.n, p) + a.d;
}

struct sphere {
  point p;
  GLfloat r;

  sphere () { };
  sphere (const point &p, GLfloat r) : p(p), r(r) { };
};

inline bool overlap (const sphere &a, const sphere &b)
{
  GLfloat d = a.r + b.r;

  return length2 (a.p - b.p) <= d * d;
}

struct cone {
  point p;
  vec3 a; // axis
  GLfloat f, fs, fc; // angle

  cone () { };
  cone (const point &p, const vec3 &a, GLfloat f)
  : p(p), a(a), f(f)
  {
    fs = sinf (f);
    fc = cosf (f);
  };
};

inline bool overlap (const cone &c, const sphere &s)
{
  vec3 u = c.p - c.a * (s.r / c.fs);
  vec3 d = s.p - u;
  return dot (c.a, d) >= length (d) * c.fc;
}

void renormalize (sector &s, point &p);

struct colour
{
  GLubyte r, g, b, a;
  colour (GLfloat r = 1., GLfloat g = 1., GLfloat b = 1., GLfloat a = 1.)
  : r(GLubyte (r * 255.F + .5F))
  , g(GLubyte (g * 255.F + .5F))
  , b(GLubyte (b * 255.F + .5F))
  , a(GLubyte (a * 255.F + .5F))
  {
  }
};

struct tex2
{
  GLfloat s, t;
  tex2 () { };
  tex2 (GLfloat s, GLfloat t) : s(s), t(t) { };
};

struct box
{
  point a, b;

  box() { };

  void reset ()
  {
    a = point ( FLT_MAX,  FLT_MAX,  FLT_MAX);
    b = point (-FLT_MAX, -FLT_MAX, -FLT_MAX);
  }

  void add (const box &o);
  void add (const point &p);
};

struct entity;
struct geometry;
struct view;
struct octant;

extern struct timer
{
  static double now;
  static double diff;
  static double fps;

  static void frame ();
  timer ();
} timer;

/*
#define MAX_EVENT_TYPES 10
enum event_type { TIMER_EV };
struct event
{
   event_type type;
};

typedef callback1<void, event&> event_cb;

class skedjuhlar
{
   public:
// only 10 types for now
   private:
      vector <list<event_cb> > event_lists;
   
   public:
      skedjuhlar () {
         event_lists.resize (MAX_EVENT_TYPES, list<event_cb>());
      }
      
      void register_event_cb (const event_type &t, const event_cb &e) {
         event_lists[t].push_back (e);
      };
      void send_event (event &e) {
         list<event_cb> &l = event_lists[e.type];
         for (list<event_cb>::iterator it = l.begin (); it != l.end (); it++) {
            (*it)(e);
         }
      };
      void check_events () {
         while (!events.empty ()) {
            event &e = events.pop_front ();
            list<event_cb> &l = event_lists[e->name];
            for (list<event_cb>::iterator it = l.begin (); it !+ l.end (); it++) {
               (*it)(e);
            }
            delete e; // ugly slow? hai hai..... 183G
         }
      }
};

extern skedjuhlar main_scheduler;
*/

namespace gl
{
#ifdef DEBUG
  extern int nesting;
  void errchk (const char *name, const char *args, const char *file, int line);
#endif

  struct vertex_v3f
  {
    point p; // vertex

    vertex_v3f () { };
    vertex_v3f (point p) : p(p) { };

    GLenum gl_format () const { return GL_V3F; }
  };

  struct vertex_t2f_n3f_v3f
  {
    tex2  t; // texture
    vec3  n; // normal
    point p; // vertex

    vertex_t2f_n3f_v3f () { };
    vertex_t2f_n3f_v3f (point p, vec3 n, tex2 t = tex2()) : p(p), n(n), t(t) { };

    GLenum gl_format () const { return GL_T2F_N3F_V3F; }
  };

  struct matrix
  {
    GLfloat data[4][4];

    const GLfloat operator ()(int i, int j) const { return data[j][i]; };
    GLfloat &operator ()(int i, int j) { return data[j][i]; };

    operator GLfloat *() { return &data[0][0]; }

    void diagonal (GLfloat v);
    void clear () { diagonal (0.F); };
    void identity () { diagonal (1.F); };

    void print (); // ugly

    static const matrix translation (const vec3 &v);
    static const matrix rotation (GLfloat degrees, const vec3 &axis);

    matrix () { };
    matrix (GLfloat diag) { diagonal (diag); };
  };

  const matrix operator *(const matrix &a, const matrix &b);
  const vec3 operator *(const matrix &a, const vec3 &v);

  struct vertex_buffer_object {
    GLuint buffer;
    GLenum format;

    void alloc ()
    {
      if (!buffer)
        glGenBuffersARB (1, &buffer);
    }

    void bind (GLenum target = GL_ARRAY_BUFFER_ARB)
    {
      glBindBufferARB (target, buffer);
      glInterleavedArrays (format, 0, 0);
    }

    void draw (GLenum mode, GLint first, GLsizei count)
    {
      bind ();
      glDrawArrays (mode, first, count);
    }

    template<class vertex>
    void set (const vertex *v, GLsizei count, GLenum usage = GL_STATIC_DRAW_ARB, GLenum target = GL_ARRAY_BUFFER_ARB)
    {
      alloc ();
      format = v->gl_format ();
      glBindBufferARB (target, buffer);
      glBufferDataARB (target, count * sizeof (vertex), v, usage);
    }

    template<class vertex>
    void set (const vector<vertex> &v, GLenum usage = GL_STATIC_DRAW_ARB, GLenum target = GL_ARRAY_BUFFER_ARB)
    {
      set (&v[0], v.size (), usage, target);
    }

    vertex_buffer_object ()
    : buffer(0)
    {
    }
    
    ~vertex_buffer_object ()
    {
      glDeleteBuffersARB (1, &buffer);
    }

    operator GLint ()
    {
      return buffer;
    }
  };

  void draw_bbox (vertex_buffer_object &vb, const sector &a, const sector &b);
}

GLuint SDL_GL_LoadTexture (SDL_Surface *surface, GLfloat *texcoord);

#endif

Revision:	1.45
Committed:	Mon Oct 18 12:01:14 2004 UTC (19 years, 7 months ago) by root
Content type:	text/plain
Branch:	MAIN
Changes since 1.44:	+49 -7 lines
Log Message:	* empty log message *
#	User	Rev	Content
1	root	1.1	#ifndef UTIL_H
2			#define UTIL_H
3
4	root	1.7	#include <cmath>
5	root	1.27	#include <cfloat>
6	root	1.39	#include <cstdlib>
7	root	1.1	#include <vector>
8	root	1.22	#include <string>
9	root	1.1
10	root	1.37	#include "opengl.h"
11	root	1.18
12	root	1.37	#include <SDL/SDL_image.h>
13	root	1.30
14	root	1.1	using namespace std;
15	root	1.18
16			extern CGcontext cgc;
17			extern CGprogram vsh;
18			extern CGprogram fsh;
19	root	1.19	extern CGparameter mv, mvp, lightpos;
20	root	1.21	extern CGprofile vsh_profile, fsh_profile;
21	root	1.1
22	root	1.36	typedef long long soffs; // 32 bit
23			typedef unsigned long long uoffs;
24	root	1.39	#define OFFS_BITS 63
25	root	1.36	#define SOFFS_MIN (soffs)-(1LL << (OFFS_BITS - 2))
26			#define SOFFS_MAX (soffs)+(1LL << (OFFS_BITS - 2))
27			#define MAXEXTENT (1ULL << (OFFS_BITS - 1))
28	root	1.1
29	root	1.39	#define ABS(n) ((n) < 0 ? -(n) : (n))
30
31	root	1.37	struct sector
32			{
33	root	1.32	soffs x, y, z;
34
35			sector (soffs x, soffs y, soffs z) : x(x), y(y), z(z) { };
36			sector (soffs xyz = 0) : x(xyz), y(xyz), z(xyz) { };
37
38			void offset (int subindex, uoffs extent)
39			{
40			if (subindex & 1) x += extent;
41			if (subindex & 2) y += extent;
42			if (subindex & 4) z += extent;
43			}
44			};
45
46			inline sector operator +(const sector &a, const sector &b)
47			{
48			return sector (a.x + b.x, a.y + b.y, a.z + b.z);
49			}
50
51			inline sector operator -(const sector &a, const sector &b)
52			{
53			return sector (a.x - b.x, a.y - b.y, a.z - b.z);
54			}
55
56			inline sector operator /(const sector &a, soffs b)
57			{
58			return sector (a.x / b, a.y / b, a.z / b);
59			}
60
61	root	1.40	inline sector operator >>(const sector &a, unsigned int b)
62			{
63			return sector (a.x >> b, a.y >> b, a.z >> b);
64			}
65
66	root	1.32	inline bool operator <=(const sector &a, const sector &b)
67			{
68			return a.x <= b.x && a.y <= b.y && a.z <= b.z;
69			}
70
71			inline soffs max (const sector &a)
72			{
73			return max (a.x, max (a.y, a.z));
74			}
75
76			inline sector translate (const sector &p, const sector &src, const sector &dst)
77			{
78			return p + (dst - src);
79			}
80
81			inline sector abs (const sector &s)
82			{
83	root	1.39	return sector (ABS (s.x), ABS (s.y), ABS (s.z));
84	root	1.32	}
85	root	1.1
86	root	1.37	struct vec3
87			{
88	root	1.1	GLfloat x, y, z;
89			vec3 () { };
90	root	1.32	vec3 (GLfloat s) : x(s), y(s), z(s) { };
91	root	1.1	vec3 (GLfloat x, GLfloat y, GLfloat z) : x(x), y(y), z(z) { };
92	root	1.32	vec3 (const sector &s) : x(s.x), y(s.y), z(s.z) { };
93	root	1.1
94			const vec3 operator -() const
95			{ return vec3 (-x, -y, -z); }
96			};
97
98			const vec3 normalize (const vec3 &v);
99			const vec3 cross (const vec3 &a, const vec3 &b);
100	root	1.6
101			inline const vec3 operator *(const vec3 &a, GLfloat s)
102			{
103			return vec3 (a.x * s, a.y * s, a.z * s);
104			}
105
106	root	1.27	inline const vec3 operator +(const vec3 &a, const vec3 &b)
107			{
108			return vec3 (a.x + b.x, a.y + b.y, a.z + b.z);
109			}
110
111			inline const vec3 operator -(const vec3 &a, const vec3 &b)
112			{
113			return vec3 (a.x - b.x, a.y - b.y, a.z - b.z);
114			}
115
116	root	1.6	inline GLfloat dot (const vec3 &a, const vec3 &b)
117			{
118			return a.x * b.x + a.y * b.y + a.z * b.z;
119			}
120
121	root	1.45	// squared length
122			inline const GLfloat length2 (const vec3 &v)
123	root	1.6	{
124	root	1.45	return dot (v, v);
125			}
126
127			inline const GLfloat length (const vec3 &v)
128			{
129			return sqrtf (length2 (v));
130	root	1.6	}
131	root	1.1
132	root	1.10	typedef vec3 point;
133
134			// a generic plane
135	root	1.37	struct plane
136			{
137	root	1.10	vec3 n;
138			GLfloat d;
139
140			plane () { };
141			plane (GLfloat a, GLfloat b, GLfloat c, GLfloat d);
142			};
143	root	1.45
144			inline GLfloat distance (const plane &a, const point &p)
145			{
146			return dot (a.n, p) + a.d;
147			}
148
149			struct sphere {
150			point p;
151			GLfloat r;
152
153			sphere () { };
154			sphere (const point &p, GLfloat r) : p(p), r(r) { };
155			};
156
157			inline bool overlap (const sphere &a, const sphere &b)
158			{
159			GLfloat d = a.r + b.r;
160
161			return length2 (a.p - b.p) <= d * d;
162			}
163
164			struct cone {
165			point p;
166			vec3 a; // axis
167			GLfloat f, fs, fc; // angle
168
169			cone () { };
170			cone (const point &p, const vec3 &a, GLfloat f)
171			: p(p), a(a), f(f)
172			{
173			fs = sinf (f);
174			fc = cosf (f);
175			};
176			};
177
178			inline bool overlap (const cone &c, const sphere &s)
179			{
180			vec3 u = c.p - c.a * (s.r / c.fs);
181			vec3 d = s.p - u;
182			return dot (c.a, d) >= length (d) * c.fc;
183			}
184	root	1.10
185			void renormalize (sector &s, point &p);
186
187	root	1.37	struct colour
188			{
189	root	1.40	GLubyte r, g, b, a;
190			colour (GLfloat r = 1., GLfloat g = 1., GLfloat b = 1., GLfloat a = 1.)
191			: r(GLubyte (r * 255.F + .5F))
192			, g(GLubyte (g * 255.F + .5F))
193			, b(GLubyte (b * 255.F + .5F))
194			, a(GLubyte (a * 255.F + .5F))
195			{
196			}
197	root	1.10	};
198	root	1.5
199	root	1.40	struct tex2
200	root	1.37	{
201	root	1.1	GLfloat s, t;
202	root	1.40	tex2 () { };
203			tex2 (GLfloat s, GLfloat t) : s(s), t(t) { };
204	root	1.1	};
205
206	root	1.37	struct box
207			{
208	root	1.27	point a, b;
209	root	1.8
210			box() { };
211	root	1.1
212			void reset ()
213			{
214	root	1.27	a = point ( FLT_MAX, FLT_MAX, FLT_MAX);
215			b = point (-FLT_MAX, -FLT_MAX, -FLT_MAX);
216	root	1.1	}
217
218			void add (const box &o);
219			void add (const point &p);
220			};
221	root	1.5
222	root	1.27	struct entity;
223			struct geometry;
224	root	1.20	struct view;
225	root	1.32	struct octant;
226	root	1.7
227	root	1.37	extern struct timer
228			{
229	root	1.7	static double now;
230			static double diff;
231	root	1.44	static double fps;
232	root	1.7
233			static void frame ();
234			timer ();
235			} timer;
236	root	1.11
237	root	1.13	/*
238	root	1.11	#define MAX_EVENT_TYPES 10
239			enum event_type { TIMER_EV };
240	root	1.37	struct event
241			{
242	root	1.11	event_type type;
243			};
244
245			typedef callback1<void, event&> event_cb;
246
247	root	1.37	class skedjuhlar
248			{
249	root	1.11	public:
250			// only 10 types for now
251			private:
252			vector <list<event_cb> > event_lists;
253
254			public:
255			skedjuhlar () {
256			event_lists.resize (MAX_EVENT_TYPES, list<event_cb>());
257			}
258
259			void register_event_cb (const event_type &t, const event_cb &e) {
260			event_lists[t].push_back (e);
261			};
262			void send_event (event &e) {
263			list<event_cb> &l = event_lists[e.type];
264			for (list<event_cb>::iterator it = l.begin (); it != l.end (); it++) {
265			(*it)(e);
266			}
267			};
268			void check_events () {
269			while (!events.empty ()) {
270			event &e = events.pop_front ();
271			list<event_cb> &l = event_lists[e->name];
272			for (list<event_cb>::iterator it = l.begin (); it !+ l.end (); it++) {
273			(*it)(e);
274			}
275			delete e; // ugly slow? hai hai..... 183G
276			}
277			}
278			};
279
280			extern skedjuhlar main_scheduler;
281	root	1.13	*/
282	root	1.22
283	root	1.37	namespace gl
284			{
285	root	1.41	#ifdef DEBUG
286			extern int nesting;
287			void errchk (const char name, const char args, const char *file, int line);
288			#endif
289
290	root	1.40	struct vertex_v3f
291			{
292			point p; // vertex
293
294			vertex_v3f () { };
295			vertex_v3f (point p) : p(p) { };
296
297			GLenum gl_format () const { return GL_V3F; }
298			};
299
300			struct vertex_t2f_n3f_v3f
301			{
302			tex2 t; // texture
303			vec3 n; // normal
304			point p; // vertex
305
306			vertex_t2f_n3f_v3f () { };
307			vertex_t2f_n3f_v3f (point p, vec3 n, tex2 t = tex2()) : p(p), n(n), t(t) { };
308
309			GLenum gl_format () const { return GL_T2F_N3F_V3F; }
310			};
311	root	1.34
312	root	1.41	struct matrix
313			{
314			GLfloat data[4][4];
315
316			const GLfloat operator ()(int i, int j) const { return data[j][i]; };
317			GLfloat &operator ()(int i, int j) { return data[j][i]; };
318
319			operator GLfloat *() { return &data[0][0]; }
320
321			void diagonal (GLfloat v);
322			void clear () { diagonal (0.F); };
323			void identity () { diagonal (1.F); };
324
325			void print (); // ugly
326
327			static const matrix translation (const vec3 &v);
328			static const matrix rotation (GLfloat degrees, const vec3 &axis);
329
330			matrix () { };
331			matrix (GLfloat diag) { diagonal (diag); };
332			};
333
334			const matrix operator *(const matrix &a, const matrix &b);
335			const vec3 operator *(const matrix &a, const vec3 &v);
336	root	1.40
337			struct vertex_buffer_object {
338			GLuint buffer;
339			GLenum format;
340
341			void alloc ()
342			{
343			if (!buffer)
344			glGenBuffersARB (1, &buffer);
345			}
346
347			void bind (GLenum target = GL_ARRAY_BUFFER_ARB)
348			{
349			glBindBufferARB (target, buffer);
350			glInterleavedArrays (format, 0, 0);
351			}
352
353			void draw (GLenum mode, GLint first, GLsizei count)
354			{
355			bind ();
356			glDrawArrays (mode, first, count);
357			}
358
359			template<class vertex>
360			void set (const vertex *v, GLsizei count, GLenum usage = GL_STATIC_DRAW_ARB, GLenum target = GL_ARRAY_BUFFER_ARB)
361			{
362			alloc ();
363			format = v->gl_format ();
364			glBindBufferARB (target, buffer);
365			glBufferDataARB (target, count * sizeof (vertex), v, usage);
366			}
367
368			template<class vertex>
369			void set (const vector<vertex> &v, GLenum usage = GL_STATIC_DRAW_ARB, GLenum target = GL_ARRAY_BUFFER_ARB)
370			{
371			set (&v[0], v.size (), usage, target);
372			}
373
374			vertex_buffer_object ()
375			: buffer(0)
376			{
377			}
378
379			~vertex_buffer_object ()
380			{
381			glDeleteBuffersARB (1, &buffer);
382			}
383
384			operator GLint ()
385			{
386			return buffer;
387			}
388			};
389
390	root	1.44	void draw_bbox (vertex_buffer_object &vb, const sector &a, const sector &b);
391	root	1.31	}
392
393	root	1.40	GLuint SDL_GL_LoadTexture (SDL_Surface surface, GLfloat texcoord);
394	root	1.27
395	root	1.1	#endif
396