X-Git-Url: http://git.mutantstargoat.com/user/nuclear/?p=dosdemo;a=blobdiff_plain;f=tools%2Fropesim%2Fsrc%2Fmain.c;h=3c616f67f8f9ece4379dfb682e6586f89fc9dcde;hp=e6221bcb1d74deb70604d7b88efa4c3de27c2d20;hb=1ae209a77434dcdf7266d2cf95d3f20c0de7a0e9;hpb=4688bc5b987bdcc392620f8d699070dda8813359

diff --git a/tools/ropesim/src/main.c b/tools/ropesim/src/main.c
index e6221bc..3c616f6 100644
--- a/tools/ropesim/src/main.c
+++ b/tools/ropesim/src/main.c
@@ -3,6 +3,7 @@
 #include <GL/glut.h>
 #include "cmesh.h"
 #include "cgmath/cgmath.h"
+#include "ropesim.h"
 
 int init(void);
 void cleanup(void);
@@ -19,6 +20,7 @@ void sball_button(int bn, int st);
 float cam_theta, cam_phi, cam_dist = 10;
 int prev_mx, prev_my;
 int bnstate[8];
+int modkeys;
 
 long start_msec;
 
@@ -29,10 +31,12 @@ cgm_vec3 gmove;
 /*cgm_quat grot = {0, 0, 0, 1};*/
 float grot_theta, grot_phi;
 float ginner_xform[16], gouter_xform[16];
-cgm_vec3 ganchor[4];
+cgm_vec3 ganchor[4], manchor[4];
 
 cgm_vec3 dbgvec[4];
 
+struct rsim_world rsim;
+
 int main(int argc, char **argv)
 {
 	glutInit(&argc, argv);
@@ -60,13 +64,19 @@ int main(int argc, char **argv)
 	return 0;
 }
 
+#define ROPE_MASSES			5
+#define ROPE_SPRINGS		(ROPE_MASSES - 1)
+#define ROPE_LEN			0.6f
+#define ROPE_MASSES_MASS	0.1f
+#define ROPE_K				80.0f
 
 int init(void)
 {
 	static const char *meshnames[] = {"suzanne", "gimbal_outer", "gimbal_inner"};
 	static struct cmesh **meshes[] = {&mesh_suz, &mesh_gout, &mesh_gin};
 	static const float amb[] = {0.05, 0.05, 0.08, 1};
-	int i;
+	int i, j;
+	struct rsim_rope *rope, *ropes_tail;
 
 	glEnable(GL_CULL_FACE);
 	glEnable(GL_DEPTH_TEST);
@@ -95,11 +105,49 @@ int init(void)
 		cmesh_remove_submesh(scn, idx);
 	}
 
+	rsim_init(&rsim);
+	ropes_tail = 0;
+
 	/* anchor points on the inner gimbal */
 	for(i=0; i<4; i++) {
 		ganchor[i].x = (float)(((i & 1) << 1) - 1) * 1.5f;
 		ganchor[i].y = (float)((i & 2) - 1) * 1.5f;
 		ganchor[i].z = 0;
+
+		manchor[i] = ganchor[i];
+		cgm_vscale(manchor + i, 0.32);
+
+
+
+		manchor[i].y += 0.15;
+
+		/* create a rope hanging from the anchor point */
+		if(!(rope = rsim_alloc_rope(ROPE_MASSES, ROPE_SPRINGS))) {
+			fprintf(stderr, "failed to allocate rope\n");
+			return -1;
+		}
+		for(j=0; j<ROPE_MASSES; j++) {
+			float t = (float)j / (float)(ROPE_MASSES - 1.0f);
+			cgm_vlerp(&rope->masses[j].p, ganchor + i, manchor + i, t);
+			rope->masses[j].m = 0.1f;
+
+			if(j < ROPE_SPRINGS) {
+				rope->springs[j].rest_len = ROPE_LEN / ROPE_SPRINGS;
+				rope->springs[j].k = ROPE_K;
+				rope->springs[j].mass[0] = rope->masses + j;
+				rope->springs[j].mass[1] = rope->masses + j + 1;
+			}
+		}
+		rsim_freeze_rope_mass(rope, rope->masses);	/* freeze first mass */
+		rsim_freeze_rope_mass(rope, rope->masses + j - 1);	/* freeze last mass */
+
+		if(!ropes_tail) {
+			rsim.ropes = ropes_tail = rope;
+		} else {
+			ropes_tail->next = rope;
+			ropes_tail = rope;
+		}
+		rope->next = 0;
 	}
 
 	return 0;
