nuclear@0: #include <stdio.h>
nuclear@0: #include "meshgen.h"
nuclear@0: #include "mesh.h"
nuclear@0: 
nuclear@0: // -------- sphere --------
nuclear@0: 
nuclear@0: #define SURAD(u)	((u) * 2.0 * M_PI)
nuclear@0: #define SVRAD(v)	((v) * M_PI)
nuclear@0: 
nuclear@0: static Vector3 sphvec(float theta, float phi)
nuclear@0: {
nuclear@0: 	return Vector3(sin(theta) * sin(phi),
nuclear@0: 			cos(phi),
nuclear@0: 			cos(theta) * sin(phi));
nuclear@0: }
nuclear@0: 
nuclear@0: void gen_sphere(Mesh *mesh, float rad, int usub, int vsub, float urange, float vrange)
nuclear@0: {
nuclear@0: 	if(usub < 4) usub = 4;
nuclear@0: 	if(vsub < 2) vsub = 2;
nuclear@0: 
nuclear@0: 	int uverts = usub + 1;
nuclear@0: 	int vverts = vsub + 1;
nuclear@0: 
nuclear@0: 	int num_verts = uverts * vverts;
nuclear@0: 	int num_quads = usub * vsub;
nuclear@0: 	int num_tri = num_quads * 2;
nuclear@0: 
nuclear@0: 	mesh->clear();
nuclear@0: 	Vector3 *varr = (Vector3*)mesh->set_attrib_data(MESH_ATTR_VERTEX, 3, num_verts, 0);
nuclear@0: 	Vector3 *narr = (Vector3*)mesh->set_attrib_data(MESH_ATTR_NORMAL, 3, num_verts, 0);
nuclear@0: 	Vector3 *tarr = (Vector3*)mesh->set_attrib_data(MESH_ATTR_TANGENT, 3, num_verts, 0);
nuclear@0: 	Vector2 *uvarr = (Vector2*)mesh->set_attrib_data(MESH_ATTR_TEXCOORD, 2, num_verts, 0);
nuclear@0: 	unsigned int *idxarr = mesh->set_index_data(num_tri * 3, 0);
nuclear@0: 
nuclear@0: 	float du = urange / (float)(uverts - 1);
nuclear@0: 	float dv = vrange / (float)(vverts - 1);
nuclear@0: 
nuclear@0: 	float u = 0.0;
nuclear@0: 	for(int i=0; i<uverts; i++) {
nuclear@0: 		float theta = SURAD(u * urange);
nuclear@0: 
nuclear@0: 		float v = 0.0;
nuclear@0: 		for(int j=0; j<vverts; j++) {
nuclear@0: 			float phi = SVRAD(v * vrange);
nuclear@0: 
nuclear@0: 			Vector3 pos = sphvec(theta, phi);
nuclear@0: 
nuclear@0: 			*varr++ = pos * rad;
nuclear@0: 			*narr++ = pos;
nuclear@0: 			*tarr++ = (sphvec(theta + 0.1f, (float)M_PI / 2.0f) - sphvec(theta - 0.1f, (float)M_PI / 2.0f)).normalized();
nuclear@0: 			*uvarr++ = Vector2(u * urange, v * vrange);
nuclear@0: 
nuclear@0: 			if(i < usub && j < vsub) {
nuclear@0: 				int idx = i * vverts + j;
nuclear@0: 				*idxarr++ = idx;
nuclear@0: 				*idxarr++ = idx + 1;
nuclear@0: 				*idxarr++ = idx + vverts + 1;
nuclear@0: 
nuclear@0: 				*idxarr++ = idx;
nuclear@0: 				*idxarr++ = idx + vverts + 1;
nuclear@0: 				*idxarr++ = idx + vverts;
nuclear@0: 			}
nuclear@0: 
nuclear@0: 			v += dv;
nuclear@0: 		}
nuclear@0: 		u += du;
nuclear@0: 	}
nuclear@0: }
nuclear@0: 
nuclear@0: 
nuclear@0: // -------- cylinder --------
nuclear@0: 
nuclear@0: static Vector3 cylvec(float theta, float height)
nuclear@0: {
nuclear@0: 	return Vector3(sin(theta), height, cos(theta));
nuclear@0: }
nuclear@0: 
nuclear@0: void gen_cylinder(Mesh *mesh, float rad, float height, int usub, int vsub, int capsub, float urange, float vrange)
nuclear@0: {
nuclear@0: 	if(usub < 4) usub = 4;
nuclear@0: 	if(vsub < 1) vsub = 1;
nuclear@0: 
nuclear@0: 	int uverts = usub + 1;
nuclear@0: 	int vverts = vsub + 1;
nuclear@0: 
nuclear@0: 	int num_body_verts = uverts * vverts;
nuclear@0: 	int num_body_quads = usub * vsub;
nuclear@0: 	int num_body_tri = num_body_quads * 2;
nuclear@0: 
nuclear@0: 	int capvverts = capsub ? capsub + 1 : 0;
nuclear@0: 	int num_cap_verts = uverts * capvverts;
nuclear@0: 	int num_cap_quads = usub * capsub;
nuclear@0: 	int num_cap_tri = num_cap_quads * 2;
nuclear@0: 
nuclear@0: 	int num_verts = num_body_verts + num_cap_verts * 2;
nuclear@0: 	int num_tri = num_body_tri + num_cap_tri * 2;
nuclear@0: 
nuclear@0: 	mesh->clear();
nuclear@0: 	Vector3 *varr = (Vector3*)mesh->set_attrib_data(MESH_ATTR_VERTEX, 3, num_verts, 0);
nuclear@0: 	Vector3 *narr = (Vector3*)mesh->set_attrib_data(MESH_ATTR_NORMAL, 3, num_verts, 0);
nuclear@0: 	Vector3 *tarr = (Vector3*)mesh->set_attrib_data(MESH_ATTR_TANGENT, 3, num_verts, 0);
nuclear@0: 	Vector2 *uvarr = (Vector2*)mesh->set_attrib_data(MESH_ATTR_TEXCOORD, 2, num_verts, 0);
nuclear@0: 	unsigned int *idxarr = mesh->set_index_data(num_tri * 3, 0);
nuclear@0: 
nuclear@0: 	float du = urange / (float)(uverts - 1);
nuclear@0: 	float dv = vrange / (float)(vverts - 1);
nuclear@0: 
nuclear@0: 	float u = 0.0;
nuclear@0: 	for(int i=0; i<uverts; i++) {
nuclear@0: 		float theta = SURAD(u);
nuclear@0: 
nuclear@0: 		float v = 0.0;
nuclear@0: 		for(int j=0; j<vverts; j++) {
nuclear@0: 			float y = (v - 0.5) * height;
nuclear@0: 			Vector3 pos = cylvec(theta, y);
nuclear@0: 
nuclear@0: 			*varr++ = Vector3(pos.x * rad, pos.y, pos.z * rad);
nuclear@0: 			*narr++ = Vector3(pos.x, 0.0, pos.z);
nuclear@0: 			*tarr++ = (cylvec(theta + 0.1, 0.0) - cylvec(theta - 0.1, 0.0)).normalized();
nuclear@0: 			*uvarr++ = Vector2(u * urange, v * vrange);
nuclear@0: 
nuclear@0: 			if(i < usub && j < vsub) {
nuclear@0: 				int idx = i * vverts + j;
nuclear@0: 
nuclear@0: 				*idxarr++ = idx;
nuclear@0: 				*idxarr++ = idx + vverts + 1;
nuclear@0: 				*idxarr++ = idx + 1;
nuclear@0: 
nuclear@0: 				*idxarr++ = idx;
nuclear@0: 				*idxarr++ = idx + vverts;
nuclear@0: 				*idxarr++ = idx + vverts + 1;
nuclear@0: 			}
nuclear@0: 
nuclear@0: 			v += dv;
nuclear@0: 		}
nuclear@0: 		u += du;
nuclear@0: 	}
nuclear@0: 
nuclear@0: 
nuclear@0: 	// now the cap!
nuclear@0: 	if(!capsub) {
nuclear@0: 		return;
nuclear@0: 	}
nuclear@0: 
nuclear@0: 	dv = 1.0 / (float)(capvverts - 1);
nuclear@0: 
nuclear@0: 	u = 0.0;
nuclear@0: 	for(int i=0; i<uverts; i++) {
nuclear@0: 		float theta = SURAD(u);
nuclear@0: 
nuclear@0: 		float v = 0.0;
nuclear@0: 		for(int j=0; j<capvverts; j++) {
nuclear@0: 			float r = v * rad;
nuclear@0: 
nuclear@0: 			Vector3 pos = cylvec(theta, height / 2.0) * r;
nuclear@0: 			pos.y = height / 2.0;
nuclear@0: 			Vector3 tang = (cylvec(theta + 0.1, 0.0) - cylvec(theta - 0.1, 0.0)).normalized();
nuclear@0: 
nuclear@0: 			*varr++ = pos;
nuclear@0: 			*narr++ = Vector3(0, 1, 0);
nuclear@0: 			*tarr++ = tang;
nuclear@0: 			*uvarr++ = Vector2(u * urange, v);
nuclear@0: 
nuclear@0: 			pos.y = -height / 2.0;
nuclear@0: 			*varr++ = pos;
nuclear@0: 			*narr++ = Vector3(0, -1, 0);
nuclear@0: 			*tarr++ = -tang;
nuclear@0: 			*uvarr++ = Vector2(u * urange, v);
nuclear@0: 
nuclear@0: 			if(i < usub && j < capsub) {
nuclear@0: 				unsigned int idx = num_body_verts + (i * capvverts + j) * 2;
nuclear@0: 
nuclear@0: 				unsigned int vidx[4] = {
nuclear@0: 					idx,
nuclear@0: 					idx + capvverts * 2,
nuclear@0: 					idx + (capvverts + 1) * 2,
nuclear@0: 					idx + 2
nuclear@0: 				};
nuclear@0: 
nuclear@0: 				*idxarr++ = vidx[0];
nuclear@0: 				*idxarr++ = vidx[2];
nuclear@0: 				*idxarr++ = vidx[1];
nuclear@0: 				*idxarr++ = vidx[0];
nuclear@0: 				*idxarr++ = vidx[3];
nuclear@0: 				*idxarr++ = vidx[2];
nuclear@0: 
nuclear@0: 				*idxarr++ = vidx[0] + 1;
nuclear@0: 				*idxarr++ = vidx[1] + 1;
nuclear@0: 				*idxarr++ = vidx[2] + 1;
nuclear@0: 				*idxarr++ = vidx[0] + 1;
nuclear@0: 				*idxarr++ = vidx[2] + 1;
nuclear@0: 				*idxarr++ = vidx[3] + 1;
nuclear@0: 			}
nuclear@0: 
nuclear@0: 			v += dv;
nuclear@0: 		}
nuclear@0: 		u += du;
nuclear@0: 	}
nuclear@0: }
nuclear@0: 
nuclear@0: // -------- cone --------
nuclear@0: 
nuclear@0: static Vector3 conevec(float theta, float y, float height)
nuclear@0: {
nuclear@0: 	float scale = 1.0 - y / height;
nuclear@0: 	return Vector3(sin(theta) * scale, y, cos(theta) * scale);
nuclear@0: }
nuclear@0: 
nuclear@0: void gen_cone(Mesh *mesh, float rad, float height, int usub, int vsub, int capsub, float urange, float vrange)
nuclear@0: {
nuclear@0: 	if(usub < 4) usub = 4;
nuclear@0: 	if(vsub < 1) vsub = 1;
nuclear@0: 
nuclear@0: 	int uverts = usub + 1;
nuclear@0: 	int vverts = vsub + 1;
nuclear@0: 
nuclear@0: 	int num_body_verts = uverts * vverts;
nuclear@0: 	int num_body_quads = usub * vsub;
nuclear@0: 	int num_body_tri = num_body_quads * 2;
nuclear@0: 
nuclear@0: 	int capvverts = capsub ? capsub + 1 : 0;
nuclear@0: 	int num_cap_verts = uverts * capvverts;
nuclear@0: 	int num_cap_quads = usub * capsub;
nuclear@0: 	int num_cap_tri = num_cap_quads * 2;
nuclear@0: 
nuclear@0: 	int num_verts = num_body_verts + num_cap_verts;
nuclear@0: 	int num_tri = num_body_tri + num_cap_tri;
nuclear@0: 
nuclear@0: 	mesh->clear();
nuclear@0: 	Vector3 *varr = (Vector3*)mesh->set_attrib_data(MESH_ATTR_VERTEX, 3, num_verts, 0);
nuclear@0: 	Vector3 *narr = (Vector3*)mesh->set_attrib_data(MESH_ATTR_NORMAL, 3, num_verts, 0);
nuclear@0: 	Vector3 *tarr = (Vector3*)mesh->set_attrib_data(MESH_ATTR_TANGENT, 3, num_verts, 0);
nuclear@0: 	Vector2 *uvarr = (Vector2*)mesh->set_attrib_data(MESH_ATTR_TEXCOORD, 2, num_verts, 0);
nuclear@0: 	unsigned int *idxarr = mesh->set_index_data(num_tri * 3, 0);
nuclear@0: 
nuclear@0: 	float du = urange / (float)(uverts - 1);
nuclear@0: 	float dv = vrange / (float)(vverts - 1);
nuclear@0: 
nuclear@0: 	float u = 0.0;
nuclear@0: 	for(int i=0; i<uverts; i++) {
nuclear@0: 		float theta = SURAD(u);
nuclear@0: 
nuclear@0: 		float v = 0.0;
nuclear@0: 		for(int j=0; j<vverts; j++) {
nuclear@0: 			float y = v * height;
nuclear@0: 			Vector3 pos = conevec(theta, y, height);
nuclear@0: 
nuclear@0: 			Vector3 tang = (conevec(theta + 0.1, 0.0, height) - conevec(theta - 0.1, 0.0, height)).normalized();
nuclear@0: 			Vector3 bitang = (conevec(theta, y + 0.1, height) - pos).normalized();
nuclear@0: 
nuclear@0: 			*varr++ = Vector3(pos.x * rad, pos.y, pos.z * rad);
nuclear@0: 			*narr++ = cross_product(tang, bitang);
nuclear@0: 			*tarr++ = tang;
nuclear@0: 			*uvarr++ = Vector2(u * urange, v * vrange);
nuclear@0: 
nuclear@0: 			if(i < usub && j < vsub) {
nuclear@0: 				int idx = i * vverts + j;
nuclear@0: 
nuclear@0: 				*idxarr++ = idx;
nuclear@0: 				*idxarr++ = idx + vverts + 1;
nuclear@0: 				*idxarr++ = idx + 1;
nuclear@0: 
nuclear@0: 				*idxarr++ = idx;
nuclear@0: 				*idxarr++ = idx + vverts;
nuclear@0: 				*idxarr++ = idx + vverts + 1;
nuclear@0: 			}
nuclear@0: 
nuclear@0: 			v += dv;
nuclear@0: 		}
nuclear@0: 		u += du;
nuclear@0: 	}
nuclear@0: 
nuclear@0: 
nuclear@0: 	// now the bottom cap!
nuclear@0: 	if(!capsub) {
nuclear@0: 		return;
nuclear@0: 	}
nuclear@0: 
nuclear@0: 	dv = 1.0 / (float)(capvverts - 1);
nuclear@0: 
nuclear@0: 	u = 0.0;
nuclear@0: 	for(int i=0; i<uverts; i++) {
nuclear@0: 		float theta = SURAD(u);
nuclear@0: 
nuclear@0: 		float v = 0.0;
nuclear@0: 		for(int j=0; j<capvverts; j++) {
nuclear@0: 			float r = v * rad;
nuclear@0: 
nuclear@0: 			Vector3 pos = conevec(theta, 0.0, height) * r;
nuclear@0: 			Vector3 tang = (cylvec(theta + 0.1, 0.0) - cylvec(theta - 0.1, 0.0)).normalized();
nuclear@0: 
nuclear@0: 			*varr++ = pos;
nuclear@0: 			*narr++ = Vector3(0, -1, 0);
nuclear@0: 			*tarr++ = tang;
nuclear@0: 			*uvarr++ = Vector2(u * urange, v);
nuclear@0: 
nuclear@0: 			if(i < usub && j < capsub) {
nuclear@0: 				unsigned int idx = num_body_verts + i * capvverts + j;
nuclear@0: 
nuclear@0: 				unsigned int vidx[4] = {
nuclear@0: 					idx,
nuclear@0: 					idx + capvverts,
nuclear@0: 					idx + (capvverts + 1),
nuclear@0: 					idx + 1
nuclear@0: 				};
nuclear@0: 
nuclear@0: 				*idxarr++ = vidx[0];
nuclear@0: 				*idxarr++ = vidx[1];
nuclear@0: 				*idxarr++ = vidx[2];
nuclear@0: 				*idxarr++ = vidx[0];
nuclear@0: 				*idxarr++ = vidx[2];
nuclear@0: 				*idxarr++ = vidx[3];
nuclear@0: 			}
nuclear@0: 
nuclear@0: 			v += dv;
nuclear@0: 		}
nuclear@0: 		u += du;
nuclear@0: 	}
nuclear@0: }
nuclear@0: 
nuclear@0: 
nuclear@0: // -------- plane --------
nuclear@0: 
nuclear@0: void gen_plane(Mesh *mesh, float width, float height, int usub, int vsub)
nuclear@0: {
nuclear@0: 	gen_heightmap(mesh, width, height, usub, vsub, 0);
nuclear@0: }
nuclear@0: 
nuclear@0: 
nuclear@0: // ----- heightmap ------
nuclear@0: 
nuclear@0: void gen_heightmap(Mesh *mesh, float width, float height, int usub, int vsub, float (*hf)(float, float, void*), void *hfdata)
nuclear@0: {
nuclear@0: 	if(usub < 1) usub = 1;
nuclear@0: 	if(vsub < 1) vsub = 1;
nuclear@0: 
nuclear@0: 	mesh->clear();
nuclear@0: 
nuclear@0: 	int uverts = usub + 1;
nuclear@0: 	int vverts = vsub + 1;
nuclear@0: 	int num_verts = uverts * vverts;
nuclear@0: 
nuclear@0: 	int num_quads = usub * vsub;
nuclear@0: 	int num_tri = num_quads * 2;
nuclear@0: 
nuclear@0: 	Vector3 *varr = (Vector3*)mesh->set_attrib_data(MESH_ATTR_VERTEX, 3, num_verts, 0);
nuclear@0: 	Vector3 *narr = (Vector3*)mesh->set_attrib_data(MESH_ATTR_NORMAL, 3, num_verts, 0);
nuclear@0: 	Vector3 *tarr = (Vector3*)mesh->set_attrib_data(MESH_ATTR_TANGENT, 3, num_verts, 0);
nuclear@0: 	Vector2 *uvarr = (Vector2*)mesh->set_attrib_data(MESH_ATTR_TEXCOORD, 2, num_verts, 0);
nuclear@0: 	unsigned int *idxarr = mesh->set_index_data(num_tri * 3, 0);
nuclear@0: 
nuclear@0: 	float du = 1.0 / (float)usub;
nuclear@0: 	float dv = 1.0 / (float)vsub;
nuclear@0: 
nuclear@0: 	float u = 0.0;
nuclear@0: 	for(int i=0; i<uverts; i++) {
nuclear@0: 		float v = 0.0;
nuclear@0: 		for(int j=0; j<vverts; j++) {
nuclear@0: 			float x = (u - 0.5) * width;
nuclear@0: 			float y = (v - 0.5) * height;
nuclear@0: 			float z = hf ? hf(u, v, hfdata) : 0.0;
nuclear@0: 
nuclear@0: 			Vector3 normal = Vector3(0, 0, 1);
nuclear@0: 			if(hf) {
nuclear@0: 				float u1z = hf(u + du, v, hfdata);
nuclear@0: 				float v1z = hf(u, v + dv, hfdata);
nuclear@0: 
nuclear@0: 				Vector3 tang = Vector3(du * width, 0, u1z - z);
nuclear@0: 				Vector3 bitan = Vector3(0, dv * height, v1z - z);
nuclear@0: 				normal = cross_product(tang, bitan).normalized();
nuclear@0: 			}
nuclear@0: 
nuclear@0: 			*varr++ = Vector3(x, y, z);
nuclear@0: 			*narr++ = normal;
nuclear@0: 			*tarr++ = Vector3(1, 0, 0);
nuclear@0: 			*uvarr++ = Vector2(u, v);
nuclear@0: 
nuclear@0: 			if(i < usub && j < vsub) {
nuclear@0: 				int idx = i * vverts + j;
nuclear@0: 
nuclear@0: 				*idxarr++ = idx;
nuclear@0: 				*idxarr++ = idx + vverts + 1;
nuclear@0: 				*idxarr++ = idx + 1;
nuclear@0: 
nuclear@0: 				*idxarr++ = idx;
nuclear@0: 				*idxarr++ = idx + vverts;
nuclear@0: 				*idxarr++ = idx + vverts + 1;
nuclear@0: 			}
nuclear@0: 
nuclear@0: 			v += dv;
nuclear@0: 		}
nuclear@0: 		u += du;
nuclear@0: 	}
nuclear@0: }
nuclear@0: 
nuclear@1: // ----- box ------
nuclear@1: void gen_box(Mesh *mesh, float xsz, float ysz, float zsz, int usub, int vsub)
nuclear@1: {
nuclear@1: 	static const float face_angles[][2] = {
nuclear@1: 		{0, 0},
nuclear@1: 		{M_PI / 2.0, 0},
nuclear@1: 		{M_PI, 0},
nuclear@1: 		{3.0 * M_PI / 2.0, 0},
nuclear@1: 		{0, M_PI / 2.0},
nuclear@1: 		{0, -M_PI / 2.0}
nuclear@1: 	};
nuclear@0: 
nuclear@1: 	if(usub < 1) usub = 1;
nuclear@1: 	if(vsub < 1) vsub = 1;
nuclear@1: 
nuclear@1: 	mesh->clear();
nuclear@1: 
nuclear@1: 	for(int i=0; i<6; i++) {
nuclear@2: 		Matrix4x4 xform, dir_xform;
nuclear@1: 		Mesh m;
nuclear@1: 
nuclear@1: 		gen_plane(&m, 1, 1, usub, vsub);
nuclear@1: 		xform.rotate(Vector3(face_angles[i][1], face_angles[i][0], 0));
nuclear@2: 		dir_xform = xform;
nuclear@1: 		xform.translate(Vector3(0, 0, 0.5));
nuclear@2: 		m.apply_xform(xform, dir_xform);
nuclear@1: 
nuclear@1: 		mesh->append(m);
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	Matrix4x4 scale;
nuclear@1: 	scale.set_scaling(Vector3(xsz, ysz, zsz));
nuclear@1: 	mesh->apply_xform(scale, Matrix4x4::identity);
nuclear@1: }
nuclear@1: 
nuclear@1: /*
nuclear@0: void gen_box(Mesh *mesh, float xsz, float ysz, float zsz)
nuclear@0: {
nuclear@0: 	mesh->clear();
nuclear@0: 
nuclear@0: 	const int num_faces = 6;
nuclear@0: 	int num_verts = num_faces * 4;
nuclear@0: 	int num_tri = num_faces * 2;
nuclear@0: 
nuclear@0: 	float x = xsz / 2.0;
nuclear@0: 	float y = ysz / 2.0;
nuclear@0: 	float z = zsz / 2.0;
nuclear@0: 
nuclear@0: 	Vector3 *varr = (Vector3*)mesh->set_attrib_data(MESH_ATTR_VERTEX, 3, num_verts, 0);
nuclear@0: 	Vector3 *narr = (Vector3*)mesh->set_attrib_data(MESH_ATTR_NORMAL, 3, num_verts, 0);
nuclear@0: 	Vector3 *tarr = (Vector3*)mesh->set_attrib_data(MESH_ATTR_TANGENT, 3, num_verts, 0);
nuclear@0: 	Vector2 *uvarr = (Vector2*)mesh->set_attrib_data(MESH_ATTR_TEXCOORD, 2, num_verts, 0);
nuclear@0: 	unsigned int *idxarr = mesh->set_index_data(num_tri * 3, 0);
nuclear@0: 
nuclear@0: 	static const Vector2 uv[] = { Vector2(0, 0), Vector2(1, 0), Vector2(1, 1), Vector2(0, 1) };
nuclear@0: 
nuclear@0: 	// front
nuclear@0: 	for(int i=0; i<4; i++) {
nuclear@0: 		*narr++ = Vector3(0, 0, 1);
nuclear@0: 		*tarr++ = Vector3(1, 0, 0);
nuclear@0: 		*uvarr++ = uv[i];
nuclear@0: 	}
nuclear@0: 	*varr++ = Vector3(-x, -y, z);
nuclear@0: 	*varr++ = Vector3(x, -y, z);
nuclear@0: 	*varr++ = Vector3(x, y, z);
nuclear@0: 	*varr++ = Vector3(-x, y, z);
nuclear@0: 	// right
nuclear@0: 	for(int i=0; i<4; i++) {
nuclear@0: 		*narr++ = Vector3(1, 0, 0);
nuclear@0: 		*tarr++ = Vector3(0, 0, -1);
nuclear@0: 		*uvarr++ = uv[i];
nuclear@0: 	}
nuclear@0: 	*varr++ = Vector3(x, -y, z);
nuclear@0: 	*varr++ = Vector3(x, -y, -z);
nuclear@0: 	*varr++ = Vector3(x, y, -z);
nuclear@0: 	*varr++ = Vector3(x, y, z);
nuclear@0: 	// back
nuclear@0: 	for(int i=0; i<4; i++) {
nuclear@0: 		*narr++ = Vector3(0, 0, -1);
nuclear@0: 		*tarr++ = Vector3(-1, 0, 0);
nuclear@0: 		*uvarr++ = uv[i];
nuclear@0: 	}
nuclear@0: 	*varr++ = Vector3(x, -y, -z);
nuclear@0: 	*varr++ = Vector3(-x, -y, -z);
nuclear@0: 	*varr++ = Vector3(-x, y, -z);
nuclear@0: 	*varr++ = Vector3(x, y, -z);
nuclear@0: 	// left
nuclear@0: 	for(int i=0; i<4; i++) {
nuclear@0: 		*narr++ = Vector3(-1, 0, 0);
nuclear@0: 		*tarr++ = Vector3(0, 0, 1);
nuclear@0: 		*uvarr++ = uv[i];
nuclear@0: 	}
nuclear@0: 	*varr++ = Vector3(-x, -y, -z);
nuclear@0: 	*varr++ = Vector3(-x, -y, z);
nuclear@0: 	*varr++ = Vector3(-x, y, z);
nuclear@0: 	*varr++ = Vector3(-x, y, -z);
nuclear@0: 	// top
nuclear@0: 	for(int i=0; i<4; i++) {
nuclear@0: 		*narr++ = Vector3(0, 1, 0);
nuclear@0: 		*tarr++ = Vector3(1, 0, 0);
nuclear@0: 		*uvarr++ = uv[i];
nuclear@0: 	}
nuclear@0: 	*varr++ = Vector3(-x, y, z);
nuclear@0: 	*varr++ = Vector3(x, y, z);
nuclear@0: 	*varr++ = Vector3(x, y, -z);
nuclear@0: 	*varr++ = Vector3(-x, y, -z);
nuclear@0: 	// bottom
nuclear@0: 	for(int i=0; i<4; i++) {
nuclear@0: 		*narr++ = Vector3(0, -1, 0);
nuclear@0: 		*tarr++ = Vector3(1, 0, 0);
nuclear@0: 		*uvarr++ = uv[i];
nuclear@0: 	}
nuclear@0: 	*varr++ = Vector3(-x, -y, -z);
nuclear@0: 	*varr++ = Vector3(x, -y, -z);
nuclear@0: 	*varr++ = Vector3(x, -y, z);
nuclear@0: 	*varr++ = Vector3(-x, -y, z);
nuclear@0: 
nuclear@0: 	// index array
nuclear@0: 	static const int faceidx[] = {0, 1, 2, 0, 2, 3};
nuclear@0: 	for(int i=0; i<num_faces; i++) {
nuclear@0: 		for(int j=0; j<6; j++) {
nuclear@0: 			*idxarr++ = faceidx[j] + i * 4;
nuclear@0: 		}
nuclear@0: 	}
nuclear@0: }
nuclear@1: */
nuclear@0: 
nuclear@0: static inline Vector3 rev_vert(float u, float v, Vector2 (*rf)(float, float, void*), void *cls)
nuclear@0: {
nuclear@0: 	Vector2 pos = rf(u, v, cls);
nuclear@0: 
nuclear@0: 	float angle = u * 2.0 * M_PI;
nuclear@0: 	float x = pos.x * cos(angle);
nuclear@0: 	float y = pos.y;
nuclear@0: 	float z = pos.x * sin(angle);
nuclear@0: 
nuclear@0: 	return Vector3(x, y, z);
nuclear@0: }
nuclear@0: 
nuclear@0: // ------ surface of revolution -------
nuclear@0: void gen_revol(Mesh *mesh, int usub, int vsub, Vector2 (*rfunc)(float, float, void*), void *cls)
nuclear@0: {
nuclear@0: 	gen_revol(mesh, usub, vsub, rfunc, 0, cls);
nuclear@0: }
nuclear@0: 
nuclear@0: void gen_revol(Mesh *mesh, int usub, int vsub, Vector2 (*rfunc)(float, float, void*),
nuclear@0: 		Vector2 (*nfunc)(float, float, void*), void *cls)
nuclear@0: {
nuclear@0: 	if(!rfunc) return;
nuclear@0: 	if(usub < 3) usub = 3;
nuclear@0: 	if(vsub < 1) vsub = 1;
nuclear@0: 
nuclear@0: 	mesh->clear();
nuclear@0: 
nuclear@0: 	int uverts = usub + 1;
nuclear@0: 	int vverts = vsub + 1;
nuclear@0: 	int num_verts = uverts * vverts;
nuclear@0: 
nuclear@0: 	int num_quads = usub * vsub;
nuclear@0: 	int num_tri = num_quads * 2;
nuclear@0: 
nuclear@0: 	Vector3 *varr = (Vector3*)mesh->set_attrib_data(MESH_ATTR_VERTEX, 3, num_verts, 0);
nuclear@0: 	Vector3 *narr = (Vector3*)mesh->set_attrib_data(MESH_ATTR_NORMAL, 3, num_verts, 0);
nuclear@0: 	Vector3 *tarr = (Vector3*)mesh->set_attrib_data(MESH_ATTR_TANGENT, 3, num_verts, 0);
nuclear@0: 	Vector2 *uvarr = (Vector2*)mesh->set_attrib_data(MESH_ATTR_TEXCOORD, 2, num_verts, 0);
nuclear@0: 	unsigned int *idxarr = mesh->set_index_data(num_tri * 3, 0);
nuclear@0: 
nuclear@0: 	float du = 1.0 / (float)(uverts - 1);
nuclear@0: 	float dv = 1.0 / (float)(vverts - 1);
nuclear@0: 
nuclear@0: 	float u = 0.0;
nuclear@0: 	for(int i=0; i<uverts; i++) {
nuclear@0: 		float v = 0.0;
nuclear@0: 		for(int j=0; j<vverts; j++) {
nuclear@0: 			Vector3 pos = rev_vert(u, v, rfunc, cls);
nuclear@0: 
nuclear@0: 			Vector3 nextu = rev_vert(fmod(u + du, 1.0), v, rfunc, cls);
nuclear@0: 			Vector3 tang = nextu - pos;
nuclear@0: 			if(tang.length_sq() < 1e-6) {
nuclear@0: 				float new_v = v > 0.5 ? v - dv * 0.25 : v + dv * 0.25;
nuclear@0: 				nextu = rev_vert(fmod(u + du, 1.0), new_v, rfunc, cls);
nuclear@0: 				tang = nextu - pos;
nuclear@0: 			}
nuclear@0: 
nuclear@0: 			Vector3 normal;
nuclear@0: 			if(nfunc) {
nuclear@0: 				normal = rev_vert(u, v, nfunc, cls);
nuclear@0: 			} else {
nuclear@0: 				Vector3 nextv = rev_vert(u, v + dv, rfunc, cls);
nuclear@0: 				Vector3 bitan = nextv - pos;
nuclear@0: 				if(bitan.length_sq() < 1e-6) {
nuclear@0: 					nextv = rev_vert(u, v - dv, rfunc, cls);
nuclear@0: 					bitan = pos - nextv;
nuclear@0: 				}
nuclear@0: 
nuclear@0: 				normal = cross_product(tang, bitan);
nuclear@0: 			}
nuclear@0: 
nuclear@0: 			*varr++ = pos;
nuclear@0: 			*narr++ = normal.normalized();
nuclear@0: 			*tarr++ = tang.normalized();
nuclear@0: 			*uvarr++ = Vector2(u, v);
nuclear@0: 
nuclear@0: 			if(i < usub && j < vsub) {
nuclear@0: 				int idx = i * vverts + j;
nuclear@0: 
nuclear@0: 				*idxarr++ = idx;
nuclear@0: 				*idxarr++ = idx + vverts + 1;
nuclear@0: 				*idxarr++ = idx + 1;
nuclear@0: 
nuclear@0: 				*idxarr++ = idx;
nuclear@0: 				*idxarr++ = idx + vverts;
nuclear@0: 				*idxarr++ = idx + vverts + 1;
nuclear@0: 			}
nuclear@0: 
nuclear@0: 			v += dv;
nuclear@0: 		}
nuclear@0: 		u += du;
nuclear@0: 	}
nuclear@0: }