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?