We can use the 3x3 rotation matrix R to determine the effect of rotations about the principle axes:
Vafter = R.Vbefore
That's straightforward.
Does anyone know what the strategy should be for solving for unknown rotations r={rx,ry,rx} radians, given we know Vbefore and Vafter?