diff --git a/stepmania/src/mathlib.c b/stepmania/src/mathlib.c
index c8193253d6..33793f70fc 100644
--- a/stepmania/src/mathlib.c
+++ b/stepmania/src/mathlib.c
@@ -69,28 +69,7 @@ void VectorInverse (vec3_t v)
 	v[2] = -v[2];
 }
 
-void ClearBounds (vec3_t mins, vec3_t maxs)
-{
-	mins[0] = mins[1] = mins[2] = 99999;
-	maxs[0] = maxs[1] = maxs[2] = -99999;
-}
-
-void AddPointToBounds (vec3_t v, vec3_t mins, vec3_t maxs)
-{
-	int		i;
-	vec_t	val;
-
-	for (i=0 ; i<3 ; i++)
-	{
-		val = v[i];
-		if (val < mins[i])
-			mins[i] = val;
-		if (val > maxs[i])
-			maxs[i] = val;
-	}
-}
-
-void AngleMatrix (const vec3_t angles, float (*matrix)[4] )
+void AngleMatrix (const vec3_t angles, float matrix[3][4] )
 {
 	float		angle;
 	float		sr, sp, sy, cr, cp, cy;
@@ -120,36 +99,6 @@ void AngleMatrix (const vec3_t angles, float (*matrix)[4] )
 	matrix[2][3] = 0.0;
 }
 
-void AngleIMatrix (const vec3_t angles, float matrix[3][4] )
-{
-	float		angle;
-	float		sr, sp, sy, cr, cp, cy;
-	
-	angle = angles[2] * (Q_PI*2 / 360);
-	sy = sin(angle);
-	cy = cos(angle);
-	angle = angles[1] * (Q_PI*2 / 360);
-	sp = sin(angle);
-	cp = cos(angle);
-	angle = angles[0] * (Q_PI*2 / 360);
-	sr = sin(angle);
-	cr = cos(angle);
-
-	// matrix = (Z * Y) * X
-	matrix[0][0] = cp*cy;
-	matrix[0][1] = cp*sy;
-	matrix[0][2] = -sp;
-	matrix[1][0] = sr*sp*cy+cr*-sy;
-	matrix[1][1] = sr*sp*sy+cr*cy;
-	matrix[1][2] = sr*cp;
-	matrix[2][0] = (cr*sp*cy+-sr*-sy);
-	matrix[2][1] = (cr*sp*sy+-sr*cy);
-	matrix[2][2] = cr*cp;
-	matrix[0][3] = 0.0;
-	matrix[1][3] = 0.0;
-	matrix[2][3] = 0.0;
-}
-
 void R_ConcatTransforms (const float in1[3][4], const float in2[3][4], float out[3][4])
 {
 	out[0][0] = in1[0][0] * in2[0][0] + in1[0][1] * in2[1][0] +
@@ -203,93 +152,3 @@ void VectorTransform (const vec3_t in1, const float in2[3][4], vec3_t out)
 	out[1] = DotProduct(in1, in2[1]) +	in2[1][3];
 	out[2] = DotProduct(in1, in2[2]) +	in2[2][3];
 }