@@ -111,6 +159,8 @@ void cleanup(void)
 	cmesh_free(mesh_gout);
 	cmesh_free(mesh_gin);
 	cmesh_free(scn);
+
+	rsim_destroy(&rsim);
 }
 
 void update(long tmsec, float dt)
@@ -118,6 +168,7 @@ void update(long tmsec, float dt)
 	int i;
 	cgm_vec3 apt0, apt1;
 	float theta, phi, brot;
+	struct rsim_rope *rope;
 
 	/*
 	cgm_mrotation_quat(ginner_xform, &grot);
@@ -130,12 +181,18 @@ void update(long tmsec, float dt)
 	cgm_mrotation_euler(ginner_xform, phi, theta, 0, CGM_EULER_XYZ);
 	cgm_mrotation_euler(gouter_xform, phi, 0, 0, CGM_EULER_XYZ);
 
+	rope = rsim.ropes;
 	for(i=0; i<4; i++) {
 		apt0 = ganchor[i];
 		cgm_vmul_m4v3(&apt0, ginner_xform);
 
 		dbgvec[i] = apt0;
+		rope->masses[0].p = apt0;
+
+		rope = rope->next;
 	}
+
+	rsim_step(&rsim, dt);
 }
 
 void display(void)
@@ -153,6 +210,7 @@ void display(void)
 	int i, count;
 	long tmsec = glutGet(GLUT_ELAPSED_TIME) - start_msec;
 	static long prev_tmsec;
+	struct rsim_rope *rope;
 
 	update(tmsec, (float)(tmsec - prev_tmsec) / 1000.0f);
 	prev_tmsec = tmsec;
@@ -194,6 +252,30 @@ void display(void)
 	}
 	glEnd();
 
+	glPushAttrib(GL_ENABLE_BIT);
+	glDisable(GL_LIGHTING);
+	glLineWidth(2);
+	glPointSize(5);
+
+	rope = rsim.ropes;
+	while(rope) {
+		glBegin(GL_LINE_STRIP);
+		glColor3f(0.2, 1, 0.2);
+		for(i=0; i<rope->num_masses; i++) {
+			glVertex3f(rope->masses[i].p.x, rope->masses[i].p.y, rope->masses[i].p.z);
+		}
+		glEnd();
+
+		glBegin(GL_POINTS);
+		glColor3f(1, 0.2, 0.2);
+		for(i=0; i<rope->num_masses; i++) {
+			glVertex3f(rope->masses[i].p.x, rope->masses[i].p.y, rope->masses[i].p.z);
+		}
+		glEnd();
+		rope = rope->next;
+	}
+	glPopAttrib();
+
 	glutSwapBuffers();
 }
 void idle(void)
@@ -222,6 +304,7 @@ void mouse(int bn, int st, int x, int y)
 	prev_mx = x;
 	prev_my = y;
 	bnstate[bn - GLUT_LEFT_BUTTON] = st == GLUT_DOWN;
+	modkeys = glutGetModifiers();
 }
 
 void motion(int x, int y)
@@ -233,21 +316,23 @@ void motion(int x, int y)
 
 	if(!(dx | dy)) return;
 
-	if(bnstate[0]) {
-		cam_theta += dx * 0.5;
-		cam_phi += dy * 0.5;
-		if(cam_phi < -90) cam_phi = -90;
-		if(cam_phi > 90) cam_phi = 90;
-	}
-
-	if(bnstate[1]) {
-		grot_theta += dx * 0.5;
-		grot_phi += dy * 0.5;
-	}
+	if(modkeys) {
+		if(bnstate[0]) {
+			grot_theta += dx * 0.5;
+			grot_phi += dy * 0.5;
+		}
+	} else {
+		if(bnstate[0]) {
+			cam_theta += dx * 0.5;
+			cam_phi += dy * 0.5;
+			if(cam_phi < -90) cam_phi = -90;
+			if(cam_phi > 90) cam_phi = 90;
+		}
 
-	if(bnstate[2]) {
-		cam_dist += dy * 0.1;
-		if(cam_dist < 0.0f) cam_dist = 0.0f;
+		if(bnstate[2]) {
+			cam_dist += dy * 0.1;
+			if(cam_dist < 0.0f) cam_dist = 0.0f;
+		}
 	}
 }