I am working with rotation matrices R[3,3] from SO(3), and anti-symmetric matrices from the map R^3->so(3).
These latter come from Omega=(x,y,z) so:
0 -z y
z 0 -x
-y x 0
I would like to be able to use these in expressions where I sometimes use the vector, and sometimes use the matrix.
Is there a convenient function to do this?
It needs to be compatible with NDSolve.