public void qmult( ref Quat c, ref Quat a, ref Quat b )
{
c.r = a.r * b.r - a.i * b.i - a.j * b.j - a.k * b.k;
c.i = a.r * b.i + a.i * b.r + a.j * b.k - a.k * b.j;
c.j = a.r * b.j + a.j * b.r + a.k * b.i - a.i * b.k;
c.k = a.r * b.k + a.k * b.r + a.i * b.j - a.j * b.i;
}