vector vb = (inverse of matrix dm) * vector va (double precision) (n.b. The matrix must be unitary, as this routine assumes that the inverse and transpose are identical)
dm double[3][3] matrix va double[3] vector
vb double[3] result vector P.T.Wallace Starlink 31 October 1993