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, fs1, fc; // angle

  cone () { };
  cone (const point &p, const vec3 &a, GLfloat f)
  : p(p), a(a), f(f)
  {
    fs1 = 1.F / sinf (f);
    fc = cosf (f);
  };
};

inline bool overlap (const cone &c, const sphere &s)
{
  vec3 d = s.p - c.p + c.a * (s.r * c.fs1);
  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.46
Committed:	Mon Oct 18 12:06:10 2004 UTC (19 years, 7 months ago) by root
Content type:	text/plain
Branch:	MAIN
Changes since 1.45:	+3 -4 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	root	1.46	GLfloat f, fs1, fc; // angle
168	root	1.45
169			cone () { };
170			cone (const point &p, const vec3 &a, GLfloat f)
171			: p(p), a(a), f(f)
172			{
173	root	1.46	fs1 = 1.F / sinf (f);
174	root	1.45	fc = cosf (f);
175			};
176			};
177
178			inline bool overlap (const cone &c, const sphere &s)
179			{
180	root	1.46	vec3 d = s.p - c.p + c.a * (s.r * c.fs1);
181	root	1.45	return dot (c.a, d) >= length (d) * c.fc;
182			}
183	root	1.10
184			void renormalize (sector &s, point &p);
185
186	root	1.37	struct colour
187			{
188	root	1.40	GLubyte r, g, b, a;
189			colour (GLfloat r = 1., GLfloat g = 1., GLfloat b = 1., GLfloat a = 1.)
190			: r(GLubyte (r * 255.F + .5F))
191			, g(GLubyte (g * 255.F + .5F))
192			, b(GLubyte (b * 255.F + .5F))
193			, a(GLubyte (a * 255.F + .5F))
194			{
195			}
196	root	1.10	};
197	root	1.5
198	root	1.40	struct tex2
199	root	1.37	{
200	root	1.1	GLfloat s, t;
201	root	1.40	tex2 () { };
202			tex2 (GLfloat s, GLfloat t) : s(s), t(t) { };
203	root	1.1	};
204
205	root	1.37	struct box
206			{
207	root	1.27	point a, b;
208	root	1.8
209			box() { };
210	root	1.1
211			void reset ()
212			{
213	root	1.27	a = point ( FLT_MAX, FLT_MAX, FLT_MAX);
214			b = point (-FLT_MAX, -FLT_MAX, -FLT_MAX);
215	root	1.1	}
216
217			void add (const box &o);
218			void add (const point &p);
219			};
220	root	1.5
221	root	1.27	struct entity;
222			struct geometry;
223	root	1.20	struct view;
224	root	1.32	struct octant;
225	root	1.7
226	root	1.37	extern struct timer
227			{
228	root	1.7	static double now;
229			static double diff;
230	root	1.44	static double fps;
231	root	1.7
232			static void frame ();
233			timer ();
234			} timer;
235	root	1.11
236	root	1.13	/*
237	root	1.11	#define MAX_EVENT_TYPES 10
238			enum event_type { TIMER_EV };
239	root	1.37	struct event
240			{
241	root	1.11	event_type type;
242			};
243
244			typedef callback1<void, event&> event_cb;
245
246	root	1.37	class skedjuhlar
247			{
248	root	1.11	public:
249			// only 10 types for now
250			private:
251			vector <list<event_cb> > event_lists;
252
253			public:
254			skedjuhlar () {
255			event_lists.resize (MAX_EVENT_TYPES, list<event_cb>());
256			}
257
258			void register_event_cb (const event_type &t, const event_cb &e) {
259			event_lists[t].push_back (e);
260			};
261			void send_event (event &e) {
262			list<event_cb> &l = event_lists[e.type];
263			for (list<event_cb>::iterator it = l.begin (); it != l.end (); it++) {
264			(*it)(e);
265			}
266			};
267			void check_events () {
268			while (!events.empty ()) {
269			event &e = events.pop_front ();
270			list<event_cb> &l = event_lists[e->name];
271			for (list<event_cb>::iterator it = l.begin (); it !+ l.end (); it++) {
272			(*it)(e);
273			}
274			delete e; // ugly slow? hai hai..... 183G
275			}
276			}
277			};
278
279			extern skedjuhlar main_scheduler;
280	root	1.13	*/
281	root	1.22
282	root	1.37	namespace gl
283			{
284	root	1.41	#ifdef DEBUG
285			extern int nesting;
286			void errchk (const char name, const char args, const char *file, int line);
287			#endif
288
289	root	1.40	struct vertex_v3f
290			{
291			point p; // vertex
292
293			vertex_v3f () { };
294			vertex_v3f (point p) : p(p) { };
295
296			GLenum gl_format () const { return GL_V3F; }
297			};
298
299			struct vertex_t2f_n3f_v3f
300			{
301			tex2 t; // texture
302			vec3 n; // normal
303			point p; // vertex
304
305			vertex_t2f_n3f_v3f () { };
306			vertex_t2f_n3f_v3f (point p, vec3 n, tex2 t = tex2()) : p(p), n(n), t(t) { };
307
308			GLenum gl_format () const { return GL_T2F_N3F_V3F; }
309			};
310	root	1.34
311	root	1.41	struct matrix
312			{
313			GLfloat data[4][4];
314
315			const GLfloat operator ()(int i, int j) const { return data[j][i]; };
316			GLfloat &operator ()(int i, int j) { return data[j][i]; };
317
318			operator GLfloat *() { return &data[0][0]; }
319
320			void diagonal (GLfloat v);
321			void clear () { diagonal (0.F); };
322			void identity () { diagonal (1.F); };
323
324			void print (); // ugly
325
326			static const matrix translation (const vec3 &v);
327			static const matrix rotation (GLfloat degrees, const vec3 &axis);
328
329			matrix () { };
330			matrix (GLfloat diag) { diagonal (diag); };
331			};
332
333			const matrix operator *(const matrix &a, const matrix &b);
334			const vec3 operator *(const matrix &a, const vec3 &v);
335	root	1.40
336			struct vertex_buffer_object {
337			GLuint buffer;
338			GLenum format;
339
340			void alloc ()
341			{
342			if (!buffer)
343			glGenBuffersARB (1, &buffer);
344			}
345
346			void bind (GLenum target = GL_ARRAY_BUFFER_ARB)
347			{
348			glBindBufferARB (target, buffer);
349			glInterleavedArrays (format, 0, 0);
350			}
351
352			void draw (GLenum mode, GLint first, GLsizei count)
353			{
354			bind ();
355			glDrawArrays (mode, first, count);
356			}
357
358			template<class vertex>
359			void set (const vertex *v, GLsizei count, GLenum usage = GL_STATIC_DRAW_ARB, GLenum target = GL_ARRAY_BUFFER_ARB)
360			{
361			alloc ();
362			format = v->gl_format ();
363			glBindBufferARB (target, buffer);
364			glBufferDataARB (target, count * sizeof (vertex), v, usage);
365			}
366
367			template<class vertex>
368			void set (const vector<vertex> &v, GLenum usage = GL_STATIC_DRAW_ARB, GLenum target = GL_ARRAY_BUFFER_ARB)
369			{
370			set (&v[0], v.size (), usage, target);
371			}
372
373			vertex_buffer_object ()
374			: buffer(0)
375			{
376			}
377
378			~vertex_buffer_object ()
379			{
380			glDeleteBuffersARB (1, &buffer);
381			}
382
383			operator GLint ()
384			{
385			return buffer;
386			}
387			};
388
389	root	1.44	void draw_bbox (vertex_buffer_object &vb, const sector &a, const sector &b);
390	root	1.31	}
391
392	root	1.40	GLuint SDL_GL_LoadTexture (SDL_Surface surface, GLfloat texcoord);
393	root	1.27
394	root	1.1	#endif
395