Multiplication of two quaternions
r_times_q = quaternion_multiply(r,q)
quaternion_multiply multiplies two quaternions.
quaternion_multiply
q = rand(4,1); p = rand(4,1); // these two routines should give the same result: quaternion_multiply(p,q) quaternion_matrix(p) * q
http://siptoolbox.sf.net