nuclear@1: #include <stdio.h>
nuclear@1: #include <stdlib.h>
nuclear@1: #include <float.h>
nuclear@1: #include <assert.h>
nuclear@1: #include "opengl.h"
nuclear@1: #include "mesh.h"
nuclear@1: #include "xform_node.h"
nuclear@1: //#include "logger.h"
nuclear@1: 
nuclear@1: int Mesh::global_sdr_loc[NUM_MESH_ATTR] = { 0, 1, 2, 3, 4, 5 };
nuclear@1: unsigned int Mesh::intersect_mode = ISECT_DEFAULT;
nuclear@1: float Mesh::vertex_sel_dist = 0.01;
nuclear@1: float Mesh::vis_vecsize = 1.0;
nuclear@1: 
nuclear@1: Mesh::Mesh()
nuclear@1: {
nuclear@1: 	clear();
nuclear@1: 
nuclear@1: 	glGenBuffers(NUM_MESH_ATTR + 1, buffer_objects);
nuclear@1: 
nuclear@1: 	for(int i=0; i<NUM_MESH_ATTR; i++) {
nuclear@1: 		vattr[i].vbo = buffer_objects[i];
nuclear@1: 	}
nuclear@1: 	ibo = buffer_objects[NUM_MESH_ATTR];
nuclear@1: 	wire_ibo = 0;
nuclear@1: }
nuclear@1: 
nuclear@1: Mesh::~Mesh()
nuclear@1: {
nuclear@1: 	glDeleteBuffers(NUM_MESH_ATTR + 1, buffer_objects);
nuclear@1: 
nuclear@1: 	if(wire_ibo) {
nuclear@1: 		glDeleteBuffers(1, &wire_ibo);
nuclear@1: 	}
nuclear@1: }
nuclear@1: 
nuclear@1: void Mesh::set_name(const char *name)
nuclear@1: {
nuclear@1: 	this->name = name;
nuclear@1: }
nuclear@1: 
nuclear@1: const char *Mesh::get_name() const
nuclear@1: {
nuclear@1: 	return name.c_str();
nuclear@1: }
nuclear@1: 
nuclear@1: bool Mesh::has_attrib(int attr) const
nuclear@1: {
nuclear@1: 	if(attr < 0 || attr >= NUM_MESH_ATTR) {
nuclear@1: 		return false;
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	// if neither of these is valid, then nobody has set this attribute
nuclear@1: 	return vattr[attr].vbo_valid || vattr[attr].data_valid;
nuclear@1: }
nuclear@1: 
nuclear@1: void Mesh::clear()
nuclear@1: {
nuclear@1: 	bones.clear();
nuclear@1: 
nuclear@1: 	for(int i=0; i<NUM_MESH_ATTR; i++) {
nuclear@1: 		vattr[i].nelem = 0;
nuclear@1: 		vattr[i].vbo_valid = false;
nuclear@1: 		vattr[i].data_valid = false;
nuclear@1: 		//vattr[i].sdr_loc = -1;
nuclear@1: 		vattr[i].data.clear();
nuclear@1: 	}
nuclear@1: 	ibo_valid = false;
nuclear@1: 	idata.clear();
nuclear@1: 
nuclear@1: 	wire_ibo_valid = false;
nuclear@1: 
nuclear@1: 	nverts = nfaces = 0;
nuclear@1: 
nuclear@1: 	/*bsph_valid = false;
nuclear@1: 	aabb_valid = false;*/
nuclear@1: }
nuclear@1: 
nuclear@1: float *Mesh::set_attrib_data(int attrib, int nelem, unsigned int num, const float *data)
nuclear@1: {
nuclear@1: 	if(attrib < 0 || attrib >= NUM_MESH_ATTR) {
nuclear@1: 		fprintf(stderr, "%s: invalid attrib: %d\n", __FUNCTION__, attrib);
nuclear@1: 		return 0;
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	if(nverts && num != nverts) {
nuclear@1: 		fprintf(stderr, "%s: attribute count missmatch (%d instead of %d)\n", __FUNCTION__, num, nverts);
nuclear@1: 		return 0;
nuclear@1: 	}
nuclear@1: 	nverts = num;
nuclear@1: 
nuclear@1: 	vattr[attrib].data.clear();
nuclear@1: 	vattr[attrib].nelem = nelem;
nuclear@1: 	vattr[attrib].data.resize(num * nelem);
nuclear@1: 
nuclear@1: 	if(data) {
nuclear@1: 		memcpy(&vattr[attrib].data[0], data, num * nelem * sizeof *data);
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	vattr[attrib].data_valid = true;
nuclear@1: 	vattr[attrib].vbo_valid = false;
nuclear@1: 	return &vattr[attrib].data[0];
nuclear@1: }
nuclear@1: 
nuclear@1: float *Mesh::get_attrib_data(int attrib)
nuclear@1: {
nuclear@1: 	if(attrib < 0 || attrib >= NUM_MESH_ATTR) {
nuclear@1: 		fprintf(stderr, "%s: invalid attrib: %d\n", __FUNCTION__, attrib);
nuclear@1: 		return 0;
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	vattr[attrib].vbo_valid = false;
nuclear@1: 	return (float*)((const Mesh*)this)->get_attrib_data(attrib);
nuclear@1: }
nuclear@1: 
nuclear@1: const float *Mesh::get_attrib_data(int attrib) const
nuclear@1: {
nuclear@1: 	if(attrib < 0 || attrib >= NUM_MESH_ATTR) {
nuclear@1: 		fprintf(stderr, "%s: invalid attrib: %d\n", __FUNCTION__, attrib);
nuclear@1: 		return 0;
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	if(!vattr[attrib].data_valid) {
nuclear@1: #if GL_ES_VERSION_2_0
nuclear@1: 		fprintf(stderr, "%s: can't read back attrib data on CrippledGL ES\n", __FUNCTION__);
nuclear@1: 		return 0;
nuclear@1: #else
nuclear@1: 		if(!vattr[attrib].vbo_valid) {
nuclear@1: 			fprintf(stderr, "%s: unavailable attrib: %d\n", __FUNCTION__, attrib);
nuclear@1: 			return 0;
nuclear@1: 		}
nuclear@1: 
nuclear@1: 		// local data copy is unavailable, grab the data from the vbo
nuclear@1: 		Mesh *m = (Mesh*)this;
nuclear@1: 		m->vattr[attrib].data.resize(nverts * vattr[attrib].nelem);
nuclear@1: 
nuclear@1: 		glBindBuffer(GL_ARRAY_BUFFER, vattr[attrib].vbo);
nuclear@1: 		void *data = glMapBuffer(GL_ARRAY_BUFFER, GL_READ_ONLY);
nuclear@1: 		memcpy(&m->vattr[attrib].data[0], data, nverts * vattr[attrib].nelem * sizeof(float));
nuclear@1: 		glUnmapBuffer(GL_ARRAY_BUFFER);
nuclear@1: 
nuclear@1: 		vattr[attrib].data_valid = true;
nuclear@1: #endif
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	return &vattr[attrib].data[0];
nuclear@1: }
nuclear@1: 
nuclear@1: void Mesh::set_attrib(int attrib, int idx, const Vector4 &v)
nuclear@1: {
nuclear@1: 	float *data = get_attrib_data(attrib);
nuclear@1: 	if(data) {
nuclear@1: 		data += idx * vattr[attrib].nelem;
nuclear@1: 		for(int i=0; i<vattr[attrib].nelem; i++) {
nuclear@1: 			data[i] = v[i];
nuclear@1: 		}
nuclear@1: 	}
nuclear@1: }
nuclear@1: 
nuclear@1: Vector4 Mesh::get_attrib(int attrib, int idx) const
nuclear@1: {
nuclear@1: 	Vector4 v(0.0, 0.0, 0.0, 1.0);
nuclear@1: 	const float *data = get_attrib_data(attrib);
nuclear@1: 	if(data) {
nuclear@1: 		data += idx * vattr[attrib].nelem;
nuclear@1: 		for(int i=0; i<vattr[attrib].nelem; i++) {
nuclear@1: 			v[i] = data[i];
nuclear@1: 		}
nuclear@1: 	}
nuclear@1: 	return v;
nuclear@1: }
nuclear@1: 
nuclear@1: unsigned int *Mesh::set_index_data(int num, const unsigned int *indices)
nuclear@1: {
nuclear@1: 	int nidx = nfaces * 3;
nuclear@1: 	if(nidx && num != nidx) {
nuclear@1: 		fprintf(stderr, "%s: index count missmatch (%d instead of %d)\n", __FUNCTION__, num, nidx);
nuclear@1: 		return 0;
nuclear@1: 	}
nuclear@1: 	nfaces = num / 3;
nuclear@1: 
nuclear@1: 	idata.clear();
nuclear@1: 	idata.resize(num);
nuclear@1: 
nuclear@1: 	if(indices) {
nuclear@1: 		memcpy(&idata[0], indices, num * sizeof *indices);
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	idata_valid = true;
nuclear@1: 	ibo_valid = false;
nuclear@1: 
nuclear@1: 	return &idata[0];
nuclear@1: }
nuclear@1: 
nuclear@1: unsigned int *Mesh::get_index_data()
nuclear@1: {
nuclear@1: 	ibo_valid = false;
nuclear@1: 	return (unsigned int*)((const Mesh*)this)->get_index_data();
nuclear@1: }
nuclear@1: 
nuclear@1: const unsigned int *Mesh::get_index_data() const
nuclear@1: {
nuclear@1: 	if(!idata_valid) {
nuclear@1: #if GL_ES_VERSION_2_0
nuclear@1: 		fprintf(stderr, "%s: can't read back index data in CrippledGL ES\n", __FUNCTION__);
nuclear@1: 		return 0;
nuclear@1: #else
nuclear@1: 		if(!ibo_valid) {
nuclear@1: 			fprintf(stderr, "%s: indices unavailable\n", __FUNCTION__);
nuclear@1: 			return 0;
nuclear@1: 		}
nuclear@1: 
nuclear@1: 		// local data copy is unavailable, gram the data from the ibo
nuclear@1: 		Mesh *m = (Mesh*)this;
nuclear@1: 		int nidx = nfaces * 3;
nuclear@1: 		m->idata.resize(nidx);
nuclear@1: 
nuclear@1: 		glBindBuffer(GL_ELEMENT_ARRAY_BUFFER, ibo);
nuclear@1: 		void *data = glMapBuffer(GL_ELEMENT_ARRAY_BUFFER, GL_READ_ONLY);
nuclear@1: 		memcpy(&m->idata[0], data, nidx * sizeof(unsigned int));
nuclear@1: 		glUnmapBuffer(GL_ELEMENT_ARRAY_BUFFER);
nuclear@1: 
nuclear@1: 		idata_valid = true;
nuclear@1: #endif
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	return &idata[0];
nuclear@1: }
nuclear@1: 
nuclear@1: void Mesh::append(const Mesh &mesh)
nuclear@1: {
nuclear@1: 	unsigned int idxoffs = nverts;
nuclear@1: 
nuclear@1: 	nverts += mesh.nverts;
nuclear@1: 	nfaces += mesh.nfaces;
nuclear@1: 
nuclear@1: 	for(int i=0; i<NUM_MESH_ATTR; i++) {
nuclear@1: 		if(has_attrib(i) && mesh.has_attrib(i)) {
nuclear@1: 			// force validating the data arrays
nuclear@1: 			get_attrib_data(i);
nuclear@1: 			mesh.get_attrib_data(i);
nuclear@1: 
nuclear@1: 			// append the mesh data
nuclear@1: 			vattr[i].data.insert(vattr[i].data.end(), mesh.vattr[i].data.begin(), mesh.vattr[i].data.end());
nuclear@1: 		}
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	if(ibo_valid || idata_valid) {
nuclear@1: 		// make index arrays valid
nuclear@1: 		get_index_data();
nuclear@1: 		mesh.get_index_data();
nuclear@1: 
nuclear@1: 		size_t orig_sz = idata.size();
nuclear@1: 
nuclear@1: 		idata.insert(idata.end(), mesh.idata.begin(), mesh.idata.end());
nuclear@1: 
nuclear@1: 		// fixup all the new indices
nuclear@1: 		for(size_t i=orig_sz; i<idata.size(); i++) {
nuclear@1: 			idata[i] += idxoffs;
nuclear@1: 		}
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	// fuck everything
nuclear@1: 	wire_ibo_valid = false;
nuclear@1: 	/*aabb_valid = false;
nuclear@1: 	bsph_valid = false;*/
nuclear@1: }
nuclear@1: 
nuclear@1: // assemble a complete vertex by adding all the useful attributes
nuclear@1: void Mesh::vertex(float x, float y, float z)
nuclear@1: {
nuclear@1: 	cur_val[MESH_ATTR_VERTEX] = Vector4(x, y, z, 1.0f);
nuclear@1: 	vattr[MESH_ATTR_VERTEX].data_valid = true;
nuclear@1: 	vattr[MESH_ATTR_VERTEX].nelem = 3;
nuclear@1: 
nuclear@1: 	for(int i=0; i<NUM_MESH_ATTR; i++) {
nuclear@1: 		if(vattr[i].data_valid) {
nuclear@1: 			for(int j=0; j<vattr[MESH_ATTR_VERTEX].nelem; j++) {
nuclear@1: 				vattr[i].data.push_back(cur_val[i][j]);
nuclear@1: 			}
nuclear@1: 		}
nuclear@1: 		vattr[i].vbo_valid = false;
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	if(idata_valid) {
nuclear@1: 		idata.clear();
nuclear@1: 	}
nuclear@1: 	ibo_valid = idata_valid = false;
nuclear@1: }
nuclear@1: 
nuclear@1: void Mesh::normal(float nx, float ny, float nz)
nuclear@1: {
nuclear@1: 	cur_val[MESH_ATTR_NORMAL] = Vector4(nx, ny, nz, 1.0f);
nuclear@1: 	vattr[MESH_ATTR_NORMAL].data_valid = true;
nuclear@1: 	vattr[MESH_ATTR_NORMAL].nelem = 3;
nuclear@1: }
nuclear@1: 
nuclear@1: void Mesh::tangent(float tx, float ty, float tz)
nuclear@1: {
nuclear@1: 	cur_val[MESH_ATTR_TANGENT] = Vector4(tx, ty, tz, 1.0f);
nuclear@1: 	vattr[MESH_ATTR_TANGENT].data_valid = true;
nuclear@1: 	vattr[MESH_ATTR_TANGENT].nelem = 3;
nuclear@1: }
nuclear@1: 
nuclear@1: void Mesh::texcoord(float u, float v, float w)
nuclear@1: {
nuclear@1: 	cur_val[MESH_ATTR_TEXCOORD] = Vector4(u, v, w, 1.0f);
nuclear@1: 	vattr[MESH_ATTR_TEXCOORD].data_valid = true;
nuclear@1: 	vattr[MESH_ATTR_TEXCOORD].nelem = 3;
nuclear@1: }
nuclear@1: 
nuclear@1: void Mesh::boneweights(float w1, float w2, float w3, float w4)
nuclear@1: {
nuclear@1: 	cur_val[MESH_ATTR_BONEWEIGHTS] = Vector4(w1, w2, w3, w4);
nuclear@1: 	vattr[MESH_ATTR_BONEWEIGHTS].data_valid = true;
nuclear@1: 	vattr[MESH_ATTR_BONEWEIGHTS].nelem = 4;
nuclear@1: }
nuclear@1: 
nuclear@1: void Mesh::boneidx(int idx1, int idx2, int idx3, int idx4)
nuclear@1: {
nuclear@1: 	cur_val[MESH_ATTR_BONEIDX] = Vector4(idx1, idx2, idx3, idx4);
nuclear@1: 	vattr[MESH_ATTR_BONEIDX].data_valid = true;
nuclear@1: 	vattr[MESH_ATTR_BONEIDX].nelem = 4;
nuclear@1: }
nuclear@1: 
nuclear@1: /// static function
nuclear@1: void Mesh::set_attrib_location(int attr, int loc)
nuclear@1: {
nuclear@1: 	if(attr < 0 || attr >= NUM_MESH_ATTR) {
nuclear@1: 		return;
nuclear@1: 	}
nuclear@1: 	Mesh::global_sdr_loc[attr] = loc;
nuclear@1: }
nuclear@1: 
nuclear@1: /// static function
nuclear@1: int Mesh::get_attrib_location(int attr)
nuclear@1: {
nuclear@1: 	if(attr < 0 || attr >= NUM_MESH_ATTR) {
nuclear@1: 		return -1;
nuclear@1: 	}
nuclear@1: 	return Mesh::global_sdr_loc[attr];
nuclear@1: }
nuclear@1: 
nuclear@1: /// static function
nuclear@1: void Mesh::clear_attrib_locations()
nuclear@1: {
nuclear@1: 	for(int i=0; i<NUM_MESH_ATTR; i++) {
nuclear@1: 		Mesh::global_sdr_loc[i] = -1;
nuclear@1: 	}
nuclear@1: }
nuclear@1: 
nuclear@1: /// static function
nuclear@1: void Mesh::set_vis_vecsize(float sz)
nuclear@1: {
nuclear@1: 	Mesh::vis_vecsize = sz;
nuclear@1: }
nuclear@1: 
nuclear@1: float Mesh::get_vis_vecsize()
nuclear@1: {
nuclear@1: 	return Mesh::vis_vecsize;
nuclear@1: }
nuclear@1: 
nuclear@1: void Mesh::apply_xform(const Matrix4x4 &xform)
nuclear@1: {
nuclear@1: 	Matrix4x4 dir_xform = xform;
nuclear@1: 	dir_xform[0][3] = dir_xform[1][3] = dir_xform[2][3] = 0.0f;
nuclear@1: 	dir_xform[3][0] = dir_xform[3][1] = dir_xform[3][2] = 0.0f;
nuclear@1: 	dir_xform[3][3] = 1.0f;
nuclear@1: 
nuclear@1: 	apply_xform(xform, dir_xform);
nuclear@1: }
nuclear@1: 
nuclear@1: void Mesh::apply_xform(const Matrix4x4 &xform, const Matrix4x4 &dir_xform)
nuclear@1: {
nuclear@1: 	for(unsigned int i=0; i<nverts; i++) {
nuclear@1: 		Vector4 v = get_attrib(MESH_ATTR_VERTEX, i);
nuclear@1: 		set_attrib(MESH_ATTR_VERTEX, i, v.transformed(xform));
nuclear@1: 
nuclear@1: 		if(has_attrib(MESH_ATTR_NORMAL)) {
nuclear@1: 			Vector3 n = get_attrib(MESH_ATTR_NORMAL, i);
nuclear@1: 			set_attrib(MESH_ATTR_NORMAL, i, n.transformed(dir_xform));
nuclear@1: 		}
nuclear@1: 		if(has_attrib(MESH_ATTR_TANGENT)) {
nuclear@1: 			Vector3 t = get_attrib(MESH_ATTR_TANGENT, i);
nuclear@1: 			set_attrib(MESH_ATTR_TANGENT, i, t.transformed(dir_xform));
nuclear@1: 		}
nuclear@1: 	}
nuclear@1: }
nuclear@1: 
nuclear@1: int Mesh::add_bone(XFormNode *bone)
nuclear@1: {
nuclear@1: 	int idx = bones.size();
nuclear@1: 	bones.push_back(bone);
nuclear@1: 	return idx;
nuclear@1: }
nuclear@1: 
nuclear@1: const XFormNode *Mesh::get_bone(int idx) const
nuclear@1: {
nuclear@1: 	if(idx < 0 || idx >= (int)bones.size()) {
nuclear@1: 		return 0;
nuclear@1: 	}
nuclear@1: 	return bones[idx];
nuclear@1: }
nuclear@1: 
nuclear@1: int Mesh::get_bones_count() const
nuclear@1: {
nuclear@1: 	return (int)bones.size();
nuclear@1: }
nuclear@1: 
nuclear@1: void Mesh::draw() const
nuclear@1: {
nuclear@1: 	((Mesh*)this)->update_buffers();
nuclear@1: 
nuclear@1: 	if(!vattr[MESH_ATTR_VERTEX].vbo_valid) {
nuclear@1: 		fprintf(stderr, "%s: invalid vertex buffer\n", __FUNCTION__);
nuclear@1: 		return;
nuclear@1: 	}
nuclear@1: 	if(global_sdr_loc[MESH_ATTR_VERTEX] == -1) {
nuclear@1: 		fprintf(stderr, "%s: shader attribute location for vertices unset\n", __FUNCTION__);
nuclear@1: 		return;
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	for(int i=0; i<NUM_MESH_ATTR; i++) {
nuclear@1: 		int loc = global_sdr_loc[i];
nuclear@1: 		if(loc >= 0 && vattr[i].vbo_valid) {
nuclear@1: 			glBindBuffer(GL_ARRAY_BUFFER, vattr[i].vbo);
nuclear@1: 			glVertexAttribPointer(loc, vattr[i].nelem, GL_FLOAT, GL_FALSE, 0, 0);
nuclear@1: 			glEnableVertexAttribArray(loc);
nuclear@1: 		}
nuclear@1: 	}
nuclear@1: 	glBindBuffer(GL_ARRAY_BUFFER, 0);
nuclear@1: 
nuclear@1: 	if(ibo_valid) {
nuclear@1: 		glBindBuffer(GL_ELEMENT_ARRAY_BUFFER, ibo);
nuclear@1: 		glDrawElements(GL_TRIANGLES, nfaces * 3, GL_UNSIGNED_INT, 0);
nuclear@1: 		glBindBuffer(GL_ELEMENT_ARRAY_BUFFER, 0);
nuclear@1: 	} else {
nuclear@1: 		glDrawArrays(GL_TRIANGLES, 0, nverts);
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	for(int i=0; i<NUM_MESH_ATTR; i++) {
nuclear@1: 		int loc = global_sdr_loc[i];
nuclear@1: 		if(loc >= 0 && vattr[i].vbo_valid) {
nuclear@1: 			glDisableVertexAttribArray(loc);
nuclear@1: 		}
nuclear@1: 	}
nuclear@1: }
nuclear@1: 
nuclear@1: void Mesh::draw_wire() const
nuclear@1: {
nuclear@1: 	((Mesh*)this)->update_wire_ibo();
nuclear@1: 
nuclear@1: 	if(!vattr[MESH_ATTR_VERTEX].vbo_valid || !wire_ibo_valid) {
nuclear@1: 		fprintf(stderr, "%s: invalid vertex buffer\n", __FUNCTION__);
nuclear@1: 		return;
nuclear@1: 	}
nuclear@1: 	if(global_sdr_loc[MESH_ATTR_VERTEX] == -1) {
nuclear@1: 		fprintf(stderr, "%s: shader attribute location for vertices unset\n", __FUNCTION__);
nuclear@1: 		return;
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	for(int i=0; i<NUM_MESH_ATTR; i++) {
nuclear@1: 		int loc = global_sdr_loc[i];
nuclear@1: 		if(loc >= 0 && vattr[i].vbo_valid) {
nuclear@1: 			glBindBuffer(GL_ARRAY_BUFFER, vattr[i].vbo);
nuclear@1: 			glVertexAttribPointer(loc, vattr[i].nelem, GL_FLOAT, GL_FALSE, 0, 0);
nuclear@1: 			glEnableVertexAttribArray(loc);
nuclear@1: 		}
nuclear@1: 	}
nuclear@1: 	glBindBuffer(GL_ARRAY_BUFFER, 0);
nuclear@1: 
nuclear@1: 	glBindBuffer(GL_ELEMENT_ARRAY_BUFFER, wire_ibo);
nuclear@1: 	glDrawElements(GL_LINES, nfaces * 6, GL_UNSIGNED_INT, 0);
nuclear@1: 	glBindBuffer(GL_ELEMENT_ARRAY_BUFFER, 0);
nuclear@1: 
nuclear@1: 	for(int i=0; i<NUM_MESH_ATTR; i++) {
nuclear@1: 		int loc = global_sdr_loc[i];
nuclear@1: 		if(loc >= 0 && vattr[i].vbo_valid) {
nuclear@1: 			glDisableVertexAttribArray(loc);
nuclear@1: 		}
nuclear@1: 	}
nuclear@1: }
nuclear@1: 
nuclear@1: void Mesh::draw_vertices() const
nuclear@1: {
nuclear@1: 	((Mesh*)this)->update_buffers();
nuclear@1: 
nuclear@1: 	if(!vattr[MESH_ATTR_VERTEX].vbo_valid) {
nuclear@1: 		fprintf(stderr, "%s: invalid vertex buffer\n", __FUNCTION__);
nuclear@1: 		return;
nuclear@1: 	}
nuclear@1: 	if(global_sdr_loc[MESH_ATTR_VERTEX] == -1) {
nuclear@1: 		fprintf(stderr, "%s: shader attribute location for vertices unset\n", __FUNCTION__);
nuclear@1: 		return;
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	for(int i=0; i<NUM_MESH_ATTR; i++) {
nuclear@1: 		int loc = global_sdr_loc[i];
nuclear@1: 		if(loc >= 0 && vattr[i].vbo_valid) {
nuclear@1: 			glBindBuffer(GL_ARRAY_BUFFER, vattr[i].vbo);
nuclear@1: 			glVertexAttribPointer(loc, vattr[i].nelem, GL_FLOAT, GL_FALSE, 0, 0);
nuclear@1: 			glEnableVertexAttribArray(loc);
nuclear@1: 		}
nuclear@1: 	}
nuclear@1: 	glBindBuffer(GL_ARRAY_BUFFER, 0);
nuclear@1: 
nuclear@1: 	glDrawArrays(GL_POINTS, 0, nverts);
nuclear@1: 
nuclear@1: 	for(int i=0; i<NUM_MESH_ATTR; i++) {
nuclear@1: 		int loc = global_sdr_loc[i];
nuclear@1: 		if(loc >= 0 && vattr[i].vbo_valid) {
nuclear@1: 			glDisableVertexAttribArray(loc);
nuclear@1: 		}
nuclear@1: 	}
nuclear@1: }
nuclear@1: 
nuclear@1: void Mesh::draw_normals() const
nuclear@1: {
nuclear@1: #ifdef USE_OLDGL
nuclear@1: 	int vert_loc = global_sdr_loc[MESH_ATTR_VERTEX];
nuclear@1: 	Vector3 *varr = (Vector3*)get_attrib_data(MESH_ATTR_VERTEX);
nuclear@1: 	Vector3 *norm = (Vector3*)get_attrib_data(MESH_ATTR_NORMAL);
nuclear@1: 
nuclear@1: 	if(!varr || !norm || vert_loc < 0) {
nuclear@1: 		return;
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	glBegin(GL_LINES);
nuclear@1: 	for(size_t i=0; i<nverts; i++) {
nuclear@1: 		glVertexAttrib3f(vert_loc, varr[i].x, varr[i].y, varr[i].z);
nuclear@1: 		Vector3 end = varr[i] + norm[i] * vis_vecsize;
nuclear@1: 		glVertexAttrib3f(vert_loc, end.x, end.y, end.z);
nuclear@1: 	}
nuclear@1: 	glEnd();
nuclear@1: 
nuclear@1: #endif	// USE_OLDGL
nuclear@1: }
nuclear@1: 
nuclear@1: void Mesh::draw_tangents() const
nuclear@1: {
nuclear@1: #ifdef USE_OLDGL
nuclear@1: 	int vert_loc = global_sdr_loc[MESH_ATTR_VERTEX];
nuclear@1: 	Vector3 *varr = (Vector3*)get_attrib_data(MESH_ATTR_VERTEX);
nuclear@1: 	Vector3 *tang = (Vector3*)get_attrib_data(MESH_ATTR_TANGENT);
nuclear@1: 
nuclear@1: 	if(!varr || !tang || vert_loc < 0) {
nuclear@1: 		return;
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	glBegin(GL_LINES);
nuclear@1: 	for(size_t i=0; i<nverts; i++) {
nuclear@1: 		glVertexAttrib3f(vert_loc, varr[i].x, varr[i].y, varr[i].z);
nuclear@1: 		Vector3 end = varr[i] + tang[i] * vis_vecsize;
nuclear@1: 		glVertexAttrib3f(vert_loc, end.x, end.y, end.z);
nuclear@1: 	}
nuclear@1: 	glEnd();
nuclear@1: 
nuclear@1: #endif	// USE_OLDGL
nuclear@1: }
nuclear@1: 
nuclear@1: #if 0
nuclear@1: void Mesh::get_aabbox(Vector3 *vmin, Vector3 *vmax) const
nuclear@1: {
nuclear@1: 	if(!aabb_valid) {
nuclear@1: 		((Mesh*)this)->calc_aabb();
nuclear@1: 	}
nuclear@1: 	*vmin = aabb.min;
nuclear@1: 	*vmax = aabb.max;
nuclear@1: }
nuclear@1: 
nuclear@1: const AABox &Mesh::get_aabbox() const
nuclear@1: {
nuclear@1: 	if(!aabb_valid) {
nuclear@1: 		((Mesh*)this)->calc_aabb();
nuclear@1: 	}
nuclear@1: 	return aabb;
nuclear@1: }
nuclear@1: 
nuclear@1: float Mesh::get_bsphere(Vector3 *center, float *rad) const
nuclear@1: {
nuclear@1: 	if(!bsph_valid) {
nuclear@1: 		((Mesh*)this)->calc_bsph();
nuclear@1: 	}
nuclear@1: 	*center = bsph.center;
nuclear@1: 	*rad = bsph.radius;
nuclear@1: 	return bsph.radius;
nuclear@1: }
nuclear@1: 
nuclear@1: const Sphere &Mesh::get_bsphere() const
nuclear@1: {
nuclear@1: 	if(!bsph_valid) {
nuclear@1: 		((Mesh*)this)->calc_bsph();
nuclear@1: 	}
nuclear@1: 	return bsph;
nuclear@1: }
nuclear@1: 
nuclear@1: /// static function
nuclear@1: void Mesh::set_intersect_mode(unsigned int mode)
nuclear@1: {
nuclear@1: 	Mesh::intersect_mode = mode;
nuclear@1: }
nuclear@1: 
nuclear@1: /// static function
nuclear@1: unsigned int Mesh::get_intersect_mode()
nuclear@1: {
nuclear@1: 	return Mesh::intersect_mode;
nuclear@1: }
nuclear@1: 
nuclear@1: /// static function
nuclear@1: void Mesh::set_vertex_select_distance(float dist)
nuclear@1: {
nuclear@1: 	Mesh::vertex_sel_dist = dist;
nuclear@1: }
nuclear@1: 
nuclear@1: /// static function
nuclear@1: float Mesh::get_vertex_select_distance()
nuclear@1: {
nuclear@1: 	return Mesh::vertex_sel_dist;
nuclear@1: }
nuclear@1: 
nuclear@1: /*bool Mesh::intersect(const Ray &ray, HitPoint *hit) const
nuclear@1: {
nuclear@1: 	assert((Mesh::intersect_mode & (ISECT_VERTICES | ISECT_FACE)) != (ISECT_VERTICES | ISECT_FACE));
nuclear@1: 
nuclear@1: 	const Vector3 *varr = (Vector3*)get_attrib_data(MESH_ATTR_VERTEX);
nuclear@1: 	const Vector3 *narr = (Vector3*)get_attrib_data(MESH_ATTR_NORMAL);
nuclear@1: 	if(!varr) {
nuclear@1: 		return false;
nuclear@1: 	}
nuclear@1: 	const unsigned int *idxarr = get_index_data();
nuclear@1: 
nuclear@1: 	// first test with the bounding box
nuclear@1: 	AABox box;
nuclear@1: 	get_aabbox(&box.min, &box.max);
nuclear@1: 	if(!box.intersect(ray)) {
nuclear@1: 		return false;
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	HitPoint nearest_hit;
nuclear@1: 	nearest_hit.dist = FLT_MAX;
nuclear@1: 	nearest_hit.obj = 0;
nuclear@1: 
nuclear@1: 	if(Mesh::intersect_mode & ISECT_VERTICES) {
nuclear@1: 		// we asked for "intersections" with the vertices of the mesh
nuclear@1: 		long nearest_vidx = -1;
nuclear@1: 		float thres_sq = Mesh::vertex_sel_dist * Mesh::vertex_sel_dist;
nuclear@1: 
nuclear@1: 		for(unsigned int i=0; i<nverts; i++) {
nuclear@1: 
nuclear@1: 			if((Mesh::intersect_mode & ISECT_FRONT) && dot_product(narr[i], ray.dir) > 0) {
nuclear@1: 				continue;
nuclear@1: 			}
nuclear@1: 
nuclear@1: 			// project the vertex onto the ray line
nuclear@1: 			float t = dot_product(varr[i] - ray.origin, ray.dir);
nuclear@1: 			Vector3 vproj = ray.origin + ray.dir * t;
nuclear@1: 
nuclear@1: 			float dist_sq = (vproj - varr[i]).length_sq();
nuclear@1: 			if(dist_sq < thres_sq) {
nuclear@1: 				if(!hit) {
nuclear@1: 					return true;
nuclear@1: 				}
nuclear@1: 				if(t < nearest_hit.dist) {
nuclear@1: 					nearest_hit.dist = t;
nuclear@1: 					nearest_vidx = i;
nuclear@1: 				}
nuclear@1: 			}
nuclear@1: 		}
nuclear@1: 
nuclear@1: 		if(nearest_vidx != -1) {
nuclear@1: 			hitvert = varr[nearest_vidx];
nuclear@1: 			nearest_hit.obj = &hitvert;
nuclear@1: 		}
nuclear@1: 
nuclear@1: 	} else {
nuclear@1: 		// regular intersection test with polygons
nuclear@1: 
nuclear@1: 		for(unsigned int i=0; i<nfaces; i++) {
nuclear@1: 			Triangle face(i, varr, idxarr);
nuclear@1: 
nuclear@1: 			// ignore back-facing polygons if the mode flags include ISECT_FRONT
nuclear@1: 			if((Mesh::intersect_mode & ISECT_FRONT) && dot_product(face.get_normal(), ray.dir) > 0) {
nuclear@1: 				continue;
nuclear@1: 			}
nuclear@1: 
nuclear@1: 			HitPoint fhit;
nuclear@1: 			if(face.intersect(ray, hit ? &fhit : 0)) {
nuclear@1: 				if(!hit) {
nuclear@1: 					return true;
nuclear@1: 				}
nuclear@1: 				if(fhit.dist < nearest_hit.dist) {
nuclear@1: 					nearest_hit = fhit;
nuclear@1: 					hitface = face;
nuclear@1: 				}
nuclear@1: 			}
nuclear@1: 		}
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	if(nearest_hit.obj) {
nuclear@1: 		if(hit) {
nuclear@1: 			*hit = nearest_hit;
nuclear@1: 
nuclear@1: 			// if we are interested in the mesh and not the faces set obj to this
nuclear@1: 			if(Mesh::intersect_mode & ISECT_FACE) {
nuclear@1: 				hit->obj = &hitface;
nuclear@1: 			} else if(Mesh::intersect_mode & ISECT_VERTICES) {
nuclear@1: 				hit->obj = &hitvert;
nuclear@1: 			} else {
nuclear@1: 				hit->obj = this;
nuclear@1: 			}
nuclear@1: 		}
nuclear@1: 		return true;
nuclear@1: 	}
nuclear@1: 	return false;
nuclear@1: }*/
nuclear@1: 
nuclear@1: 
nuclear@1: // ------ private member functions ------
nuclear@1: 
nuclear@1: void Mesh::calc_aabb()
nuclear@1: {
nuclear@1: 	// the cast is to force calling the const version which doesn't invalidate
nuclear@1: 	if(!((const Mesh*)this)->get_attrib_data(MESH_ATTR_VERTEX)) {
nuclear@1: 		return;
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	aabb.min = Vector3(FLT_MAX, FLT_MAX, FLT_MAX);
nuclear@1: 	aabb.max = -aabb.min;
nuclear@1: 
nuclear@1: 	for(unsigned int i=0; i<nverts; i++) {
nuclear@1: 		Vector4 v = get_attrib(MESH_ATTR_VERTEX, i);
nuclear@1: 		for(int j=0; j<3; j++) {
nuclear@1: 			if(v[j] < aabb.min[j]) {
nuclear@1: 				aabb.min[j] = v[j];
nuclear@1: 			}
nuclear@1: 			if(v[j] > aabb.max[j]) {
nuclear@1: 				aabb.max[j] = v[j];
nuclear@1: 			}
nuclear@1: 		}
nuclear@1: 	}
nuclear@1: 	aabb_valid = true;
nuclear@1: }
nuclear@1: 
nuclear@1: void Mesh::calc_bsph()
nuclear@1: {
nuclear@1: 	// the cast is to force calling the const version which doesn't invalidate
nuclear@1: 	if(!((const Mesh*)this)->get_attrib_data(MESH_ATTR_VERTEX)) {
nuclear@1: 		return;
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	Vector3 v;
nuclear@1: 	bsph.center = Vector3(0, 0, 0);
nuclear@1: 
nuclear@1: 	// first find the center
nuclear@1: 	for(unsigned int i=0; i<nverts; i++) {
nuclear@1: 		v = get_attrib(MESH_ATTR_VERTEX, i);
nuclear@1: 		bsph.center += v;
nuclear@1: 	}
nuclear@1: 	bsph.center /= (float)nverts;
nuclear@1: 
nuclear@1: 	bsph.radius = 0.0f;
nuclear@1: 	for(unsigned int i=0; i<nverts; i++) {
nuclear@1: 		v = get_attrib(MESH_ATTR_VERTEX, i);
nuclear@1: 		float dist_sq = (v - bsph.center).length_sq();
nuclear@1: 		if(dist_sq > bsph.radius) {
nuclear@1: 			bsph.radius = dist_sq;
nuclear@1: 		}
nuclear@1: 	}
nuclear@1: 	bsph.radius = sqrt(bsph.radius);
nuclear@1: 
nuclear@1: 	bsph_valid = true;
nuclear@1: }
nuclear@1: #endif
nuclear@1: 
nuclear@1: void Mesh::update_buffers()
nuclear@1: {
nuclear@1: 	for(int i=0; i<NUM_MESH_ATTR; i++) {
nuclear@1: 		if(has_attrib(i) && !vattr[i].vbo_valid) {
nuclear@1: 			glBindBuffer(GL_ARRAY_BUFFER, vattr[i].vbo);
nuclear@1: 			glBufferData(GL_ARRAY_BUFFER, nverts * vattr[i].nelem * sizeof(float), &vattr[i].data[0], GL_STATIC_DRAW);
nuclear@1: 			vattr[i].vbo_valid = true;
nuclear@1: 		}
nuclear@1: 	}
nuclear@1: 	glBindBuffer(GL_ARRAY_BUFFER, 0);
nuclear@1: 
nuclear@1: 	if(idata_valid && !ibo_valid) {
nuclear@1: 		glBindBuffer(GL_ELEMENT_ARRAY_BUFFER, ibo);
nuclear@1: 		glBufferData(GL_ELEMENT_ARRAY_BUFFER, nfaces * 3 * sizeof(unsigned int), &idata[0], GL_STATIC_DRAW);
nuclear@1: 		ibo_valid = true;
nuclear@1: 	}
nuclear@1: 	glBindBuffer(GL_ELEMENT_ARRAY_BUFFER, 0);
nuclear@1: }
nuclear@1: 
nuclear@1: void Mesh::update_wire_ibo()
nuclear@1: {
nuclear@1: 	update_buffers();
nuclear@1: 
nuclear@1: 	if(wire_ibo_valid) {
nuclear@1: 		return;
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	if(!wire_ibo) {
nuclear@1: 		glGenBuffers(1, &wire_ibo);
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	unsigned int *wire_idxarr = new unsigned int[nfaces * 6];
nuclear@1: 	unsigned int *dest = wire_idxarr;
nuclear@1: 
nuclear@1: 	if(ibo_valid) {
nuclear@1: 		// we're dealing with an indexed mesh
nuclear@1: 		const unsigned int *idxarr = ((const Mesh*)this)->get_index_data();
nuclear@1: 
nuclear@1: 		for(unsigned int i=0; i<nfaces; i++) {
nuclear@1: 			*dest++ = idxarr[0];
nuclear@1: 			*dest++ = idxarr[1];
nuclear@1: 			*dest++ = idxarr[1];
nuclear@1: 			*dest++ = idxarr[2];
nuclear@1: 			*dest++ = idxarr[2];
nuclear@1: 			*dest++ = idxarr[0];
nuclear@1: 			idxarr += 3;
nuclear@1: 		}
nuclear@1: 	} else {
nuclear@1: 		// not an indexed mesh ...
nuclear@1: 		for(unsigned int i=0; i<nfaces; i++) {
nuclear@1: 			int vidx = i * 3;
nuclear@1: 			*dest++ = vidx;
nuclear@1: 			*dest++ = vidx + 1;
nuclear@1: 			*dest++ = vidx + 1;
nuclear@1: 			*dest++ = vidx + 2;
nuclear@1: 			*dest++ = vidx + 2;
nuclear@1: 			*dest++ = vidx;
nuclear@1: 		}
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	glBindBuffer(GL_ELEMENT_ARRAY_BUFFER, wire_ibo);
nuclear@1: 	glBufferData(GL_ELEMENT_ARRAY_BUFFER, nfaces * 6 * sizeof(unsigned int), wire_idxarr, GL_STATIC_DRAW);
nuclear@1: 	delete [] wire_idxarr;
nuclear@1: 	wire_ibo_valid = true;
nuclear@1: 	glBindBuffer(GL_ELEMENT_ARRAY_BUFFER, 0);
nuclear@1: }
nuclear@1: 
nuclear@1: 
nuclear@1: // ------ class Triangle ------
nuclear@1: Triangle::Triangle()
nuclear@1: {
nuclear@1: 	normal_valid = false;
nuclear@1: 	id = -1;
nuclear@1: }
nuclear@1: 
nuclear@1: Triangle::Triangle(const Vector3 &v0, const Vector3 &v1, const Vector3 &v2)
nuclear@1: {
nuclear@1: 	v[0] = v0;
nuclear@1: 	v[1] = v1;
nuclear@1: 	v[2] = v2;
nuclear@1: 	normal_valid = false;
nuclear@1: 	id = -1;
nuclear@1: }
nuclear@1: 
nuclear@1: Triangle::Triangle(int n, const Vector3 *varr, const unsigned int *idxarr)
nuclear@1: {
nuclear@1: 	if(idxarr) {
nuclear@1: 		v[0] = varr[idxarr[n * 3]];
nuclear@1: 		v[1] = varr[idxarr[n * 3 + 1]];
nuclear@1: 		v[2] = varr[idxarr[n * 3 + 2]];
nuclear@1: 	} else {
nuclear@1: 		v[0] = varr[n * 3];
nuclear@1: 		v[1] = varr[n * 3 + 1];
nuclear@1: 		v[2] = varr[n * 3 + 2];
nuclear@1: 	}
nuclear@1: 	normal_valid = false;
nuclear@1: 	id = n;
nuclear@1: }
nuclear@1: 
nuclear@1: void Triangle::calc_normal()
nuclear@1: {
nuclear@1: 	normal = cross_product(v[1] - v[0], v[2] - v[0]).normalized();
nuclear@1: 	normal_valid = true;
nuclear@1: }
nuclear@1: 
nuclear@1: const Vector3 &Triangle::get_normal() const
nuclear@1: {
nuclear@1: 	if(!normal_valid) {
nuclear@1: 		((Triangle*)this)->calc_normal();
nuclear@1: 	}
nuclear@1: 	return normal;
nuclear@1: }
nuclear@1: 
nuclear@1: void Triangle::transform(const Matrix4x4 &xform)
nuclear@1: {
nuclear@1: 	v[0].transform(xform);
nuclear@1: 	v[1].transform(xform);
nuclear@1: 	v[2].transform(xform);
nuclear@1: 	normal_valid = false;
nuclear@1: }
nuclear@1: 
nuclear@1: void Triangle::draw() const
nuclear@1: {
nuclear@1: 	Vector3 n[3];
nuclear@1: 	n[0] = get_normal();
nuclear@1: 	n[1] = get_normal();
nuclear@1: 	n[2] = get_normal();
nuclear@1: 
nuclear@1: 	int vloc = Mesh::get_attrib_location(MESH_ATTR_VERTEX);
nuclear@1: 	int nloc = Mesh::get_attrib_location(MESH_ATTR_NORMAL);
nuclear@1: 
nuclear@1: 	glEnableVertexAttribArray(vloc);
nuclear@1: 	glVertexAttribPointer(vloc, 3, GL_FLOAT, GL_FALSE, 0, &v[0].x);
nuclear@1: 	glVertexAttribPointer(nloc, 3, GL_FLOAT, GL_FALSE, 0, &n[0].x);
nuclear@1: 
nuclear@1: 	glDrawArrays(GL_TRIANGLES, 0, 3);
nuclear@1: 
nuclear@1: 	glDisableVertexAttribArray(vloc);
nuclear@1: 	glDisableVertexAttribArray(nloc);
nuclear@1: 	CHECKGLERR;
nuclear@1: }
nuclear@1: 
nuclear@1: void Triangle::draw_wire() const
nuclear@1: {
nuclear@1: 	static const int idxarr[] = {0, 1, 1, 2, 2, 0};
nuclear@1: 	int vloc = Mesh::get_attrib_location(MESH_ATTR_VERTEX);
nuclear@1: 
nuclear@1: 	glEnableVertexAttribArray(vloc);
nuclear@1: 	glVertexAttribPointer(vloc, 3, GL_FLOAT, GL_FALSE, 0, &v[0].x);
nuclear@1: 
nuclear@1: 	glDrawElements(GL_LINES, 6, GL_UNSIGNED_INT, idxarr);
nuclear@1: 
nuclear@1: 	glDisableVertexAttribArray(vloc);
nuclear@1: 	CHECKGLERR;
nuclear@1: }
nuclear@1: 
nuclear@1: Vector3 Triangle::calc_barycentric(const Vector3 &pos) const
nuclear@1: {
nuclear@1: 	Vector3 norm = get_normal();
nuclear@1: 
nuclear@1: 	float area_sq = fabs(dot_product(cross_product(v[1] - v[0], v[2] - v[0]), norm));
nuclear@1: 	if(area_sq < 1e-5) {
nuclear@1: 		return Vector3(0, 0, 0);
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	float asq0 = fabs(dot_product(cross_product(v[1] - pos, v[2] - pos), norm));
nuclear@1: 	float asq1 = fabs(dot_product(cross_product(v[2] - pos, v[0] - pos), norm));
nuclear@1: 	float asq2 = fabs(dot_product(cross_product(v[0] - pos, v[1] - pos), norm));
nuclear@1: 
nuclear@1: 	return Vector3(asq0 / area_sq, asq1 / area_sq, asq2 / area_sq);
nuclear@1: }
nuclear@1: 
nuclear@1: /*bool Triangle::intersect(const Ray &ray, HitPoint *hit) const
nuclear@1: {
nuclear@1: 	Vector3 normal = get_normal();
nuclear@1: 
nuclear@1: 	float ndotdir = dot_product(ray.dir, normal);
nuclear@1: 	if(fabs(ndotdir) < 1e-4) {
nuclear@1: 		return false;
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	Vector3 vertdir = v[0] - ray.origin;
nuclear@1: 	float t = dot_product(normal, vertdir) / ndotdir;
nuclear@1: 
nuclear@1: 	Vector3 pos = ray.origin + ray.dir * t;
nuclear@1: 	Vector3 bary = calc_barycentric(pos);
nuclear@1: 
nuclear@1: 	if(bary.x + bary.y + bary.z > 1.00001) {
nuclear@1: 		return false;
nuclear@1: 	}
nuclear@1: 
nuclear@1: 	if(hit) {
nuclear@1: 		hit->dist = t;
nuclear@1: 		hit->pos = ray.origin + ray.dir * t;
nuclear@1: 		hit->normal = normal;
nuclear@1: 		hit->obj = this;
nuclear@1: 	}
nuclear@1: 	return true;
nuclear@1: }*/