-
-
-
-void AngleQuaternion( const vec3_t angles, vec4_t quaternion )
-{
-	float		angle;
-	float		sr, sp, sy, cr, cp, cy;
-
-	// FIXME: rescale the inputs to 1/2 angle
-	angle = angles[2] * 0.5;
-	sy = sin(angle);
-	cy = cos(angle);
-	angle = angles[1] * 0.5;
-	sp = sin(angle);
-	cp = cos(angle);
-	angle = angles[0] * 0.5;
-	sr = sin(angle);
-	cr = cos(angle);
-
-	quaternion[0] = sr*cp*cy-cr*sp*sy; // X
-	quaternion[1] = cr*sp*cy+sr*cp*sy; // Y
-	quaternion[2] = cr*cp*sy-sr*sp*cy; // Z
-	quaternion[3] = cr*cp*cy+sr*sp*sy; // W
-}
-
-void QuaternionMatrix( const vec4_t quaternion, float (*matrix)[4] )
-{
-
-	matrix[0][0] = 1.0 - 2.0 * quaternion[1] * quaternion[1] - 2.0 * quaternion[2] * quaternion[2];
-	matrix[1][0] = 2.0 * quaternion[0] * quaternion[1] + 2.0 * quaternion[3] * quaternion[2];
-	matrix[2][0] = 2.0 * quaternion[0] * quaternion[2] - 2.0 * quaternion[3] * quaternion[1];
-
-	matrix[0][1] = 2.0 * quaternion[0] * quaternion[1] - 2.0 * quaternion[3] * quaternion[2];
-	matrix[1][1] = 1.0 - 2.0 * quaternion[0] * quaternion[0] - 2.0 * quaternion[2] * quaternion[2];
-	matrix[2][1] = 2.0 * quaternion[1] * quaternion[2] + 2.0 * quaternion[3] * quaternion[0];
-
-	matrix[0][2] = 2.0 * quaternion[0] * quaternion[2] + 2.0 * quaternion[3] * quaternion[1];
-	matrix[1][2] = 2.0 * quaternion[1] * quaternion[2] - 2.0 * quaternion[3] * quaternion[0];
-	matrix[2][2] = 1.0 - 2.0 * quaternion[0] * quaternion[0] - 2.0 * quaternion[1] * quaternion[1];
-}
-
-void QuaternionSlerp( const vec4_t p, vec4_t q, float t, vec4_t qt )
-{
-	int i;
-	float omega, cosom, sinom, sclp, sclq;
-
-	// decide if one of the quaternions is backwards
-	float a = 0;
-	float b = 0;
-	for (i = 0; i < 4; i++) {
-		a += (p[i]-q[i])*(p[i]-q[i]);
-		b += (p[i]+q[i])*(p[i]+q[i]);
-	}
-	if (a > b) {
-		for (i = 0; i < 4; i++) {
-			q[i] = -q[i];
-		}
-	}
-
-	cosom = p[0]*q[0] + p[1]*q[1] + p[2]*q[2] + p[3]*q[3];
-
-	if ((1.0 + cosom) > 0.00000001) {
-		if ((1.0 - cosom) > 0.00000001) {
-			omega = acos( cosom );
-			sinom = sin( omega );
-			sclp = sin( (1.0 - t)*omega) / sinom;
-			sclq = sin( t*omega ) / sinom;
-		}
-		else {
-			sclp = 1.0 - t;
-			sclq = t;
-		}
-		for (i = 0; i < 4; i++) {
-			qt[i] = sclp * p[i] + sclq * q[i];
-		}
-	}
-	else {
-		qt[0] = -p[1];
-		qt[1] = p[0];
-		qt[2] = -p[3];
-		qt[3] = p[2];
-		sclp = sin( (1.0 - t) * 0.5 * Q_PI);
-		sclq = sin( t * 0.5 * Q_PI);
-		for (i = 0; i < 3; i++) {
-			qt[i] = sclp * p[i] + sclq * qt[i];
-		}
-	}
-}
-
-
diff --git a/stepmania/src/mathlib.h b/stepmania/src/mathlib.h
index c2b9288e09..de751586b9 100644
--- a/stepmania/src/mathlib.h
+++ b/stepmania/src/mathlib.h
@@ -39,11 +39,7 @@ int VectorCompare (vec3_t v1, vec3_t v2);
 vec_t VectorNormalize (vec3_t v);
 void VectorInverse (vec3_t v);
 
-void ClearBounds (vec3_t mins, vec3_t maxs);
-void AddPointToBounds (vec3_t v, vec3_t mins, vec3_t maxs);
-
 void AngleMatrix (const vec3_t angles, float matrix[3][4] );
-void AngleIMatrix (const vec3_t angles, float matrix[3][4] );
 void R_ConcatTransforms (const float in1[3][4], const float in2[3][4], float out[3][4]);
 
 void VectorIRotate (const vec3_t in1, const float in2[3][4], vec3_t out);
@@ -51,11 +47,6 @@ void VectorRotate (const vec3_t in1, const float in2[3][4], vec3_t out);
 
 void VectorTransform (const vec3_t in1, const float in2[3][4], vec3_t out);
 
-void AngleQuaternion( const vec3_t angles, vec4_t quaternion );
-void QuaternionMatrix( const vec4_t quaternion, float (*matrix)[4] );
-void QuaternionSlerp( const vec4_t p, vec4_t q, float t, vec4_t qt );
-
-
 #ifdef __cplusplus
 }
 #endif