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

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

inline bool overlap (const cone &c, const sphere &s)
{
  vec3 d = s.p - c.p + c.a * (s.r * c.fs1);
  GLfloat dd = dot (c.a, d);
  return dd > 0.F && dd * dd >= length2 (d) * c.fc2;
}

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.47
Committed:	Mon Oct 18 12:24:29 2004 UTC (19 years, 7 months ago) by root
Content type:	text/plain
Branch:	MAIN
Changes since 1.46:	+4 -3 lines
Log Message:	* empty log message *
#	Content
1	#ifndef UTIL_H
2	#define UTIL_H
3
4	#include <cmath>
5	#include <cfloat>
6	#include <cstdlib>
7	#include <vector>
8	#include <string>
9
10	#include "opengl.h"
11
12	#include <SDL/SDL_image.h>
13
14	using namespace std;
15
16	extern CGcontext cgc;
17	extern CGprogram vsh;
18	extern CGprogram fsh;
19	extern CGparameter mv, mvp, lightpos;
20	extern CGprofile vsh_profile, fsh_profile;
21
22	typedef long long soffs; // 32 bit
23	typedef unsigned long long uoffs;
24	#define OFFS_BITS 63
25	#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
29	#define ABS(n) ((n) < 0 ? -(n) : (n))
30
31	struct sector
32	{
33	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	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	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	return sector (ABS (s.x), ABS (s.y), ABS (s.z));
84	}
85
86	struct vec3
87	{
88	GLfloat x, y, z;
89	vec3 () { };
90	vec3 (GLfloat s) : x(s), y(s), z(s) { };
91	vec3 (GLfloat x, GLfloat y, GLfloat z) : x(x), y(y), z(z) { };
92	vec3 (const sector &s) : x(s.x), y(s.y), z(s.z) { };
93
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
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	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	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	// squared length
122	inline const GLfloat length2 (const vec3 &v)
123	{
124	return dot (v, v);
125	}
126
127	inline const GLfloat length (const vec3 &v)
128	{
129	return sqrtf (length2 (v));
130	}
131
132	typedef vec3 point;
133
134	// a generic plane
135	struct plane
136	{
137	vec3 n;
138	GLfloat d;
139
140	plane () { };
141	plane (GLfloat a, GLfloat b, GLfloat c, GLfloat d);
142	};
143
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, fs1, fc2; // angle
168
169	cone () { };
170	cone (const point &p, const vec3 &a, GLfloat f)
171	: p(p), a(a), f(f)
172	{
173	fs1 = 1.F / sinf (f);
174	fc2 = cosf (f); fc2 *= fc2;
175	};
176	};
177
178	inline bool overlap (const cone &c, const sphere &s)
179	{
180	vec3 d = s.p - c.p + c.a * (s.r * c.fs1);
181	GLfloat dd = dot (c.a, d);
182	return dd > 0.F && dd * dd >= length2 (d) * c.fc2;
183	}
184
185	void renormalize (sector &s, point &p);
186
187	struct colour
188	{
189	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	};
198
199	struct tex2
200	{
201	GLfloat s, t;
202	tex2 () { };
203	tex2 (GLfloat s, GLfloat t) : s(s), t(t) { };
204	};
205
206	struct box
207	{
208	point a, b;
209
210	box() { };
211
212	void reset ()
213	{
214	a = point ( FLT_MAX, FLT_MAX, FLT_MAX);
215	b = point (-FLT_MAX, -FLT_MAX, -FLT_MAX);
216	}
217
218	void add (const box &o);
219	void add (const point &p);
220	};
221
222	struct entity;
223	struct geometry;
224	struct view;
225	struct octant;
226
227	extern struct timer
228	{
229	static double now;
230	static double diff;
231	static double fps;
232
233	static void frame ();
234	timer ();
235	} timer;
236
237	/*
238	#define MAX_EVENT_TYPES 10
239	enum event_type { TIMER_EV };
240	struct event
241	{
242	event_type type;
243	};
244
245	typedef callback1<void, event&> event_cb;
246
247	class skedjuhlar
248	{
249	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	*/
282
283	namespace gl
284	{
285	#ifdef DEBUG
286	extern int nesting;
287	void errchk (const char name, const char args, const char *file, int line);
288	#endif
289
290	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
312	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
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	void draw_bbox (vertex_buffer_object &vb, const sector &a, const sector &b);
391	}
392
393	GLuint SDL_GL_LoadTexture (SDL_Surface surface, GLfloat texcoord);
394
395	#endif
396